search

lanl / SEPIA

Simulation-Enabled Prediction, Inference, and Analysis: physics-informed statistical learning.
Other
29 stars 6 forks source link

Theta Constraints #18

Closed jgattiker closed 3 years ago

jgattiker commented 4 years ago

A method for the user to implement (reasonably) arbitrary constraints on thetas, e.g. theta1>(theta2)**2.

natalieklein229 commented 3 years ago

users can now pass theta_init (initial theta values following the constraint) and theta_fcon (function that evaluates the constraint on a theta array) to SepiaModel, and it will be evaluated during MCMC. A test case is in test/test_SepiaMCMC.py that confirms it does not break the code and that the sampled thetas do obey the constraint.