Function metropolis_hastings(vector<MCMCRecord>&, vector<PhaseDgm>&, double, double, double, double, vector<double>, vector<double>, double, bool, int, int)

Function Documentation

void metropolis_hastings(vector<MCMCRecord> &chain, vector<PhaseDgm> &Comp, double MassPrior, double MUncPrior, double RadPrior, double RUncPrior, vector<double> Tgap, vector<double> ave_rho, double P0, bool isothermal, int steps, int numlayers)