Changeset 3828 in Sophya for trunk/Cosmo


Ignore:
Timestamp:
Aug 3, 2010, 2:07:19 PM (15 years ago)
Author:
cmv
Message:

print pas utile, cmv 03/08/2010

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Cosmo/SimLSS/cmvginit3d.cc

    r3821 r3828  
    369369 cout<<"For Histo2DErr: d="<<dktmin<<","<<dkzmin
    370370     <<" max="<<ktnyqmax<<","<<kznyqmax<<" n="<<nherrt<<","<<nherrz<<endl;
    371 
    372  //-----------------------------------------------------------------
    373  cout<<"\n--- Computing spectra variance up to Kmax at z="<<pkz->GetZ()<<endl;
    374  // En fait on travaille sur un cube inscrit dans la sphere de rayon kmax:
    375  // sphere: Vs = 4Pi/3 k^3 , cube inscrit (cote k*sqrt(2)): Vc = (k*sqrt(2))^3
    376  // Vc/Vs = 0.675   ->  keff = kmax * (0.675)^(1/3) = kmax * 0.877
    377  double knyqmax_mod = 0.877*knyqmax;
    378  ldlkint = (log10(knyqmax_mod)-log10(k1int))/npt;
    379  varpk_int.SetInteg(0.01,ldlkint,-1.,4);
    380  double sr2int_kmax = varpk_int.Variance(k1int,knyqmax_mod);
    381  cout<<"varpk_int(<"<<knyqmax_mod<<")="<<sr2int_kmax<<"  ->  sigma="<<sqrt(sr2int_kmax)<<endl;
    382 
    383371 PrtTim(">>>> End Initialisation of GeneFluct3D");
    384372
Note: See TracChangeset for help on using the changeset viewer.