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/MultiDistribution.h // This is really -*- C++ -*- #ifndef MultiDistribution_h #define MultiDistribution_h #include "Distribution.h" #include "cliDistribution.h" #include "clivector.h" #include "Serializable.h" namespace BIE { //+ CLICLASS MultiDistribution SUPER BinnedDistribution class MultiDistribution : public BinnedDistribution { private: RecordType * uniontype; int N; int Nbins; int Ndim; vector<BinnedDistribution*> dists; vector<int> bin_counts; vector<int> bin_offsets; vector<int> dim_counts; vector<int> dim_offsets; inline int ND(int i) { for (int indx=1; indx<N; indx++) if (i<bin_offsets[indx]) return indx-1; return N-1; } void initialize(); bool AccumulateData(double w, RecordBuffer * datapoint); public: //+ CLICONSTR int BinnedDistribution** MultiDistribution(int number, BinnedDistribution** list); //+ CLICONSTR cliDistribution* MultiDistribution(vector<BinnedDistribution*> *list); //+ CLICONSTR MultiDistribution(); ~MultiDistribution(); MultiDistribution* New(); void ComputeDistribution(); double PDF(State& z); double logPDF(State& z); vector<double> lower(); vector<double> upper(); vector<double> Mean(); vector<double> StdDev(); vector<double> Moments(int i); vector<double> Sample(); int numberData() {return bin_offsets[N];} double getValue(int i) { int j = ND(i); return dists[j]->getValue(i-bin_offsets[j]); } vector<double> getLow(int i) { int j = ND(i); return dists[j]->getLow(i-bin_offsets[j]); } vector<double> getHigh(int i) { int j = ND(i); return dists[j]->getHigh(i-bin_offsets[j]); } int getdim(int i) { int j = ND(i); return dists[j]->getdim(i-bin_offsets[j]); } int Ndist(void) {return N;} // 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(BinnedDistribution); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(uniontype); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(N); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Nbins); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Ndim); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(dists); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(bin_counts); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(bin_offsets); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(dim_counts); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(dim_offsets); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; } // namespace BIE BIE_CLASS_TYPE_INFO(BIE::MultiDistribution) BIE_CLASS_EXPORT_KEY(BIE::MultiDistribution) #endif