Temperature Issues in MCMC Calculation

[mzb123]

In my model, I have fixed the temperature of the primary star at 6125 K, while iterating over the temperature ratio parameter teffratio .

I have successfully run the MCMC calculation, but I encountered issues when attempting to extract the temperature of the secondary star. When running the following code, the secondary star temperature and teffratio obtained are comparable. This presents itself as an error. Could you provide guidance on resolving this issue?

The attachment includes my bundle file. OO Leo-LC9-fixed-spot on primary.zip

b.uncertainties_from_distribution_collection(solution='emcee_sol', parameters=['teff@secondary'], tex=True)

[mzb123]

b.get_parameter('teff', component='secondary', context='component').get_distribution('emcee_sol')
Out:<distl.composite {teffratio@binary} * <distl.delta loc=1.0611573111573112 unit=> unit=>
Does this mean just doing the panning？

[kecnry]

Thanks for reporting this and attaching the file - I do think this appears to be a bug, I'm guessing with the way the units are handled when propagating the constraints, but that will take a little more investigation before we can fix the bug. In the meantime, you can get the expected result by sampling manually from the teffratio distribution and multiplying by the (fixed) primary temperature to get a distribution for the secondary temperature.

import phoebe

b = phoebe.open('OO Leo-LC9-fixed-spot on primary.bundle')
b.plot_distribution_collection('emcee_sol', parameters=['teffratio', 'teff@secondary', 'teff@primary'], show=True)

>>> print(b.filter(qualifier='teff'))

ParameterSet: 3 parameters
           teff@primary@component: 6125.0 K
C        teff@secondary@component: 6390.493770110515 K
                  teff@constraint: {teffratio@binary@component} * {teff@primary@component}