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/GalaxyModelOneD.h // This is really -*- C++ -*- #ifndef GalaxyModelOneD_h #define GalaxyModelOneD_h #include <deque> #include <gaussQ.h> #include <Model.h> #include <Histogram1D.h> #include <CacheGalaxyModel.h> #include "Serializable.h" namespace BIE { //+ CLICLASS GalaxyModelOneD SUPER Model class GalaxyModelOneD : public Model { public: //+ CLICONSTR int int SampleDistribution* GalaxyModelOneD(int ndim, int mdim, SampleDistribution *_dist); GalaxyModelOneD(); ~GalaxyModelOneD(); //+ CLIMETHOD void SetKnots int void SetKnots(int num) { NUM = num; delete intgr; intgr = new JacoQuad(NUM, ALPHA, BETA); } //+ CLIMETHOD void LogLineOfSight double void LogLineOfSight(double Smin) { logs = true; smin = Smin; } //+ CLIMETHOD void LinearLineOfSight void LinearLineOfSight() { logs = false; } //+ CLIMETHOD void SetExtinction double double void SetExtinction(double A, double Z) {A1 = A; Z1 = Z;} //+ CLIMETHOD void ExtinctionCoefficient double void ExtinctionCoefficient(double ak) { AK = ak; } //+ CLIMETHOD void CacheLimit int void CacheLimit(int num) { cache_limit=num; } //+ CLIMETHOD void ModelLimits double double double double void ModelLimits(double amin, double amax, double hmin, double hmax) { AMIN = amin; AMAX = amax; HMIN = hmin; HMAX = hmax; } //+ CLIMETHOD void MaximumDistance double void MaximumDistance(double rmax) { RMAX=rmax; } void Initialize(State&); void Initialize(vector<double>& w, vector< vector<double> >& p); void ResetCache() {} virtual double NormEval(double x, double y, SampleDistribution *d=NULL); 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=NULL); string ParameterDescription(int i); int DataDimension() { return 1; } static double LMAG; static double HMAG; static double A1; static double Z1; static double K0; static double SIGK; 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; // 1200.0 static bool logs; // Default: false static double smin; // Default: 0.01 protected: static double Log10; static double onedeg; bool cache_full; unsigned int cache_limit; virtual void manageCache(coordPair&); int M; int Mcur; int Ndim; static const char* LENGTH_FIELDNAME; static const char* HEIGHT_FIELDNAME; static const char* PARAM_NAMES[]; int nbins; vector<double> lowb; vector<double> highb; JacoQuad *intgr; vector<double> wt; vector< vector<double> > pt; 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; EvalType type; // AUTO GENERATED BY ../persistence/autopersist.py 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(SIGK); 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(logs); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(smin); 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(cache_limit); 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(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(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(good_bounds); 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::GalaxyModelOneD) BIE_CLASS_EXPORT_KEY(BIE::GalaxyModelOneD) #endif