Changeset 3155 in Sophya for trunk/Cosmo/SimLSS/cmvobserv3d.cc

Timestamp:

Jan 19, 2007, 6:18:38 PM (19 years ago)

Author:

cmv

Message:

PrtTim deplaces de genefluct3d.cc -> cmvobserv3d.cc cmv 19/01/2007

File:

• trunk/Cosmo/SimLSS/cmvobserv3d.cc (modified) (22 diffs)

trunk/Cosmo/SimLSS/cmvobserv3d.cc

@@ -135 +135 @@
      <<", npt="<<npt<<endl;
  cout<<"R="<<R<<" Rg="<<Rg<<" Mpc, sigmaR="<<sigmaR<<endl;
- cout<<"nstar= "<<nstar<<"  mstar="<<mstar<<"  alpha="<<alpha<<endl;
+ cout<<"nstar= "<<nstar<<"  mstar="<<mstar<<"  alpa="<<alpha<<endl;
  cout<<"schmin="<<schmin<<" ("<<lschmin
      <<"), schmax="<<schmax<<" ("<<lschmax<<") Msol"
@@ -204 +204 @@
  cout<<"varpk_int="<<sr2int<<"  ->  sigma="<<sqrt(sr2int)<<endl;
 
- PrtTim("End of definition");
-
  //-----------------------------------------------------------------
  cout<<endl<<"\n--- Compute galaxy density number"<<endl;
@@ -219 +217 @@
  cout<<"Galaxy mass density= "<<mass_by_mpc3<<" Msol/Mpc^3 between limits"<<endl;
 
+ PrtTim(">>>> End of definition");
+
  //-----------------------------------------------------------------
  // FFTW3 (p26): faster if sizes 2^a 3^b 5^c 7^d 11^e 13^f  with e+f=0 ou 1
- cout<<endl<<"\n--- Compute Spectrum Realization"<<endl;
+ cout<<endl<<"\n--- Initialisation de GeneFluct3D"<<endl;
 
  pkz.SetZ(zref);
  TArray< complex<r_8> > pkgen;
  GeneFluct3D fluct3d(pkgen);
- fluct3d.SetPrtLevel(3);
+ fluct3d.SetPrtLevel(2);
  fluct3d.SetNThread(nthread);
  fluct3d.SetSize(nx,ny,nz,dx,dy,dz);
@@ -256 +256 @@
  cout<<"varpk_int(<"<<knyqmax_mod<<")="<<sr2int_kmax<<"  ->  sigma="<<sqrt(sr2int_kmax)<<endl;
 
+ PrtTim(">>>> End Initialisation de GeneFluct3D");
+
  cout<<"\n--- Computing a realization in Fourier space"<<endl;
  if(computefourier0) fluct3d.ComputeFourier0(pkz);
    else              fluct3d.ComputeFourier(pkz);
+ PrtTim(">>>> End Computing a realization in Fourier space");
 
  if(1) {
@@ -269 +272 @@
    tagobs = "hpkgen"; posobs.PutObject(hpkgen,tagobs);
    }
+   PrtTim(">>>> End Checking realization spectra");
  }
 
@@ -280 +284 @@
    tagobs = "hpkgen2"; posobs.PutObject(hpkgen2,tagobs);
    }
+   PrtTim(">>>> End Checking realization 2D spectra");
  }
 
@@ -285 +290 @@
    cout<<"\n--- Computing convolution by pixel shape"<<endl;
    fluct3d.FilterByPixel();
- }
-
- if(wfits) fluct3d.WriteFits("!cmvobserv3d_k0.fits");
- if(wppf) fluct3d.WritePPF("cmvobserv3d_k0.ppf",false);
+   PrtTim(">>>> End Computing convolution by pixel shape");
+ }
+
+ if(wfits) {
+   fluct3d.WriteFits("!cmvobserv3d_k0.fits");
+   PrtTim(">>>> End WriteFits");
+ }
+ if(wppf) {
+   fluct3d.WritePPF("cmvobserv3d_k0.ppf",false);
+   PrtTim(">>>> End WritePPF");
+ }
 
  if(1) {
@@ -299 +311 @@
    tagobs = "hpkgenf"; posobs.PutObject(hpkgenf,tagobs);
    }
+   PrtTim(">>>> End Checking realization spectra");
  }
 
@@ -310 +323 @@
    tagobs = "hpkgenf2"; posobs.PutObject(hpkgenf2,tagobs);
    }
+   PrtTim(">>>> End Checking realization 2D spectra");
  }
 
@@ -316 +330 @@
  double rmin,rmax; rgen.MinMax(rmin,rmax);
  cout<<"rgen.Min = "<<rmin<<" , Max="<<rmax<<endl;
-
- if(wfits) fluct3d.WriteFits("!cmvobserv3d_r0.fits");
- if(wppf) fluct3d.WritePPF("cmvobserv3d_r0.ppf",true);
+   PrtTim(">>>> End Computing a realization in real space");
+
+ if(wfits) {
+   fluct3d.WriteFits("!cmvobserv3d_r0.fits");
+   PrtTim(">>>> End WriteFits");
+ }
+ if(wppf) {
+   fluct3d.WritePPF("cmvobserv3d_r0.ppf",true);
+   PrtTim(">>>> End WritePPF");
+ }
 
  int_8 nm;
@@ -328 +349 @@
    cout<<"rgen:("<<nm<<") Mean = "<<rm<<", Sigma^2 = "
        <<rs2<<" -> "<<sqrt(rs2)<<",  n(<-1)="<<nlowone<<endl;
+   PrtTim(">>>> End Check mean and variance in real space");
  }
 
@@ -335 +357 @@
    int_8 nvarr = fluct3d.VarianceFrReal(R,varr);
    cout<<"R="<<R<<" : sigmaR^2="<<varr<<" -> "<<sqrt(varr)<<",   n="<<nvarr<<endl;
+   PrtTim(">>>> End Check variance sigmaR in real space");
  }
 
@@ -343 +366 @@
  cout<<"1+rgen: ("<<nm<<") Mean = "<<rm<<", Sigma^2 = "
      <<rs2<<" -> "<<sqrt(rs2)<<endl;
+ PrtTim(">>>> End Converting fluctuations into mass");
 
  cout<<"\n--- Converting mass into galaxy number"<<endl;
@@ -350 +374 @@
  cout<<"galaxy: ("<<nm<<") Mean = "<<rm<<", Sigma^2 = "
      <<rs2<<" -> "<<sqrt(rs2)<<endl;
+ PrtTim(">>>> End Converting mass into galaxy number");
 
  cout<<"\n--- Set negative pixels to BAD"<<endl;
@@ -357 +382 @@
  cout<<"galaxy: ("<<nm<<") Mean = "<<rm<<", Sigma^2 = "
      <<rs2<<" -> "<<sqrt(rs2)<<endl;
+ PrtTim(">>>> End Set negative pixels to BAD etc...");
 
  cout<<"\n--- Apply poisson on galaxy number"<<endl;
@@ -364 +390 @@
  cout<<"galaxy: ("<<nm<<") Mean = "<<rm<<", Sigma^2 = "
      <<rs2<<" -> "<<sqrt(rs2)<<endl;
+ PrtTim(">>>> End Apply poisson on galaxy number");
 
  cout<<"\n--- Convert Galaxy number to HI mass"<<endl;
@@ -382 +409 @@
      <<rs2<<" -> "<<sqrt(rs2)<<endl;
  cout<<"Equivalent: "<<rm*nm/fluct3d.NPix()<<" Msol / pixels"<<endl;
+ PrtTim(">>>> End Convert Galaxy number to HI mass");
 
  cout<<"\n--- Set BAD pixels to Zero"<<endl;
@@ -389 +417 @@
  cout<<"galaxy: ("<<nm<<") Mean = "<<rm<<", Sigma^2 = "
      <<rs2<<" -> "<<sqrt(rs2)<<endl;
-
- if(wfits) fluct3d.WriteFits("!cmvobserv3d_r.fits");
- if(wppf) fluct3d.WritePPF("cmvobserv3d_r.ppf",true);
+ PrtTim(">>>> End Set BAD pixels to Zero etc...");
+
+ if(wfits) {
+   fluct3d.WriteFits("!cmvobserv3d_r.fits");
+   PrtTim(">>>> End WriteFits");
+ }
+ if(wppf) {
+   fluct3d.WritePPF("cmvobserv3d_r.ppf",true);
+   PrtTim(">>>> End WritePPF");
+ }
 
  cout<<"\n--- Add noise to HI Flux snoise="<<snoise<<endl;
@@ -398 +433 @@
  cout<<"HI mass with noise: ("<<nm<<") Mean = "<<rm<<", Sigma^2 = "
      <<rs2<<" -> "<<sqrt(rs2)<<endl;
+ PrtTim(">>>> End Add noise");
 
  //-----------------------------------------------------------------
@@ -405 +441 @@
    cout<<endl<<"\n--- ReComputing spectrum from real space"<<endl;
    fluct3d.ReComputeFourier();
- }
-
- if(wfits) fluct3d.WriteFits("!cmvobserv3d_k.fits");
- if(wppf) fluct3d.WritePPF("cmvobserv3d_k.ppf",false);
+   PrtTim(">>>> End ReComputing spectrum");
+ }
+
+ if(wfits) {
+   fluct3d.WriteFits("!cmvobserv3d_k.fits");
+   PrtTim(">>>> End WriteFits");
+ }
+ if(wppf) {
+   fluct3d.WritePPF("cmvobserv3d_k.ppf",false);
+   PrtTim(">>>> End WritePPF");
+ }
 
  if(1) {
@@ -417 +460 @@
    fluct3d.ComputeSpectrum(hpkrec);
    tagobs = "hpkrec"; posobs.PutObject(hpkrec,tagobs);
+   PrtTim(">>>> End Computing final spectrum");
  }
 
@@ -428 +472 @@
    tagobs = "hpkrec2"; posobs.PutObject(hpkrec2,tagobs);
    }
- }
-
+   PrtTim(">>>> End Computing final 2D spectrum");
+ }
+
+ PrtTim(">>>> End Of Job");
  return 0;
 }