1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
|
#include "strmodel1d.hpp"
using namespace mu;
using namespace std;
using namespace opt_utilities;
strmodel1d* strmodel1d::do_clone()const
{
return new strmodel1d(*this);
}
strmodel1d::strmodel1d()
{
set_buildin_fun();
}
strmodel1d::strmodel1d(const strmodel1d& rhs)
:par_vec(rhs.par_vec),
par_names(rhs.par_names),
expr(rhs.expr),
var_name(rhs.var_name)
{
set_buildin_fun();
set_expr(expr,par_names,var_name);
}
strmodel1d& strmodel1d::operator=(const strmodel1d& rhs)
{
set_buildin_fun();
expr=rhs.expr;
par_names=rhs.par_names;
var_name=rhs.var_name;
par_vec=rhs.par_vec;
set_expr(expr,par_names,var_name);
return *this;
}
void strmodel1d::set_buildin_fun()
{
mp.DefineFun("sin",sin);
mp.DefineFun("cos",cos);
mp.DefineFun("exp",exp);
mp.DefineFun("tan",tan);
mp.DefineFun("log",log);
}
void strmodel1d::set_expr(const string& _expr,
const std::vector<std::string>& _par_names,
const std::string& _var_name)
{
expr=_expr;
par_names=_par_names;
var_name=_var_name;
this->clear_param_info();
par_vec.resize(par_names.size());
mp.ClearVar();
// mp.ClearFun();
for(unsigned int i=0;i<par_vec.size();++i)
{
mp.DefineVar(par_names[i].c_str(),&par_vec[i]);
this->push_param_info(param_info<std::vector<double> >(par_names[i],0));
}
mp.DefineVar(var_name.c_str(),&x);
mp.SetExpr(expr.c_str());
}
double strmodel1d::do_eval(const double& _x,const vector<double>& p)
{
for(unsigned int i=0;i<par_vec.size();++i)
{
//get_element(par_vec,i)=get_element(p,i);
set_element(par_vec,i,get_element(p,i));
// cout<<par_vec[i]<<" ";
}
x=_x;
// cout<<x<<" ";
double result(mp.Eval());
// cout<<result<<endl;
return result;
}
|
