Changeset 473 in Sophya for trunk/SophyaLib/Samba/spheregorski.h


Ignore:
Timestamp:
Oct 18, 1999, 4:37:44 PM (26 years ago)
Author:
ansari
Message:

modifs francois : passage en double etc. 18-OCT-99

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/SophyaLib/Samba/spheregorski.h

    r470 r473  
    3232
    3333SphereGorski();
    34 SphereGorski(int_4 m);
     34SphereGorski(int m);
    3535SphereGorski(const SphereGorski<T>& s);
    3636virtual ~SphereGorski();
     
    3939
    4040/* Nombre de pixels du decoupage */
    41 virtual int_4 NbPixels() const;
    42 inline void setNbPixels(int_4 n) {nPix_= n;}
     41virtual int NbPixels() const;
     42inline void setNbPixels(int n) {nPix_= n;}
    4343
    4444/* Valeur du contenu du pixel d'indice "RING" k  */
    45 virtual T& PixVal(int_4 k);
    46 virtual T const& PixVal(int_4 k) const;
     45virtual T& PixVal(int k);
     46virtual T const& PixVal(int k) const;
    4747
    4848/* Nombre de tranches en theta */
    49 int_4 NbThetaSlices() const;
    50 void GetThetaSlice(int_4 index,r_4& theta,TVector<float>& phi,TVector<T>& value) const;
     49int NbThetaSlices() const;
     50void GetThetaSlice(int index,double& theta,TVector<double>& phi,TVector<T>& value) const;
    5151
    5252/* Indice "RING" du pixel vers lequel pointe une direction definie par
    5353ses  coordonnees spheriques */                                   
    54 virtual int_4 PixIndexSph(float theta, float phi) const;
     54virtual int PixIndexSph(double theta,double phi) const;
    5555
    5656/* Coordonnees spheriques du milieu du pixel d'indice "RING" k   */
    57 virtual void PixThetaPhi(int_4 k, float& teta, float& phi) const;
     57virtual void PixThetaPhi(int k,double& theta,double& phi) const;
    5858
    5959/* Pixel Solid angle  (steradians) */
    60 virtual r_8 PixSolAngle(int_4 dummy) const;
    61 inline void setPixSolAngle(r_8 x) {omeg_= x;}
     60virtual double PixSolAngle(int dummy) const;
     61inline void setPixSolAngle(double x) {omeg_= x;}
    6262
    6363// --------------- Specific methods
    6464
    65 virtual void Resize(int_4 m);
     65virtual void Resize(int m);
    6666
    6767inline virtual char* TypeOfMap() const {return "RING";};
    6868
    69 /* Valeur du contenu du pixel d'indice "NEST" k                                 */
    70 virtual T& PixValNest(int_4 k);
    71 virtual T const& PixValNest(int_4 k) const;
     69/* Valeur du contenu du pixel d'indice "NEST" k */
     70virtual T& PixValNest(int k);
     71virtual T const& PixValNest(int k) const;
    7272
    7373/* Indice "NEST" du pixel vers lequel pointe une direction definie par
    7474ses  coordonnees spheriques */                                   
    75 virtual int_4 PixIndexSphNest(float theta, float phi) const;
     75virtual int PixIndexSphNest(double theta,double phi) const;
    7676
    7777/* Coordonnees spheriques du milieu du pixel d'indice "NEST" k       */
    78 virtual void PixThetaPhiNest(int_4 k, float& teta, float& phi) const;
     78virtual void PixThetaPhiNest(int k,double& theta,double& phi) const;
    7979
    8080/* algorithme de pixelisation */
    81 void Pixelize(int_4);
     81void Pixelize(int);
    8282
    8383/* convertit index nested en ring  */
    84 int_4 NestToRing(int_4 ) const;
     84int NestToRing(int) const;
    8585
    8686/* convertit index ring en nested" */
    87 int_4 RingToNest(int_4 ) const;
     87int RingToNest(int) const;
    8888
    89 /*    analyse en harmoniques spheriques des valeurs des pixels de la
     89/* analyse en harmoniques spheriques des valeurs des pixels de la
    9090   sphere : appel du module anafast (Gorski-Hivon) */
    9191//void anharm(int, float, float*);
     
    9797
    9898/* retourne/fixe la valeur du parametre Gorski */
    99 inline virtual int_4 SizeIndex() const {return(nSide_);}
    100 inline void setSizeIndex(int_4 n) {nSide_= n;}
     99inline virtual int SizeIndex() const {return(nSide_);}
     100inline void setSizeIndex(int n) {nSide_= n;}
    101101
    102102/* retourne les pointeurs /remplit les tableaux */
    103103inline const NDataBlock<T>* getDataBlock() const { return (&pixels_); }
    104 inline void setDataBlock(T* data, int_4 m) { pixels_.FillFrom(m,data); }
     104inline void setDataBlock(T* data,int m) { pixels_.FillFrom(m,data); }
    105105
    106106/* retourne/fixe les parametres des modules anafast et synfast */
    107 void getParafast(int_4& nlmax,int_4& nmmax,int_4& iseed,float& fwhm,float& quadr,float& cut) const;
    108 void setParafast(int_4 nlmax,int_4 nmmax,int_4 iseed,float fwhm,float quadr,float cut,char* filename);
     107void getParafast(int& nlmax,int& nmmax,int& iseed,float& fwhm,float& quadr,float& cut) const;
     108void setParafast(int nlmax,int nmmax,int iseed,float fwhm,float quadr,float cut,char* filename);
    109109
    110110/* retourne/fixe le nom du fichier qui contient le spectre de puissance */
     
    119119void InitNul();
    120120
    121 int  nest2ring(int nside, int ipnest) const;
    122 int  ring2nest(int nside, int ipring) const;
     121int  nest2ring(int nside,int ipnest) const;
     122int  ring2nest(int nside,int ipring) const;
    123123
    124 int  ang2pix_ring(int nside, double theta, double phi) const;
    125 int  ang2pix_nest(int nside, double theta, double phi) const;
    126 void pix2ang_ring(int nside, int ipix, double& theta,  double& phi) const;
    127 void pix2ang_nest(int nside, int ipix, double& theta, double& phi) const;
     124int  ang2pix_ring(int nside,double theta,double phi) const;
     125int  ang2pix_nest(int nside,double theta,double phi) const;
     126void pix2ang_ring(int nside,int ipix,double& theta,double& phi) const;
     127void pix2ang_nest(int nside,int ipix,double& theta,double& phi) const;
    128128
    129129// ------------- variables internes -----------------------
    130 int_4 nSide_;
    131 int_4 nPix_;
    132 r_8 omeg_;
     130int nSide_;
     131int nPix_;
     132double omeg_;
    133133
    134134NDataBlock<T> pixels_;
Note: See TracChangeset for help on using the changeset viewer.