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/GalaxyModelOneDHashed.h // This is really -*- C++ -*- #ifndef GalaxyModelOneDHashed_h #define GalaxyModelOneDHashed_h #include "hash_map.h" #include <gaussQ.h> #include <Model.h> #include <Histogram1D.h> #include <CacheGalaxyModel.h> #include "Serializable.h" namespace BIE { //+ CLICLASS GalaxyModelOneDHashed SUPER Model class GalaxyModelOneDHashed : public Model { public: //+ CLICONSTR int int SampleDistribution* GalaxyModelOneDHashed(int ndim, int mdim, SampleDistribution *_dist); ~GalaxyModelOneDHashed(); //+ 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;} void Initialize(State&); void Initialize(vector<double>& w, vector< vector<double> >& p); void ResetCache() {} virtual double NormEval(double x, double y, SampleDistribution *d=NULL); vector<double> Evaluate(double x, double y, SampleDistribution *d=NULL); string ParameterDescription(int i); int DataDimension() { return 1; } static double A1; static double Z1; static double K0; static double SIGK; static int NUM; static double ALPHA; static double BETA; static double R0; 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; 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 CacheGalaxyModel* generate(const coordPair &, SampleDistribution *sd=NULL); virtual CacheGalaxyModel* compute_bins(const coordPair &, SampleDistribution *sd); bool good_bounds; void check_bounds(); mmapGalCM cache; coordPair lastP; int missed; EvalType type; // AUTO GENERATED BY ../persistence/autopersist.py protected: GalaxyModelOneDHashed() {} 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(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(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(cache); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(lastP); 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::GalaxyModelOneDHashed) BIE_CLASS_EXPORT_KEY(BIE::GalaxyModelOneDHashed) #endif