From 412f4f074e758b4889a41a32af571e6f6b0a3701 Mon Sep 17 00:00:00 2001
From: astrojhgu <astrojhgu@ed2142bd-67ad-457f-ba7c-d818d4011675>
Date: Wed, 13 Apr 2011 16:42:47 +0000
Subject: git-svn-id: file:///home/svn/opt_utilities@187
 ed2142bd-67ad-457f-ba7c-d818d4011675

---
 statistics/leastsq.hpp | 64 ++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 64 insertions(+)

diff --git a/statistics/leastsq.hpp b/statistics/leastsq.hpp
index 3b5ab37..f23f0d1 100644
--- a/statistics/leastsq.hpp
+++ b/statistics/leastsq.hpp
@@ -8,6 +8,7 @@
 #define LEAST_SQ_HPP
 #define OPT_HEADER
 #include <core/fitter.hpp>
+#include <misc/optvec.hpp>
 #include <iostream>
 #include <vector>
 #include <cmath>
@@ -84,9 +85,72 @@ namespace opt_utilities
       return result;
     }
   };
+
+  template<typename T,typename Ts,typename Tstr>
+  class leastsq<optvec<T>,optvec<T>,optvec<T>,Ts,Tstr>
+    :public statistic<optvec<T>,optvec<T>,optvec<T>,Ts,Tstr>
+  {
+  private:
+    bool verb;
+    int n;
+    
+    typedef optvec<T> Tx;
+    typedef optvec<T> Ty;
+    typedef optvec<T> Tp;
+    statistic<Ty,Tx,Tp,Ts,Tstr>* do_clone()const
+    {
+      // return const_cast<statistic<Ty,Tx,Tp>*>(this);
+      return new leastsq<Ty,Tx,Tp,Ts,Tstr>(*this);
+    }
+
+    const char* do_get_type_name()const
+    {
+      return "least square statistic";
+    }
+    
+  public:
+    void verbose(bool v)
+    {
+      verb=v;
+    }
+  public:
+    leastsq()
+      :verb(false)
+    {}
+    
+    
+
+    Ts do_eval(const Tp& p)
+    {
+      Ts result(0);
+      for(int i=(this->get_data_set()).size()-1;i>=0;--i)
+	{
+	  Ty chi(this->get_data_set().get_data(0).get_y().size());
+	  for(int j=0;j<chi.size();++j)
+	    {
+	      Ty model_y(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]);
+		}
+	      else
+		{
+		  chi[j]=(this->get_data_set().get_data(i).get_y()[j]-model_y[j]);
+		}
+	    }
+	  result+=sum(chi*chi);
+	  
+	}
+      return result;
+    }
+  };
+  
+
   
 }
 
+
+
 #endif
 //EOF
 
-- 
cgit v1.2.2