Multi-temperature FitΒΆ

A method to fit an experimental spectrum directly from SpectrumFactory, with fit_spectrum()

Typical output is similar to the radis-examples Multi-temperature fit :

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 radis.tools.fitting : Tvib12Tvib3Trot_NonLTEModel(). Other models can be used, such as in the one-temperature fit example

More advanced tools for interactive fitting of multi-dimensional, multi-slabs spectra can be found in fitroom. Finally, the GPU-accelerated example shows how to obtain real-time interactive spectra.

from os.path import join

from radis import Spectrum, SpectrumFactory

Get Fitted Data

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",
)
/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()`

Calculate

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
    },
)
  • plot multi temperature fit
  • plot multi temperature fit
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

Total running time of the script: ( 1 minutes 3.959 seconds)