Changeset 3805 in Sophya for trunk/Cosmo/SimLSS/cmvginit3d.cc
- Timestamp:
- Jul 24, 2010, 12:01:34 AM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/Cosmo/SimLSS/cmvginit3d.cc
r3800 r3805 218 218 cout<<endl<<"\n--- Create Spectrum"<<endl; 219 219 220 Initial Spectrumpkini(ns,as);220 InitialPowerLaw pkini(ns,as); 221 221 222 222 TransfertEisenstein tf(h100,om0-ob0,ob0,T_CMB_Par,false); … … 224 224 tfnos.SetNoOscEnv(2); 225 225 226 Growth Factorgrowth(om0,ol0);227 // Growth Factorgrowth(1.,0.); // D(z) = 1/(1+z)226 GrowthEisenstein growth(om0,ol0); 227 // GrowthEisenstein growth(1.,0.); // D(z) = 1/(1+z) 228 228 double growth_at_z = growth(zref); 229 229 cout<<"...Growth factor at z="<<zref<<" = "<<growth_at_z<<endl; 230 230 231 PkSpectrum0 pk0(pkini,tf); 232 PkSpectrum0 pk0nos(pkini,tfnos); 233 234 PkSpectrumZ pkz(pk0,growth,zref); 235 PkSpectrumZ pkznos(pk0nos,growth,zref); 231 PkSpecCalc pkz(pkini,tf,growth,zref); 232 PkSpecCalc pkznos(pkini,tfnos,growth,zref); 236 233 237 234 //-----------------------------------------------------------------
Note:
See TracChangeset
for help on using the changeset viewer.