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

Timestamp:

Sep 11, 2008, 12:37:39 PM (17 years ago)

Author:

cmv

Message:

possibilite de travailler en float cmv 11/09/2008

File:

• trunk/Cosmo/SimLSS/genefluct3d.h (modified) (4 diffs)

trunk/Cosmo/SimLSS/genefluct3d.h

@@ -18 +18 @@
 #include "pkspectrum.h"
 
+#define WITH_FFTW_THREAD
+//#define GEN3D_FLOAT
+
+#if defined(GEN3D_FLOAT)
+#define GEN3D_TYPE r_4
+#define GEN3D_FFTW_PLAN fftwf_plan
+#define GEN3D_FFTW_COMPLEX fftwf_complex
+#else
+#define GEN3D_TYPE r_8
+#define GEN3D_FFTW_PLAN fftw_plan
+#define GEN3D_FFTW_COMPLEX fftw_complex
+#endif
 
 namespace SOPHYA {
@@ -41 +53 @@
   long LosComRedshift(double zinc=0.001,long npoints=-1);
 
-  TArray< complex<r_8> >& GetComplexArray(void) {return T_;}
-  fftw_complex* GetComplexPointer(void) {return fdata_;}
-  TArray<r_8>& GetRealArray(void) {return R_;}
-  r_8* GetRealPointer(void) {return data_;}
+  TArray< complex<GEN3D_TYPE> >& GetComplexArray(void) {return T_;}
+  GEN3D_FFTW_COMPLEX * GetComplexPointer(void) {return fdata_;}
+  TArray<GEN3D_TYPE>& GetRealArray(void) {return R_;}
+  GEN3D_TYPE* GetRealPointer(void) {return data_;}
 
   // Pour adressage data_[ip]
@@ -157 +169 @@
 
   // la gestion de la FFT
-  bool is_set_fftw_plan;
-  fftw_plan pf_,pb_;
+  bool is_set_fft_plan;
+  GEN3D_FFTW_PLAN pf_,pb_;
   unsigned short nthread_;
   int lp_;
@@ -164 +176 @@
   // le stockage du Cube de donnees et les pointeurs
   bool array_allocated_;  // true if array has been allocated
-  TArray< complex<r_8> > T_;
-  fftw_complex *fdata_;
-  TArray<r_8> R_;
-  double *data_;
+  TArray< complex<GEN3D_TYPE> > T_;
+  GEN3D_FFTW_COMPLEX *fdata_;
+  TArray<GEN3D_TYPE> R_;
+  GEN3D_TYPE *data_;
 
   // l'observateur