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/TemperedSimulation.h // This is really -*- C++ -*- #ifndef TemperedSimulation_h #define TemperedSimulation_h #include "Serializable.h" #include "LikelihoodComputation.h" #include "MHWidth.h" #include "MHWidthOne.h" #include "Simulation.h" namespace BIE { class MixturePrior; //+ CLICLASS TemperedSimulation SUPER Simulation class TemperedSimulation : public Simulation { public: static int Ninter; static bool use_tempering; static double tpow; static int minmc; static int ITMAX; //+ CLICONSTR int int int double MHWidth* BaseDataTree* Model* Integration* Converge* MixturePrior* LikelihoodComputation* MCAlgorithm* TemperedSimulation(int max, int ndim, int minmc_p, double maxT, MHWidth* width, BaseDataTree* d, Model* m, Integration* i, Converge* c, MixturePrior* mp, LikelihoodComputation* l, MCAlgorithm* mca); //+ CLICONSTR int int int double MHWidth* BaseDataTree* Model* Integration* Converge* MixturePrior* LikelihoodComputation* MCAlgorithm* Simulation* TemperedSimulation(int max, int ndim, int minmc_p, double maxT, MHWidth* width, BaseDataTree* d, Model* m, Integration* i, Converge* c, MixturePrior* mp, LikelihoodComputation* l, MCAlgorithm* mca, Simulation *last); TemperedSimulation() {} virtual ~TemperedSimulation(); virtual void Reinitialize(MHWidth* width, MixturePrior *mp); //+ CLIMETHOD void SetNinter int void SetNinter(int N) { Ninter = N; } virtual void PrintStepDiagnostic(); int GetAcceptance(void) { return accnum; } int GetSwap(void) { return swapnum; } void GetLastCycle(double& u, double& d) { u=up; d=down; } virtual void PrintStateDiagnostic(); vector<double> GetMixstat(void); protected: virtual void MCMethod(); void compute_proposal_state(int k); double MaxT; MHWidth* width; vector<double> beta, factor; vector<int> mstat0, mstat1; int Ntot; int accnum, swapnum; double up, down; double upM, downM; int updownNum; Uniform * unit; private: // Save the RNG // template<class Archive> void post_save(Archive &ar, const unsigned int file_version) const { vector<Serializable*> s; s.push_back(BIEgen); ar << BOOST_SERIALIZATION_NVP(s); } // Restore the RNG // template<class Archive> void post_load(Archive &ar, const unsigned int file_version) { vector<Serializable *>s; ar >> BOOST_SERIALIZATION_NVP(s); BIEgen = dynamic_cast<BIEACG*>(s[0]); } // AUTO GENERATED BY ../persistence/autopersist.py private: friend class boost::serialization::access; BOOST_SERIALIZATION_SPLIT_MEMBER(); template<class Archive> void save(Archive & ar, const unsigned int version) const { this->pre_save(ar, version); try { ar << BOOST_SERIALIZATION_BASE_OBJECT_NVP(Simulation); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(Ninter); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(use_tempering); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(tpow); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(minmc); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(ITMAX); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(MaxT); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(width); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(beta); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(factor); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(mstat0); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(mstat1); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(Ntot); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(accnum); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(swapnum); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(up); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(down); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(upM); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(downM); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(updownNum); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(unit); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_save(ar, version); } template<class Archive> void load(Archive & ar, const unsigned int version) { this->pre_load(ar, version); try { ar >> BOOST_SERIALIZATION_BASE_OBJECT_NVP(Simulation); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(Ninter); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(use_tempering); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(tpow); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(minmc); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(ITMAX); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(MaxT); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(width); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(beta); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(factor); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(mstat0); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(mstat1); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(Ntot); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(accnum); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(swapnum); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(up); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(down); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(upM); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(downM); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(updownNum); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(unit); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_load(ar, version); } }; } BIE_CLASS_TYPE_INFO(BIE::TemperedSimulation) BIE_CLASS_EXPORT_KEY(BIE::TemperedSimulation) #endif