From 87ede9b87ebb7323d167d09f0c7b8efde03c736c Mon Sep 17 00:00:00 2001 From: astrojhgu Date: Sun, 11 Mar 2012 18:00:22 +0000 Subject: git-svn-id: file:///home/svn/opt_utilities@229 ed2142bd-67ad-457f-ba7c-d818d4011675 --- interface/pymodel.hpp | 2 +- statistics/chisq.hpp | 16 ++++++++-------- statistics/cstat.hpp | 4 ++-- statistics/leastsq.hpp | 4 ++-- 4 files changed, 13 insertions(+), 13 deletions(-) diff --git a/interface/pymodel.hpp b/interface/pymodel.hpp index 8165d6d..3a8b495 100644 --- a/interface/pymodel.hpp +++ b/interface/pymodel.hpp @@ -123,7 +123,7 @@ namespace opt_utilities typename element_type_trait::element_type pvalue= boost::python::extract::element_type>(args_values[i]); - push_param_info(param_info(pname,pvalue)); + this->push_param_info(param_info(pname,pvalue)); } } private: diff --git a/statistics/chisq.hpp b/statistics/chisq.hpp index 8a3881a..b9a32f7 100644 --- a/statistics/chisq.hpp +++ b/statistics/chisq.hpp @@ -62,7 +62,7 @@ namespace opt_utilities Ts result(0); for(int i=(this->get_data_set()).size()-1;i>=0;--i) { - Ty chi=(this->get_data_set().get_data(i).get_y()-eval_model(this->get_data_set().get_data(i).get_x(),p))/this->get_data_set().get_data(i).get_y_upper_err(); + Ty chi=(this->get_data_set().get_data(i).get_y()-this->eval_model(this->get_data_set().get_data(i).get_x(),p))/this->get_data_set().get_data(i).get_y_upper_err(); result+=chi*chi; } @@ -134,15 +134,15 @@ namespace opt_utilities Tx x1=this->get_data_set().get_data(i).get_x()-this->get_data_set().get_data(i).get_x_lower_err(); Tx x2=this->get_data_set().get_data(i).get_x()+this->get_data_set().get_data(i).get_x_upper_err(); Tx x=this->get_data_set().get_data(i).get_x(); - Ty errx1=(eval_model(x1,p)-eval_model(x,p)); - Ty errx2=(eval_model(x2,p)-eval_model(x,p)); + Ty errx1=(this->eval_model(x1,p)-this->eval_model(x,p)); + Ty errx2=(this->eval_model(x2,p)-this->eval_model(x,p)); //Ty errx=0; #else Ty errx1=0; Ty errx2=0; #endif - Ty y_model=eval_model(this->get_data_set().get_data(i).get_x(),p); + Ty y_model=this->eval_model(this->get_data_set().get_data(i).get_x(),p); Ty y_obs=this->get_data_set().get_data(i).get_y(); Ty y_err; @@ -182,13 +182,13 @@ namespace opt_utilities Ty chi=(y_obs-y_model)/std::sqrt(y_err*y_err+errx*errx); - // Ty chi=(this->get_data_set().get_data(i).get_y()-eval_model(this->get_data_set().get_data(i).get_x(),p)); + // Ty chi=(this->get_data_set().get_data(i).get_y()-this->eval_model(this->get_data_set().get_data(i).get_x(),p)); // cerr<get_data_set()[i].x,p)<eval_model(this->get_data_set()[i].x,p)<get_data_set()[i].y_upper_err<get_data_set()[i].x<<"\t"<get_data_set()[i].y<<"\t"<get_data_set()[i].x,p)<get_data_set()[i].x<<"\t"<get_data_set()[i].y<<"\t"<eval_model(this->get_data_set()[i].x,p)<get_data_set().get_data(0).get_y().size()); for(int j=0;jget_data_set().get_data(i).get_x(),p)); + Ty model_y(this->eval_model(this->get_data_set().get_data(i).get_x(),p)); if(model_y[j]>this->get_data_set().get_data(i).get_y()[j]) { chi[j]=(this->get_data_set().get_data(i).get_y()[j]-model_y[j])/this->get_data_set().get_data(i).get_y_upper_err()[j]; diff --git a/statistics/cstat.hpp b/statistics/cstat.hpp index 1d845e3..f1c1d8f 100644 --- a/statistics/cstat.hpp +++ b/statistics/cstat.hpp @@ -60,7 +60,7 @@ namespace opt_utilities Ts result(0); for(int i=(this->get_data_set()).size()-1;i>=0;--i) { - Ty model_y=eval_model(this->get_data_set().get_data(i).get_x(),p); + Ty model_y=this->eval_model(this->get_data_set().get_data(i).get_x(),p); result-=contract(this->get_data_set().get_data(i).get_y(),std::log(model_y),result); } @@ -125,7 +125,7 @@ namespace opt_utilities Ts result(0); for(int i=(this->get_data_set()).size()-1;i>=0;--i) { - Ty model_y=eval_model(this->get_data_set().get_data(i).get_x(),p); + Ty model_y=this->eval_model(this->get_data_set().get_data(i).get_x(),p); result-=contract1(this->get_data_set().get_data(i).get_y(),std::log(model_y),result); } diff --git a/statistics/leastsq.hpp b/statistics/leastsq.hpp index f23f0d1..1ecd2dc 100644 --- a/statistics/leastsq.hpp +++ b/statistics/leastsq.hpp @@ -62,7 +62,7 @@ namespace opt_utilities Ts result(0); for(int i=(this->get_data_set()).size()-1;i>=0;--i) { - Ty chi=(this->get_data_set().get_data(i).get_y()-eval_model(this->get_data_set().get_data(i).get_x(),p)); + Ty chi=(this->get_data_set().get_data(i).get_y()-this->eval_model(this->get_data_set().get_data(i).get_x(),p)); result+=chi*chi; } @@ -128,7 +128,7 @@ namespace opt_utilities Ty chi(this->get_data_set().get_data(0).get_y().size()); for(int j=0;jget_data_set().get_data(i).get_x(),p)); + Ty model_y(this->eval_model(this->get_data_set().get_data(i).get_x(),p)); if(model_y[j]>this->get_data_set().get_data(i).get_y()[j]) { chi[j]=(this->get_data_set().get_data(i).get_y()[j]-model_y[j]); -- cgit v1.2.2