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/PFSimulation.h // This is really -*- C++ -*- #ifndef PFSimulation_h #define PFSimulation_h #include "RunSimulation.h" #include "Distribution.h" #include "EnsembleStat.h" #include "Tile.h" #include "Model.h" #include "Integration.h" #include "PostMixturePrior.h" #include "InitialMixturePrior.h" #include "DataTree.h" #include "CountConverge.h" #include "SubsampleConverge.h" #include "MappedGrid.h" #include "KdTessellation.h" #include "QuadGrid.h" #include "CountConverge.h" #include "SubsampleConverge.h" #include "LikelihoodComputation.h" #include "BIEmpi.h" #include "TemperedSimulation.h" #include "ReversibleJump.h" #include "MetropolisHastings.h" #include "gvariable.h" #include "Serializable.h" namespace BIE { //+ CLICLASS PFSimulation SUPER RunSimulation class PFSimulation : public RunSimulation { public: static double ESSprob; static int maxTries; //+ CLICONSTR int int Simulation* Ensemble* PFSimulation(int nsteps, int npart, Simulation *sim, Ensemble *factory); ~PFSimulation(); //+ CLIMETHOD void Run void Run(); //+ CLIMETHOD void Resume void Resume(); //+ CLIMETHOD void ResumeLog void ResumeLog() { resumelog = true; } //+ CLIMETHOD void Nsample int void Nsample(int n) { nMCMC = n; } //+ CLIMETHOD void Restart void Restart() { restart = true; } //+ CLIMETHOD void SetTessTool TessToolSender* void SetTessTool(TessToolSender *tt); //+ CLIMETHOD void SampleNext void SampleNext(); //+ CLIMETHOD void SwitchOnCLI void SwitchOnCLI(); //+ CLIMETHOD void SuspendSimulation virtual void SuspendSimulation(); //+ CLIMETHOD void ResumeSimulation virtual void ResumeSimulation(); //+ CLIMETHOD void SetAutoTessTool int virtual void SetAutoTessTool(int freq); private: vector<Chain*> particles; Frontier * frontier; vector<int> dataSoFar; Simulation * sim; double curESS; int ttfreq; string prefixed_initfile; string prefixed_metfile; int nsteps; // Total number of steps unsigned npart; // Number of particles int nmix; // Two component mixture double alpha; // Dirichlet shape int ndim; // Two model parameters int Mcur; // Current number in mixture // (fixed to nmix for now) int nMCMC; // Number of MCMC steps per resample int swapstat[2]; // For swap statistics void MCMCStep(void); void UpdateStep(void); bool Resample(void); void AccumulateStats(void); Ensemble * sstat, * factory; int curLevel; void ExchangeParticles(void); void WeightNormalize(void); void IncreaseResolution(void); void PrintParticles(string filename); void WeightedStats(); bool resumelog; void resumeLog(void); Uniform * unit; // AUTO GENERATED BY ../persistence/autopersist.py protected: PFSimulation() {} 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(RunSimulation); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ESSprob); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(maxTries); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(particles); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(frontier); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(dataSoFar); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(sim); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(curESS); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ttfreq); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(prefixed_initfile); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(prefixed_metfile); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(nsteps); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(npart); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(nmix); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(alpha); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ndim); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Mcur); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(nMCMC); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(swapstat); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(sstat); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(factory); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(curLevel); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(resumelog); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(unit); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; } BIE_CLASS_TYPE_INFO(BIE::PFSimulation) BIE_CLASS_EXPORT_KEY(BIE::PFSimulation) #endif