Changeset 857 in Sophya for trunk/SophyaProg/Tests/tspm2.cc
- Timestamp:
- Apr 10, 2000, 3:04:05 PM (25 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaProg/Tests/tspm2.cc
r819 r857 1 1 #include <iostream.h> 2 2 #include <math.h> 3 #include "skymapinit.h" 4 #include "spheregorski.h" 5 #include "spherethetaphi.h" 3 #include "sambainit.h" 4 #include "skymap.h" 6 5 7 6 #include "tod.h" … … 87 86 InitTim(); // Initializing the CPU timer 88 87 if ((narg > 1) && (strcmp(arg[1],"-h") == 0) ) { 89 cout << " tspm [ Gorski_M=32] : GorskiSpherical Map Test " << endl;88 cout << " tspm [HEALPix_M=32] : HEALPix Spherical Map Test " << endl; 90 89 exit(0); 91 90 } … … 93 92 int m=32; 94 93 if (narg >1) m = atoi(arg[1]); 95 cout << " ===== GorskiSpherical Map Test M= " << m << endl;94 cout << " ===== HEALPix Spherical Map Test M= " << m << endl; 96 95 97 96 POutPersist s("spheres.ppf"); … … 99 98 100 99 { 101 Sphere Gorski<double> sph(m);100 SphereHEALPix<double> sph(m); 102 101 103 102 cout << "Filling spherical map NPixels= " << sph.NbPixels() << endl; … … 110 109 111 110 { 112 FIO_Sphere Gorski<double> fiog(&sph) ;111 FIO_SphereHEALPix<double> fiog(&sph) ; 113 112 nomobj = "sphg1"; 114 113 fiog.Write(s, nomobj); 115 cout << "SphMap Sphere Gorski<double> written to POutPersist with name " << nomobj << endl;114 cout << "SphMap SphereHEALPix<double> written to POutPersist with name " << nomobj << endl; 116 115 } 117 116 … … 149 148 150 149 { 151 Sphere Gorski<float> sph(m);150 SphereHEALPix<float> sph(m); 152 151 153 152 cout << "Filling spherical map2 NPixels= " << sph.NbPixels() << endl; … … 166 165 167 166 { 168 FIO_Sphere Gorski<float> fiog(&sph) ;167 FIO_SphereHEALPix<float> fiog(&sph) ; 169 168 nomobj = "sphg2"; 170 169 fiog.Write(s, nomobj); 171 cout << "SphMap Sphere Gorski<float> written to POutPersist with name " << nomobj << endl;170 cout << "SphMap SphereHEALPix<float> written to POutPersist with name " << nomobj << endl; 172 171 } 173 172 … … 176 175 FitsIoServer fios; 177 176 { 178 cout << "Test of Write/Read Sphere Gorski<float> to FITS (sphg_r4.fits) " << endl;177 cout << "Test of Write/Read SphereHEALPix<float> to FITS (sphg_r4.fits) " << endl; 179 178 fios.save(sph, "sphg_r4.fits"); 180 179 181 Sphere Gorski<float> sphr(4);180 SphereHEALPix<float> sphr(4); 182 181 fios.load(sphr, "sphg_r4.fits", 2); 183 cout << " Read from file - Sphere Gorski<float> NPixels= " << sphr.NbPixels() << endl;182 cout << " Read from file - SphereHEALPix<float> NPixels= " << sphr.NbPixels() << endl; 184 183 int ndiff = 0; 185 184 for(int k=0; k<sphr.NbPixels(); k++) {
Note:
See TracChangeset
for help on using the changeset viewer.