aboutsummaryrefslogtreecommitdiffstats
path: root/models
diff options
context:
space:
mode:
authorastrojhgu <astrojhgu@ed2142bd-67ad-457f-ba7c-d818d4011675>2012-01-19 13:06:54 +0000
committerastrojhgu <astrojhgu@ed2142bd-67ad-457f-ba7c-d818d4011675>2012-01-19 13:06:54 +0000
commit5855acc5ef57ec0fcc2b43447df939ece970c40c (patch)
tree10c13374f2668f34ae69aa37f62929b0581c51d4 /models
parent87c4c388b6848c6ec57e905509b876dd265588f1 (diff)
downloadopt-utilities-5855acc5ef57ec0fcc2b43447df939ece970c40c.tar.bz2
git-svn-id: file:///home/svn/opt_utilities@218 ed2142bd-67ad-457f-ba7c-d818d4011675
Diffstat (limited to 'models')
-rw-r--r--models/quad_pl.hpp52
1 files changed, 52 insertions, 0 deletions
diff --git a/models/quad_pl.hpp b/models/quad_pl.hpp
new file mode 100644
index 0000000..0619d2b
--- /dev/null
+++ b/models/quad_pl.hpp
@@ -0,0 +1,52 @@
+/**
+ \file quad_pl.hpp
+ \brief power law model
+ \author Junhua Gu
+ */
+
+
+#ifndef QUAD_PL_MODEL_H_
+#define QUAD_PL_MODEL_H_
+#define OPT_HEADER
+#include <core/fitter.hpp>
+#include <cmath>
+#include <misc/optvec.hpp>
+
+namespace opt_utilities
+{
+ template <typename T>
+ class quad_pl
+ :public model<T,T,std::vector<T>,std::string>
+ {
+ private:
+ quad_pl* do_clone()const
+ {
+ return new quad_pl<T>(*this);
+ }
+
+ const char* do_get_type_name()const
+ {
+ return "quad power law";
+ }
+ public:
+ quad_pl()
+ {
+ this->push_param_info(param_info<std::vector<T> >("norm",1));
+ this->push_param_info(param_info<std::vector<T> >("gamma",1));
+ this->push_param_info(param_info<std::vector<T> >("corr",0));
+ }
+
+ T do_eval(const T& x,const std::vector<T>& param)
+ {
+ T norm=get_element(param,0);
+ T gamma=get_element(param,1);
+ T corr=get_element(param,2);
+ return norm*pow(x,gamma)*exp(corr*log(x)*log(x));
+ }
+ };
+}
+
+
+
+#endif
+//EOF