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/include/GalaxyModelND.h // This is really -*- C++ -*- #ifndef GalaxyModelND_h #define GalaxyModelND_h #include <gaussQ.h> #include <Model.h> #include <Histogram1D.h> #include <BIEconfig.h> #include <cmath> #include <deque> #include <utility> #include "CacheGalaxyModel.h" #include "Serializable.h" namespace BIE { //+ CLICLASS GalaxyModelND SUPER Model class GalaxyModelND : public Model { public: //+ CLICONSTR int int SampleDistribution* GalaxyModelND(int ndim, int mdim, SampleDistribution *histo); ~GalaxyModelND(); //+ CLIMETHOD void SetKnots int void SetKnots(int num) { NUM = num; delete intgr; intgr = new JacoQuad(NUM, ALPHA, BETA); } //+ CLIMETHOD void CacheLimit int void CacheLimit(int num) { cache_limit=num; } //+ CLIMETHOD void SetExtinction double double void SetExtinction(double A, double Z) {A1 = A; Z1 = Z;} //@! From state vector void Initialize(State&); //@! From component weights and component parameter vectors void Initialize(vector<double>& w, vector< vector<double> >& p); void ResetCache() {} virtual double NormEval(double x, double y, SampleDistribution *d); virtual double NormEval(double x, double y); virtual double NormEvalMeasure(double x, double y) { return cos(y); } vector<double> Evaluate(double x, double y, SampleDistribution *d); string ParameterDescription(int i); int DataDimension() { return 1; } static double LMAG; static double HMAG; static double A1; static double Z1; static double K0; static int NUM; static double ALPHA; static double BETA; static double R0; // Maximum radius (kpc) (default 20.0) static double RMAX; static double AK; static double AMIN; // Defaults: 0.2 static double AMAX; // 8.0 static double HMIN; // 50.0 static double HMAX; // 2000.0 static string BASISDATA; // Defaults: "basisdata.2d" protected: static double Log10; static double onedeg; int M; int Mcur; int Ndim; unsigned int cache_limit; virtual void manageCache(coordPair&); static const char* LENGTH_FIELDNAME; static const char* HEIGHT_FIELDNAME; static const char* PARAM_NAMES[]; int nbins; vector<dvector> lowb; vector<dvector> highb; vector<double> flux; int Nflux; vector<double> lolim; vector<double> hilim; int nparam; vector<double> w; vector<dvector> pos; vector<dvector> sig; JacoQuad *intgr; vector<double> wt; vector< vector<double> > pt; vector<double> work; virtual void generate(double L, double B, SampleDistribution *sd=NULL); virtual void compute_bins(); bool good_bounds; void check_bounds(); mmapGalCM cache; mmapGalCM::iterator mit; deque<coordPair> cacheList; CacheGalaxyModel *current; int missed; EvalType type; // AUTO GENERATED BY ../persistence/autopersist.py protected: GalaxyModelND() {} 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(LMAG); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(HMAG); 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(K0); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(NUM); 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(AK); 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(BASISDATA); 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(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(cache_limit); 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(flux); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Nflux); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(lolim); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(hilim); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(nparam); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(w); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(pos); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(sig); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(intgr); 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(current); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(missed); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(type); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; } BIE_CLASS_TYPE_INFO(BIE::GalaxyModelND) BIE_CLASS_EXPORT_KEY(BIE::GalaxyModelND) #endif