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/PopModelCacheSF.h // This is really -*- C++ -*- #ifndef PopModelCacheSF_h #define PopModelCacheSF_h #include "Serializable.h" #include "hash_map.h" #include <utility> #include <list> using namespace std; #include <gaussQ.h> #include <Model.h> #include <Histogram1D.h> #include <BIEconfig.h> #include <PopulationsApp.h> #include <coord.h> #include <map> namespace BIE { typedef vector<double> dvector; typedef vector<real> rvector; class PopCacheSF; typedef hashers<coordPair, PopCacheSF*, hash<coordPair>, eqcoord>::MapType mmapSF; typedef list<PopCacheSF*> CacheListSF; //+ CLICLASS PopModelCacheSF SUPER Model class PopModelCacheSF : public Model { public: //+ CLICONSTR int int int SampleDistribution* PopulationsApp* string PopModelCacheSF(int ndim, int mdim, int num, SampleDistribution *histo, PopulationsApp *Pop, string basisfile ); ~PopModelCacheSF(); //+ 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: // Enforce normalization in population weights void NormalizePopWeights(vector<double>*& p); void NormalizePopWeights(State& s); 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; // 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; vector<double> ptt; double *ss; virtual PopCacheSF* generate(const coordPair& P, SampleDistribution *d); virtual PopCacheSF* compute_bins(const coordPair& P, SampleDistribution *d); bool good_bounds; void check_bounds(); mmapSF cache; CacheListSF cachelist; coordPair lastP; int nages; int nbasis; vector<double> ages; vector<double> popweight; vector<rvector> basis; 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; private: template<class Archive> void post_load(Archive & ar, const unsigned int version) { ss = &ptt[0]; } // AUTO GENERATED BY ../persistence/autopersist.py protected: PopModelCacheSF() {} 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(ptt); 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(nbasis); 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(basis); 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 PopCacheSF: 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; int N; vector<real> zero; map<int, real, less<int> > bins_map; map<int, vector<real>, less<int> > bins1_map; vector< vector<real> > bins1; vector<real> bins; list<PopCacheSF*>::iterator it; void reset_bin1(void); void reset_bin(void); real* bin1(int j); double bin(int j); PopCacheSF(double ll, double bb, int nages) { L = ll; B = bb; N = nages; zero = vector<real>(N, 0.0); } ~PopCacheSF(); // AUTO GENERATED BY ../persistence/autopersist.py protected: PopCacheSF() {} 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(N); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(zero); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(bins_map); 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; } try { ar & BOOST_SERIALIZATION_NVP(bins); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; } BIE_CLASS_TYPE_INFO(BIE::PopModelCacheSF) BIE_CLASS_TYPE_INFO(BIE::PopCacheSF) BIE_CLASS_EXPORT_KEY(BIE::PopModelCacheSF) BIE_CLASS_EXPORT_KEY(BIE::PopCacheSF) #endif