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/PopModelND.h // This is really -*- C++ -*- #ifndef PopModelND_h #define PopModelND_h #include "Serializable.h" #include <gaussQ.h> #include <Model.h> #include <Histogram1D.h> #include <BIEconfig.h> #include <PopulationsApp.h> namespace BIE { typedef vector<double> dvector; typedef vector<real> rvector; class NDPopCachedModel; //+ CLICLASS PopModelND SUPER Model class PopModelND : public Model { public: //+ CLICONSTR int int SampleDistribution* PopModelND(int ndim, int mdim, SampleDistribution *dist); ~PopModelND(); //+ CLIMETHOD void SetKnots int void SetKnots(int num) { NUM = num; delete intgr; intgr = new JacoQuad(NUM, ALPHA, BETA); } void Initialize(State&); void Initialize(vector<double>& w, vector<double>*& p); void ResetCache() {listCntr = 0;} 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 1; } static int NUM; static double METAL; static double ALPHA; static double BETA; static double R0; 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 double NU; // Defaults: 1.0 static string DataDir; // Defaults: "/proj/weinberg/3Mtab"; protected: PopulationsApp *pops; static double Log10; static double onedeg; 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, *pt; vector<double> work; virtual void generate(double L, double B, SampleDistribution *sd); virtual void compute_bins(); bool good_bounds; void check_bounds(); vector<NDPopCachedModel*> cache; int listCntr; int missed; int nages; EvalType type; // AUTO GENERATED BY ../persistence/autopersist.py protected: PopModelND() {} 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(NUM); 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(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(NU); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(DataDir); 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(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(cache); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(listCntr); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(missed); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(nages); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(type); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; class NDPopCachedModel: public Serializable { public: double L; double B; vector<real> bins; vector<real> R; vector<real> z; vector<real> fac; vector<rvector>* work; NDPopCachedModel(double ll, double bb) { L = ll; B = bb; work = NULL; } ~NDPopCachedModel () { delete [] work; } void SetNumberComp(int nages, int nlos) { delete [] work; work = new vector<rvector> [nages]; for (int i=0; i<nages; i++) work[i] = vector<rvector>(nlos); } // AUTO GENERATED BY ../persistence/autopersist.py protected: NDPopCachedModel() {} 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(bins); 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(fac); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(work); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; } BIE_CLASS_TYPE_INFO(BIE::PopModelND) BIE_CLASS_TYPE_INFO(BIE::NDPopCachedModel) BIE_CLASS_EXPORT_KEY(BIE::PopModelND) BIE_CLASS_EXPORT_KEY(BIE::NDPopCachedModel) #endif