Note: This page was generated from a Jupyter notebook which can be found at docs/tutorial_notebooks/3_model_iterations_and_plots.ipynb in your DYNAMITE directory. —-

3. Output Overview

By the end of the notebook, you will have run a grid of Schwarzschild models. This will involve the following steps: 1. run a small grid of orbit-based models with DYNAMITE 2. understand the structure of the output 3. plot the output 4. look at the iteration process while running the grid of models

You should run this from the directory docs/tutorial_notebooks.

Setup

We will run DYNAMITE on CALIFA data of NGC 6278. To prepare the input data files, you should first run the tutorial “Data Preparation for Gauss Hermite kinematics” (data_prep_for_gauss_hermites.ipynb). The relevant files you need for this tutorial are:

| tutorial_notebooks
| ├── NGC6278_input
| │   ├── dynamite_input
| │   │   ├── gauss_hermite_kins.ecsv
| │   │   ├── aperture.dat
| │   │   ├── bins.dat
| │   │   ├── mge.ecsv
| │   │   └── ...
| │   └── ...
| │   └── ...
| └── NGC6278_config.yaml
| └── *.ipynb
|

Read the configuration file

[1]:

import dynamite as dyn

print('DYNAMITE')
print('    version', dyn.__version__)
print('    installed at ', dyn.__path__)

fname = 'NGC6278_config.yaml'
c = dyn.config_reader.Configuration(fname, reset_logging=True)

[INFO] 11:56:41 - dynamite.config_reader.Configuration - Config file NGC6278_config.yaml read.
[INFO] 11:56:41 - dynamite.config_reader.Configuration - io_settings...
[INFO] 11:56:41 - dynamite.config_reader.Configuration - Output directory tree: NGC6278_output/.
[INFO] 11:56:41 - dynamite.config_reader.Configuration - system_attributes...
[INFO] 11:56:41 - dynamite.config_reader.Configuration - model_components...

DYNAMITE
    version 2.0.0
    installed at  ['/Users/leahskusa/Desktop/uni/dynamics_assistant/dynamite/dynamite_test/DYN_TEST/lib/python3.9/site-packages/dynamite-2.0.0-py3.9.egg/dynamite']

[INFO] 11:56:42 - dynamite.config_reader.Configuration - system_parameters...
[INFO] 11:56:42 - dynamite.config_reader.Configuration - orblib_settings...
[INFO] 11:56:42 - dynamite.config_reader.Configuration - weight_solver_settings...
[INFO] 11:56:42 - dynamite.config_reader.Configuration - Will attempt to recover partially run models.
[INFO] 11:56:42 - dynamite.config_reader.Configuration - parameter_space_settings...
[INFO] 11:56:42 - dynamite.config_reader.Configuration - multiprocessing_settings...
[INFO] 11:56:42 - dynamite.config_reader.Configuration - ... using 4 CPUs for orbit integration.
[INFO] 11:56:42 - dynamite.config_reader.Configuration - ... using 4 CPUs for weight solving.
[INFO] 11:56:42 - dynamite.config_reader.Configuration - legacy_settings...
[INFO] 11:56:42 - dynamite.config_reader.Configuration - System assembled
[INFO] 11:56:42 - dynamite.config_reader.Configuration - Configuration validated
[INFO] 11:56:42 - dynamite.config_reader.Configuration - Instantiated parameter space
[INFO] 11:56:42 - dynamite.model.AllModels - Previous models have been found: Reading NGC6278_output/all_models.ecsv into AllModels.table
[INFO] 11:56:42 - dynamite.config_reader.Configuration - Instantiated AllModels object

All the options in the configuration file are held in the object c. For example, let’s look at the io_settings. Output from this tutorial will be saved in the output_directory.

[2]:

# delete previous output if available
c.remove_existing_orblibs()
c.remove_existing_all_models_file(wipe_other_files=False)
c.backup_config_file(keep=3, delete_other=True)

[INFO] 11:56:42 - dynamite.config_reader.Configuration - Model output tree NGC6278_output/models/ removed.
[INFO] 11:56:42 - dynamite.config_reader.Configuration - Deleted existing NGC6278_output/all_models.ecsv.
[INFO] 11:56:42 - dynamite.model.AllModels - No previous models (file NGC6278_output/all_models.ecsv) have been found: Making an empty table in AllModels.table
[INFO] 11:56:42 - dynamite.config_reader.Configuration - Instantiated empty AllModels object
[INFO] 11:56:42 - dynamite.config_reader.Configuration - Config file backup: NGC6278_output/NGC6278_config_001.yaml.

[3]:

c.settings.io_settings

[3]:

{'input_directory': 'NGC6278_input/dynamite_input/',
 'output_directory': 'NGC6278_output/',
 'all_models_file': 'all_models.ecsv',
 'model_directory': 'NGC6278_output/models/',
 'plot_directory': 'NGC6278_output/plots/'}

In fact, by creating the configuration object c, we have also created the output_directory and copied a version of the configuration file there,

[4]:

ls NGC6278_output

NGC6278_config_000.yaml  NGC6278_config_001.yaml  plots/

Run the models

Making the ModelIterator object will start running a grid of orbit-based models. This next step will take about 5-10 minutes using 4 cpus

[5]:

import time

t = time.perf_counter()

smi = dyn.model_iterator.ModelIterator(config=c)

delta_t = time.perf_counter()-t
print(f'Computation time: {delta_t} seconds = {delta_t/60} minutes')

[INFO] 11:56:53 - dynamite.model_iterator.ModelIterator - LegacyGridSearch: "iteration 0"
[INFO] 11:56:53 - dynamite.parameter_space.LegacyGridSearch - LegacyGridSearch added 1 new model(s) out of 1
[INFO] 11:56:53 - dynamite.model_iterator.ModelInnerIterator - ... running model 1 out of 1
[INFO] 11:56:53 - dynamite.orblib.LegacyOrbitLibrary - Calculating initial conditions
[INFO] 11:58:34 - dynamite.orblib.LegacyOrbitLibrary - ...done - cmd_orb_start exit code 0. Logfile: NGC6278_output/models/orblib_000_000/datfil/orbstart.log.
[INFO] 11:58:34 - dynamite.orblib.LegacyOrbitLibrary - Integrating orbit library tube orbits
[INFO] 11:58:49 - dynamite.orblib.LegacyOrbitLibrary - ...done - cmd_tube_orbs exit code 0. Logfiles: NGC6278_output/models/orblib_000_000/datfil/orblib.log, NGC6278_output/models/orblib_000_000/datfil/triaxmass.log, NGC6278_output/models/orblib_000_000/datfil/triaxmassbin.log.
[INFO] 11:58:49 - dynamite.orblib.LegacyOrbitLibrary - Integrating orbit library box orbits
[INFO] 11:59:02 - dynamite.orblib.LegacyOrbitLibrary - ...done - cmd_box_orbs exit code 0. Logfile: NGC6278_output/models/orblib_000_000/datfil/orblibbox.log.
[INFO] 11:59:02 - dynamite.weight_solvers.LegacyWeightSolver - Using WeightSolver: LegacyWeightSolver
[INFO] 11:59:02 - dynamite.weight_solvers.LegacyWeightSolver - Fitting orbit library to the kinematic data: ml5.00
[INFO] 11:59:06 - dynamite.weight_solvers.LegacyWeightSolver - ...done, NNLS problem solved -  cmd_nnls_5.0 exit code 0. Logfile: NGC6278_output/models/orblib_000_000/ml5.00/nnls.log
[INFO] 11:59:06 - dynamite.config_reader.Configuration - Config file copied to NGC6278_output/models/orblib_000_000/ml5.00/NGC6278_config.yaml.
[INFO] 11:59:08 - dynamite.model_iterator.ModelInnerIterator - Model 1: NGC6278_output/models/orblib_000_000/ml5.00/model_done_staging.ecsv written.
[INFO] 11:59:08 - dynamite.model_iterator.ModelInnerIterator - Iteration done, staging files deleted.
[INFO] 11:59:08 - dynamite.plotter.Plotter - kinchi2 vs. model id plot created (1 models).
[INFO] 11:59:08 - dynamite.plotter.Plotter - Plot NGC6278_output/plots/kinchi2_progress_plot.png saved in NGC6278_output/plots/
[INFO] 11:59:08 - dynamite.plotter.Plotter - Making chi2 plot scaled according to kinchi2
[INFO] 11:59:09 - dynamite.plotter.Plotter - Plot NGC6278_output/plots/kinchi2_plot.png saved in NGC6278_output/plots/
[INFO] 11:59:09 - dynamite.plotter.Plotter - Plotting kinematic maps for 1 kin_sets.
[INFO] 11:59:09 - dynamite.plotter.Plotter - Plotting kinematic maps for kin_set no 0: califa
/Users/leahskusa/Desktop/uni/dynamics_assistant/dynamite/dynamite_test/DYN_TEST/lib/python3.9/site-packages/plotbin-3.1.3-py3.9.egg/plotbin/sauron_colormap.py:105: UserWarning: Trying to register the cmap 'sauron' which already exists.
/Users/leahskusa/Desktop/uni/dynamics_assistant/dynamite/dynamite_test/DYN_TEST/lib/python3.9/site-packages/plotbin-3.1.3-py3.9.egg/plotbin/sauron_colormap.py:106: UserWarning: Trying to register the cmap 'sauron_r' which already exists.
[INFO] 11:59:19 - dynamite.model_iterator.ModelIterator - LegacyGridSearch: "iteration 1"
[INFO] 11:59:19 - dynamite.parameter_space.LegacyGridSearch - LegacyGridSearch added 4 new model(s) out of 4
[INFO] 11:59:19 - dynamite.model_iterator.ModelInnerIterator - ... running model 1 out of 4
[INFO] 11:59:19 - dynamite.model_iterator.ModelInnerIterator - ... running model 2 out of 4
[INFO] 11:59:19 - dynamite.orblib.LegacyOrbitLibrary - Calculating initial conditions
[INFO] 11:59:19 - dynamite.orblib.LegacyOrbitLibrary - Calculating initial conditions
[INFO] 12:01:07 - dynamite.orblib.LegacyOrbitLibrary - ...done - cmd_orb_start exit code 0. Logfile: NGC6278_output/models/orblib_001_001/datfil/orbstart.log.
[INFO] 12:01:07 - dynamite.orblib.LegacyOrbitLibrary - Integrating orbit library tube orbits
[INFO] 12:01:08 - dynamite.orblib.LegacyOrbitLibrary - ...done - cmd_orb_start exit code 0. Logfile: NGC6278_output/models/orblib_001_000/datfil/orbstart.log.
[INFO] 12:01:08 - dynamite.orblib.LegacyOrbitLibrary - Integrating orbit library tube orbits
[INFO] 12:01:25 - dynamite.orblib.LegacyOrbitLibrary - ...done - cmd_tube_orbs exit code 0. Logfiles: NGC6278_output/models/orblib_001_001/datfil/orblib.log, NGC6278_output/models/orblib_001_001/datfil/triaxmass.log, NGC6278_output/models/orblib_001_001/datfil/triaxmassbin.log.
[INFO] 12:01:25 - dynamite.orblib.LegacyOrbitLibrary - Integrating orbit library box orbits
[INFO] 12:01:25 - dynamite.orblib.LegacyOrbitLibrary - ...done - cmd_tube_orbs exit code 0. Logfiles: NGC6278_output/models/orblib_001_000/datfil/orblib.log, NGC6278_output/models/orblib_001_000/datfil/triaxmass.log, NGC6278_output/models/orblib_001_000/datfil/triaxmassbin.log.
[INFO] 12:01:25 - dynamite.orblib.LegacyOrbitLibrary - Integrating orbit library box orbits
[INFO] 12:01:40 - dynamite.orblib.LegacyOrbitLibrary - ...done - cmd_box_orbs exit code 0. Logfile: NGC6278_output/models/orblib_001_000/datfil/orblibbox.log.
[INFO] 12:01:40 - dynamite.weight_solvers.LegacyWeightSolver - Using WeightSolver: LegacyWeightSolver
[INFO] 12:01:40 - dynamite.weight_solvers.LegacyWeightSolver - Fitting orbit library to the kinematic data: ml5.00
[INFO] 12:01:40 - dynamite.orblib.LegacyOrbitLibrary - ...done - cmd_box_orbs exit code 0. Logfile: NGC6278_output/models/orblib_001_001/datfil/orblibbox.log.
[INFO] 12:01:40 - dynamite.weight_solvers.LegacyWeightSolver - Using WeightSolver: LegacyWeightSolver
[INFO] 12:01:40 - dynamite.weight_solvers.LegacyWeightSolver - Fitting orbit library to the kinematic data: ml5.00
[INFO] 12:01:46 - dynamite.weight_solvers.LegacyWeightSolver - ...done, NNLS problem solved -  cmd_nnls_5.0 exit code 0. Logfile: NGC6278_output/models/orblib_001_000/ml5.00/nnls.log
[INFO] 12:01:46 - dynamite.config_reader.Configuration - Config file copied to NGC6278_output/models/orblib_001_000/ml5.00/NGC6278_config.yaml.
[INFO] 12:01:46 - dynamite.weight_solvers.LegacyWeightSolver - ...done, NNLS problem solved -  cmd_nnls_5.0 exit code 0. Logfile: NGC6278_output/models/orblib_001_001/ml5.00/nnls.log
[INFO] 12:01:46 - dynamite.config_reader.Configuration - Config file copied to NGC6278_output/models/orblib_001_001/ml5.00/NGC6278_config.yaml.
[INFO] 12:01:47 - dynamite.model_iterator.ModelInnerIterator - Model 1: NGC6278_output/models/orblib_001_000/ml5.00/model_done_staging.ecsv written.
[INFO] 12:01:47 - dynamite.model_iterator.ModelInnerIterator - Model 2: NGC6278_output/models/orblib_001_001/ml5.00/model_done_staging.ecsv written.
[INFO] 12:01:47 - dynamite.model_iterator.ModelInnerIterator - ... running model 3 out of 4
[INFO] 12:01:47 - dynamite.model_iterator.ModelInnerIterator - ... running model 4 out of 4
[INFO] 12:01:47 - dynamite.weight_solvers.LegacyWeightSolver - Using WeightSolver: LegacyWeightSolver
[INFO] 12:01:47 - dynamite.weight_solvers.LegacyWeightSolver - Fitting orbit library to the kinematic data: ml1.00
[INFO] 12:01:48 - dynamite.weight_solvers.LegacyWeightSolver - Using WeightSolver: LegacyWeightSolver
[INFO] 12:01:48 - dynamite.weight_solvers.LegacyWeightSolver - Fitting orbit library to the kinematic data: ml9.00
[INFO] 12:01:54 - dynamite.weight_solvers.LegacyWeightSolver - ...done, NNLS problem solved -  cmd_nnls_1.0 exit code 0. Logfile: NGC6278_output/models/orblib_000_000/ml1.00/nnls.log
[INFO] 12:01:54 - dynamite.config_reader.Configuration - Config file copied to NGC6278_output/models/orblib_000_000/ml1.00/NGC6278_config.yaml.
[INFO] 12:01:54 - dynamite.weight_solvers.LegacyWeightSolver - ...done, NNLS problem solved -  cmd_nnls_9.0 exit code 0. Logfile: NGC6278_output/models/orblib_000_000/ml9.00/nnls.log
[INFO] 12:01:54 - dynamite.config_reader.Configuration - Config file copied to NGC6278_output/models/orblib_000_000/ml9.00/NGC6278_config.yaml.
[INFO] 12:01:56 - dynamite.model_iterator.ModelInnerIterator - Model 3: NGC6278_output/models/orblib_000_000/ml1.00/model_done_staging.ecsv written.
[INFO] 12:01:56 - dynamite.model_iterator.ModelInnerIterator - Model 4: NGC6278_output/models/orblib_000_000/ml9.00/model_done_staging.ecsv written.
[INFO] 12:01:56 - dynamite.model_iterator.ModelInnerIterator - Iteration done, staging files deleted.
[INFO] 12:01:56 - dynamite.plotter.Plotter - kinchi2 vs. model id plot created (5 models).
[INFO] 12:01:56 - dynamite.plotter.Plotter - Plot NGC6278_output/plots/kinchi2_progress_plot.png saved in NGC6278_output/plots/
[INFO] 12:01:56 - dynamite.plotter.Plotter - Making chi2 plot scaled according to kinchi2
[INFO] 12:01:57 - dynamite.plotter.Plotter - Plot NGC6278_output/plots/kinchi2_plot.png saved in NGC6278_output/plots/
[INFO] 12:01:57 - dynamite.plotter.Plotter - Plotting kinematic maps for 1 kin_sets.
[INFO] 12:01:57 - dynamite.plotter.Plotter - Plotting kinematic maps for kin_set no 0: califa
/Users/leahskusa/Desktop/uni/dynamics_assistant/dynamite/dynamite_test/DYN_TEST/lib/python3.9/site-packages/plotbin-3.1.3-py3.9.egg/plotbin/sauron_colormap.py:105: UserWarning: Trying to register the cmap 'sauron' which already exists.
/Users/leahskusa/Desktop/uni/dynamics_assistant/dynamite/dynamite_test/DYN_TEST/lib/python3.9/site-packages/plotbin-3.1.3-py3.9.egg/plotbin/sauron_colormap.py:106: UserWarning: Trying to register the cmap 'sauron_r' which already exists.
[INFO] 12:02:09 - dynamite.model_iterator.ModelIterator - LegacyGridSearch: "iteration 2"
[INFO] 12:02:09 - dynamite.parameter_space.LegacyGridSearch - LegacyGridSearch added 2 new model(s) out of 2
[INFO] 12:02:09 - dynamite.model_iterator.ModelInnerIterator - ... running model 1 out of 2
[INFO] 12:02:09 - dynamite.model_iterator.ModelInnerIterator - ... running model 2 out of 2
[INFO] 12:02:09 - dynamite.weight_solvers.LegacyWeightSolver - Using WeightSolver: LegacyWeightSolver
[INFO] 12:02:09 - dynamite.weight_solvers.LegacyWeightSolver - Fitting orbit library to the kinematic data: ml1.00
[INFO] 12:02:09 - dynamite.weight_solvers.LegacyWeightSolver - Using WeightSolver: LegacyWeightSolver
[INFO] 12:02:09 - dynamite.weight_solvers.LegacyWeightSolver - Fitting orbit library to the kinematic data: ml9.00
[INFO] 12:02:15 - dynamite.weight_solvers.LegacyWeightSolver - ...done, NNLS problem solved -  cmd_nnls_1.0 exit code 0. Logfile: NGC6278_output/models/orblib_001_001/ml1.00/nnls.log
[INFO] 12:02:15 - dynamite.config_reader.Configuration - Config file copied to NGC6278_output/models/orblib_001_001/ml1.00/NGC6278_config.yaml.
[INFO] 12:02:15 - dynamite.weight_solvers.LegacyWeightSolver - ...done, NNLS problem solved -  cmd_nnls_9.0 exit code 0. Logfile: NGC6278_output/models/orblib_001_001/ml9.00/nnls.log
[INFO] 12:02:15 - dynamite.config_reader.Configuration - Config file copied to NGC6278_output/models/orblib_001_001/ml9.00/NGC6278_config.yaml.
[INFO] 12:02:16 - dynamite.model_iterator.ModelInnerIterator - Model 1: NGC6278_output/models/orblib_001_001/ml1.00/model_done_staging.ecsv written.
[INFO] 12:02:17 - dynamite.model_iterator.ModelInnerIterator - Model 2: NGC6278_output/models/orblib_001_001/ml9.00/model_done_staging.ecsv written.
[INFO] 12:02:17 - dynamite.model_iterator.ModelInnerIterator - Iteration done, staging files deleted.
[INFO] 12:02:17 - dynamite.plotter.Plotter - kinchi2 vs. model id plot created (7 models).
[INFO] 12:02:17 - dynamite.plotter.Plotter - Plot NGC6278_output/plots/kinchi2_progress_plot.png saved in NGC6278_output/plots/
[INFO] 12:02:17 - dynamite.plotter.Plotter - Making chi2 plot scaled according to kinchi2
[INFO] 12:02:17 - dynamite.plotter.Plotter - Plot NGC6278_output/plots/kinchi2_plot.png saved in NGC6278_output/plots/
[INFO] 12:02:17 - dynamite.plotter.Plotter - Plotting kinematic maps for 1 kin_sets.
[INFO] 12:02:17 - dynamite.plotter.Plotter - Plotting kinematic maps for kin_set no 0: califa
/Users/leahskusa/Desktop/uni/dynamics_assistant/dynamite/dynamite_test/DYN_TEST/lib/python3.9/site-packages/plotbin-3.1.3-py3.9.egg/plotbin/sauron_colormap.py:105: UserWarning: Trying to register the cmap 'sauron' which already exists.
/Users/leahskusa/Desktop/uni/dynamics_assistant/dynamite/dynamite_test/DYN_TEST/lib/python3.9/site-packages/plotbin-3.1.3-py3.9.egg/plotbin/sauron_colormap.py:106: UserWarning: Trying to register the cmap 'sauron_r' which already exists.
[INFO] 12:02:29 - dynamite.model_iterator.ModelIterator - LegacyGridSearch: "iteration 3"
[INFO] 12:02:29 - dynamite.parameter_space.LegacyGridSearch - LegacyGridSearch added 3 new model(s) out of 3
[INFO] 12:02:29 - dynamite.model_iterator.ModelInnerIterator - ... running model 1 out of 3
[INFO] 12:02:29 - dynamite.orblib.LegacyOrbitLibrary - Calculating initial conditions
[INFO] 12:04:21 - dynamite.orblib.LegacyOrbitLibrary - ...done - cmd_orb_start exit code 0. Logfile: NGC6278_output/models/orblib_003_000/datfil/orbstart.log.
[INFO] 12:04:21 - dynamite.orblib.LegacyOrbitLibrary - Integrating orbit library tube orbits
[INFO] 12:04:40 - dynamite.orblib.LegacyOrbitLibrary - ...done - cmd_tube_orbs exit code 0. Logfiles: NGC6278_output/models/orblib_003_000/datfil/orblib.log, NGC6278_output/models/orblib_003_000/datfil/triaxmass.log, NGC6278_output/models/orblib_003_000/datfil/triaxmassbin.log.
[INFO] 12:04:40 - dynamite.orblib.LegacyOrbitLibrary - Integrating orbit library box orbits
[INFO] 12:04:55 - dynamite.orblib.LegacyOrbitLibrary - ...done - cmd_box_orbs exit code 0. Logfile: NGC6278_output/models/orblib_003_000/datfil/orblibbox.log.
[INFO] 12:04:55 - dynamite.weight_solvers.LegacyWeightSolver - Using WeightSolver: LegacyWeightSolver
[INFO] 12:04:55 - dynamite.weight_solvers.LegacyWeightSolver - Fitting orbit library to the kinematic data: ml5.00
[INFO] 12:05:01 - dynamite.weight_solvers.LegacyWeightSolver - ...done, NNLS problem solved -  cmd_nnls_5.0 exit code 0. Logfile: NGC6278_output/models/orblib_003_000/ml5.00/nnls.log
[INFO] 12:05:01 - dynamite.config_reader.Configuration - Config file copied to NGC6278_output/models/orblib_003_000/ml5.00/NGC6278_config.yaml.
[INFO] 12:05:03 - dynamite.model_iterator.ModelInnerIterator - Model 1: NGC6278_output/models/orblib_003_000/ml5.00/model_done_staging.ecsv written.
[INFO] 12:05:03 - dynamite.model_iterator.ModelInnerIterator - ... running model 2 out of 3
[INFO] 12:05:03 - dynamite.model_iterator.ModelInnerIterator - ... running model 3 out of 3
[INFO] 12:05:03 - dynamite.weight_solvers.LegacyWeightSolver - Using WeightSolver: LegacyWeightSolver
[INFO] 12:05:03 - dynamite.weight_solvers.LegacyWeightSolver - Fitting orbit library to the kinematic data: ml3.00
[INFO] 12:05:03 - dynamite.weight_solvers.LegacyWeightSolver - Using WeightSolver: LegacyWeightSolver
[INFO] 12:05:03 - dynamite.weight_solvers.LegacyWeightSolver - Fitting orbit library to the kinematic data: ml7.00
[INFO] 12:05:10 - dynamite.weight_solvers.LegacyWeightSolver - ...done, NNLS problem solved -  cmd_nnls_7.0 exit code 0. Logfile: NGC6278_output/models/orblib_001_001/ml7.00/nnls.log
[INFO] 12:05:10 - dynamite.config_reader.Configuration - Config file copied to NGC6278_output/models/orblib_001_001/ml7.00/NGC6278_config.yaml.
[INFO] 12:05:10 - dynamite.weight_solvers.LegacyWeightSolver - ...done, NNLS problem solved -  cmd_nnls_3.0 exit code 0. Logfile: NGC6278_output/models/orblib_001_001/ml3.00/nnls.log
[INFO] 12:05:10 - dynamite.config_reader.Configuration - Config file copied to NGC6278_output/models/orblib_001_001/ml3.00/NGC6278_config.yaml.
[INFO] 12:05:11 - dynamite.model_iterator.ModelInnerIterator - Model 3: NGC6278_output/models/orblib_001_001/ml7.00/model_done_staging.ecsv written.
[INFO] 12:05:12 - dynamite.model_iterator.ModelInnerIterator - Model 2: NGC6278_output/models/orblib_001_001/ml3.00/model_done_staging.ecsv written.
[INFO] 12:05:12 - dynamite.model_iterator.ModelInnerIterator - Iteration done, staging files deleted.
[INFO] 12:05:12 - dynamite.plotter.Plotter - kinchi2 vs. model id plot created (10 models).
[INFO] 12:05:12 - dynamite.plotter.Plotter - Plot NGC6278_output/plots/kinchi2_progress_plot.png saved in NGC6278_output/plots/
[INFO] 12:05:12 - dynamite.plotter.Plotter - Making chi2 plot scaled according to kinchi2
[INFO] 12:05:12 - dynamite.plotter.Plotter - Plot NGC6278_output/plots/kinchi2_plot.png saved in NGC6278_output/plots/
[INFO] 12:05:12 - dynamite.plotter.Plotter - Plotting kinematic maps for 1 kin_sets.
[INFO] 12:05:12 - dynamite.plotter.Plotter - Plotting kinematic maps for kin_set no 0: califa
/Users/leahskusa/Desktop/uni/dynamics_assistant/dynamite/dynamite_test/DYN_TEST/lib/python3.9/site-packages/plotbin-3.1.3-py3.9.egg/plotbin/sauron_colormap.py:105: UserWarning: Trying to register the cmap 'sauron' which already exists.
/Users/leahskusa/Desktop/uni/dynamics_assistant/dynamite/dynamite_test/DYN_TEST/lib/python3.9/site-packages/plotbin-3.1.3-py3.9.egg/plotbin/sauron_colormap.py:106: UserWarning: Trying to register the cmap 'sauron_r' which already exists.

Computation time: 512.4313896360001 seconds = 8.540523160600001 minutes

The following files have been created in the models directory,

[6]:

ls NGC6278_output/models

orblib_000_000/ orblib_001_000/ orblib_001_001/ orblib_003_000/

Each directory holds a different orbit library

orblib_XXX_YYY

where XXX labels the iteration when it was created, and YYY labels the position within that iteration. Looking inside one of these directories, we see the following files:

[7]:

ls NGC6278_output/models/orblib_000_000

cmd_box_orbs   cmd_nnls_9.0   datfil/        interpolgrid   ml9.00/
cmd_nnls_1.0   cmd_orb_start  fort.30        ml1.00/
cmd_nnls_5.0   cmd_tube_orbs  infil/         ml5.00/

which are:

  • cmd_*: bash scripts for running Fortran programs

  • datfil/: directory holding the orbit library for the reference potential

  • infil/: input files for running Fortran programs

  • ml*/: directories containing output orbital weights (and other results) for different values of ml

Each ml* directory hold outputs for a re-scaled version of the same potential, where the value of ml is a mass scaling applied to a reference potential. The reference potential uses the the first value of ml encountered in the parameter search.

Some plots are automatically created,

[8]:

ls NGC6278_output/plots

anisotropy_var.png         kinchi2_progress_plot.png
betaz_var.png              kinematic_map_califa.png
chi2_plot.png              orbit_linear_only.png
enclosedmassm_linear.png   triaxial_qpt.png
kinchi2_plot.png

and they represent the following quantities:

  1. kinchi2_progress_plot : chi2 values vs model ID fbc249911dc042b9afffd1a9402696ec

  2. kinchi2_progress_plot : model parameters vs chi2 values. If more than 2 paramters were left free, this would be a traingle plot of chi2 values, 336c4c273bef4f2c8a66b91921960e97

  3. kinematic_map_califa.png : the kinematic maps for the current minimum-chi2 model 48ade32c5a2a462487eefc64b59353a8

A summary of all the models run so far is saved in the file NGC6278_output/all_models.ecsv. This is an Astropy ECSV file. A table holding this data is stored in c,

[9]:

c.all_models.table

[9]:
Table length=10
m-bha-bhc-dhf-dhq-starsp-starsu-starsmlchi2kinchi2time_modifiedorblib_doneweights_doneall_donewhich_iterdirectory
float64float64float64float64float64float64float64float64float64float64objectboolboolboolint64object
100000.00.0018.010.00.540.990.99995.036462.86802261109616324.4636723920962022-05-05T09:59:08.000TrueTrueTrue0orblib_000_000/ml5.00/
100000.00.0018.03.16227766016837950.540.990.99995.035574.41551156808417097.6742735452282022-05-05T10:01:47.000TrueTrueTrue1orblib_001_000/ml5.00/
100000.00.0018.031.6227766016837930.540.990.99995.037813.29817202782616009.4274496465442022-05-05T10:01:47.000TrueTrueTrue1orblib_001_001/ml5.00/
100000.00.0018.010.00.540.990.99991.0496485.33015648247361077.5145082952022-05-05T10:01:56.000TrueTrueTrue1orblib_000_000/ml1.00/
100000.00.0018.010.00.540.990.99999.057060.05166216243435834.601580949852022-05-05T10:01:56.000TrueTrueTrue1orblib_000_000/ml9.00/
100000.00.0018.031.6227766016837930.540.990.99991.0480137.15078561887345603.0541964542022-05-05T10:02:16.000TrueTrueTrue2orblib_001_001/ml1.00/
100000.00.0018.031.6227766016837930.540.990.99999.060138.7603498813839254.2978415724142022-05-05T10:02:16.000TrueTrueTrue2orblib_001_001/ml9.00/
100000.00.0018.017.782794100389230.540.990.99995.055173.9578457800717192.771321220812022-05-05T10:05:03.000TrueTrueTrue3orblib_003_000/ml5.00/
100000.00.0018.031.6227766016837930.540.990.99993.077632.8579272309150532.391745789882022-05-05T10:05:12.000TrueTrueTrue3orblib_001_001/ml3.00/
100000.00.0018.031.6227766016837930.540.990.99997.041558.80246617462421283.7540067106242022-05-05T10:05:11.000TrueTrueTrue3orblib_001_001/ml7.00/

At this stage, you could:

  • run more models, perhaps first adjusting settings in the configuration file, as explained in 2_quickstart.ipynb,

    • increasing the n_max_mods and/or n_max_iter

    • adjust parameter bounds and/or which parameters are kept free

  • plot other visualisations

Plotting

DYNAMITE provides other plotting methods in the Plotter:

[11]:

plotter = dyn.plotter.Plotter(config=c)

In all the functions that require to use the values of the \(\chi^2\) to make the plots, the user can choose which \(\chi^2\) to use, by specifying the value of the parameter which_chi2. The recommended value to use is which_chi2='kinchi2'.

The plots produced by the functions introduced below are all saved in the plot directory specified in the directory plots within the output directory specified in the configuration file. After running all the cells in this notebook, you will find your plots in the directory NGC6278_output/plots.

The mass_plot function generates a cumulative mass plot, showing the enclosed mass profiles for the mass-follows-light component (red), for the dark matter (blue), and for the sum of the two (black). The solid lines correspond to the best-fit model, the shaded areas represent 1 sigma uncertainties. You can specify the radial extent of the plot and the type of file you want to be saved with the figure (e.g., '.png', '.pdf', … if figtype=None, the default is used and a '.png' figure is created).

[12]:

fig1 = plotter.mass_plot(which_chi2='kinchi2', Rmax_arcs=50, figtype=None)

Selecting  10  models

[INFO] 12:08:37 - dynamite.plotter.Plotter - Plot NGC6278_output/plots/enclosedmassm_linear.png saved in NGC6278_output/plots/
../_images/tutorial_notebooks_3_model_iterations_and_plots_21_2.png

The orbit_plot function generates a plot showing the stellar orbit distribution, described as probability density of orbits; circularity (\(\lambda_z\)) is represented here as a function of the distance from the galactic centre r (in arcsec). You can specify the type of file you want to be saved with the figure (e.g., '.png', '.pdf', …). In this case, Rmax_arcs represents the upper radial limit for orbit selection, in arcsec, meaning that only orbits extending up to Rmax_arcs are plotted.

[13]:

fig2 = plotter.orbit_plot(Rmax_arcs=50)

[INFO] 12:08:38 - dynamite.plotter.Plotter - Plot NGC6278_output/plots/orbit_linear_only.png saved in NGC6278_output/plots/
../_images/tutorial_notebooks_3_model_iterations_and_plots_23_1.png

The beta_plot function generates two plots, showing the intrinsic and projected anisotropy profiles.

[14]:

fig3, fig4 = plotter.beta_plot(which_chi2='kinchi2', Rmax_arcs=50)

[INFO] 12:08:42 - dynamite.plotter.Plotter - Figure NGC6278_output/plots/anisotropy_var.png saved in NGC6278_output/plots/
[INFO] 12:08:42 - dynamite.plotter.Plotter - Figure NGC6278_output/plots/betaz_var.png saved in NGC6278_output/plots/
../_images/tutorial_notebooks_3_model_iterations_and_plots_25_1.png
../_images/tutorial_notebooks_3_model_iterations_and_plots_25_2.png

The qpu_plot function creates a plot showing the intrinsic flattenings \(q\) and \(p\), with the blue and black lines respectively, as a function of the distance from the galactic centre (in arcsec). The value of \(T = (1-p^2)/(1-q^2)\) is also shown (red line).

[15]:

fig5 = plotter.qpu_plot(which_chi2='kinchi2', Rmax_arcs=50,figtype =None)

[INFO] 12:08:43 - dynamite.plotter.Plotter - Plot NGC6278_output/plots/triaxial_qpt.png saved in NGC6278_output/plots/
../_images/tutorial_notebooks_3_model_iterations_and_plots_27_1.png

The Iteration Process

We will now have a look at the iteration process by examining the \(\chi^2\)-value.

In the table above (c.all_models.table), the model parameters are given for each step of the iteration. Most parameters are fixed and therefore constant, but the mass-to-light ratio ml and the dark matter parameter f-dh are free. We will examine how these parameters are changed with each iteration.

First, we get the upper and lower limits of the parameters from the configuration file (for more details on this step have a look at the notebook parameter_space.ipynb):

[16]:

import numpy as np
import matplotlib.pyplot as plt

# extract the lo/hi limits for the two free parameters f and ml

f = c.parspace.get_parameter_from_name('f-dh')
f_lims_raw = [f.par_generator_settings['lo'], f.par_generator_settings['hi']]
f_lims = [f.get_par_value_from_raw_value(lim0) for lim0 in f_lims_raw]

ml = c.parspace.get_parameter_from_name('ml')
ml_lims_raw = [ml.par_generator_settings['lo'], ml.par_generator_settings['hi']]
ml_lims = [ml.get_par_value_from_raw_value(lim0) for lim0 in ml_lims_raw]

From the table, we can get a list of all iterations that were made:

[17]:

# get list of iterations
iterations = np.unique(c.all_models.table['which_iter'])

Now plot the \(\chi^2\)-value of the models for every iteration:

[18]:

# plot chi2 of models vs iterations
for iter0 in iterations:
    table = c.all_models.table
    table = table[table['which_iter']==iter0]
    plt.figure(figsize=(6,5))
    plt.scatter(table['f-dh'],
                table['ml'],
                c=table['chi2'],
                cmap=plt.cm.viridis_r,
                s=200)


    cbar = plt.colorbar()
    cbar.set_label('$\chi^2$')
    plt.gca().set_title(f'iteration {iter0}')

    plt.gca().set_xlim(*f_lims)
    plt.gca().set_ylim(*ml_lims)
    plt.gca().set_xlabel(f.LaTeX)
    plt.gca().set_ylabel(ml.LaTeX)
    plt.gca().set_xscale('log')
    plt.tight_layout()
    plt.show()
../_images/tutorial_notebooks_3_model_iterations_and_plots_33_0.png
../_images/tutorial_notebooks_3_model_iterations_and_plots_33_1.png
../_images/tutorial_notebooks_3_model_iterations_and_plots_33_2.png
../_images/tutorial_notebooks_3_model_iterations_and_plots_33_3.png

We can see how the parameter values get closer together for each iteration. The next plot gives an overview over all iterations:

[19]:

# plot the models: f-dh vs ml altogether
plt.scatter(c.all_models.table['f-dh'],
            c.all_models.table['ml'],
            c=c.all_models.table['chi2'] - np.min(c.all_models.table['chi2']),
            cmap=plt.cm.viridis_r,
            s=200)
cbar = plt.colorbar()
cbar.set_label('$\chi^2$')
plt.gca().set_title(f'all iterations')
plt.gca().set_xlim(*f_lims)
plt.gca().set_ylim(*ml_lims)
plt.gca().set_xlabel(f.LaTeX)
plt.gca().set_ylabel(ml.LaTeX)
plt.gca().set_xscale('log')
plt.tight_layout()
plt.show()
../_images/tutorial_notebooks_3_model_iterations_and_plots_35_0.png

We can see how with every iteration, the \(\chi^2\)-value becomes smaller, as does the scatter between the values.

The plot kinchi2_progress_plot already shown above is another representation of the iteration process. The different models are plotted on a grid in the parameter space, where the best fit model is marked with a cross.

4a4c2310ea244081881ab630ed10eab4

[ ]: