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/Ensemble.h // This is really -*- C++ -*- #ifndef Ensemble_h #define Ensemble_h #include <deque> #include <map> #include "BIEconfig.h" #include "ACG.h" #include "Uniform.h" #include "Normal.h" #include "DiscUnif.h" #include "Distribution.h" #include "Vector.h" #include "Serializable.h" using namespace std; namespace BIE { class StateData: public Serializable { public: double prob; State p; StateData() {}; StateData(double, State&); StateData(const StateData&); void Broadcast(); // 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(Serializable); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(prob); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(p); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; class PairIndex: public Serializable { public: PairIndex() { }; int index; double value; bool operator==(PairIndex t) { return (fabs(value - t.value)/(fabs(value) + fabs(t.value)) < 1.0e-12); } // 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(Serializable); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(index); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(value); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; class StateCache: public Serializable { private: int begin, key_pos; int count, cur_index; deque<StateData> * states; map<int, vector<PairIndex> > pairs; map<int, DiscreteUniform *> disc; public: StateCache(int nburn, deque<StateData>& states, int keypos=0); StateCache(const StateCache&); ~StateCache(); vector<double> Mean(int m); vector<double> Sample(int m); double logPDF(State &P); unsigned size() { return states->size(); } // AUTO GENERATED BY ../persistence/autopersist.py protected: StateCache() {} 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(Serializable); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(begin); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(key_pos); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(count); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(cur_index); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(states); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(pairs); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(disc); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; //+ CLICLASS Ensemble SUPER SampleDistribution [1] class Ensemble : public SampleDistribution { protected: int M; vector<int> cnts, mcnt; int minc; int maxc; int N; int count; int offset; int burnIn; int cur_index; int ibeg; int clev; vector<double> ccum; vector< pair<double, int> > fraction; deque<StateData> states; Uniform * unit; DiscreteUniform * disc; bool dimensions_set; int sampleSubspace(); void makeSampleFraction(vector<int>&); bool verbose; static int instance; int _id; public: static int keep; static double thresh; static bool continuous; static int key_pos; //+ CLICONSTR Ensemble(); //+ CLICONSTR int int Ensemble(int nmix, int ndim); //+ CLICONSTR int int int int string int Ensemble(int nmix, int ndim, int level, int nburn, string filename, int keypos); Ensemble(const Ensemble&); ~Ensemble(); //+ CLIMETHOD void Reset int int virtual void Reset(int nmix, int ndim); //+ CLIMETHOD void Reset int int int int string int virtual void Reset(int nmix, int ndim, int level, int nburn, string filename, int keypos); //+ CLIMETHOD void setDimensions int int virtual void setDimensions(int nmix, int ndim); //+ CLIMETHOD void setContinuous virtual void setContinuous(); //+ CLIMETHOD void setMaxRange virtual void setMaxRange(); //+ CLIMETHOD void setNKeep int virtual void setNKeep(int n) { keep = n; } //+ CLIMETHOD void Order int virtual void Order(int n); //+ CLIMETHOD void setVerboseOn virtual void setVerboseOn() { verbose = true; } //+ CLIMETHOD void setVerboseOff virtual void setVerboseOff() { verbose = false; } virtual bool AccumData(double, State&); virtual bool AccumData(StateData &s) { return AccumData(s.prob, s.p); } virtual bool AccumulateData(double v, RecordBuffer* datapoint) { cerr << "Not implemented for an Ensemble . . \n"; return false; } virtual void ComputeDistribution() {} virtual void ComputeDistribution(int n) { burnIn=n; ComputeDistribution(); } virtual void stats(int n, double& m, double& s); virtual int Nstates() { return states.size() + offset; } virtual double PDFSubspace(int m); virtual double logPDF(State &p) { return -INFINITY; } virtual double logPDFMarginal(int m, int n, vector<double>& V) { return -INFINITY; } virtual Ensemble* New() { return new Ensemble(M, N); } virtual vector<double> Mean(int m) { return vector<double>(N, 0.0); } virtual vector<double> StdDev(int m) { return vector<double>(N, 0.0); } virtual vector<double> StdDev(void) { return Distribution::StdDev(); } virtual vector<double> Mean(void) { return Distribution::Mean(); } virtual vector<double> Sample(int m) { return Distribution::Sample(); } virtual vector<double> Sample() { return Sample(sampleSubspace()); } virtual vector<double> Width(int m) {return vector<double>(N, 0.0); } virtual vector<double> StdDevMarginal(int m, int n) { return Distribution::StdDev(); } virtual vector<double> MeanMarginal(int m, int n) { return Distribution::Mean(); } virtual vector<double> SampleMarginal(int m, int n) { return Distribution::Sample(); } virtual vector<double> WidthMarginal(int m, int n) {return vector<double>(N, 0.0); } virtual void PrintDiag(ostream& out) {} virtual void PrintDiag(string& outfile) {} virtual int MinC() { return minc; } virtual int MaxC() { return maxc; } virtual void Broadcast(); // 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(SampleDistribution); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(M); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(cnts); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(mcnt); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(minc); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(maxc); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(N); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(count); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(offset); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(burnIn); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(cur_index); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ibeg); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(clev); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ccum); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(fraction); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(states); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(unit); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(disc); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(dimensions_set); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(verbose); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(instance); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(_id); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(keep); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(thresh); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(continuous); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(key_pos); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; } // namespace BIE BIE_CLASS_TYPE_INFO(BIE::StateData) BIE_CLASS_TYPE_INFO(BIE::PairIndex) BIE_CLASS_TYPE_INFO(BIE::StateCache) BIE_CLASS_TYPE_INFO(BIE::Ensemble) BIE_CLASS_EXPORT_KEY(BIE::StateData) BIE_CLASS_EXPORT_KEY(BIE::PairIndex) BIE_CLASS_EXPORT_KEY(BIE::StateCache) BIE_CLASS_EXPORT_KEY(BIE::Ensemble) #endif