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/PopulationsApp.h // This is really -*- C++ -*- #ifndef PopulationsApp_h #define PopulationsApp_h #include "Serializable.h" #include <cmath> #include <ostream> #include <string> #include <vector> #include <algorithm> using namespace std; #include <BIEException.h> #include <Populations.h> namespace BIE { enum BirthRate {Burst, Exponential}; //+ CLICLASS PopulationsApp class PopulationsApp: public Serializable { typedef vector<double> dvector; typedef vector<int> ivector; typedef hashers<PopKey, double, hash<PopKey>, eqpopkey>::MapType HMapKey; private: bool cmd; size_t nmetal; vector<double> metal; vector<string> smetal; size_t nages; vector<double> ages, age_boundary; vector<string> sages; size_t ndistm; vector<double> distm; vector<string> sdistm; size_t nflux0, nflux, nMarge, nSize, nSize2; vector<double> minmag; vector<double> maxmag; vector<double> wmag; vector<int> nmag; vector<bool> use_flux; vector<double> flux; vector<int> ii; vector<int> iflx, iflx1, jflx, jflx1; vector<int> iflxmin; vector<int> iflxmax; vector<int> iflxoff; vector<double> flxwght, wght; vector<dvector> f2; vector<ivector> i2; vector<double> ff; vector<double> vvector; vector<double> fflux; int ncorners; string datadir, isoname; int AgeIndex(double age); int MetalIndex(double z); int DistModIndex(double m); int get_index(vector<int>&); int get_index2(vector<int>&); double array_interp(double *array, double *flux); double cumul_interp(double *array, double *flux); void Cumulate(int iage, int idistm, int imetal); void debug_accumulation(int iage, int idistm, int imetal); void debug_norm(); PopKey key; void read_data(void); vector< vector<HMapKey> > Hash; vector< vector< vector< vector<double> > > > Array, Cumul; public: static bool verbose; static bool use_array; //+ CLICONSTR string clivectori* bool PopulationsApp(string directory, vector<int>* Marginalize=NULL, bool CMD=false); //+ CLICONSTR string string clivectori* PopulationsApp(string directory, string isoname, vector<int>* Marginalize=NULL); ~PopulationsApp(); double GetValue(vector<double>& flux, double age, double metal, double distmod); double GetValue(vector<int>& indx, int age, int metal, int distmod); double Interpolate(vector<double>& flux, double age, double metal, double distmod); double Bin(vector<double>& fmin, vector<double>& fmax, double age, double metal, double distmod); void SetParameters(void) {}; vector<bool> get_marge(void) const { return use_flux; } vector<double> Metal() { return metal; } vector<double> Ages() { return ages; } vector<double> DistMod() { return distm; } vector<double> AgeBoundaries() { return age_boundary; } vector<double> MinMag() { return minmag; } vector<double> MaxMag() { return maxmag; } vector<double> MagDelta() { return wmag; } vector<int> NumMag() { return nmag; } bool CMDflag(void) { return cmd; } void VerifyDouble(); void VerifyInt(); //+ CLIMETHOD void UseArray void UseArray() { use_array = true; } //+ CLIMETHOD void UseHash void UseHash() { use_array = false; } // AUTO GENERATED BY ../persistence/autopersist.py protected: PopulationsApp() {} 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(cmd); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(nmetal); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(metal); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(smetal); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(nages); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ages); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(age_boundary); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(sages); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ndistm); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(distm); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(sdistm); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(nflux0); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(nflux); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(nMarge); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(nSize); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(nSize2); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(minmag); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(maxmag); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(wmag); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(nmag); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(use_flux); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(flux); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ii); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(iflx); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(iflx1); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(jflx); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(jflx1); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(iflxmin); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(iflxmax); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(iflxoff); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(flxwght); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(wght); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(f2); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(i2); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ff); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(vvector); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(fflux); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ncorners); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(datadir); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(isoname); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(key); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Hash); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Array); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(Cumul); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; } BIE_CLASS_TYPE_INFO(BIE::PopulationsApp) BIE_CLASS_EXPORT_KEY(BIE::PopulationsApp) #endif