Changeset 857 in Sophya for trunk/SophyaProg/Tests/tspm.cc


Ignore:
Timestamp:
Apr 10, 2000, 3:04:05 PM (25 years ago)
Author:
ansari
Message:

Adapatation a SphereHEALPix, correction tsttmatrix.cc tsttvector.cc - Reza 10/4/2000

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/SophyaProg/Tests/tspm.cc

    r768 r857  
    11#include <iostream.h>
    22#include "skymapinit.h"
    3 #include "spheregorski.h"
    4 #include "spherethetaphi.h"
     3#include "skymap.h"
    54
    65#include "tod.h"
     
    3029  InitTim();   // Initializing the CPU timer
    3130  if ((narg > 1) && (strcmp(arg[1],"-h") == 0) )  {
    32     cout << " tspm [Gorski_M=32] [M_TetaPhi=64] : Spherical Map Test " << endl;
     31    cout << " tspm [HEALPix_M=32] [M_TetaPhi=64] : Spherical Map Test " << endl;
    3332    exit(0);
    3433    }
     
    3635  int  m=32;
    3736  if (narg >1)  m = atoi(arg[1]);
    38   cout <<  " ===== Gorski Spherical Map Test M= " << m << endl;
     37  cout <<  " ===== HEALPix Spherical Map Test M= " << m << endl;
    3938
    40   SphereGorski<double> sph(m);
     39  SphereHEALPix<double> sph(m);
    4140
    4241  cout << "Filling spherical map NPixels= " << sph.NbPixels() << endl;
     
    5049  // Computing mean and sigma on the sphere
    5150  MeanSig(sph, gmoy, gsig);
    52   cout << "SphereGorski<double> Mean= " << gmoy << "  Sigma = " << gsig << endl;
     51  cout << "SphereHEALPix<double> Mean= " << gmoy << "  Sigma = " << gsig << endl;
    5352  PrtTim("End of Mean-Sig ");
    5453
     
    5655  {
    5756  POutPersist s("sphg.ppf");
    58   FIO_SphereGorski<double> fiog(&sph) ;
     57  FIO_SphereHEALPix<double> fiog(&sph) ;
    5958  fiog.Write(s);
    60   cout << "SphereGorski<double>  written to sphg.ppf "  << endl;
     59  cout << "SphereHEALPix<double>  written to sphg.ppf "  << endl;
    6160  }
    6261
    6362  // Reading from the file
    6463  {
    65   FIO_SphereGorski<double> fiog("sphg.ppf");
    66   SphereGorski<double>  sph2 = fiog;
     64  FIO_SphereHEALPix<double> fiog("sphg.ppf");
     65  SphereHEALPix<double>  sph2 = fiog;
    6766  PrtTim("End of Write/Read ");
    6867
Note: See TracChangeset for help on using the changeset viewer.