Changeset 2958 in Sophya for trunk/SophyaLib/Samba/sphericaltransformserver.cc

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/sphericaltransformserver.cc (modified) (4 diffs)

trunk/SophyaLib/Samba/sphericaltransformserver.cc

@@ -232 +232 @@
       //              b(m,theta) = sum_over_l>m (lambda_l_m(theta) * a_l_m) 
       //              ------------------------------------------------------
+      // ===> Optimisation Reza, Mai 2006 
+      /*---  Le bout de code suivant est remplace par l'appel a la nouvelle fonction
+        qui calcule la somme au vol 
       LambdaLMBuilder lb(theta,nlmax,nmmax);
       //  somme sur m de 0 a l'infini
-      int m;
-      for (m = 0; m <= nmmax; m++)
+      for (int_4 m = 0; m <= nmmax; m++)
         {
           b_m_theta(m) = (T)( lb.lamlm(m,m) ) * alm(m,m);
@@ -243 +245 @@
             }
         }
+      ------- Fin version PRE-Mai2006 */
+      LambdaLMBuilder::ComputeBmFrAlm(theta,nlmax,nmmax, alm, b_m_theta);
+      //Fin Optimisation Reza, Mai 2006 <==== 
+
     //        obtains the negative m of b(m,theta) (= complex conjugate)
 
-      for (m=1;m<=nmmax;m++)
+      for (int_4 m=1;m<=nmmax;m++)
         {
           b_m_theta(-m) = conj(b_m_theta(m));
@@ -555 +561 @@
 //      for each m and l
 //      -----------------------------------------------------------------------
+
+      // ===> Optimisation Reza, Mai 2006 
+      /*---  Le bout de code suivant est remplace par l'appel a la nouvelle fonction
+        qui calcule la somme au vol 
       //        PrtTim("avant instanciation LM ");
       LambdaLMBuilder lb(theta,nlmax,nmmax);
@@ -576 +586 @@
             }
         }
+      ------- Fin version PRE-Mai2006 */
+      r_8 domega=map.PixSolAngle(map.PixIndexSph(theta,phi0));
+      phase *= complex<T>((T)domega, 0.);
+      LambdaLMBuilder::ComputeAlmFrPhase(theta,nlmax,nmmax, phase, alm);
+      //Fin Optimisation Reza, Mai 2006 <====