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/GaussTestLikelihoodFunction.h // This is really -*- C++ -*- #ifndef GaussTestLikelihoodFunction_h #define GaussTestLikelihoodFunction_h #include "LikelihoodFunction.h" #include "Serializable.h" namespace BIE { //+ CLICLASS GaussTestLikelihoodFunction SUPER LikelihoodFunction class GaussTestLikelihoodFunction : public LikelihoodFunction { private: vector<double> fdata, pdata; double xmin, xmax, dx; double *centers, *variance, *weights; uint32 number; int dim; int nbins; int N; bool point, cauchyData, cauchyModel; public: //+ CLICONSTR GaussTestLikelihoodFunction(); //+ CLICONSTR int int GaussTestLikelihoodFunction(int Nbins, int N0); //+ CLICONSTR int GaussTestLikelihoodFunction(int N0); ~GaussTestLikelihoodFunction(); //+ CLIMETHOD void useCauchyModel bool void useCauchyModel(bool b) { cauchyModel = b; } //+ CLIMETHOD void useCauchyData bool void useCauchyData(bool b) { cauchyData = b; } //+ CLIMETHOD void newModel string void newModel(string file); //+ CLIMETHOD void SetDim int void SetDim(int n); double LikeProb(std::vector<double> &z, SampleDistribution* sd, double norm, Tile *t, State *s, int indx) { return LocalLikelihood(s); } string ParameterDescription(int i); protected: void defaultModel(); void makeSyntheticData(); void makeSyntheticPointData(); double LocalLikelihood(State* s); // AUTO GENERATED BY ../persistence/autopersist.py private: friend class boost::serialization::access; BOOST_SERIALIZATION_SPLIT_MEMBER(); template<class Archive> void save(Archive & ar, const unsigned int version) const { this->pre_save(ar, version); try { ar << BOOST_SERIALIZATION_BASE_OBJECT_NVP(LikelihoodFunction); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(fdata); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(pdata); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(xmin); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(xmax); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(dx); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(number); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(dim); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(nbins); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(N); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(point); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(cauchyData); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar << BOOST_SERIALIZATION_NVP(cauchyModel); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } char __buf[128]; for(uint32 __i=0; __i<number; __i++) { sprintf(__buf, "centers_%d", __i); try { ar << boost::serialization::make_nvp(__buf, centers[__i]); } catch (::boost::archive::archive_exception & e) { throw new BoostSerializationException(string(e.what()),__buf, __FILE__,__LINE__); } catch (BoostSerializationException * e) { throw new BoostSerializationException(e, __buf,__FILE__,__LINE__); } } for(uint32 __i=0; __i<number; __i++) { sprintf(__buf, "variance_%d", __i); try { ar << boost::serialization::make_nvp(__buf, variance[__i]); } catch (::boost::archive::archive_exception & e) { throw new BoostSerializationException(string(e.what()),__buf, __FILE__,__LINE__); } catch (BoostSerializationException * e) { throw new BoostSerializationException(e, __buf,__FILE__,__LINE__); } } for(uint32 __i=0; __i<number; __i++) { sprintf(__buf, "weights_%d", __i); try { ar << boost::serialization::make_nvp(__buf, weights[__i]); } catch (::boost::archive::archive_exception & e) { throw new BoostSerializationException(string(e.what()),__buf, __FILE__,__LINE__); } catch (BoostSerializationException * e) { throw new BoostSerializationException(e, __buf,__FILE__,__LINE__); } } this->post_save(ar, version); } template<class Archive> void load(Archive & ar, const unsigned int version) { this->pre_load(ar, version); try { ar >> BOOST_SERIALIZATION_BASE_OBJECT_NVP(LikelihoodFunction); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(fdata); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(pdata); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(xmin); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(xmax); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(dx); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(number); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(dim); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(nbins); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(N); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(point); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(cauchyData); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } try { ar >> BOOST_SERIALIZATION_NVP(cauchyModel); BIE_CATCH_BOOST_SERIALIZATION_EXCEPTION; } char __buf[128]; centers = new double[number]; for(uint32 __i=0; __i<number; __i++) { sprintf(__buf, "centers_%d", __i); try { ar >> boost::serialization::make_nvp(__buf, centers[__i]); } catch (::boost::archive::archive_exception & e) { throw new BoostSerializationException(string(e.what()),__buf, __FILE__,__LINE__); } catch (BoostSerializationException * e) { throw new BoostSerializationException(e, __buf,__FILE__,__LINE__); } } variance = new double[number]; for(uint32 __i=0; __i<number; __i++) { sprintf(__buf, "variance_%d", __i); try { ar >> boost::serialization::make_nvp(__buf, variance[__i]); } catch (::boost::archive::archive_exception & e) { throw new BoostSerializationException(string(e.what()),__buf, __FILE__,__LINE__); } catch (BoostSerializationException * e) { throw new BoostSerializationException(e, __buf,__FILE__,__LINE__); } } weights = new double[number]; for(uint32 __i=0; __i<number; __i++) { sprintf(__buf, "weights_%d", __i); try { ar >> boost::serialization::make_nvp(__buf, weights[__i]); } catch (::boost::archive::archive_exception & e) { throw new BoostSerializationException(string(e.what()),__buf, __FILE__,__LINE__); } catch (BoostSerializationException * e) { throw new BoostSerializationException(e, __buf,__FILE__,__LINE__); } } this->post_load(ar, version); } }; } //namespace BIE BIE_CLASS_TYPE_INFO(BIE::GaussTestLikelihoodFunction) BIE_CLASS_EXPORT_KEY(BIE::GaussTestLikelihoodFunction) #endif