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/EnsembleStat.h // This is really -*- C++ -*- #ifndef EnsembleStat_h #define EnsembleStat_h #include "Ensemble.h" #include "Serializable.h" using namespace std; namespace BIE { //+ CLICLASS EnsembleStat SUPER Ensemble [1] class EnsembleStat : public Ensemble { protected: vector<int> icut; vector< vector<int> > mcut; vector<int> mcum; bool var_computed; vector< vector<double> > vmean; vector< vector< vector<double> > > mmean; vector< vector<double> > vvar; vector< vector< vector<double> > > mvar; Normal * normal; double factor; public: static int minsub; //+ CLICONSTR EnsembleStat(); //+ CLICONSTR int int EnsembleStat(int nmix, int ndim); //+ CLICONSTR int int int int string int EnsembleStat(int nmix, int ndim, int level, int nburn, string filename, int keypos); EnsembleStat(const EnsembleStat&); ~EnsembleStat(); //+ CLIMETHOD void Reset int int void Reset(int nmix, int ndim); //+ CLIMETHOD void Reset int int int int string int void Reset(int nmix, int ndim, int level, int nburn, string filename, int keypos); //+ CLIMETHOD void setDimensions int int void setDimensions(int nmix, int ndim); vector<Matrix> covar; vector< vector<Matrix> > covarM; vector<Matrix> Teigen; vector< vector<Matrix> > TeigeM; vector<Matrix> eigen; vector< vector<Matrix> > eigenM; vector<Vector> mean; vector< vector<Vector> > meanM; vector<Vector> evals; vector< vector<Vector> > evalsM; vector<Vector> Tmean; vector< vector<Vector> > TmeanM; //+ CLIMETHOD void ComputeDistribution //+ CLIMETHOD void ComputeDistribution int void ComputeDistribution(); void ComputeDistribution(int n) { burnIn=n; ComputeDistribution(); } void ComputeDistributionMarginal(ostream& out); int trim_location(int m) { return icut[m]; } int trim_location_marginal(int m, int n) { return mcut[m][n]; } EnsembleStat* New() { return new EnsembleStat(M, N); } double PDF(State&); double logPDF(State&); double logPDFMarginal(int m, int n, vector<double>& V); double PDF(int, State&); vector<double> lower(void); vector<double> upper(void); vector<double> Mean(int m); vector<double> StdDev(int m); vector<double> Moments(int m, int k); vector<double> Width(int m) { return StdDev(m); } vector<double> Sample(int m); double SampleOne(int m, int j); vector<double> StdDevMarginal(int m, int n); vector<double> MeanMarginal(int m, int n); vector<double> SampleMarginal(int m, int n); vector<double> WidthMarginal(int m, int n) { return StdDevMarginal(m, n); } void PrintDiag(ostream& out); void PrintDiag(string& outfile); // 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(Ensemble); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(icut); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(mcut); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(mcum); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(var_computed); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(vmean); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(mmean); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(vvar); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(mvar); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(normal); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(factor); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(minsub); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(covar); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(covarM); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Teigen); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(TeigeM); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(eigen); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(eigenM); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(mean); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(meanM); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(evals); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(evalsM); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Tmean); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(TmeanM); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; } // namespace BIE BIE_CLASS_TYPE_INFO(BIE::EnsembleStat) BIE_CLASS_EXPORT_KEY(BIE::EnsembleStat) #endif