Changeset 3157 in Sophya for trunk/Cosmo/SimLSS/geneutils.h

Timestamp:

Jan 25, 2007, 6:04:48 PM (19 years ago)

Author:

cmv

Message:

intro du facteur de croissance dans la simul cmv 25/01/2007

File:

• trunk/Cosmo/SimLSS/geneutils.h (modified) (5 diffs)

trunk/Cosmo/SimLSS/geneutils.h

@@ -7 +7 @@
 #include "histos.h"
 #include "tvector.h"
+#include "cspline.h"
 
 #include <vector>
@@ -12 +13 @@
 namespace SOPHYA {
 
+//----------------------------------------------------
 class InterpFunc {
 public:
@@ -19 +21 @@
   double XMin(void) {return _xmin;}
   double XMax(void) {return _xmax;}
+  inline double X(long i) {return _xmin + i*_dx;}
 
   //! Retourne l'element le plus proche de f donnant y=f(x)
@@ -24 +27 @@
          {
          x -= _xmin;
-         long i = long(x/_xlarg*_nm1+0.5);  // On prend le "i" le plus proche
+         long i = long(x/_dx+0.5);  // On prend le "i" le plus proche
          if(i<0) i=0; else if(i>=_nm1) i=_nm1-1;
          return _y[i];
@@ -43 +46 @@
 
 protected:
-  double _xmin,_xmax,_xlarg,_dx;
+  double _xmin,_xmax,_dx;
+  long _nm1;  // n-1
   vector<double>& _y;
-  long _nm1;  // n-1
+};
+
+//----------------------------------------------------
+class InverseFunc {
+public:
+  InverseFunc(vector<double>& x,vector<double>& y);
+  virtual ~InverseFunc(void);
+  int ComputeLinear(long n,vector<double>& xfcty);
+  int ComputeParab(long n,vector<double>& xfcty);
+  double YMin(void) {return _ymin;}
+  double YMax(void) {return _ymax;}
+protected:
+  inline void find_in_y(double x,long& klo,long& khi)
+    {
+      long k;
+      klo=0, khi=_y.size()-1;
+      while (khi-klo > 1) {
+        k = (khi+klo) >> 1;
+        if (_y[k] > x) khi=k; else klo=k;
+      }
+    }
+
+  vector<double>& _x;
+  vector<double>& _y;
+  double _ymin,_ymax,_dy;
 };
 
 }
 
+//----------------------------------------------------
 int FuncToHisto(GenericFunc& func,Histo& h,bool logaxex=false);
 int FuncToVec(GenericFunc& func,TVector<r_8>& h,double xmin,double xmax,bool logaxex=false);
 
+//----------------------------------------------------
 double AngSol(double dtheta,double dphi,double theta0=M_PI/2.);
 double AngSol(double dtheta);
 
+//----------------------------------------------------
 unsigned long PoissRandLimit(double mu,double mumax=10.);