.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "auto_examples/plot_multi_temperature_fit.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note Run this example online : - Click :ref:`here ` to download the full example code - Then start `Radis-Lab `__, upload the Jupyter notebook, and run it from there. .. rst-class:: sphx-glr-example-title .. _sphx_glr_auto_examples_plot_multi_temperature_fit.py: .. _example_multi_temperature_fit: ===================== Multi-temperature Fit ===================== A method to fit an experimental spectrum directly from :py:class:`~radis.lbl.factory.SpectrumFactory`, with :py:meth:`~radis.lbl.factory.SpectrumFactory.fit_spectrum` Typical output is similar to the `radis-examples Multi-temperature fit `__ : .. image:: https://raw.githubusercontent.com/radis/radis-examples/master/docs/multi-temperature-fit.gif The method requires a fitting model. An example model is provided in :py:mod:`radis.tools.fitting` : :py:func:`~radis.tools.fitting.Tvib12Tvib3Trot_NonLTEModel`. Other models can be used, such as in the :ref:`one-temperature fit example ` More advanced tools for interactive fitting of multi-dimensional, multi-slabs spectra can be found in :py:mod:`fitroom`. Finally, the :ref:`GPU-accelerated example` shows how to obtain real-time interactive spectra. .. GENERATED FROM PYTHON SOURCE LINES 27-32 .. code-block:: default from os.path import join from radis import Spectrum, SpectrumFactory .. GENERATED FROM PYTHON SOURCE LINES 33-34 Get Fitted Data .. GENERATED FROM PYTHON SOURCE LINES 34-50 .. code-block:: default from radis.test.utils import getValidationCase, setup_test_line_databases from radis.tools.fitting import Tvib12Tvib3Trot_NonLTEModel setup_test_line_databases() # Data from Dang, adapted by Klarenaar, digitized by us s_exp = Spectrum.from_txt( getValidationCase( join("test_CO2_3Tvib_vs_klarenaar_data", "klarenaar_2017_digitized_data.csv") ), "transmittance_noslit", wunit="cm-1", unit="", delimiter=",", name="Klarenaar 2017", ) .. rst-class:: sphx-glr-script-out .. code-block:: none /home/docs/checkouts/readthedocs.org/user_builds/radis/envs/master/lib/python3.8/site-packages/radis/spectrum/spectrum.py:4400: UserWarning: Wavespace is not evenly spaced (0.000%) for transmittance_noslit. This may create problems if later convolving with slit function (`s.apply_slit()`). You can use `s.resample_even()` .. GENERATED FROM PYTHON SOURCE LINES 51-52 Calculate .. GENERATED FROM PYTHON SOURCE LINES 52-89 .. code-block:: default sf = SpectrumFactory( 2284.2, 2284.6, wstep=0.001, # cm-1 pressure=20 * 1e-3, # bar db_use_cached=True, lvl_use_cached=True, cutoff=1e-25, isotope="1,2", path_length=10, # cm-1 mole_fraction=0.1 * 28.97 / 44.07, truncation=1, # cm-1 medium="vacuum", export_populations=None, # 'vib', # parsum_mode="tabulation" ) sf.warnings["MissingSelfBroadeningWarning"] = "ignore" sf.warnings["PerformanceWarning"] = "ignore" sf.load_databank("HITEMP-CO2-TEST") s_best, best = sf.fit_spectrum( s_exp.take("transmittance_noslit"), model=Tvib12Tvib3Trot_NonLTEModel, fit_parameters={ "T12": 517, "T3": 2641, "Trot": 491, }, bounds={"T12": [300, 2000], "T3": [300, 5000], "Trot": [300, 2000]}, fixed_parameters={"vib_distribution": "treanor"}, plot=True, solver_options={ "method": "TNC", "maxiter": 80, # 👈 increase to let the fit converge }, ) .. rst-class:: sphx-glr-horizontal * .. image-sg:: /auto_examples/images/sphx_glr_plot_multi_temperature_fit_001.png :alt: plot multi temperature fit :srcset: /auto_examples/images/sphx_glr_plot_multi_temperature_fit_001.png :class: sphx-glr-multi-img * .. image-sg:: /auto_examples/images/sphx_glr_plot_multi_temperature_fit_002.png :alt: plot multi temperature fit :srcset: /auto_examples/images/sphx_glr_plot_multi_temperature_fit_002.png :class: sphx-glr-multi-img .. rst-class:: sphx-glr-script-out .. code-block:: none Using database: HITEMP-CO2-TEST 'HITEMP-CO2-TEST': {'info': 'HITEMP-2010, CO2, 3 main isotope (CO2-626, 636, 628), 2283.7-2285.1 cm-1', 'path': ['/home/docs/checkouts/readthedocs.org/user_builds/radis/envs/master/lib/python3.8/site-packages/radis/test/files/cdsd_hitemp_09_fragment.txt'], 'format': 'cdsd-hitemp', 'parfuncfmt': 'hapi', 'levelsfmt': 'radis'} /home/docs/checkouts/readthedocs.org/user_builds/radis/envs/master/lib/python3.8/site-packages/radis/misc/warning.py:354: MissingReferenceWarning: Missing doi for CDSD-HITEMP. Use HITEMP-2010? Calculating Non-Equilibrium Spectrum Physical Conditions ---------------------------------------- Tgas 491 K Trot 491 K Tvib 517,517,2641 K isotope 1,2 mole_fraction 0.06573632856818698 molecule CO2 path_length 10 cm pressure_mbar 20.0 mbar rot_distribution boltzmann self_absorption True state X vib_distribution treanor wavenum_max 2284.6000 cm-1 wavenum_min 2284.2000 cm-1 Computation Parameters ---------------------------------------- Tref 296 K add_at_used broadening_method voigt cutoff 1e-25 cm-1/(#.cm-2) dbformat cdsd-hitemp dbpath /home/docs/checkouts/readthedocs.org/user_builds/radis/envs/master/lib/python3.8/site-packages/radis... folding_thresh 1e-06 include_neighbouring_lines True memory_mapping_engine auto neighbour_lines 0 cm-1 optimization simple parfuncfmt hapi parsum_mode full summation pseudo_continuum_threshold 0 sparse_ldm auto truncation 1 cm-1 waveunit cm-1 wstep 0.001 cm-1 zero_padding -1 ---------------------------------------- Calculating energy levels with Dunham expansion for CO2(X1Σu+)(iso1) (0s) 0.0% (1s) 2.6% (1s) 5.3% (2s) 7.9% (2s) 10.5% (3s) 13.2% (4s) 15.8% (4s) 18.4% (5s) 21.1% (5s) 23.7% (5s) 26.3% (5s) 28.9% (6s) 31.6% (6s) 34.2% (6s) 36.8% (7s) 39.5% (7s) 42.1% (7s) 44.7% (7s) 47.4% (7s) 50.0% (7s) 52.6% (7s) 55.3% (7s) 57.9% (7s) 60.5% (8s) 63.2% (8s) 65.8% (8s) 68.4% (8s) 71.1% (8s) 73.7% (8s) 76.3% (8s) 78.9% (8s) 81.6% (8s) 84.2% (8s) 86.8% (8s) 89.5% (8s) 92.1% (8s) 94.7% (8s) 97.4% (8s) 100.0% Database generated up to v1=37, v2=65, v3=22, J=344 Calculating energy levels with Dunham expansion for CO2(X1Σu+)(iso2) (0s) 0.0% (1s) 2.6% (1s) 5.3% (2s) 7.9% (3s) 10.5% (3s) 13.2% (4s) 15.8% (4s) 18.4% (5s) 21.1% (5s) 23.7% (5s) 26.3% (6s) 28.9% (6s) 31.6% (6s) 34.2% (6s) 36.8% (7s) 39.5% (7s) 42.1% (7s) 44.7% (7s) 47.4% (7s) 50.0% (8s) 52.6% (8s) 55.3% (8s) 57.9% (8s) 60.5% (8s) 63.2% (8s) 65.8% (8s) 68.4% (8s) 71.1% (8s) 73.7% (8s) 76.3% (8s) 78.9% (8s) 81.6% (8s) 84.2% (8s) 86.8% (8s) 89.5% (8s) 92.1% (8s) 94.7% (8s) 97.4% (8s) 100.0% Database generated up to v1=37, v2=67, v3=23, J=344 Fetching Evib & Erot. If using this code several times you should consider updating the database directly. See functions in factory.py /home/docs/checkouts/readthedocs.org/user_builds/radis/envs/master/lib/python3.8/site-packages/radis/misc/warning.py:354: NegativeEnergiesWarning: There are negative rotational energies in the database 27.57s - Spectrum calculated /home/docs/checkouts/readthedocs.org/user_builds/radis/envs/master/lib/python3.8/site-packages/radis/misc/warning.py:354: NegativeEnergiesWarning: There are negative rotational energies in the database ------------------------------ TYPICAL FIT CALCULATION TIME: Fit (in progress) profiler : spectrum_calculation 0.126s ████████████████ check_line_databank 0.001s check_non_eq_param 0.000s reinitialize 0.003s copy_database 0.001s memory_usage_warning 0.003s reset_population 0.000s calc_noneq_population_multiTvib 0.090s ███████████ part_function 0.084s ██████████ others 0.006s scaled_non_eq_linestrength 0.003s calc_emission_integral 0.004s applied_linestrength_cutoff 0.002s calc_lineshift 0.001s calc_hwhm 0.005s generate_wavenumber_arrays 0.000s calc_line_broadening 0.016s ██ precompute_LDM_lineshapes 0.002s LDM_Initialized_vectors 0.000s LDM_closest_matching_line 0.000s LDM_Distribute_lines 0.010s █ LDM_convolve 0.003s calc_other_spectral_quan 0.000s generate_spectrum_obj 0.000s others 0.006s ------------------------------ /home/docs/checkouts/readthedocs.org/user_builds/radis/envs/master/lib/python3.8/site-packages/radis/misc/warning.py:354: NegativeEnergiesWarning: There are negative rotational energies in the database /home/docs/checkouts/readthedocs.org/user_builds/radis/envs/master/lib/python3.8/site-packages/radis/misc/warning.py:354: NegativeEnergiesWarning: There are negative rotational energies in the database Now starting the fitting process: --------------------------------- T12=1150.0,T3=2650.0,Trot=1150.0, Residual: 0.0143 🏆 T12=1170.0,T3=2650.0,Trot=1150.0, Residual: 0.0144 T12=1150.0,T3=2670.0,Trot=1150.0, Residual: 0.0143 🏆 T12=1150.0,T3=2650.0,Trot=1170.0, Residual: 0.0144 T12=1150.0,T3=2650.0,Trot=1150.0, Residual: 0.0143 T12=1170.0,T3=2650.0,Trot=1150.0, Residual: 0.0144 T12=1150.0,T3=2670.0,Trot=1150.0, Residual: 0.0143 🏆 T12=1150.0,T3=2650.0,Trot=1170.0, Residual: 0.0144 T12=789.0,T3=5000.0,Trot=462.0, Residual: 0.0129 🏆 T12=809.0,T3=5000.0,Trot=462.0, Residual: 0.0130 T12=789.0,T3=4980.0,Trot=462.0, Residual: 0.0127 🏆 T12=789.0,T3=5000.0,Trot=482.0, Residual: 0.0106 🏆 T12=969.0,T3=3825.0,Trot=806.0, Residual: 0.0078 🏆 T12=989.0,T3=3825.0,Trot=806.0, Residual: 0.0080 T12=969.0,T3=3845.0,Trot=806.0, Residual: 0.0078 🏆 T12=969.0,T3=3825.0,Trot=826.0, Residual: 0.0080 T12=879.0,T3=4412.0,Trot=634.0, Residual: 0.0059 🏆 T12=899.0,T3=4412.0,Trot=634.0, Residual: 0.0061 T12=879.0,T3=4432.0,Trot=634.0, Residual: 0.0059 T12=879.0,T3=4412.0,Trot=654.0, Residual: 0.0058 🏆 T12=879.0,T3=4413.0,Trot=634.0, Residual: 0.0059 T12=899.0,T3=4413.0,Trot=634.0, Residual: 0.0061 T12=879.0,T3=4433.0,Trot=634.0, Residual: 0.0059 T12=879.0,T3=4413.0,Trot=654.0, Residual: 0.0058 🏆 T12=600.0,T3=4940.0,Trot=545.0, Residual: 0.0054 🏆 T12=620.0,T3=4940.0,Trot=545.0, Residual: 0.0055 T12=600.0,T3=4960.0,Trot=545.0, Residual: 0.0055 T12=600.0,T3=4940.0,Trot=565.0, Residual: 0.0050 🏆 T12=740.0,T3=4676.0,Trot=590.0, Residual: 0.0050 T12=760.0,T3=4676.0,Trot=590.0, Residual: 0.0052 T12=740.0,T3=4696.0,Trot=590.0, Residual: 0.0050 T12=740.0,T3=4676.0,Trot=610.0, Residual: 0.0048 🏆 T12=670.0,T3=4808.0,Trot=567.0, Residual: 0.0050 T12=690.0,T3=4808.0,Trot=567.0, Residual: 0.0051 T12=670.0,T3=4828.0,Trot=567.0, Residual: 0.0050 T12=670.0,T3=4808.0,Trot=587.0, Residual: 0.0047 🏆 T12=670.0,T3=4808.0,Trot=567.0, Residual: 0.0050 T12=690.0,T3=4808.0,Trot=567.0, Residual: 0.0051 T12=670.0,T3=4828.0,Trot=567.0, Residual: 0.0050 T12=670.0,T3=4808.0,Trot=587.0, Residual: 0.0047 🏆 T12=405.0,T3=4881.0,Trot=782.0, Residual: 0.0048 T12=425.0,T3=4881.0,Trot=782.0, Residual: 0.0046 🏆 T12=405.0,T3=4901.0,Trot=782.0, Residual: 0.0048 T12=405.0,T3=4881.0,Trot=802.0, Residual: 0.0049 T12=537.0,T3=4845.0,Trot=675.0, Residual: 0.0038 🏆 T12=557.0,T3=4845.0,Trot=675.0, Residual: 0.0038 🏆 T12=537.0,T3=4865.0,Trot=675.0, Residual: 0.0038 T12=537.0,T3=4845.0,Trot=695.0, Residual: 0.0037 🏆 T12=537.0,T3=4845.0,Trot=675.0, Residual: 0.0038 T12=557.0,T3=4845.0,Trot=675.0, Residual: 0.0038 T12=537.0,T3=4865.0,Trot=675.0, Residual: 0.0038 T12=537.0,T3=4845.0,Trot=695.0, Residual: 0.0037 🏆 T12=576.0,T3=4833.0,Trot=708.0, Residual: 0.0036 🏆 T12=596.0,T3=4833.0,Trot=708.0, Residual: 0.0036 T12=576.0,T3=4853.0,Trot=708.0, Residual: 0.0037 T12=576.0,T3=4833.0,Trot=728.0, Residual: 0.0036 T12=576.0,T3=4833.0,Trot=708.0, Residual: 0.0036 🏆 T12=596.0,T3=4833.0,Trot=708.0, Residual: 0.0036 T12=576.0,T3=4853.0,Trot=708.0, Residual: 0.0037 T12=576.0,T3=4833.0,Trot=728.0, Residual: 0.0036 T12=615.0,T3=3528.0,Trot=665.0, Residual: 0.0031 🏆 T12=635.0,T3=3528.0,Trot=665.0, Residual: 0.0032 T12=615.0,T3=3548.0,Trot=665.0, Residual: 0.0031 🏆 T12=615.0,T3=3528.0,Trot=685.0, Residual: 0.0033 T12=595.0,T3=4180.0,Trot=687.0, Residual: 0.0031 T12=615.0,T3=4180.0,Trot=687.0, Residual: 0.0031 T12=595.0,T3=4200.0,Trot=687.0, Residual: 0.0031 T12=595.0,T3=4180.0,Trot=707.0, Residual: 0.0032 T12=605.0,T3=3854.0,Trot=676.0, Residual: 0.0030 🏆 T12=625.0,T3=3854.0,Trot=676.0, Residual: 0.0030 T12=605.0,T3=3874.0,Trot=676.0, Residual: 0.0030 T12=605.0,T3=3854.0,Trot=696.0, Residual: 0.0031 T12=605.0,T3=3854.0,Trot=676.0, Residual: 0.0030 🏆 T12=625.0,T3=3854.0,Trot=676.0, Residual: 0.0030 T12=605.0,T3=3874.0,Trot=676.0, Residual: 0.0030 T12=605.0,T3=3854.0,Trot=696.0, Residual: 0.0031 T12=588.0,T3=2852.0,Trot=508.0, Residual: 0.0031 T12=608.0,T3=2852.0,Trot=508.0, Residual: 0.0034 T12=588.0,T3=2872.0,Trot=508.0, Residual: 0.0032 T12=588.0,T3=2852.0,Trot=528.0, Residual: 0.0029 🏆 T12=597.0,T3=3353.0,Trot=592.0, Residual: 0.0028 🏆 T12=617.0,T3=3353.0,Trot=592.0, Residual: 0.0029 T12=597.0,T3=3373.0,Trot=592.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=612.0, Residual: 0.0028 🏆 T12=593.0,T3=3103.0,Trot=550.0, Residual: 0.0029 T12=613.0,T3=3103.0,Trot=550.0, Residual: 0.0031 T12=593.0,T3=3123.0,Trot=550.0, Residual: 0.0029 T12=593.0,T3=3103.0,Trot=570.0, Residual: 0.0028 🏆 T12=595.0,T3=3228.0,Trot=571.0, Residual: 0.0028 T12=615.0,T3=3228.0,Trot=571.0, Residual: 0.0030 T12=595.0,T3=3248.0,Trot=571.0, Residual: 0.0029 T12=595.0,T3=3228.0,Trot=591.0, Residual: 0.0028 T12=596.0,T3=3291.0,Trot=582.0, Residual: 0.0028 T12=616.0,T3=3291.0,Trot=582.0, Residual: 0.0030 T12=596.0,T3=3311.0,Trot=582.0, Residual: 0.0028 T12=596.0,T3=3291.0,Trot=602.0, Residual: 0.0028 T12=596.0,T3=3322.0,Trot=587.0, Residual: 0.0028 T12=616.0,T3=3322.0,Trot=587.0, Residual: 0.0030 T12=596.0,T3=3342.0,Trot=587.0, Residual: 0.0028 T12=596.0,T3=3322.0,Trot=607.0, Residual: 0.0028 T12=597.0,T3=3338.0,Trot=589.0, Residual: 0.0028 T12=617.0,T3=3338.0,Trot=589.0, Residual: 0.0030 T12=597.0,T3=3358.0,Trot=589.0, Residual: 0.0028 T12=597.0,T3=3338.0,Trot=609.0, Residual: 0.0028 T12=597.0,T3=3346.0,Trot=591.0, Residual: 0.0028 T12=617.0,T3=3346.0,Trot=591.0, Residual: 0.0030 T12=597.0,T3=3366.0,Trot=591.0, Residual: 0.0028 T12=597.0,T3=3346.0,Trot=611.0, Residual: 0.0028 T12=597.0,T3=3349.0,Trot=591.0, Residual: 0.0028 T12=617.0,T3=3349.0,Trot=591.0, Residual: 0.0030 T12=597.0,T3=3369.0,Trot=591.0, Residual: 0.0028 T12=597.0,T3=3349.0,Trot=611.0, Residual: 0.0028 T12=597.0,T3=3351.0,Trot=592.0, Residual: 0.0028 T12=617.0,T3=3351.0,Trot=592.0, Residual: 0.0030 T12=597.0,T3=3371.0,Trot=592.0, Residual: 0.0028 T12=597.0,T3=3351.0,Trot=612.0, Residual: 0.0028 T12=597.0,T3=3352.0,Trot=592.0, Residual: 0.0028 T12=617.0,T3=3352.0,Trot=592.0, Residual: 0.0029 T12=597.0,T3=3372.0,Trot=592.0, Residual: 0.0028 T12=597.0,T3=3352.0,Trot=612.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=592.0, Residual: 0.0028 T12=617.0,T3=3353.0,Trot=592.0, Residual: 0.0029 T12=597.0,T3=3373.0,Trot=592.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=612.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=592.0, Residual: 0.0028 T12=617.0,T3=3353.0,Trot=592.0, Residual: 0.0029 T12=597.0,T3=3373.0,Trot=592.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=612.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=592.0, Residual: 0.0028 T12=617.0,T3=3353.0,Trot=592.0, Residual: 0.0029 T12=597.0,T3=3373.0,Trot=592.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=612.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=592.0, Residual: 0.0028 T12=617.0,T3=3353.0,Trot=592.0, Residual: 0.0029 T12=597.0,T3=3373.0,Trot=592.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=612.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=592.0, Residual: 0.0028 T12=617.0,T3=3353.0,Trot=592.0, Residual: 0.0029 T12=597.0,T3=3373.0,Trot=592.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=612.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=592.0, Residual: 0.0028 T12=617.0,T3=3353.0,Trot=592.0, Residual: 0.0029 T12=597.0,T3=3373.0,Trot=592.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=612.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=592.0, Residual: 0.0028 T12=617.0,T3=3353.0,Trot=592.0, Residual: 0.0029 T12=597.0,T3=3373.0,Trot=592.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=612.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=592.0, Residual: 0.0028 T12=617.0,T3=3353.0,Trot=592.0, Residual: 0.0029 T12=597.0,T3=3373.0,Trot=592.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=612.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=592.0, Residual: 0.0028 T12=617.0,T3=3353.0,Trot=592.0, Residual: 0.0029 T12=597.0,T3=3373.0,Trot=592.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=612.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=592.0, Residual: 0.0028 T12=617.0,T3=3353.0,Trot=592.0, Residual: 0.0029 T12=597.0,T3=3373.0,Trot=592.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=612.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=592.0, Residual: 0.0028 T12=617.0,T3=3353.0,Trot=592.0, Residual: 0.0029 T12=597.0,T3=3373.0,Trot=592.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=612.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=592.0, Residual: 0.0028 T12=617.0,T3=3353.0,Trot=592.0, Residual: 0.0029 T12=597.0,T3=3373.0,Trot=592.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=612.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=592.0, Residual: 0.0028 T12=617.0,T3=3353.0,Trot=592.0, Residual: 0.0029 T12=597.0,T3=3373.0,Trot=592.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=612.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=592.0, Residual: 0.0028 T12=617.0,T3=3353.0,Trot=592.0, Residual: 0.0029 T12=597.0,T3=3373.0,Trot=592.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=612.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=592.0, Residual: 0.0028 T12=617.0,T3=3353.0,Trot=592.0, Residual: 0.0029 T12=597.0,T3=3373.0,Trot=592.0, Residual: 0.0028 T12=597.0,T3=3353.0,Trot=612.0, Residual: 0.0028 T12=429.0,T3=1266.0,Trot=300.0, Residual: 0.0132 T12=449.0,T3=1266.0,Trot=300.0, Residual: 0.0132 T12=429.0,T3=1286.0,Trot=300.0, Residual: 0.0129 T12=429.0,T3=1266.0,Trot=320.0, Residual: 0.0134 T12=513.0,T3=2310.0,Trot=446.0, Residual: 0.0029 T12=533.0,T3=2310.0,Trot=446.0, Residual: 0.0031 T12=513.0,T3=2330.0,Trot=446.0, Residual: 0.0029 T12=513.0,T3=2310.0,Trot=466.0, Residual: 0.0030 T12=555.0,T3=2831.0,Trot=519.0, Residual: 0.0026 🏆 T12=575.0,T3=2831.0,Trot=519.0, Residual: 0.0028 T12=555.0,T3=2851.0,Trot=519.0, Residual: 0.0026 T12=555.0,T3=2831.0,Trot=539.0, Residual: 0.0025 🏆 T12=534.0,T3=2571.0,Trot=483.0, Residual: 0.0025 T12=554.0,T3=2571.0,Trot=483.0, Residual: 0.0028 T12=534.0,T3=2591.0,Trot=483.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=503.0, Residual: 0.0025 T12=523.0,T3=2440.0,Trot=464.0, Residual: 0.0027 T12=543.0,T3=2440.0,Trot=464.0, Residual: 0.0029 T12=523.0,T3=2460.0,Trot=464.0, Residual: 0.0026 T12=523.0,T3=2440.0,Trot=484.0, Residual: 0.0027 T12=529.0,T3=2505.0,Trot=473.0, Residual: 0.0026 T12=549.0,T3=2505.0,Trot=473.0, Residual: 0.0028 T12=529.0,T3=2525.0,Trot=473.0, Residual: 0.0026 T12=529.0,T3=2505.0,Trot=493.0, Residual: 0.0026 T12=531.0,T3=2538.0,Trot=478.0, Residual: 0.0026 T12=551.0,T3=2538.0,Trot=478.0, Residual: 0.0028 T12=531.0,T3=2558.0,Trot=478.0, Residual: 0.0026 T12=531.0,T3=2538.0,Trot=498.0, Residual: 0.0026 T12=532.0,T3=2554.0,Trot=480.0, Residual: 0.0025 T12=552.0,T3=2554.0,Trot=480.0, Residual: 0.0028 T12=532.0,T3=2574.0,Trot=480.0, Residual: 0.0025 T12=532.0,T3=2554.0,Trot=500.0, Residual: 0.0026 T12=533.0,T3=2562.0,Trot=481.0, Residual: 0.0025 T12=553.0,T3=2562.0,Trot=481.0, Residual: 0.0028 T12=533.0,T3=2582.0,Trot=481.0, Residual: 0.0025 T12=533.0,T3=2562.0,Trot=501.0, Residual: 0.0026 T12=533.0,T3=2567.0,Trot=482.0, Residual: 0.0025 T12=553.0,T3=2567.0,Trot=482.0, Residual: 0.0028 T12=533.0,T3=2587.0,Trot=482.0, Residual: 0.0025 T12=533.0,T3=2567.0,Trot=502.0, Residual: 0.0026 T12=534.0,T3=2569.0,Trot=482.0, Residual: 0.0025 T12=554.0,T3=2569.0,Trot=482.0, Residual: 0.0028 T12=534.0,T3=2589.0,Trot=482.0, Residual: 0.0025 T12=534.0,T3=2569.0,Trot=502.0, Residual: 0.0025 T12=534.0,T3=2570.0,Trot=482.0, Residual: 0.0025 T12=554.0,T3=2570.0,Trot=482.0, Residual: 0.0028 T12=534.0,T3=2590.0,Trot=482.0, Residual: 0.0025 T12=534.0,T3=2570.0,Trot=502.0, Residual: 0.0025 T12=534.0,T3=2570.0,Trot=482.0, Residual: 0.0025 T12=554.0,T3=2570.0,Trot=482.0, Residual: 0.0028 T12=534.0,T3=2590.0,Trot=482.0, Residual: 0.0025 T12=534.0,T3=2570.0,Trot=502.0, Residual: 0.0025 T12=534.0,T3=2570.0,Trot=483.0, Residual: 0.0025 T12=554.0,T3=2570.0,Trot=483.0, Residual: 0.0028 T12=534.0,T3=2590.0,Trot=483.0, Residual: 0.0025 T12=534.0,T3=2570.0,Trot=503.0, Residual: 0.0025 T12=534.0,T3=2570.0,Trot=483.0, Residual: 0.0025 T12=554.0,T3=2570.0,Trot=483.0, Residual: 0.0028 T12=534.0,T3=2590.0,Trot=483.0, Residual: 0.0025 T12=534.0,T3=2570.0,Trot=503.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=483.0, Residual: 0.0025 T12=554.0,T3=2571.0,Trot=483.0, Residual: 0.0028 T12=534.0,T3=2591.0,Trot=483.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=503.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=483.0, Residual: 0.0025 T12=554.0,T3=2571.0,Trot=483.0, Residual: 0.0028 T12=534.0,T3=2591.0,Trot=483.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=503.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=483.0, Residual: 0.0025 T12=554.0,T3=2571.0,Trot=483.0, Residual: 0.0028 T12=534.0,T3=2591.0,Trot=483.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=503.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=483.0, Residual: 0.0025 T12=554.0,T3=2571.0,Trot=483.0, Residual: 0.0028 T12=534.0,T3=2591.0,Trot=483.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=503.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=483.0, Residual: 0.0025 T12=554.0,T3=2571.0,Trot=483.0, Residual: 0.0028 T12=534.0,T3=2591.0,Trot=483.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=503.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=483.0, Residual: 0.0025 T12=554.0,T3=2571.0,Trot=483.0, Residual: 0.0028 T12=534.0,T3=2591.0,Trot=483.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=503.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=483.0, Residual: 0.0025 T12=554.0,T3=2571.0,Trot=483.0, Residual: 0.0028 T12=534.0,T3=2591.0,Trot=483.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=503.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=483.0, Residual: 0.0025 T12=554.0,T3=2571.0,Trot=483.0, Residual: 0.0028 T12=534.0,T3=2591.0,Trot=483.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=503.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=483.0, Residual: 0.0025 T12=554.0,T3=2571.0,Trot=483.0, Residual: 0.0028 T12=534.0,T3=2591.0,Trot=483.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=503.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=483.0, Residual: 0.0025 T12=554.0,T3=2571.0,Trot=483.0, Residual: 0.0028 T12=534.0,T3=2591.0,Trot=483.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=503.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=483.0, Residual: 0.0025 T12=554.0,T3=2571.0,Trot=483.0, Residual: 0.0028 T12=534.0,T3=2591.0,Trot=483.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=503.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=483.0, Residual: 0.0025 T12=554.0,T3=2571.0,Trot=483.0, Residual: 0.0028 T12=534.0,T3=2591.0,Trot=483.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=503.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=483.0, Residual: 0.0025 T12=554.0,T3=2571.0,Trot=483.0, Residual: 0.0028 T12=534.0,T3=2591.0,Trot=483.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=503.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=483.0, Residual: 0.0025 T12=554.0,T3=2571.0,Trot=483.0, Residual: 0.0028 T12=534.0,T3=2591.0,Trot=483.0, Residual: 0.0025 T12=534.0,T3=2571.0,Trot=503.0, Residual: 0.0025 T12=300.0,T3=2544.0,Trot=476.0, Residual: 0.0055 T12=320.0,T3=2544.0,Trot=476.0, Residual: 0.0052 T12=300.0,T3=2564.0,Trot=476.0, Residual: 0.0055 T12=300.0,T3=2544.0,Trot=496.0, Residual: 0.0056 T12=417.0,T3=2557.0,Trot=480.0, Residual: 0.0035 T12=437.0,T3=2557.0,Trot=480.0, Residual: 0.0031 T12=417.0,T3=2577.0,Trot=480.0, Residual: 0.0035 T12=417.0,T3=2557.0,Trot=500.0, Residual: 0.0036 T12=475.0,T3=2564.0,Trot=481.0, Residual: 0.0026 T12=495.0,T3=2564.0,Trot=481.0, Residual: 0.0025 🏆 T12=475.0,T3=2584.0,Trot=481.0, Residual: 0.0026 T12=475.0,T3=2564.0,Trot=501.0, Residual: 0.0027 T12=505.0,T3=2567.0,Trot=482.0, Residual: 0.0024 🏆 T12=525.0,T3=2567.0,Trot=482.0, Residual: 0.0025 T12=505.0,T3=2587.0,Trot=482.0, Residual: 0.0024 T12=505.0,T3=2567.0,Trot=502.0, Residual: 0.0025 T12=396.0,T3=2403.0,Trot=300.0, Residual: 0.0089 T12=416.0,T3=2403.0,Trot=300.0, Residual: 0.0088 T12=396.0,T3=2423.0,Trot=300.0, Residual: 0.0090 T12=396.0,T3=2403.0,Trot=320.0, Residual: 0.0080 T12=450.0,T3=2485.0,Trot=391.0, Residual: 0.0048 T12=470.0,T3=2485.0,Trot=391.0, Residual: 0.0048 T12=450.0,T3=2505.0,Trot=391.0, Residual: 0.0049 T12=450.0,T3=2485.0,Trot=411.0, Residual: 0.0040 T12=505.0,T3=2567.0,Trot=482.0, Residual: 0.0024 🏆 T12=525.0,T3=2567.0,Trot=482.0, Residual: 0.0025 T12=505.0,T3=2587.0,Trot=482.0, Residual: 0.0024 T12=505.0,T3=2567.0,Trot=502.0, Residual: 0.0025 Best ['T12', 'T3', 'Trot'] = [ 504.58098239 2567.21733341 481.80164708]['', '', ''] reached at iteration 310/324 .. rst-class:: sphx-glr-timing **Total running time of the script:** ( 1 minutes 3.959 seconds) .. _sphx_glr_download_auto_examples_plot_multi_temperature_fit.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: plot_multi_temperature_fit.py ` .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: plot_multi_temperature_fit.ipynb `