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/PopModelCacheF.h // This is really -*- C++ -*- #ifndef PopModelCacheF_h #define PopModelCacheF_h #include "Serializable.h" #include <unistd.h> #include <sys/resource.h> #include "hash_map.h" #include <utility> using namespace std; #include <gaussQ.h> #include <Model.h> #include <Histogram1D.h> #include <BIEconfig.h> #include <PopulationsApp.h> #include <coord.h> #include <list> #include <map> namespace BIE { typedef vector<double> dvector; typedef vector<real> rvector; class PopCacheF; typedef hashers<coordPair, PopCacheF*, hash<coordPair>, eqcoord>::MapType mmap; typedef list<PopCacheF*> CacheList; //+ CLICLASS PopModelCacheF SUPER Model class PopModelCacheF : public Model { public: //+ CLICONSTR int int int SampleDistribution* PopulationsApp* PopModelCacheF(int ndim, int mdim, int num, SampleDistribution *histo, PopulationsApp *Pop); ~PopModelCacheF(); //+ CLIMETHOD void SetLineOfSight double double bool void SetLineOfSight(double rmin, double rmax, bool rlog) { RMIN = rmin; RMAX = rmax; RLOG = rlog; } //+ CLIMETHOD void UseVectorCache void UseVectorCache() { use_vector = true; } //+ CLIMETHOD void UseHashCache void UseHashCache() { use_vector = false; } //+ CLIMETHOD void SetMagNorm double double void SetMagNorm(double x, double y) { mnorm = x * pow(M_PI/(180.0*60.0*60.0), 2.0); mback = y; mm_use = true; } //+ CLIMETHOD double GetEMag0 double GetEMag0() { return m0flux; } //+ CLIMETHOD double GetEMag1 double GetEMag1() { return m1flux; } //+ CLIMETHOD double GetEMag2 double GetEMag2() { return m2flux; } void Initialize(State&); void Initialize(vector<double>& w, vector<double>*& p); void ResetCache() {} virtual double NormEval(double x, double y, SampleDistribution *d); virtual vector<double> Evaluate(double x, double y, SampleDistribution *d); string ParameterDescription(int i); int DataDimension() { return Nflux; } void PrintCache(); //+ METHOD void SetExtinctionFactor double void SetExtinctionFactor(double efac) { EFAC = efac; } static double METAL; static double ALPHA; static double BETA; static double R0; static double A1; static double Z1; static double Z0; static double EFAC; static double EMAX; static bool DELTA; static double RT0; static double DRT0; static int NEBV; static bool RING; static double RINGWIDTH; static double RINGAMP; static double RINGLOC; static double CENTERWIDTH; static double CENTERAMP; static double NEARBYWIDTH; static double NEARBYAMP; static bool RLOG; static double RMIN; static double RMAX; static double AMIN; // Defaults: 0.2 static double AMAX; // 8.0 static double HMIN; // 50.0 static double HMAX; // 1200.0 static BirthRate BRmodel; // Default: Exponential static double BURST; // Defaults: 0.0 static double NU; // Defaults: 1.0 static string DataDir; // Defaults: "/proj/weinberg/3Mtab"; static double zerotol; // Defaults: 1.0e-12 static bool use_vector; // Defaults: true protected: int NUM; bool cmd; PopulationsApp *pops; static double Log10; static double onedeg; pid_t pid; int nbins; vector<dvector> lowb; vector<dvector> highb; int Nflux; // Stellar components // int @autopersist(nparam); // vector<double> @autopersist(w); // vector<dvector> @autopersist(pos), @autopersist(sig); JacoQuad *intgr; int M; int Mcur; int Ndim; vector<double> wt; vector<double> *pt; vector<double> work; virtual PopCacheF* generate(const coordPair& P, SampleDistribution *d); virtual PopCacheF* compute_bins(const coordPair& P, SampleDistribution *d); bool good_bounds; void check_bounds(); mmap cache; CacheList cachelist; coordPair lastP; int nages; vector<double> ages; vector<double> popweight; vector<bool> marge; unsigned maxvmsize; unsigned vmsize; unsigned max_element; EvalType type; double m0flux; double m1flux; double m2flux; double mmean; double mmean2; double mavg; double mnorm; double mback; void compute_mm(); double mm_fraction(double mm); bool mm_computed; bool mm_use; // AUTO GENERATED BY ../persistence/autopersist.py protected: PopModelCacheF() {} 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(Model); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(METAL); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ALPHA); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(BETA); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(R0); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(A1); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Z1); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Z0); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(EFAC); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(EMAX); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(DELTA); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(RT0); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(DRT0); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(NEBV); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(RING); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(RINGWIDTH); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(RINGAMP); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(RINGLOC); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(CENTERWIDTH); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(CENTERAMP); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(NEARBYWIDTH); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(NEARBYAMP); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(RLOG); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(RMIN); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(RMAX); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(AMIN); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(AMAX); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(HMIN); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(HMAX); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(BRmodel); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(BURST); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(NU); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(DataDir); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(zerotol); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(use_vector); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(NUM); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(cmd); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(pops); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Log10); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(onedeg); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(pid); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(nbins); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(lowb); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(highb); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Nflux); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(intgr); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(M); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Mcur); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Ndim); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(wt); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(pt); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(work); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(good_bounds); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(cache); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(cachelist); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(lastP); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(nages); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ages); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(popweight); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(marge); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(maxvmsize); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(vmsize); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(max_element); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(type); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(m0flux); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(m1flux); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(m2flux); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(mmean); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(mmean2); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(mavg); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(mnorm); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(mback); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(mm_computed); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(mm_use); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; class PopCacheF: public Serializable { public: double L; double B; vector<real> R; vector<real> z; vector<real> s; vector<real> mm0; vector<real> mm1; vector<real> dm; list<PopCacheF*>::iterator it; map<int, double, less<int> > bins_map; vector<real> bins; void reset_bin(void); double bin(int j); map<int, double, less<int> > bins1_map; vector<real> bins1; void reset_bin1(void); double bin1(int j); PopCacheF(double ll, double bb) { L = ll; B = bb; } // AUTO GENERATED BY ../persistence/autopersist.py protected: PopCacheF() {} 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(L); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(B); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(R); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(z); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(s); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(mm0); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(mm1); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(dm); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(bins_map); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(bins); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(bins1_map); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(bins1); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; } BIE_CLASS_TYPE_INFO(BIE::PopModelCacheF) BIE_CLASS_TYPE_INFO(BIE::PopCacheF) BIE_CLASS_EXPORT_KEY(BIE::PopModelCacheF) BIE_CLASS_EXPORT_KEY(BIE::PopCacheF) #endif