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/CMD.h // This may look like C code, but it is really -*- C++ -*- #ifndef _CMD_h #define _CMD_h #include "Serializable.h" #include <cmath> #include <iostream> #include <string> #include <vector> #include <algorithm> #include <map> using namespace std; #include "hash_map.h" #include <ACG.h> #include <Uniform.h> #include <Populations.h> size_t ivecHash(vector<int> P); HASHMAP_NAMESPACE_BEGIN template <> struct hash< vector<int> > { size_t operator()(vector<int> s) const { return ivecHash(s); } }; HASHMAP_NAMESPACE_END struct eqivector { bool operator()(vector<int> c1, vector<int> c2) const { int n1 = c1.size(), n2 = c2.size(); if (n1 != n2) return false; for (int i=0; i<n1; i++) if (c1[i]!=c2[i]) return false; return true; } }; namespace BIE { //+ CLICLASS CMD class CMD: public Serializable { typedef vector<double> dvector; typedef vector<int> ivector; typedef hashers<ivector, float, hash<ivector>, eqivector>::MapType HMap; private: int ncolors; vector<double> colormin; vector<double> colormax; vector<double> wcolor; vector<int> nbins; vector<int> ncol; vector<ivector> colors; double distmod; vector<double> colr; vector<int> ii; vector<int> icol; vector<int> icolmin; vector<int> icolmax; vector<int> icoloff; vector<double> colwght, off; vector<dvector> f2; vector<ivector> i2; void initialize_realization(); ACG *gen; Uniform *uniform; bool realize_setup; double realize_norm; map<double, int, less<double> > realize_data; string datafile; void read_data(istream&, double); HMap Hash; public: static bool verbose; static int seed; //+ CLICONSTR string CMD(string directory, double minval=1.0e-10); ~CMD(); double GetValue(vector<double>& flux); double GetValue(vector<int>& indx); double Interpolate(vector<double>& flux, double distmod); double Bin(vector<double>& fmin, vector<double>& fmax, double distmod); void SetParameters(void) {}; vector<int> get_color(int n) { int nn = min<int>(max<int>(0, n), ncolors-1); return colors[nn]; } vector<ivector> get_colors(void) { return colors; } int get_ncolors(void) { return ncolors; } vector<double> MinMag() { return colormin; } vector<double> MaxMag() { return colormax; } vector<double> ColorDelta() { return wcolor; } vector<int> NumBins() { return nbins; } double DistMod() { return distmod; } void VerifyDouble(double distmod); void VerifyInt(); vector<double> genPoint(); void HashStat(); // AUTO GENERATED BY ../persistence/autopersist.py protected: CMD() {} 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(ncolors); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(colormin); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(colormax); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(wcolor); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(nbins); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ncol); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(colors); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(distmod); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(colr); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ii); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(icol); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(icolmin); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(icolmax); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(icoloff); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(colwght); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(off); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(f2); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(i2); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(gen); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(uniform); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(realize_setup); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(realize_norm); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(realize_data); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(datafile); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Hash); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(verbose); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(seed); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; } BIE_CLASS_TYPE_INFO(BIE::CMD) BIE_CLASS_EXPORT_KEY(BIE::CMD) #endif