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

Timestamp:

Jun 23, 2006, 12:35:17 PM (19 years ago)

Author:

ansari

Message:

Ajouts methode GetThetaSliceDataPtr() dans SphericalMap, SphereThetaPhi SphereHEALPix SphereECP pour optimisation calcul transforme Ylm + passage TypeOfMap TETAFI->ECP pour SphereECP , Reza 23/6/2006

File:

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

trunk/SophyaLib/SkyMap/spherehealpix.cc

@@ -340 +340 @@
 void SphereHEALPix<T>::GetThetaSlice(int_4 index,r_8& theta,TVector<r_8>& phi,TVector<T>& value) const 
 {
-  if (index<0 || index >= NbThetaSlices())  {
-      cout << " SphereHEALPix::GetThetaSlice : Pixel index out of range" <<endl;
-      throw RangeCheckError(" SphereHEALPix::GetThetaSlice()  index out of range");
-  }
+  if (index<0 || index >= NbThetaSlices()) 
+    throw RangeCheckError(" SphereHEALPix::GetThetaSlice()  index out of range");
 
   int_4 iring= sliceBeginIndex_(index);
@@ -388 +386 @@
 
 template<class T>
-void SphereHEALPix<T>::GetThetaSlice(int_4 sliceIndex,r_8& theta, r_8& phi0, TVector<int_4>& pixelIndices,TVector<T>& value) const 
+void SphereHEALPix<T>::GetThetaSlice(int_4 sliceIndex,r_8& theta, r_8& phi0, 
+                                     TVector<int_4>& pixelIndices,TVector<T>& value) const 
 
 {
 
   if (sliceIndex<0 || sliceIndex >= NbThetaSlices()) 
-    {
-      cout << " SphereHEALPix::GetThetaSlice : Pixel index out of range" <<endl;
-      throw RangeCheckError(" SphereHEALPix::GetThetaSlice : Pixel index out of range");
-    }
+    throw RangeCheckError(" SphereHEALPix::GetThetaSlice() index out of range");
   int_4 iring= sliceBeginIndex_(sliceIndex);
   int_4 lring  =  sliceLenght_(sliceIndex);
@@ -417 +413 @@
     PixThetaPhi(ring2nest(nSide_,iring), theta, phi0);    
   }
+}
+
+//! return a pointer to the specified slice pixel data in RING ordering, NULL in NESTED
+template<class T>
+T* SphereHEALPix<T>::GetThetaSliceDataPtr(int_4 sliceIndex) 
+
+{
+  if (sliceIndex<0 || sliceIndex >= NbThetaSlices()) 
+    throw RangeCheckError(" SphereHEALPix::GetThetaSliceDataPtr(): index out of range");
+  if (fgring_) 
+    return pixels_.Begin()+sliceBeginIndex_(sliceIndex);
+  else return NULL;
 }