Changeset 2958 in Sophya for trunk/SophyaLib/Samba/alm.h

Timestamp:

Jun 1, 2006, 1:34:50 PM (19 years ago)

Author:

ansari

Message:

1/ passage en sa_size_t (au lieu de int) dans Alm<T> et Bm<T>
2/ Ajout des methodes optimisees (statiques) pour calcul transforme Ylm
ds LambdaLMBuilder et utilisation ds SphericalTransformServer

Reza , 1/06/2006

File:

• trunk/SophyaLib/Samba/alm.h (modified) (3 diffs)

trunk/SophyaLib/Samba/alm.h

@@ -21 +21 @@
 Alm() : TriangularMatrix<complex<T> >()  {;};
     /* instanciate the class from the maximum value of the index \f$l\f$ in the sequence of the \f$a_{lm}\f$'s */
-Alm(int nlmax) : TriangularMatrix<complex<T> >(nlmax+1) {;}
+Alm(sa_size_t nlmax) : TriangularMatrix<complex<T> >(nlmax+1) {;}
 Alm(const Alm<T>& a,  bool share=false)  : TriangularMatrix<complex<T> >(a, share)  {;}
 
 Alm(const TVector<T>& clin, const r_8 fwhm) ;
 /*! resize with a new lmax */
-inline void ReSizeToLmax(int_4 nlmax) {this->ReSizeRow(nlmax+1);}
-inline int_4 Lmax() const {return this->rowNumber()-1;}
+inline void ReSizeToLmax(sa_size_t nlmax) {this->ReSizeRow(nlmax+1);}
+inline sa_size_t Lmax() const {return this->rowNumber()-1;}
 TVector<T> powerSpectrum() const;
 
@@ -48 +48 @@
 Bm(): nmmax_(0) {;};
     /* instanciate from the maximum absolute value of the index */
-Bm(int mmax) : nmmax_((uint_4)mmax)   { bm_.ReSize( (uint_4)(2*mmax+1) );}
+Bm(sa_size_t mmax) : nmmax_(mmax)   { bm_.ReSize( (2*mmax+1) );}
  Bm(const Bm<T>& b,  bool share=false) : nmmax_(b.nmmax_), bm_(b.bm_, share) {;}
  /*! resize with a new value of mmax */
-inline void ReSizeToMmax(int_4 mmax) 
+inline void ReSizeToMmax(sa_size_t mmax) 
   {
-    nmmax_=(uint_4)mmax;
+    nmmax_= mmax;
     bm_.ReSize(2* nmmax_+1);
   }
-//inline int_4 Size() const {return 2*nmmax_+1;};
-inline T& operator()(int m) {return bm_( adr_i(m));};
-inline T const& operator()(int m) const {return *(bm_.Begin()+ adr_i(m));};
+//inline sa_size_t Size() const {return 2*nmmax_+1;};
+inline T& operator()(sa_size_t m) {return bm_( adr_i(m));};
+inline T const& operator()(sa_size_t m) const {return *(bm_.Begin()+ adr_i(m));};
 /*! return the current value of the maximum absolute value of the index */
-inline int_4 Mmax() const 
+inline sa_size_t Mmax() const 
    {
-     return (int_4)nmmax_;
+     return nmmax_;
    }
 inline Bm<T>& operator = (const Bm<T>& b)
@@ -76 +76 @@
   private:
 /*! return the address in the array representing the vector */
-inline uint_4 adr_i(int i) const 
+inline sa_size_t adr_i(sa_size_t i) const 
 {
-  return(uint_4)(i+nmmax_);
+  return (i+nmmax_);
 }
 
 
 
- uint_4 nmmax_;
+ sa_size_t nmmax_;
  NDataBlock<T> bm_;