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/TemperedDifferentialEvolution.h // This is really -*- C++ -*- #ifndef TemperedDifferentialEvolution_h #define TemperedDifferentialEvolution_h #include <LikelihoodComputation.h> #include <MHWidth.h> #include <MHWidthOne.h> #include <Simulation.h> #include <DifferentialEvolution.h> #include <Chain.h> #include "Serializable.h" namespace BIE { //+ CLICLASS TemperedDifferentialEvolution SUPER DifferentialEvolution class TemperedDifferentialEvolution : public DifferentialEvolution { private: int Nlevels, pk; unsigned tcount; unsigned jcount; public: static int Ninter; static double tpow; static int minmc; static unsigned tfreq; //+ CLICONSTR int int int string int double BaseDataTree* Model* Integration* Converge* MixturePrior* LikelihoodComputation* MCAlgorithm* //+ CLICONSTR int int int string int double BaseDataTree* Model* Integration* Converge* MixturePrior* LikelihoodComputation* MCAlgorithm* Simulation* TemperedDifferentialEvolution(int max, int ndim, int number, string deinit, int minmc_p, double maxT, BaseDataTree* d, Model* m, Integration* i, Converge* c, MixturePrior* mp, LikelihoodComputation* l, MCAlgorithm* mca, Simulation* last=0); //+ CLIMETHOD void SetTempFreq int void SetTempFreq(int n) { tfreq = n; } //+ CLIMETHOD void SetEquilSteps int void SetEquilSteps(int n) { Ninter = n; } //+ CLIMETHOD void SetTpow double void SetTpow(double x) { tpow = x; } //+ CLIMETHOD void SetMinMC int void SetMinMC(int n) { minmc = n; } virtual void PrintStepDiagnostic(); protected: vector<double> up; vector<double> down; vector<double> beta; vector<double> factor; int curswp; int curtry; int totswp; int tottry; double MaxT; vector<unsigned long> mstat0, stat0T; vector<unsigned long> mstat1, stat1T; bool tstat; virtual void MCMethod(); virtual void OneDifferentialEvolutionStep(); virtual void PowerUP(int); virtual void PowerDOWN(int); // AUTO GENERATED BY ../persistence/autopersist.py protected: TemperedDifferentialEvolution() {} 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(DifferentialEvolution); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Nlevels); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(tcount); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(jcount); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Ninter); 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(tfreq); 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(beta); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(factor); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(curswp); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(curtry); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(totswp); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(tottry); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(MaxT); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(mstat0); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(stat0T); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(mstat1); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(stat1T); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(tstat); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; } BIE_CLASS_TYPE_INFO(BIE::TemperedDifferentialEvolution) BIE_CLASS_EXPORT_KEY(BIE::TemperedDifferentialEvolution) #endif