Changeset 3828 in Sophya

Timestamp:

Aug 3, 2010, 2:07:19 PM (15 years ago)

Author:

cmv

Message:

print pas utile, cmv 03/08/2010

File:

• trunk/Cosmo/SimLSS/cmvginit3d.cc (modified) (1 diff)

trunk/Cosmo/SimLSS/cmvginit3d.cc

@@ -369 +369 @@
  cout<<"For Histo2DErr: d="<<dktmin<<","<<dkzmin
      <<" max="<<ktnyqmax<<","<<kznyqmax<<" n="<<nherrt<<","<<nherrz<<endl;
-
- //-----------------------------------------------------------------
- cout<<"\n--- Computing spectra variance up to Kmax at z="<<pkz->GetZ()<<endl;
- // En fait on travaille sur un cube inscrit dans la sphere de rayon kmax:
- // sphere: Vs = 4Pi/3 k^3 , cube inscrit (cote k*sqrt(2)): Vc = (k*sqrt(2))^3
- // Vc/Vs = 0.675   ->  keff = kmax * (0.675)^(1/3) = kmax * 0.877
- double knyqmax_mod = 0.877*knyqmax;
- ldlkint = (log10(knyqmax_mod)-log10(k1int))/npt;
- varpk_int.SetInteg(0.01,ldlkint,-1.,4);
- double sr2int_kmax = varpk_int.Variance(k1int,knyqmax_mod);
- cout<<"varpk_int(<"<<knyqmax_mod<<")="<<sr2int_kmax<<"  ->  sigma="<<sqrt(sr2int_kmax)<<endl;
-
  PrtTim(">>>> End Initialisation of GeneFluct3D");