diff options
author | astrojhgu <astrojhgu@ed2142bd-67ad-457f-ba7c-d818d4011675> | 2009-12-17 17:01:15 +0000 |
---|---|---|
committer | astrojhgu <astrojhgu@ed2142bd-67ad-457f-ba7c-d818d4011675> | 2009-12-17 17:01:15 +0000 |
commit | 4170b13f464439ab97efadb210dab0d749fa4b42 (patch) | |
tree | f6d31514c19a9a015c4d22477f5a3fb8bccd250b /vmodel/bl.hpp | |
parent | eeafd6e502bc60da664de3cf0eb9251f349df228 (diff) | |
download | opt-utilities-4170b13f464439ab97efadb210dab0d749fa4b42.tar.bz2 |
rename vmodel to vmodels
git-svn-id: file:///home/svn/opt_utilities@103 ed2142bd-67ad-457f-ba7c-d818d4011675
Diffstat (limited to 'vmodel/bl.hpp')
-rw-r--r-- | vmodel/bl.hpp | 67 |
1 files changed, 0 insertions, 67 deletions
diff --git a/vmodel/bl.hpp b/vmodel/bl.hpp deleted file mode 100644 index 84d6ca5..0000000 --- a/vmodel/bl.hpp +++ /dev/null @@ -1,67 +0,0 @@ -#ifndef BROKEN_LINE_MODEL_H_ -#define BROKEN_LINE_MODEL_H_ -#define OPT_HEADER -#include <core/fitter.hpp> -#include <cmath> -#include <misc/optvec.hpp> -namespace opt_utilities -{ - template <typename T> - class bl - :public model<optvec<T>,optvec<T>,optvec<T>,std::string> - { - private: - model<optvec<T>,optvec<T>,optvec<T>,std::string >* do_clone()const - { - return new bl<T>(*this); - } - - const char* do_get_type_name()const - { - return "broken linear model"; - } - public: - bl() - { - this->push_param_info(param_info<optvec<T> >("break point y value",1)); - this->push_param_info(param_info<optvec<T> >("break point x value",1)); - this->push_param_info(param_info<optvec<T> >("slop 1",1)); - this->push_param_info(param_info<optvec<T> >("slop 2",1)); - } - - public: - optvec<T> do_eval(const optvec<T>& x,const optvec<T>& param) - { - T x_b=get_element(param,0); - T f_b=get_element(param,1); - T k1=get_element(param,2); - T k2=get_element(param,3); - optvec<double> result(x.size()); - for(int i=0;i<x.size();++i) - { - if(x[i]<x_b) - { - result[i]=k1*(x[i]-x_b)+f_b; - } - else - { - result[i]=k2*(x[i]-x_b)+f_b; - } - } - return result; - } - - private: - std::string do_get_information()const - { - return "broken linear model\n" - "y=k1*(x-x_b)+y_b for x<x_b\n" - "y=k2*(x-x_b)+y_b otherwise\n"; - } - }; -} - - - -#endif -//EOF |