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/PopModel/HistogramCache.h // This is really -*- C++ -*- #include "Serializable.h" #include <algorithm> #include <string> #include <vector> using namespace std; #include "gfunction.h" #include "hash_map.h" class HistoKey: public Serializable { protected: vector<unsigned short> indx; unsigned int dim; public: HistoKey() {dim=0;} HistoKey(vector<unsigned short>& Indx) { dim=Indx.size(); indx = vector<unsigned short>(dim); for (unsigned int k=0; k<dim; k++) indx[k]=Indx[k]; } HistoKey(const HistoKey& t) {indx=t.indx; dim=t.dim;} void reset(vector<short unsigned>& Indx) { if (Indx.size() != (unsigned)dim) { dim = Indx.size(); indx = vector<unsigned short>(dim); } for (unsigned k=0; k<Indx.size(); k++) indx[k]=Indx[k]; } int get_I(int i) const { return indx[i]; } vector<unsigned short> get_I() const { return indx; } int get_dim() const { return dim; } // 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(indx); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(dim); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; size_t munge_knuth(HistoKey& x); int hash_func_knuth(HistoKey& x, int Size); int hash_func(HistoKey& x, int Size); HASHMAP_NAMESPACE_BEGIN template<> struct hash<HistoKey> { size_t operator()(HistoKey s) const { return munge_knuth(s); } }; HASHMAP_NAMESPACE_END struct eqhkey { bool operator()(const HistoKey& k1, const HistoKey& k2) const { bool ok = true; for (int k=0; k<k1.get_dim(); k++) { if (k1.get_I(k) != k2.get_I(k)) ok = false; } return ok; } }; typedef hashers<HistoKey, double, hash<HistoKey>, eqhkey>::MapType mmapHK; class HistogramCache: public Serializable { private: mmapHK *histo; HistoKey key; vector<double> low, high, width; vector<unsigned short> ndim, indx; unsigned int dim; mmapHK::iterator hitr; public: static int numBuckets; HistogramCache(vector<double>& lo, vector<double>& hi, vector<double>& wid); HistogramCache(HistogramCache& p); ~HistogramCache(); void addValue(vector<double>& z, double value=1); double getValue(vector<unsigned short>& indx); double getValue(vector<double>& z); bool addToExistingValue(vector<double>& z, double value=1); int getNumber(void) { return histo->size(); } mmapHK::iterator Reset() { return (hitr = histo->begin());} mmapHK::iterator Next() { return ++hitr; } mmapHK::iterator Current() { return hitr; } void printList(); // AUTO GENERATED BY ../persistence/autopersist.py protected: HistogramCache() {} 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(histo); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(key); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(low); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(high); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(width); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ndim); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(indx); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(dim); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(numBuckets); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } };