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/SimpleGalaxyModel.h // This is really -*- C++ -*- #ifndef SimpleGalaxyModel_h #define SimpleGalaxyModel_h #include "hash_map.h" #include "Model.h" #include "StdCandle2.h" #include "coord.h" #include "Serializable.h" namespace BIE { class CachedModel: public Serializable { public: double L; double B; StdCandle * std; CachedModel(double ll, double bb, StdCandle *p) { L = ll; B = bb; std = p; } ~CachedModel() { delete std; } // AUTO GENERATED BY ../persistence/autopersist.py protected: CachedModel() {} 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(std); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; typedef hashers<coordPair, CachedModel*, hash<coordPair>, eqcoord>::MapType mmapCM; //+ CLICLASS SimpleGalaxyModel SUPER Model class SimpleGalaxyModel : public Model { public: //+ CLICONSTR int int bool SimpleGalaxyModel(int ndim, int mdim, bool cmap=false); ~SimpleGalaxyModel(); void Initialize(State&); void Initialize(vector<double>& w, vector< vector<double> >& p); void ResetCache() {} double NormEval(double x, double y, SampleDistribution *d); double NormEval(double x, double y); double NormEvalMeasure(double x, double B) { if (mapped) B = asin(pow(fabs(B), 1.0/(1.0 - MU))) * copysign(1.0, B); return cos(B); } vector<double> EvaluateBinned(double x, double y, BinnedDistribution *d); vector<double> EvaluatePoint(double x, double y, PointDistribution *d); string ParameterDescription(int i); int DataDimension() { return 1; } //+ CLIMETHOD void SetStdCandle int int void SetStdCandle(int num, int ntab) { StdCandle::NUM = num; StdCandle::NTAB = ntab; } static double A1; static double Z1; static double K0; static double SIGK; static double KLOW; static double KLIM; static double MU; protected: private: bool mapped; int M; int Ndim; // Constants related to record type. static const char* LENGTH_FIELDNAME; static const char* HEIGHT_FIELDNAME; static const char* PARAM_NAMES[]; StdCandle std0; mmapCM cache; mmapCM::iterator mit; CachedModel* current; vector< vector<double> > pt; vector<double> wt; vector<double> norm, normoff; // AUTO GENERATED BY ../persistence/autopersist.py protected: SimpleGalaxyModel() {} 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(KLOW); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(KLIM); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(MU); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(mapped); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(M); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Ndim); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(std0); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(cache); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(current); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(pt); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(wt); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(norm); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(normoff); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; } // namespace BIE BIE_CLASS_TYPE_INFO(BIE::CachedModel) BIE_CLASS_TYPE_INFO(BIE::SimpleGalaxyModel) BIE_CLASS_EXPORT_KEY(BIE::CachedModel) BIE_CLASS_EXPORT_KEY(BIE::SimpleGalaxyModel) #endif