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/modelCMD.h // This is really -*- C++ -*- #include <time.h> #include <cstdlib> #include <cmath> #include <iostream> #ifndef modelCMD_h #define modelCMD_h using namespace std; #include "gvariable.h" #include <Normal.h> #include <Uniform.h> #include "Serializable.h" namespace BIE { class star: public Serializable { public: float h, j, k; star() {} // 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(Serializable); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(h); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(j); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(k); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; class CMD2: public Serializable { typedef vector<int> ivector; typedef vector<float> fvector; typedef vector<fvector> ffvector; private: int Nbins, iSFR; fvector SFR, t_init, t_final, metal; Normal *normal; Uniform *unit; //bool bar, disk, normalize; int ntotal, nages, initialized, imetal; int imetal1, imetal2; ivector number_ages; ffvector log_age, age, M_ini, M_act; ffvector log_L, log_Te, log_G; ffvector Mbol, Mu, Mb, Mv, Mr, Mi, Mj; ffvector Mh, Mk, CO, McTPAGB; fvector isoc_Z; ivector isoc_sizes; fvector list_ages, list_ages_log; float random_age, logage; fvector sfrh_time, sfrh_bar, sfrh_disk; fvector amr_time, amr_bar_box, amr_bar_infall_04; fvector amr_bar_infall_06, amr_bar_infall_max; fvector amr_bar_outfall_05, amr_bar_outfall_10; fvector amr_bar_inout_04, amr_bar_inout_05; fvector amr_bar_inout_06, amr_disk_box_014; fvector amr_disk_box_008; fvector amr_disk_infall_03, amr_disk_infall_05; fvector amr_disk_infall_08; fvector amr_disk_infall_max, amr_disk_outfall_06; fvector amr_disk_outfall_13, amr_disk_outfall_20; fvector amr_disk_outfall_25, amr_disk_inout_01; fvector amr_disk_inout_005, amr_disk_inout_02; float rm1, rm2, rm3, rm4, deltat; int t1, t2; int iage_t1_Z1, iage_t2_Z1; int iage_t1_Z2, iage_t2_Z2; float matched_age_t1, matched_age_t2; fvector Mh_t1_Z1, Mj_t1_Z1, Mk_t1_Z1; fvector log_age_t1_Z1, M_ini_t1_Z1, M_act_t1_Z1; fvector Mh_t1_Z2, Mj_t1_Z2, Mk_t1_Z2; fvector log_age_t1_Z2, M_ini_t1_Z2, M_act_t1_Z2; fvector Mh_t2_Z1, Mj_t2_Z1, Mk_t2_Z1; fvector log_age_t2_Z1, M_ini_t2_Z1, M_act_t2_Z1; fvector Mh_t2_Z2, Mj_t2_Z2, Mk_t2_Z2; fvector log_age_t2_Z2, M_ini_t2_Z2, M_act_t2_Z2; ffvector mass_A, mass_B, mass_C, mass_D; float massmax, minmass; float minmass_t1_Z1, minmass_t2_Z1; float minmass_t1_Z2, minmass_t2_Z2; float maxmass_t1_Z1, maxmass_t2_Z1; float maxmass_t1_Z2, maxmass_t2_Z2; float Z; float threshold; public: // CMD2() {}; CMD2(); CMD2(int); CMD2(int, bool, bool); CMD2(int, int); CMD2(int, int, int); ~CMD2(); int init_SFR(int); int init_isoc(); int sample_from_SFR_history(); int find_mass_range(); float get_gaussian_rn(float); float get_random_time(int, float); float get_random_gauss(float, float); float get_result_disk_jk(int, int); //float get_result_disk_jkk(int, int); float get_result_disk_jk(float, float, float, float); //float get_result_disk_jkk(float, float, float, float); float get_result_bar_jk(int, int); //float get_result_bar_jkk(int, int); float get_result_bar_jk(float, float, float, float); //float get_result_bar_jkk(float, float, float, float); int set_random_time(int, void **); float get_random_age(); int get_IMF_flag(float mass); //int run_model(int, void **); void run_model(); void CMD2_init(int, int, int); void setbar(); void setdisk(); void setnbinsk(int); void setnbinsj(int); void setnbinsjk(int); void setnbins(int, int, int); void setjrange(float, float); void setkrange(float, float); void setjkrange(float, float); void setranges(float, float, float, float, float, float); void dontnormalize(); void donormalize(); void setthreshold(float); // float get_jmin(); float get_jmax(); float get_kmin(); float get_kmax(); /* int istar; int idebug, ierror, ifixed; int ibuffer; //int iprocessed; fvector jkm, jhm, hm, jm, km, dhm, djm, dkm, met, mact, agem, minim; fvector hm_disk, jm_disk, km_disk, hm_bar, jm_bar, km_bar; */ vector<float> kvsjk_bar, kvsj_bar, kvsjk_disk, kvsj_disk; vector<int> ikvsjk_bar, ikvsj_bar, ikvsjk_disk, ikvsj_disk; float kmin, kmax, jmin, jmax, jkmin, jkmax; int nbins_k, nbins_j, nbins_jk; //struct star { // float h; // float j; // float k; //}; vector<star> dstar; vector<star> dstar_f; vector<int> dstar_i; vector<star> bstar; vector<star> bstar_f; vector<int> bstar_i; struct sortbyk { bool operator () (const star & lhs , const star & rhs) const { return lhs.k < rhs.k; } }; // 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(Serializable); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(kvsjk_bar); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(kvsj_bar); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(kvsjk_disk); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(kvsj_disk); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ikvsjk_bar); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ikvsj_bar); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ikvsjk_disk); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(ikvsj_disk); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(kmin); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(kmax); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(jmin); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(jmax); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(jkmin); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(jkmax); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(nbins_k); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(nbins_j); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(nbins_jk); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(dstar); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(dstar_f); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(dstar_i); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(bstar); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(bstar_f); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar & BOOST_SERIALIZATION_NVP(bstar_i); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } this->post_serialize(ar, version); } }; } BIE_CLASS_TYPE_INFO(BIE::star) BIE_CLASS_TYPE_INFO(BIE::CMD2) BIE_CLASS_EXPORT_KEY(BIE::star) BIE_CLASS_EXPORT_KEY(BIE::CMD2) #endif