Changeset 3330 in Sophya for trunk/Cosmo/SimLSS/genefluct3d.h

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/genefluct3d.h (modified) (3 diffs)

trunk/Cosmo/SimLSS/genefluct3d.h

@@ -51 +51 @@
   // Pour adressage fdata_[ip][0-1]
   inline int_8 IndexC(long i,long j,long k) {return (int_8)(k+NCz_*(j+Ny_*i));}
-  // - On peut aussi adresser:
-  // TArray< complex<r_8> >& pk = gf3d.GetComplexArray(); pk(i,j,k) = ...;
-  // TArray<r_8>& rgen = gf3d.GetRealArray(); rgen(i,j,k) = ...;
+  // On peut aussi adresser:
+  // TArray< complex<r_8> >& pk = gf3d.GetComplexArray();
+  //    pk(k,j,i) avec k=[0,NCz_[  j=[0,Ny_[   i=[0,Nx_[
+  //    pk[IndexC(i,j,k)]
+  // TArray<r_8>& rgen = gf3d.GetRealArray();
+  //    rgen(k,j,i) avec k=[0,NTz_[  j=[0,Ny_[   i=[0,Nx_[
+  //                mais seul k=[0,Nz_[ est utile
+  //    rgen[IndexR(i,j,k)]
+  // ATTENTION: TArray adresse en memoire a l'envers du C !
+  //            Tarray(n1,n2,n3) == Carray[n3][n2][n1]
 
   vector<long> GetNpix(void) {return N_;}
   int_8 NPix(void) {return NRtot_;}
+  long GetNx(void) {return Nx_;}
+  long GetNy(void) {return Ny_;}
+  long GetNz(void) {return Nz_;}
 
   // Return |K_i| module relative to pixel indices
@@ -83 +93 @@
 
   void ComputeReal(void);
+  void ApplyGrowthFactor(void);
 
   void ReComputeFourier(void);
@@ -88 +99 @@
   int  ComputeSpectrum(HistoErr& herr);
   int  ComputeSpectrum2D(Histo2DErr& herr);
-  void ApplyGrowthFactor(void);
+  int  ComputeSpectrum(HistoErr& herr,double sigma,bool pixcor);
+  int  ComputeSpectrum2D(Histo2DErr& herr,double sigma,bool pixcor);
 
   int_8 VarianceFrReal(double R,double& var);