.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "auto_examples/3_Fitting/plot5_legacyFit_Trot-Tvib.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note You can download :ref:`below ` the full example code and run it with 🔬 `Radis-Lab `__, .. rst-class:: sphx-glr-example-title .. _sphx_glr_auto_examples_3_Fitting_plot5_legacyFit_Trot-Tvib.py: .. _example_multi_temperature_fit: ===================== Legacy #2: non-equilibrium CO2 (Tvib_12, Tvib_3, Trot) ===================== A `3 temperature fitting example `__ . reproducing the validation case of Klarenaar 2017 [1]_, who calculated a transmittance spectrum from the initial data of Dang 1973 [2]_, with a 1 rotational temperature + 3 vibrational temperature (Treanor distributions) model. |CO2| Energies are calculated from Dunham developments in an uncoupled harmonic oscillator - rigid rotor model. The example is based on one of `RADIS validation cases `_. It makes use of the RADIS `Spectrum `_ class and the associated compare and load functions 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. .. [1] Klarenaar et al 2017, "Time evolution of vibrational temperatures in a CO2 glow discharge measured with infrared absorption spectroscopy" doi/10.1088/1361-6595/aa902e .. [2] Dang et al 1982, "Detailed vibrational population distributions in a CO2 laser discharge as measured with a tunable diode laser" doi/10.1007/BF00694640 .. GENERATED FROM PYTHON SOURCE LINES 40-45 .. code-block:: Python from os.path import join from radis import Spectrum, SpectrumFactory .. GENERATED FROM PYTHON SOURCE LINES 46-47 Get Fitted Data .. GENERATED FROM PYTHON SOURCE LINES 47-63 .. code-block:: Python 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/checkouts/latest/radis/spectrum/spectrum.py:4921: 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 64-65 Calculate .. GENERATED FROM PYTHON SOURCE LINES 65-102 .. code-block:: Python 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_legacy( 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/3_Fitting/images/sphx_glr_plot5_legacyFit_Trot-Tvib_001.png :alt: plot5 legacyFit Trot Tvib :srcset: /auto_examples/3_Fitting/images/sphx_glr_plot5_legacyFit_Trot-Tvib_001.png :class: sphx-glr-multi-img * .. image-sg:: /auto_examples/3_Fitting/images/sphx_glr_plot5_legacyFit_Trot-Tvib_002.png :alt: plot5 legacyFit Trot Tvib :srcset: /auto_examples/3_Fitting/images/sphx_glr_plot5_legacyFit_Trot-Tvib_002.png :class: sphx-glr-multi-img .. rst-class:: sphx-glr-script-out .. code-block:: none /home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/latest/radis/db/molparam.py:252: FutureWarning: The 'delim_whitespace' keyword in pd.read_csv is deprecated and will be removed in a future version. Use ``sep='\s+'`` instead 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/checkouts/latest/radis/test/files/cdsd_hitemp_09_fragment.txt'], 'format': 'cdsd-hitemp', 'parfuncfmt': 'hapi', 'levelsfmt': 'radis'} /home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/latest/radis/misc/warning.py:427: 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 0.02 bar 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/checkouts/latest/radis/test/files/cdsd_hitemp... diluent air 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% (0s) 2.6% (1s) 5.3% (2s) 7.9% (2s) 10.5% (3s) 13.2% (3s) 15.8% (3s) 18.4% (4s) 21.1% (4s) 23.7% (4s) 26.3% (5s) 28.9% (5s) 31.6% (5s) 34.2% (5s) 36.8% (6s) 39.5% (6s) 42.1% (6s) 44.7% (6s) 47.4% (6s) 50.0% (6s) 52.6% (6s) 55.3% (6s) 57.9% (6s) 60.5% (6s) 63.2% (6s) 65.8% (6s) 68.4% (6s) 71.1% (6s) 73.7% (7s) 76.3% (7s) 78.9% (7s) 81.6% (7s) 84.2% (7s) 86.8% (7s) 89.5% (7s) 92.1% (7s) 94.7% (7s) 97.4% (7s) 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% (2s) 10.5% (3s) 13.2% (3s) 15.8% (4s) 18.4% (4s) 21.1% (4s) 23.7% (5s) 26.3% (5s) 28.9% (5s) 31.6% (5s) 34.2% (5s) 36.8% (5s) 39.5% (6s) 42.1% (6s) 44.7% (6s) 47.4% (6s) 50.0% (6s) 52.6% (6s) 55.3% (6s) 57.9% (7s) 60.5% (7s) 63.2% (7s) 65.8% (7s) 68.4% (7s) 71.1% (7s) 73.7% (7s) 76.3% (7s) 78.9% (7s) 81.6% (7s) 84.2% (7s) 86.8% (7s) 89.5% (7s) 92.1% (7s) 94.7% (7s) 97.4% (7s) 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/checkouts/latest/radis/misc/warning.py:427: NegativeEnergiesWarning: There are negative rotational energies in the database 22.78s - Spectrum calculated /home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/latest/radis/misc/warning.py:427: NegativeEnergiesWarning: There are negative rotational energies in the database ------------------------------ TYPICAL FIT CALCULATION TIME: Fit (in progress) profiler : spectrum_calculation 0.106s ████████████████ check_line_databank 0.002s check_non_eq_param 0.000s reinitialize 0.005s copy_database 0.001s memory_usage_warning 0.003s reset_population 0.000s calc_noneq_population_multiTvib 0.071s ██████████ part_function 0.066s █████████ others 0.006s scaled_non_eq_linestrength 0.003s calc_emission_integral 0.003s applied_linestrength_cutoff 0.001s calc_lineshift 0.001s calc_hwhm 0.005s generate_wavenumber_arrays 0.000s calc_line_broadening 0.014s ██ precompute_LDM_lineshapes 0.001s LDM_Initialized_vectors 0.000s LDM_closest_matching_line 0.000s LDM_Distribute_lines 0.009s █ LDM_convolve 0.003s calc_other_spectral_quan 0.000s generate_spectrum_obj 0.000s others 0.007s ------------------------------ /home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/latest/radis/misc/warning.py:427: NegativeEnergiesWarning: There are negative rotational energies in the database /home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/latest/radis/misc/warning.py:427: NegativeEnergiesWarning: There are negative rotational energies in the database Now starting the fitting process: --------------------------------- /home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/latest/radis/tools/fitting.py:486: OptimizeWarning: Unknown solver options: maxiter 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=782.0,T3=5000.0,Trot=459.0, Residual: 0.0133 🏆 T12=802.0,T3=5000.0,Trot=459.0, Residual: 0.0134 T12=782.0,T3=4980.0,Trot=459.0, Residual: 0.0131 🏆 T12=782.0,T3=5000.0,Trot=479.0, Residual: 0.0110 🏆 T12=966.0,T3=3825.0,Trot=804.0, Residual: 0.0077 🏆 T12=986.0,T3=3825.0,Trot=804.0, Residual: 0.0080 T12=966.0,T3=3845.0,Trot=804.0, Residual: 0.0077 🏆 T12=966.0,T3=3825.0,Trot=824.0, Residual: 0.0079 T12=874.0,T3=4412.0,Trot=632.0, Residual: 0.0059 🏆 T12=894.0,T3=4412.0,Trot=632.0, Residual: 0.0061 T12=874.0,T3=4432.0,Trot=632.0, Residual: 0.0059 T12=874.0,T3=4412.0,Trot=652.0, Residual: 0.0058 🏆 T12=874.0,T3=4413.0,Trot=632.0, Residual: 0.0059 T12=894.0,T3=4413.0,Trot=632.0, Residual: 0.0061 T12=874.0,T3=4433.0,Trot=632.0, Residual: 0.0059 T12=874.0,T3=4413.0,Trot=652.0, Residual: 0.0058 🏆 T12=576.0,T3=4971.0,Trot=546.0, Residual: 0.0056 🏆 T12=596.0,T3=4971.0,Trot=546.0, Residual: 0.0056 T12=576.0,T3=4991.0,Trot=546.0, Residual: 0.0056 T12=576.0,T3=4971.0,Trot=566.0, Residual: 0.0051 🏆 T12=725.0,T3=4692.0,Trot=589.0, Residual: 0.0050 🏆 T12=745.0,T3=4692.0,Trot=589.0, Residual: 0.0051 T12=725.0,T3=4712.0,Trot=589.0, Residual: 0.0050 T12=725.0,T3=4692.0,Trot=609.0, Residual: 0.0048 🏆 T12=650.0,T3=4832.0,Trot=567.0, Residual: 0.0050 T12=670.0,T3=4832.0,Trot=567.0, Residual: 0.0051 T12=650.0,T3=4852.0,Trot=567.0, Residual: 0.0050 T12=650.0,T3=4832.0,Trot=587.0, Residual: 0.0047 🏆 T12=688.0,T3=4762.0,Trot=578.0, Residual: 0.0049 T12=708.0,T3=4762.0,Trot=578.0, Residual: 0.0051 T12=688.0,T3=4782.0,Trot=578.0, Residual: 0.0050 T12=688.0,T3=4762.0,Trot=598.0, Residual: 0.0047 🏆 T12=688.0,T3=4762.0,Trot=578.0, Residual: 0.0049 T12=708.0,T3=4762.0,Trot=578.0, Residual: 0.0051 T12=688.0,T3=4782.0,Trot=578.0, Residual: 0.0050 T12=688.0,T3=4762.0,Trot=598.0, Residual: 0.0047 🏆 T12=315.0,T3=5000.0,Trot=850.0, Residual: 0.0059 T12=335.0,T3=5000.0,Trot=850.0, Residual: 0.0057 T12=315.0,T3=4980.0,Trot=850.0, Residual: 0.0059 T12=315.0,T3=5000.0,Trot=870.0, Residual: 0.0060 T12=501.0,T3=4881.0,Trot=714.0, Residual: 0.0039 🏆 T12=521.0,T3=4881.0,Trot=714.0, Residual: 0.0038 🏆 T12=501.0,T3=4901.0,Trot=714.0, Residual: 0.0039 T12=501.0,T3=4881.0,Trot=734.0, Residual: 0.0039 T12=594.0,T3=4821.0,Trot=646.0, Residual: 0.0039 T12=614.0,T3=4821.0,Trot=646.0, Residual: 0.0040 T12=594.0,T3=4841.0,Trot=646.0, Residual: 0.0040 T12=594.0,T3=4821.0,Trot=666.0, Residual: 0.0038 🏆 T12=548.0,T3=4851.0,Trot=680.0, Residual: 0.0038 🏆 T12=568.0,T3=4851.0,Trot=680.0, Residual: 0.0038 🏆 T12=548.0,T3=4871.0,Trot=680.0, Residual: 0.0038 T12=548.0,T3=4851.0,Trot=700.0, Residual: 0.0037 🏆 T12=548.0,T3=4851.0,Trot=680.0, Residual: 0.0038 T12=568.0,T3=4851.0,Trot=680.0, Residual: 0.0038 T12=548.0,T3=4871.0,Trot=680.0, Residual: 0.0038 T12=548.0,T3=4851.0,Trot=700.0, Residual: 0.0037 🏆 T12=584.0,T3=4793.0,Trot=714.0, Residual: 0.0036 🏆 T12=604.0,T3=4793.0,Trot=714.0, Residual: 0.0037 T12=584.0,T3=4813.0,Trot=714.0, Residual: 0.0036 T12=584.0,T3=4793.0,Trot=734.0, Residual: 0.0036 T12=584.0,T3=4793.0,Trot=714.0, Residual: 0.0036 🏆 T12=604.0,T3=4793.0,Trot=714.0, Residual: 0.0037 T12=584.0,T3=4813.0,Trot=714.0, Residual: 0.0036 T12=584.0,T3=4793.0,Trot=734.0, Residual: 0.0036 T12=584.0,T3=3812.0,Trot=702.0, Residual: 0.0032 🏆 T12=604.0,T3=3812.0,Trot=702.0, Residual: 0.0032 🏆 T12=584.0,T3=3832.0,Trot=702.0, Residual: 0.0032 T12=584.0,T3=3812.0,Trot=722.0, Residual: 0.0034 T12=584.0,T3=4302.0,Trot=708.0, Residual: 0.0032 T12=604.0,T3=4302.0,Trot=708.0, Residual: 0.0032 T12=584.0,T3=4322.0,Trot=708.0, Residual: 0.0033 T12=584.0,T3=4302.0,Trot=728.0, Residual: 0.0033 T12=584.0,T3=4057.0,Trot=705.0, Residual: 0.0032 🏆 T12=604.0,T3=4057.0,Trot=705.0, Residual: 0.0032 🏆 T12=584.0,T3=4077.0,Trot=705.0, Residual: 0.0032 T12=584.0,T3=4057.0,Trot=725.0, Residual: 0.0033 T12=584.0,T3=4057.0,Trot=705.0, Residual: 0.0032 T12=604.0,T3=4057.0,Trot=705.0, Residual: 0.0032 🏆 T12=584.0,T3=4077.0,Trot=705.0, Residual: 0.0032 T12=584.0,T3=4057.0,Trot=725.0, Residual: 0.0033 T12=599.0,T3=2565.0,Trot=470.0, Residual: 0.0038 T12=619.0,T3=2565.0,Trot=470.0, Residual: 0.0042 T12=599.0,T3=2585.0,Trot=470.0, Residual: 0.0038 T12=599.0,T3=2565.0,Trot=490.0, Residual: 0.0035 T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028 🏆 T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029 T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028 🏆 T12=595.0,T3=2938.0,Trot=529.0, Residual: 0.0030 T12=615.0,T3=2938.0,Trot=529.0, Residual: 0.0033 T12=595.0,T3=2958.0,Trot=529.0, Residual: 0.0030 T12=595.0,T3=2938.0,Trot=549.0, Residual: 0.0029 T12=593.0,T3=3124.0,Trot=558.0, Residual: 0.0029 T12=613.0,T3=3124.0,Trot=558.0, Residual: 0.0031 T12=593.0,T3=3144.0,Trot=558.0, Residual: 0.0029 T12=593.0,T3=3124.0,Trot=578.0, Residual: 0.0028 🏆 T12=592.0,T3=3218.0,Trot=573.0, Residual: 0.0028 T12=612.0,T3=3218.0,Trot=573.0, Residual: 0.0030 T12=592.0,T3=3238.0,Trot=573.0, Residual: 0.0028 T12=592.0,T3=3218.0,Trot=593.0, Residual: 0.0028 🏆 T12=592.0,T3=3264.0,Trot=580.0, Residual: 0.0028 T12=612.0,T3=3264.0,Trot=580.0, Residual: 0.0030 T12=592.0,T3=3284.0,Trot=580.0, Residual: 0.0028 T12=592.0,T3=3264.0,Trot=600.0, Residual: 0.0028 T12=592.0,T3=3288.0,Trot=584.0, Residual: 0.0028 T12=612.0,T3=3288.0,Trot=584.0, Residual: 0.0029 T12=592.0,T3=3308.0,Trot=584.0, Residual: 0.0028 T12=592.0,T3=3288.0,Trot=604.0, Residual: 0.0028 T12=592.0,T3=3299.0,Trot=586.0, Residual: 0.0028 T12=612.0,T3=3299.0,Trot=586.0, Residual: 0.0029 T12=592.0,T3=3319.0,Trot=586.0, Residual: 0.0028 T12=592.0,T3=3299.0,Trot=606.0, Residual: 0.0028 T12=592.0,T3=3305.0,Trot=587.0, Residual: 0.0028 T12=612.0,T3=3305.0,Trot=587.0, Residual: 0.0029 T12=592.0,T3=3325.0,Trot=587.0, Residual: 0.0028 T12=592.0,T3=3305.0,Trot=607.0, Residual: 0.0028 T12=591.0,T3=3308.0,Trot=587.0, Residual: 0.0028 T12=611.0,T3=3308.0,Trot=587.0, Residual: 0.0029 T12=591.0,T3=3328.0,Trot=587.0, Residual: 0.0028 T12=591.0,T3=3308.0,Trot=607.0, Residual: 0.0028 T12=591.0,T3=3309.0,Trot=587.0, Residual: 0.0028 T12=611.0,T3=3309.0,Trot=587.0, Residual: 0.0029 T12=591.0,T3=3329.0,Trot=587.0, Residual: 0.0028 T12=591.0,T3=3309.0,Trot=607.0, Residual: 0.0028 T12=591.0,T3=3310.0,Trot=587.0, Residual: 0.0028 T12=611.0,T3=3310.0,Trot=587.0, Residual: 0.0029 T12=591.0,T3=3330.0,Trot=587.0, Residual: 0.0028 T12=591.0,T3=3310.0,Trot=607.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028 T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029 T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028 T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029 T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028 T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029 T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028 T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029 T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028 T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029 T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028 T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029 T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028 T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029 T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028 T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029 T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028 T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029 T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028 T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029 T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028 T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029 T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028 T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029 T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028 T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029 T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028 T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029 T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028 T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028 T12=434.0,T3=1587.0,Trot=300.0, Residual: 0.0083 T12=454.0,T3=1587.0,Trot=300.0, Residual: 0.0083 T12=434.0,T3=1607.0,Trot=300.0, Residual: 0.0080 T12=434.0,T3=1587.0,Trot=320.0, Residual: 0.0084 T12=513.0,T3=2449.0,Trot=444.0, Residual: 0.0029 T12=533.0,T3=2449.0,Trot=444.0, Residual: 0.0031 T12=513.0,T3=2469.0,Trot=444.0, Residual: 0.0030 T12=513.0,T3=2449.0,Trot=464.0, Residual: 0.0026 🏆 T12=552.0,T3=2880.0,Trot=516.0, Residual: 0.0027 T12=572.0,T3=2880.0,Trot=516.0, Residual: 0.0029 T12=552.0,T3=2900.0,Trot=516.0, Residual: 0.0027 T12=552.0,T3=2880.0,Trot=536.0, Residual: 0.0025 🏆 T12=532.0,T3=2664.0,Trot=480.0, Residual: 0.0027 T12=552.0,T3=2664.0,Trot=480.0, Residual: 0.0029 T12=532.0,T3=2684.0,Trot=480.0, Residual: 0.0028 T12=532.0,T3=2664.0,Trot=500.0, Residual: 0.0025 🏆 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=537.0,T3=2718.0,Trot=489.0, Residual: 0.0027 T12=557.0,T3=2718.0,Trot=489.0, Residual: 0.0029 T12=537.0,T3=2738.0,Trot=489.0, Residual: 0.0027 T12=537.0,T3=2718.0,Trot=509.0, Residual: 0.0025 T12=540.0,T3=2745.0,Trot=493.0, Residual: 0.0027 T12=560.0,T3=2745.0,Trot=493.0, Residual: 0.0029 T12=540.0,T3=2765.0,Trot=493.0, Residual: 0.0027 T12=540.0,T3=2745.0,Trot=513.0, Residual: 0.0025 T12=541.0,T3=2759.0,Trot=495.0, Residual: 0.0027 T12=561.0,T3=2759.0,Trot=495.0, Residual: 0.0029 T12=541.0,T3=2779.0,Trot=495.0, Residual: 0.0027 T12=541.0,T3=2759.0,Trot=515.0, Residual: 0.0025 T12=542.0,T3=2765.0,Trot=496.0, Residual: 0.0027 T12=562.0,T3=2765.0,Trot=496.0, Residual: 0.0029 T12=542.0,T3=2785.0,Trot=496.0, Residual: 0.0027 T12=542.0,T3=2765.0,Trot=516.0, Residual: 0.0025 T12=542.0,T3=2769.0,Trot=497.0, Residual: 0.0027 T12=562.0,T3=2769.0,Trot=497.0, Residual: 0.0029 T12=542.0,T3=2789.0,Trot=497.0, Residual: 0.0027 T12=542.0,T3=2769.0,Trot=517.0, Residual: 0.0025 T12=542.0,T3=2771.0,Trot=497.0, Residual: 0.0027 T12=562.0,T3=2771.0,Trot=497.0, Residual: 0.0029 T12=542.0,T3=2791.0,Trot=497.0, Residual: 0.0027 T12=542.0,T3=2771.0,Trot=517.0, Residual: 0.0025 T12=542.0,T3=2771.0,Trot=497.0, Residual: 0.0027 T12=562.0,T3=2771.0,Trot=497.0, Residual: 0.0029 T12=542.0,T3=2791.0,Trot=497.0, Residual: 0.0027 T12=542.0,T3=2771.0,Trot=517.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=490.0,T3=2299.0,Trot=407.0, Residual: 0.0035 T12=510.0,T3=2299.0,Trot=407.0, Residual: 0.0037 T12=490.0,T3=2319.0,Trot=407.0, Residual: 0.0036 T12=490.0,T3=2299.0,Trot=427.0, Residual: 0.0031 T12=516.0,T3=2536.0,Trot=452.0, Residual: 0.0029 T12=536.0,T3=2536.0,Trot=452.0, Residual: 0.0031 T12=516.0,T3=2556.0,Trot=452.0, Residual: 0.0030 T12=516.0,T3=2536.0,Trot=472.0, Residual: 0.0025 T12=529.0,T3=2654.0,Trot=475.0, Residual: 0.0028 T12=549.0,T3=2654.0,Trot=475.0, Residual: 0.0029 T12=529.0,T3=2674.0,Trot=475.0, Residual: 0.0028 T12=529.0,T3=2654.0,Trot=495.0, Residual: 0.0025 T12=536.0,T3=2713.0,Trot=486.0, Residual: 0.0027 T12=556.0,T3=2713.0,Trot=486.0, Residual: 0.0029 T12=536.0,T3=2733.0,Trot=486.0, Residual: 0.0028 T12=536.0,T3=2713.0,Trot=506.0, Residual: 0.0025 T12=539.0,T3=2743.0,Trot=492.0, Residual: 0.0027 T12=559.0,T3=2743.0,Trot=492.0, Residual: 0.0029 T12=539.0,T3=2763.0,Trot=492.0, Residual: 0.0027 T12=539.0,T3=2743.0,Trot=512.0, Residual: 0.0025 T12=541.0,T3=2757.0,Trot=495.0, Residual: 0.0027 T12=561.0,T3=2757.0,Trot=495.0, Residual: 0.0029 T12=541.0,T3=2777.0,Trot=495.0, Residual: 0.0027 T12=541.0,T3=2757.0,Trot=515.0, Residual: 0.0025 T12=541.0,T3=2765.0,Trot=496.0, Residual: 0.0027 T12=561.0,T3=2765.0,Trot=496.0, Residual: 0.0029 T12=541.0,T3=2785.0,Trot=496.0, Residual: 0.0027 T12=541.0,T3=2765.0,Trot=516.0, Residual: 0.0025 T12=542.0,T3=2768.0,Trot=497.0, Residual: 0.0027 T12=562.0,T3=2768.0,Trot=497.0, Residual: 0.0029 T12=542.0,T3=2788.0,Trot=497.0, Residual: 0.0027 T12=542.0,T3=2768.0,Trot=517.0, Residual: 0.0025 T12=542.0,T3=2770.0,Trot=497.0, Residual: 0.0027 T12=562.0,T3=2770.0,Trot=497.0, Residual: 0.0029 T12=542.0,T3=2790.0,Trot=497.0, Residual: 0.0027 T12=542.0,T3=2770.0,Trot=517.0, Residual: 0.0025 T12=542.0,T3=2771.0,Trot=497.0, Residual: 0.0027 T12=562.0,T3=2771.0,Trot=497.0, Residual: 0.0029 T12=542.0,T3=2791.0,Trot=497.0, Residual: 0.0027 T12=542.0,T3=2771.0,Trot=517.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027 T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029 T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027 T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025 Best ['T12', 'T3', 'Trot'] = [ 532.40831841 2664.44122594 499.63988523]['', '', ''] reached at iteration 209/404 .. rst-class:: sphx-glr-timing **Total running time of the script:** (1 minutes 7.698 seconds) .. _sphx_glr_download_auto_examples_3_Fitting_plot5_legacyFit_Trot-Tvib.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: plot5_legacyFit_Trot-Tvib.ipynb ` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: plot5_legacyFit_Trot-Tvib.py `