Changeset 3330 in Sophya for trunk/Cosmo/SimLSS/geneutils.cc

Timestamp:

Oct 1, 2007, 7:10:50 PM (18 years ago)

Author:

cmv

Message:

intro ComputeSpectrum avec soustraction de bruit et deconvolution par la fct de transfert du pixel cmv 01/10/2007

File:

• trunk/Cosmo/SimLSS/geneutils.cc (modified) (3 diffs)

trunk/Cosmo/SimLSS/geneutils.cc

@@ -55 +55 @@
 //-------------------------------------------------------------------
 // Classe d'inversion d'une fonction STRICTEMENT MONOTONE CROISSANTE
-// - Le vecteur y a "Nin" elements y_i tels que "y_i = f(x_i)"
-// - On a x(i) < x(i+1) et y(i) < y(i+1)
-// - La classe renvoie ymin=y(0) , ymax=y(Nin -1)
-//   et le vecteur x = f^-1(y) de "Nout" elements
-//   Les y_i sont regulierement espaces et ymin et ymax
-//   La re-interpolation inverse est faite par lineaire
+//
+// - On part de deux vecteurs x,y de "Nin" elements tels que "y_i = f[x_i]"
+//   ou la fonction "f" est strictement monotone croissante:
+//          x(i) < x(i+1)   et   y(i) < y(i+1)
+// - Le but de la classe est de remplir un vecteur X de "Nout" elements
+//   tels que:    X_j = f^-1[Y_j]    avec   j=[0,Nout[
+//   avec les Y_j regulierement espaces entre ymin=y(0) , ymax=y(Nin -1)
+//   cad:     X_j = f^-1[ ymin+j*(ymax-ymin)/(Nout-1) ]
+// - La construction du vecteur X est realisee
+//   par interpolation lineaire (ComputeLinear) ou parabolique (ComputeParab)
 InverseFunc::InverseFunc(vector<double>& x,vector<double>& y)
   : _ymin(0.) , _ymax(0.) , _x(x) , _y(y)
@@ -92 +96 @@
 }
 
-int InverseFunc::ComputeLinear(long n,vector<double>& xfcty)
+int InverseFunc::ComputeLinear(long nout,vector<double>& xfcty)
 // Compute table "xfcty" by linear interpolation of "x" versus "y"
-//   on "n" points from "ymin" to "ymax":
-// xfcty[i] = interpolation of function "x" for "ymin+i*(ymax-ymin)/(n-1.)"
-{
-  if(n<3) return -1;
-
-  xfcty.resize(n);
+//   on "nout" points from "ymin" to "ymax":
+// xfcty[i] = interpolation of function "x" for "ymin+i*(ymax-ymin)/(nout-1)"
+{
+  if(nout<3) return -1;
+
+  xfcty.resize(nout);
 
   long i1,i2;
   double x;
-  for(int_4 i=0;i<n;i++) {
-    double y = _ymin + i*(_ymax-_ymin)/(n-1.);
+  for(int_4 i=0;i<nout;i++) {
+    double y = _ymin + i*(_ymax-_ymin)/(nout-1.);
     find_in_y(y,i1,i2);
     double dy = _y[i2]-_y[i1];
@@ -118 +122 @@
 }
 
-int InverseFunc::ComputeParab(long n,vector<double>& xfcty)
-{
-  if(n<3) return -1;
-
-  xfcty.resize(n);
+int InverseFunc::ComputeParab(long nout,vector<double>& xfcty)
+{
+  if(nout<3) return -1;
+
+  xfcty.resize(nout);
 
   long i1,i2,i3;
   double x;
-  for(int_4 i=0;i<n;i++) {
-    double y = _ymin + i*(_ymax-_ymin)/(n-1.);
+  for(int_4 i=0;i<nout;i++) {
+    double y = _ymin + i*(_ymax-_ymin)/(nout-1.);
     find_in_y(y,i1,i2);
     // On cherche le 3ieme point selon la position de y / au 2 premiers