Changeset 3524 in Sophya for trunk/Cosmo/SimLSS/genefluct3d.cc

Timestamp:

Sep 22, 2008, 4:07:02 PM (17 years ago)

Author:

cmv

Message:

intro du remplissage du NTuple de debug cmv 22/09/2008

File:

• trunk/Cosmo/SimLSS/genefluct3d.cc (modified) (8 diffs)

trunk/Cosmo/SimLSS/genefluct3d.cc

@@ -15 +15 @@
 #include "fabtwriter.h"
 #include "fioarr.h"
+#include "ntuple.h"
 
 #include "arrctcast.h"
@@ -83 +84 @@
  nthread_ = 0;
  lp_ = 0;
- array_allocated_ = false;
+ array_allocated_ = false; array_type = 0;
  cosmo_ = NULL;
  growth_ = NULL;
@@ -153 +154 @@
  try {
    T_.ReSize(3,SzK_);
-   array_allocated_ = true;
+   array_allocated_ = true; array_type=0;
    if(lp_>1) cout<<"  allocating: "<<T_.Size()*sizeof(complex<GEN3D_TYPE>)/1.e6<<" Mo"<<endl;
  } catch (...) {
@@ -552 +553 @@
 
 //-------------------------------------------------------
+void GeneFluct3D::NTupleCheck(POutPersist &pos,string ntname,unsigned long nent)
+// Remplit le NTuple "ntname" avec "nent" valeurs du cube (reel ou complex) et l'ecrit dans "pos"
+{
+  if(ntname.size()<=0 || nent==0) return;
+  int nvar = 0;
+  if(array_type==1) nvar = 3;
+  else if(array_type==2) nvar = 4;
+  else return;
+  char *vname[4] = {"t","z","re","im"};
+  float xnt[4];
+  NTuple nt(nvar,vname);
+
+  if(array_type==1) {
+    unsigned long nmod = Nx_*Ny_*Nz_/nent; if(nmod==0) nmod=1;
+    unsigned long n=0;
+    for(long i=0;i<Nx_;i++) for(long j=0;j<Ny_;j++) for(long l=0;l<Nz_;l++) {
+      if(n==nmod) {
+        int_8 ip = IndexR(i,j,l);
+        xnt[0]=sqrt(i*i+j*j); xnt[1]=l; xnt[2]=data_[ip];
+        nt.Fill(xnt);
+        n=0;
+      }
+      n++;
+    }
+  } else {
+    unsigned long nmod = Nx_*Ny_*NCz_/nent; if(nmod==0) nmod=1;
+    unsigned long n=0;
+    for(long i=0;i<Nx_;i++) for(long j=0;j<Ny_;j++) for(long l=0;l<NCz_;l++) {
+      if(n==nmod) {
+        xnt[0]=sqrt(i*i+j*j); xnt[1]=l; xnt[2]=T_(l,j,i).real(); xnt[3]=T_(l,j,i).imag();
+        nt.Fill(xnt);
+        n=0;
+      }
+      n++;
+    }
+  }
+
+  pos.PutObject(nt,ntname);
+}
+
+//-------------------------------------------------------
 void GeneFluct3D::Print(void)
 {
@@ -611 +653 @@
  }
 
+ array_type = 2;
+
  if(lp_>0) cout<<"...computing power"<<endl;
  double p = compute_power_carte();
@@ -654 +698 @@
  }
 
+ array_type = 2;
  manage_coefficients();   // gros effet pour les spectres que l'on utilise !
 
@@ -851 +896 @@
  // On fait la FFT
  GEN3D_FFTW_EXECUTE(pb_);
+ array_type = 1;
 }
 
@@ -861 +907 @@
  // On fait la FFT
  GEN3D_FFTW_EXECUTE(pf_);
+ array_type = 2;
+
  // On corrige du pb de la normalisation de FFTW3
  double v = (double)NRtot_;