Changeset 880 in Sophya for trunk/SophyaProg/PMixer/extractRS.cc


Ignore:
Timestamp:
Apr 11, 2000, 3:34:55 PM (25 years ago)
Author:
ansari
Message:

Adaptation SphereGorski -> SphereHEALPix apres modifs Sophie - Reza 11/4/2000

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/SophyaProg/PMixer/extractRS.cc

    r877 r880  
    1212#include "tvector.h"
    1313#include "vector3d.h"
    14 #include "nbrandom.h"
     14#include "srandgen.h"
    1515#include "timing.h"
    1616#include "sambainit.h"
     
    8888 
    8989 
    90   cout << " extractRS/Info : NComp = " <<  nskycomp << " SphereGorski_NSide= " << hp_nside << endl;
     90  cout << " extractRS/Info : NComp = " <<  nskycomp << " SphereHEALPix_NSide= " << hp_nside << endl;
    9191  cout << "  ... MapPath = " << (string)mapPath << "  DebugLev= " << debuglev
    9292       << "  PrintLev= " << printlev << endl;
     
    9595  // We create an output persist file for writing debug objects
    9696  if (debuglev > 0) so = new POutPersist("extrRSdbg.ppf");
    97   SphereGorski<float> SourceMap(hp_nside);
    98   SphereGorski<float> OutputMap(hp_nside);
     97  SphereHEALPix<float> SourceMap(hp_nside);
     98  SphereHEALPix<float> OutputMap(hp_nside);
    9999  bool okout = true;
    100100 
     
    190190        K = dc.DParam(key, 1, 1.);
    191191        if (debuglev > 4) {  // Writing the input map to the outppf
    192           FIO_TMatrix<float> fiog;
     192          FIO_TArray<float> fiog;
    193193          fiog.Write(*so, key);
    194194        }
     
    282282      //        MeanSig(SourceMap.DataBlock(),moy,sig);
    283283      //        cout << " MeanSig for Source Map - Mean= " << moy << " Sigma= " << sig << endl;
    284       //        cout << "Source Map (SphereGorski<float>) written to FITS file "
     284      //        cout << "Source Map (SphereHEALPix<float>) written to FITS file "
    285285      //             << (string)(arg[3]) << endl;
    286286      //       }
     
    291291        MeanSig(OutputMap.DataBlock(),moy,sig);
    292292        cout << " MeanSig for Output Map - Mean= " << moy << " Sigma= " << sig << endl;
    293         cout << "Output Map (SphereGorski<float>) written to FITS file "
     293        cout << "Output Map (SphereHEALPix<float>) written to FITS file "
    294294             << (string)(arg[2]) << endl;
    295295      }
     
    298298      if (narg > 2) {
    299299        POutPersist s(arg[3]);
    300         FIO_SphereGorski<float> fiog(OutputMap);
     300        FIO_SphereHEALPix<float> fiog(OutputMap);
    301301        fiog.Write(s);
    302         cout << "Output Map (SphereGorski<float>) written to POutPersist file " 
     302        cout << "Output Map (SphereHEALPix<float>) written to POutPersist file " 
    303303             << (string)(arg[3]) << endl;
    304304        if(printlev>2) PrtTim("End of WritePPF ");
Note: See TracChangeset for help on using the changeset viewer.