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

Timestamp:

Jul 24, 2010, 6:15:07 PM (15 years ago)

Author:

cmv

Message:

suite de la mise au point pour lecture fichiers CAMB, cmv 24/07/2010

File:

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

trunk/Cosmo/SimLSS/genefluct3d.cc

@@ -978 +978 @@
 //             (tous les pixels d'un plan Z sont mis au meme redshift z que celui du milieu)
 {
- if(lp_>0) cout<<"--- ApplyGrowthFactor: evol="<<type_evol<<endl;
+  if(lp_>0) cout<<"--- ApplyGrowthFactor: evol="<<type_evol<<" (zpk="<<compute_pk_redsh_ref_<<")"<<endl;
  check_array_alloc();
 
@@ -1035 +1035 @@
 
  double zpk = compute_pk_redsh_ref_;
- double dpsd_orig = - cosmo_->H(zpk) * (1.+zpk) * growth_->DsDz(zpk,good_dzinc_) / (*growth_)(zpk);
+ double grw_orig = (*growth_)(zpk);
+ double dpsd_orig = - cosmo_->H(zpk) * (1.+zpk) * growth_->DsDz(zpk,good_dzinc_) / grw_orig;
+ if(lp_>0) cout<<"    original growth="<<grw_orig<<" dpsd="<<dpsd_orig<<" computed at z="<<zpk<<endl;
 
  InterpFunc interpinv(loscom2zred_min_,loscom2zred_max_,loscom2zred_);
@@ -1049 +1051 @@
          else dz = fabs(dz); // tous les plans Z au meme redshift
        double z = interpinv(dz);   // interpolation par morceau
+       double grw = (*growth_)(z);
        double dpsd = interdpd(z);
        int_8 ip = IndexR(i,j,l);
-       data_[ip] *= dpsd / dpsd_orig;
+       // on remet le beta au bon z
+       // on corrige du growth factor car data_ a ete calcule avec pk(zpk)
+       data_[ip] *= (dpsd / dpsd_orig) * (grw / grw_orig);
      }
    }
@@ -1062 +1067 @@
 // Calcule une realisation dans l'espace reel
 {
- if(lp_>0) cout<<"--- ComputeReal ---"<<endl;
- check_array_alloc();
-
- // On fait la FFT
- GEN3D_FFTW_EXECUTE(pb_);
- array_type = 1;
+  if(lp_>0) cout<<"--- ComputeReal --- from spectrum at z="<<compute_pk_redsh_ref_<<endl;
+  check_array_alloc();
+
+  // On fait la FFT
+  GEN3D_FFTW_EXECUTE(pb_);
+  array_type = 1;
 }