search

olivares-j / Kalkayotl

This repository contains a code to infer distances to stellar clusters and its stars from Gaia parallaxes.
GNU General Public License v3.0
11 stars 8 forks source link

initialization #25

Open olivares-j opened 1 year ago

olivares-j commented 1 year ago

Initialization is fundamental for the 6D version due to radial velocity uncertainties and possible outliers. Unfortunately in the init_nuts function of pymc convergence tolerance is hard-coded. If the sampler fails with the derivative of the RV is zero it meas the initialization is not good enough and we will need to modify those tolerances.

Solutions:

  1. Modify the pymc routine and make a pull request
  2. Work with your own pymc fork.
olivares-j commented 1 year ago

To use the initialization values with numpyro or blackjax try to do the converstion of the LKJ packed by yourself.