Changeset 2960 in Sophya for trunk/SophyaLib/SkyMap/spherehealpix.cc

Timestamp:

Jun 1, 2006, 4:35:47 PM (19 years ago)

Author:

ansari

Message:

petites corrections + doc(xygen) / Protection Resize() ds SphereHEALPix et SphereThetaPhi - Reza 1/06/2006

File:

• trunk/SophyaLib/SkyMap/spherehealpix.cc (modified) (4 diffs)

trunk/SophyaLib/SkyMap/spherehealpix.cc

@@ -28 +28 @@
 
 \verbatim
-
+    Adapted from :
     -----------------------------------------------------------------------
      version 0.8.2  Aug97 TAC  Eric Hivon, Kris Gorski
@@ -84 +84 @@
 
 {
-
-  if(m <= 0 || m > 8192) 
-    {
-      cout << "SphereHEALPix : m hors bornes [0,8192], m= " << m << endl;
-      throw RangeCheckError("SphereHEALPix<T>::SphereHEALPix() - Out of bound nside (< 8192)!");
-    }
-  // verifier que m est une puissance de deux 
-  int x= m;
-  while(x%2 == 0) x/=2;
-  if(x != 1) 
-    {  
-      cout<<"SphereHEALPix: m doit etre une puissance de deux, m= "<<m<<endl;
-      throw ParmError("SphereHEALPix<T>::SphereHEALPix() - nside != 2^n !");
-    }
   InitNul();
   Pixelize(m);
   SetThetaSlices();
 }
-//++
+//!    Copy constructor
 template<class T>
 SphereHEALPix<T>::SphereHEALPix(const SphereHEALPix<T>& s, bool share)
   : pixels_(s.pixels_, share), sliceBeginIndex_(s.sliceBeginIndex_, share), 
                                 sliceLenght_(s.sliceLenght_, share)
-//    copy constructor
 //--
 {
@@ -214 +199 @@
 
 /*!  \fn   void SOPHYA::SphereHEALPix::Resize(int_4 m)
-  \param <m>   "nside" of the Gorski algorithm
+  \param <m>   "nside" of the HEALPix algorithm
 
   The total number of pixels will be Npix =  12*nside**2
@@ -225 +210 @@
 
 {
+  if ((m <= 0 && nSide_ > 0))  {
+    cout << "SphereHEALPix<T>::Resize(m) with m<=0, NOT resized" << endl;
+    return;
+  }
+  InitNul();
+  Pixelize(m);
+  SetThetaSlices();
+}
+
+template<class T>
+void  SphereHEALPix<T>::Pixelize( int_4 m) 
+
+//    prépare la pixelisation Gorski (m a la même signification 
+//    que pour le constructeur)
+//
+//    
+{
   if (m<=0 || m> 8192) {
-    cout << "SphereHEALPix : m hors bornes [0,8192], m= " << m << endl;
-    exit(1);
+    cout << "SphereHEALPix<T>::Pixelize() m=" << m <<" out of range [0,8192]" << endl;
+    throw ParmError("SphereHEALPix<T>::Pixelize() m out of range");
   }
   // verifier que m est une puissance de deux 
   int x= m;
   while (x%2==0) x/=2;
-  if(x != 1) 
-    {  
-      cout<<"SphereHEALPix: m doit etre une puissance de deux, m= "<<m<<endl;
-      exit(1);
-    }
-  InitNul();
-  Pixelize(m);
-  SetThetaSlices();
-}
-
-template<class T>
-void  SphereHEALPix<T>::Pixelize( int_4 m) 
-
-//    prépare la pixelisation Gorski (m a la même signification 
-//    que pour le constructeur)
-//
-//    
-{
+  if(x != 1) {  
+    cout << "SphereHEALPix<T>::Pixelize() m=" << m << " != 2^n " << endl;
+    throw ParmError("SphereHEALPix<T>::Pixelize() m!=2^n");
+  }
+
   // On memorise les arguments d'appel
   nSide_= m;