Changeset 3510 in Sophya for trunk/SophyaLib/Samba/alm.cc

Timestamp:

Aug 8, 2008, 3:11:45 PM (17 years ago)

Author:

ansari

Message:

Nettoyage alm.h .cc , utilisation RandomGenerator par Alm<T> et SphericalTransformServer , Reza 08/08/2008

File:

• trunk/SophyaLib/Samba/alm.cc (modified) (6 diffs)

trunk/SophyaLib/Samba/alm.cc

@@ -1 +1 @@
 #include "sopnamsp.h"
 #include "alm.h"
+
+
+/*!
+  \class Alm 
+  \ingroup Samba 
+  Class for manipulating the coefficients \f$a_{lm}\f$ of the development
+  in spherical harmonics of  a function efined on a sphere.
+*/
+ 
+/*!
+  fwhm specifies the gaussian beam half witdh in arc.minutes
+*/
 template <class T>
 Alm<T>::Alm(const TVector<T>& clin, const r_8 fwhm) 
 
 {
+  int_4 nlmax= clin.NElts()-1;
+
+  //alm.ReSize(nlmax);
+  this->ReSizeRow(nlmax+1);
+  RandomGenerator rg(1, false);
+  GenFromCl(clin, fwhm, rg);
+}
+
+/*!
+  fwhm specifies the gaussian beam half witdh in arc.minutes
+*/
+template <class T>
+Alm<T>::Alm(const TVector<T>& clin, const r_8 fwhm, RandomGenerator & rg)
+{
+  int_4 nlmax= clin.NElts()-1;
+
+  //alm.ReSize(nlmax);
+  this->ReSizeRow(nlmax+1);
+  GenFromCl(clin, fwhm, rg);  
+}
 
 
+
+template <class T>
+void Alm<T>::GenFromCl(const TVector<T> & clin, const r_8 fwhm, RandomGenerator & rg)
+{
   /*=======================================================================
      creates the a_lm from the power spectrum, 
@@ -21 +57 @@
   int_4 nlmax= clin.NElts()-1;
 
-  //alm.ReSize(nlmax);
-  this->ReSizeRow(nlmax+1);
-
   r_8 sig_smooth = fwhm/sqrt(8.*log(2.))/(60.*180.)* M_PI;
   int_4 n_l = nlmax+1;
@@ -29 +62 @@
 
   //    --- smoothes the initial power spectrum ---
-  TVector<T> cl=clin;
+  TVector<T> cl(clin, false);
   int l;
   for (l=0;l<n_l;l++)
@@ -45 +78 @@
       T rms=sqrt(cl(l));
       //        ------ m = 0 ------
-      complex<T> zeta1(NorRand());
+      complex<T> zeta1((T)rg.Gaussian() );
       (*this)(l,0)   = zeta1 * rms;
 
@@ -52 +85 @@
         {
           complex<T> aux1(hsqrt2);
-          complex<T> aux2(NorRand(),NorRand());
+          complex<T> aux2((T)rg.Gaussian() , (T)rg.Gaussian() );
           zeta1=aux1*aux2;
           (*this)(l,m)=rms*zeta1;
@@ -79 +112 @@
   return powsp;
 }
+
+/*!
+  \class Bm
+  \ingroup Samba 
+  Class for a vector with an index running from \f$-m_{max}\f$ to \f$+m_{max}\f$ 
+  (then the size of the vector will be actually \f$2m_{max}+1)\f$.
+  This class is used by the spherical harmonics transform server 
+*/
+
  
 #ifdef __CXX_PRAGMA_TEMPLATES__