Changeset 1954 in Sophya for trunk/SophyaLib/SkyMap

Timestamp:

Mar 28, 2002, 2:10:57 PM (24 years ago)

Author:

ansari

Message:

protection/lancement exception ds HEALPix::xxx et FIO_SphereHEALPix<T>::SetDataObj() - Reza 28/3/2002

Location:

trunk/SophyaLib/SkyMap

Files:

• HEALPixUtils.cc (modified) (6 diffs)
• fiospherehealpix.cc (modified) (4 diffs)

trunk/SophyaLib/SkyMap/HEALPixUtils.cc

@@ -175 +175 @@
   
   if(  nside<1 ||  nside>ns_max ) {
-    cout << "nside out of range" << endl;
-    exit(0);
+    char buff[64];
+    sprintf(buff,"HEALPix::nest2ring nside(=%d) out of range ", nside);
+    throw RangeCheckError(PExcLongMessage(buff));
   }
   npix = 12 *  nside* nside;
   if( ipnest<0 || ipnest>npix-1 ) {
-    cout << "ipnest out of range" << endl;
-    exit(0);
+    char buff[64];
+    sprintf(buff,"HEALPix::nest2ring ipnest(=%d) out of range ", ipnest);
+    throw RangeCheckError(PExcLongMessage(buff));
   }
 
@@ -261 +263 @@
   
   if( nside<1 || nside>ns_max ) {
-    cout << "nside out of range" << endl;
-    exit(0);
+    char buff[64];
+    sprintf(buff,"HEALPix::ring2nest nside(=%d) out of range ", nside);
+    throw RangeCheckError(PExcLongMessage(buff));
   }
   npix = 12 * nside*nside;
   if( ipring<0 || ipring>npix-1 ) {
-    cout << "ipring out of range" << endl;
-    exit(0);
+    char buff[64];
+    sprintf(buff,"HEALPix::ring2nest ipring(=%d) out of range ", ipring);
+    throw RangeCheckError(PExcLongMessage(buff));
   }
   
@@ -364 +368 @@
 
   if( nside<1 || nside>ns_max ) {
-    cout << "nside out of range" << endl;
-    exit(0);
+    char buff[64];
+    sprintf(buff,"HEALPix::ang2pix_ring nside(=%d) out of range ", nside);
+    throw RangeCheckError(PExcLongMessage(buff));
   }
 
   if( theta<0. || theta>Pi) {
-    cout << "theta out of range" << endl;
-    exit(0);
+    char buff[64];
+    sprintf(buff,"HEALPix::ang2pix_ring theta(=%g) out of range ", theta);
+    throw RangeCheckError(PExcLongMessage(buff));
   }
 
@@ -447 +453 @@
   
   if( nside<1 || nside>ns_max ) {
-    cout << "nside out of range" << endl;
-    exit(0);
+    char buff[64];
+    sprintf(buff,"HEALPix:::ang2pix_nest nside(=%d) out of range ", nside);
+    throw RangeCheckError(PExcLongMessage(buff));
   }
   if( theta<0 || theta>Pi ) {
-    cout << "theta out of range" << endl;
-    exit(0);
+    char buff[64];
+    sprintf(buff,"HEALPix:::ang2pix_nest theta(=%g) out of range ", theta);
+    throw RangeCheckError(PExcLongMessage(buff));
   }
   z  = cos(theta);
@@ -531 +539 @@
   
   if( nside<1 || nside>ns_max ) {
-    cout << "nside out of range" << endl;
-    exit(0);
+    char buff[64];
+    sprintf(buff,"HEALPix::pix2ang_ring nside(=%d) out of range ", nside);
+    throw RangeCheckError(PExcLongMessage(buff));
   }
   npix = 12*nside*nside;      // ! total number of points
   if( ipix<0 || ipix>npix-1 ) {
-    cout << "ipix out of range" << endl;
-    exit(0);
+    char buff[64];
+    sprintf(buff,"HEALPix::pix2ang_ring ipix(=%d) out of range ", ipix);
+    throw RangeCheckError(PExcLongMessage(buff));
   }
   
@@ -609 +619 @@
   
   if( nside<1 || nside>ns_max ) {
-    cout << "nside out of range" << endl;
-    exit(0);
+    char buff[64];
+    sprintf(buff,"HEALPix::pix2ang_nest nside(=%d) out of range ", nside);
+    throw RangeCheckError(PExcLongMessage(buff));
   }
   npix = 12 * nside*nside;
   if( ipix<0 || ipix>npix-1 ) {
-    cout << "ipix out of range" << endl;
-    exit(0);
+    char buff[64];
+    sprintf(buff,"HEALPix::pix2ang_nest ipix(=%d) out of range ", ipix);
+    throw RangeCheckError(PExcLongMessage(buff));
   }
   

trunk/SophyaLib/SkyMap/fiospherehealpix.cc

@@ -4 +4 @@
 #include "pexceptions.h"
 #include "fiondblock.h"
+#include "datatype.h"
+#include <typeinfo>
 
 //*******************************************************************
@@ -55 +57 @@
 {
   SphereHEALPix<T> * po = dynamic_cast< SphereHEALPix<T> * >(&o);
-  if (po == NULL) return;
+  if (po == NULL) { 
+    char buff[160];
+    sprintf(buff,"FIO_SphereHEALPix<%s>::SetDataObj(%s) - Object type error ! ",
+            DataTypeInfo<T>::getTypeName().c_str(), typeid(o).name());
+    throw TypeMismatchExc(PExcLongMessage(buff));
+  }
   if (ownobj && dobj_) delete dobj_;
   dobj_ = po; ownobj = false;
@@ -93 +100 @@
   is.GetR8(Omega);
   dobj_->setParameters(nSide,nPix, Omega);
-
+  cout << " DBG-FIO_SphereHEALPix<T> nSide=" << nSide << " nPix=" << nPix << endl;
 // On lit les DataBlocks;
   FIO_NDataBlock<T> fio_nd(&dobj_->pixels_);
@@ -101 +108 @@
   FIO_NDataBlock<int_4> fio_ndi2(&dobj_->sliceLenght_);
   fio_ndi2.Read(is);
+
+  cout << " DBG-2 SizeIndex()= " << dobj_->SizeIndex() 
+       << " NbPixels()= " << dobj_->NbPixels() << endl;
 }