Changeset 880 in Sophya for trunk/SophyaProg/PMixer/extractRS.cc
- Timestamp:
- Apr 11, 2000, 3:34:55 PM (25 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaProg/PMixer/extractRS.cc
r877 r880 12 12 #include "tvector.h" 13 13 #include "vector3d.h" 14 #include " nbrandom.h"14 #include "srandgen.h" 15 15 #include "timing.h" 16 16 #include "sambainit.h" … … 88 88 89 89 90 cout << " extractRS/Info : NComp = " << nskycomp << " Sphere Gorski_NSide= " << hp_nside << endl;90 cout << " extractRS/Info : NComp = " << nskycomp << " SphereHEALPix_NSide= " << hp_nside << endl; 91 91 cout << " ... MapPath = " << (string)mapPath << " DebugLev= " << debuglev 92 92 << " PrintLev= " << printlev << endl; … … 95 95 // We create an output persist file for writing debug objects 96 96 if (debuglev > 0) so = new POutPersist("extrRSdbg.ppf"); 97 Sphere Gorski<float> SourceMap(hp_nside);98 Sphere Gorski<float> OutputMap(hp_nside);97 SphereHEALPix<float> SourceMap(hp_nside); 98 SphereHEALPix<float> OutputMap(hp_nside); 99 99 bool okout = true; 100 100 … … 190 190 K = dc.DParam(key, 1, 1.); 191 191 if (debuglev > 4) { // Writing the input map to the outppf 192 FIO_T Matrix<float> fiog;192 FIO_TArray<float> fiog; 193 193 fiog.Write(*so, key); 194 194 } … … 282 282 // MeanSig(SourceMap.DataBlock(),moy,sig); 283 283 // cout << " MeanSig for Source Map - Mean= " << moy << " Sigma= " << sig << endl; 284 // cout << "Source Map (Sphere Gorski<float>) written to FITS file "284 // cout << "Source Map (SphereHEALPix<float>) written to FITS file " 285 285 // << (string)(arg[3]) << endl; 286 286 // } … … 291 291 MeanSig(OutputMap.DataBlock(),moy,sig); 292 292 cout << " MeanSig for Output Map - Mean= " << moy << " Sigma= " << sig << endl; 293 cout << "Output Map (Sphere Gorski<float>) written to FITS file "293 cout << "Output Map (SphereHEALPix<float>) written to FITS file " 294 294 << (string)(arg[2]) << endl; 295 295 } … … 298 298 if (narg > 2) { 299 299 POutPersist s(arg[3]); 300 FIO_Sphere Gorski<float> fiog(OutputMap);300 FIO_SphereHEALPix<float> fiog(OutputMap); 301 301 fiog.Write(s); 302 cout << "Output Map (Sphere Gorski<float>) written to POutPersist file "302 cout << "Output Map (SphereHEALPix<float>) written to POutPersist file " 303 303 << (string)(arg[3]) << endl; 304 304 if(printlev>2) PrtTim("End of WritePPF ");
Note:
See TracChangeset
for help on using the changeset viewer.