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

Timestamp:

Jun 21, 2006, 5:53:06 PM (19 years ago)

Author:

ansari

Message:

Remplacement appel a RfourierSynthesisFromB (TF reel) par fourierSynthesisFromB (TF complexe) pour corriger pb transforme Ylm sur HEALPix (zones polar cap) , Reza 21/6/2006

File:

• trunk/SophyaLib/Samba/sphericaltransformserver.cc (modified) (3 diffs)

trunk/SophyaLib/Samba/sphericaltransformserver.cc

@@ -175 +175 @@
   map.Resize(pixelSizeIndex);
   string sphere_type=map.TypeOfMap();
-      int premiereTranche = 0;
-      int derniereTranche = map.NbThetaSlices()-1;
-  if (sphere_type.substr(0,4) == "RING")
-  {
-    nsmax=map.SizeIndex();
-  }
-  else
-  {
-    // pour une sphere Gorski le nombre de pixels est 12*nsmax**2
-    // on calcule une quantite equivalente a nsmax pour la sphere-theta-phi
-    // en vue de l'application du critere Healpix : nlmax<=3*nsmax-1
-    // c'est approximatif ; a raffiner.
-    if (sphere_type.substr(0,6) == "TETAFI")
-    {
-      nsmax=(int_4)sqrt(map.NbPixels()/12.);
-      premiereTranche++;
-      derniereTranche--;
-    }
-    else
-    {
-      cout << " unknown type of sphere : " << sphere_type << endl;
-      throw IOExc(" unknown type of sphere: " + (string)sphere_type );
-    }
-  //  cout << "GenerateFromAlm: the sphere is of type : " << sphere_type << endl;
-  //  cout << "GenerateFromAlm: size index (nside) of the sphere= " << nsmax << endl;
-  //  cout << "GenerateFromAlm: nlmax (from Alm) = " << nlmax << endl;
-  //  if (nlmax>3*nsmax-1) 
-  //  {
-      //     cout << "GenerateFromAlm: nlmax should be <= 3*nside-1" << endl;
-  //   if (strncmp(sphere_type,"TETAFI",6) == 0)
-  //    {
-  //      cout << "GenerateFromAlm: nlmax should be <= 3*nside-1" << endl;
-  //      cout << " (for this criterium, nsmax is computed as sqrt(nbPixels/12))" << endl;
-  //    }
-  //}
-  }
+  int premiereTranche = 0;
+  int derniereTranche = map.NbThetaSlices()-1;
+
   Bm<complex<T> > b_m_theta(nmmax);
 
-  //  map.Resize(nsmax);
 
 
@@ -227 +193 @@
       map.GetThetaSlice(ith,theta,phi0, pixNumber,datan);
       nph = pixNumber.NElts();
-
+      if (nph < 2) continue;  // On laisse tomber les tranches avec un point
       //       -----------------------------------------------------
       //              for each theta, and each m, computes
@@ -259 +225 @@
       // ---------------------------------------------------------------*/
 
-
-      if (sphere_type.substr(0,4) == "RING")
-        {
+      /* ----- Reza, Juin 2006 : 
+         En verifiant la difference entre deux cartes 
+         cl -> map -> alm -> map2 et mapdiff = map-map2 
+         je me suis apercu qu'il y avait des differences importantes - dans les
+         deux zones 'polar cap' de HEALPix - apres avoir pas mal chercher, il semble
+         que la routine RfourierSynthesisFromB en soit responsable 
+         Je fais donc tout passer dans fourierSynthesisFromB 
+      if ( (sphere_type == "RING") || (sphere_type == "NESTED") ) {
           TVector<T> Temp = RfourierSynthesisFromB(b_m_theta,nph,phi0);
-          for (int i=0;i< nph;i++) map(pixNumber(i))=Temp(i);
+          for (int i=0;i< nph;i++) map(pixNumber(i))=Temp(i); 
         }
       else 
+      */
         // pour des pixelisations quelconques (autres que HEALPix
         //  nph n'est pas toujours pair