aboutsummaryrefslogtreecommitdiffstats
path: root/vmodel/bl.hpp
diff options
context:
space:
mode:
authorastrojhgu <astrojhgu@ed2142bd-67ad-457f-ba7c-d818d4011675>2009-12-17 17:01:15 +0000
committerastrojhgu <astrojhgu@ed2142bd-67ad-457f-ba7c-d818d4011675>2009-12-17 17:01:15 +0000
commit4170b13f464439ab97efadb210dab0d749fa4b42 (patch)
treef6d31514c19a9a015c4d22477f5a3fb8bccd250b /vmodel/bl.hpp
parenteeafd6e502bc60da664de3cf0eb9251f349df228 (diff)
downloadopt-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.hpp67
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