BIE Overview Quick Start Theoretical overview Available methods Parallel chains Computation engine Command Line Interface Graphical User Interface Assigning Output Visualization tool Data handling Software technology Parallel debugging Results to date Recent developments Links User guide Future features Project goals Project Team Copyright and license

/home/weinberg/src/BIE/include/PostMixturePrior.h // This is really -*- C++ -*- #ifndef PostMixturePrior_h #define PostMixturePrior_h #include <BIEconfig.h> #include <Uniform.h> #include <Normal.h> #include <MixturePrior.h> #include <Dirichlet.h> #include <Ensemble.h> #include "Serializable.h" namespace BIE { //+ CLICLASS PostMixturePrior SUPER MixturePrior class PostMixturePrior : public MixturePrior { private: Ensemble* _dist; bool _manage_dist; vector<double> _lower; vector<double> _upper; bool _clamp; bool _sample; public: static double width_factor; static int max_iter; //+ CLICONSTR PostMixturePrior() {}; //+ CLICONSTR MixturePrior* Ensemble* PostMixturePrior(MixturePrior* old, Ensemble* dist); //+ CLICONSTR MixturePrior* Ensemble* int int int int string int PostMixturePrior(MixturePrior* old, Ensemble* dist, int nmix, int ndim, int level, int nburn, string filename, int keypos); ~PostMixturePrior(); PostMixturePrior* New(); //+ CLIMETHOD void FreeRange void FreeRange() { _clamp = false; } //+ CLIMETHOD void useSamples void useSamples() { _sample = true; } //+ CLIMETHOD void useWidths void useWidths() { _sample = false; } double PDF(State& x) { return _dist->PDF(x); } double logPDF(State& x) { return _dist->logPDF(x); } double logPDFMarginal(int M, int n, vector<double>& x) { return _dist->logPDFMarginal(M, n, x); } vector<double> lower(void) { return _lower; } vector<double> upper(void) { return _upper; } vector<double> Mean() { return _dist->Mean(); } vector<double> StdDev() { return _dist->StdDev(); } vector<double> Moments(int m) { return _dist->Moments(m); } vector<double> Sample() { return vector<double>(Ndim*Nmix, 0.0); } vector<double> Sample(int m) { return _dist->Sample(m); } virtual void SamplePrior(int& m, vector<double>& w, vector< vector<double> >& p); virtual void SamplePrior(int M, int n, vector<double>& V); /* Based on prior type, return random variates for mixture model The @param width is ignored for this sample method; the sample comes from the covariance matrix directly scaled by the @param width_factor */ virtual void SampleProposal(Chain& ch, MHWidth* width); virtual int MinC() { return _dist->MinC(); } virtual int MaxC() { return _dist->MaxC(); } // AUTO GENERATED BY ../persistence/autopersist.py private: friend class boost::serialization::access; template<class Archive> void serialize(Archive & ar, const unsigned int version) { this->pre_serialize(ar, version); try { ar & BOOST_SERIALIZATION_BASE_OBJECT_NVP(MixturePrior); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(_dist); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(_manage_dist); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(_lower); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(_upper); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(_clamp); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(_sample); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; } BIE_CLASS_TYPE_INFO(BIE::PostMixturePrior) BIE_CLASS_EXPORT_KEY(BIE::PostMixturePrior) #endif