diff options
| -rw-r--r-- | vmodels/nbeta1d.hpp | 62 |
1 files changed, 62 insertions, 0 deletions
diff --git a/vmodels/nbeta1d.hpp b/vmodels/nbeta1d.hpp new file mode 100644 index 0000000..b98d179 --- /dev/null +++ b/vmodels/nbeta1d.hpp @@ -0,0 +1,62 @@ +#ifndef VNBETA_MODEL_H_ +#define VNBETA_MODEL_H_ +#define OPT_HEADER +#include <core/fitter.hpp> +#include <misc/optvec.hpp> +#include <cmath> + +namespace opt_utilities +{ + template <typename T> + class nbeta1d + :public model<optvec<T>,optvec<T>,optvec<T>,std::string> + { + typedef optvec<T> Tv; + private: + nbeta1d<T>* do_clone()const + { + return new nbeta1d<T>(*this); + } + + const char* do_get_type_name()const + { + return "1d density beta model"; + } + public: + nbeta1d() + { + this->push_param_info(param_info<Tv>("n0",1)); + this->push_param_info(param_info<Tv>("rc",10)); + this->push_param_info(param_info<Tv>("beta",2./3.)); + } + + public: + Tv do_eval(const Tv& x,const Tv& param) + { + Tv result(x.size()); + T n0=get_element(param,0); + T r_c=get_element(param,1); + T beta=get_element(param,2); + + //return x*get_element(param,0)+get_element(param,1); + for(size_t i=0;i!=x.size();++i) + { + + result[i]=n0*pow(1+(x[i]*x[i])/(r_c*r_c),-3./2.*beta); + + } + return result; + } + + private: + std::string do_get_information()const + { + return ""; + } + }; +} + + + +#endif +//EOF |