git-svn-id: file:///home/svn/opt_utilities@1 ed2142bd-67ad-457f-ba7c-d818d4011675

1 files changed, 180 insertions, 0 deletions

diff --git a/misc/data_loaders.hpp b/misc/data_loaders.hpp
new file mode 100644
index 0000000..6d2195d
--- /dev/null
+++ b/misc/data_loaders.hpp
@@ -0,0 +1,180 @@
+#ifndef DATA_LOADERS_H
+#define DATA_LOADERS_H
+#include <core/fitter.hpp>
+#include <core/default_data_set.hpp>
+#include <iostream>
+#include <fstream>
+#include <cmath>
+
+namespace opt_utilities
+{
+  template <typename Ty,typename Tx>
+  class dl_x_y_ye;
+
+  template <typename Ty,typename Tx>
+  std::istream& operator>>(std::istream& ifs,dl_x_y_ye<Ty,Tx>& dl);
+
+  template <typename Ty,typename Tx>
+  class dl_x_xe_y_ye;
+
+  template <typename Ty,typename Tx>
+  std::istream& operator>>(std::istream& ifs,dl_x_xe_y_ye<Ty,Tx>& dl);
+
+  template <typename Ty,typename Tx>
+  class dl_x_xu_xl_y_yu_yl;
+  
+  template <typename Ty,typename Tx>
+  std::istream& operator>> (std::istream& ifs,dl_x_xu_xl_y_yu_yl<Ty,Tx>& dl);
+
+  template <typename Ty,typename Tx>
+  class dl_x_y_ye
+  {
+  private:
+    default_data_set<Ty,Tx> ds;
+  public:
+    data_set<Ty,Tx>& get_data_set()
+    {
+      return ds;
+    }
+
+    void load_from(std::istream& ifs)
+    {
+      for(;;)
+	{
+	  Tx x;
+	  Tx x_err;
+	  Ty y;
+	  Ty y_err(1);
+	  
+	  ifs>>x>>y>>y_err;
+	  
+	  if(!ifs.good())
+	    {
+	      break;
+	    }
+	  data<Ty,Tx> d(x,y,y_err,y_err,x_err,x_err);
+	  ds.push_back(d);
+	}
+      //return ifs;
+    }
+
+    void load_from(const char* name)
+    {
+      std::ifstream ifs(name);
+      load_from(ifs);
+
+    }
+    //friend std::istream& operator>> <>(std::istream& ifs,dl_x_y_ye<Ty,Tx>& dl);
+  };
+
+  template <typename Ty,typename Tx>
+  std::istream& operator>>(std::istream& ifs,dl_x_y_ye<Ty,Tx>& dl)
+  {
+    dl.load_from(ifs);
+    return ifs;
+  }
+
+
+  template <typename Ty,typename Tx>
+  class dl_x_xe_y_ye
+  {
+  private:
+    default_data_set<Ty,Tx> ds;
+  public:
+    data_set<Ty,Tx>& get_data_set()
+    {
+      return ds;
+    }
+
+    void load_from(std::istream& ifs)
+    {
+      for(;;)
+	{
+	  Tx x;
+	  Tx x_err;
+	  Ty y;
+	  Ty y_err(1);
+	  
+	  ifs>>x>>x_err>>y>>y_err;
+	  
+	  if(!ifs.good())
+	    {
+	      break;
+	    }
+	  data<Ty,Tx> d(x,y,y_err,y_err,x_err,x_err);
+	  ds.push_back(d);
+	}
+    }
+
+    void load_from(const char* name)
+    {
+      std::ifstream ifs(name);
+      load_from(ifs); 
+    }
+  };
+
+  template <typename Ty,typename Tx>
+  std::istream& operator>>(std::istream& ifs,dl_x_xe_y_ye<Ty,Tx>& dl)
+  {
+    dl.load_from(ifs);
+    return ifs;
+  }
+
+  template <typename Ty,typename Tx>
+  class dl_x_xu_xl_y_yu_yl
+  {
+  private:
+  
+    default_data_set<Ty,Tx> ds;
+  public:
+    data_set<Ty,Tx>& get_data_set()
+    {
+      return ds;
+    }
+    
+    void load_from(std::istream& ifs)
+    {
+      for(;;)
+	{
+	  Tx x;
+	  Tx xl,xu;
+	  Ty y;
+	  Ty yl(1),yu(1);
+	  
+	  ifs>>x>>xu>>xl>>y>>yu>>yl;
+	  
+	  xu=std::abs(xu);
+	  xl=std::abs(xl);
+	  yu=std::abs(yu);
+	  yl=std::abs(yl);
+	  
+	  if(!ifs.good())
+	    {
+	      break;
+	    }
+	  data<Ty,Tx> d(x,y,yl,yu,xl,xu);
+	  ds.push_back(d);
+	}
+    }
+
+    void load_from(const char* name)
+    {
+      std::ifstream ifs(name);
+      load_from(ifs);
+    }
+  };
+
+  template <typename Ty,typename Tx>
+  std::istream& operator>> (std::istream& ifs,dl_x_xu_xl_y_yu_yl<Ty,Tx>& dl)
+  {
+    dl.load_from(ifs);
+    return ifs;
+  }
+
+  
+
+}
+
+
+#endif
+//EOF