Changeset 2990 in Sophya for trunk/SophyaLib/SkyMap/spherethetaphi.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/spherethetaphi.cc (modified) (15 diffs)

trunk/SophyaLib/SkyMap/spherethetaphi.cc

@@ -192 +192 @@
   NPix_  = a.NPix_;
   Omega_ = a.Omega_;
-  int k;
+  int_4 k;
   for (k=0; k< NPix_; k++) pixels_(k) = a.pixels_(k);
   for (k=0; k< a.NPhi_.Size(); k++) NPhi_(k) = a.NPhi_(k);
@@ -246 +246 @@
 {
   double dphi;
-  int i,k;
+  int_4 i,k;
   bool fgzn = false;
 
@@ -270 +270 @@
 void SphereThetaPhi<T>::PixThetaPhi(int_4 k,double& theta,double& phi) const
 {
-  int i; 
+  int_4 i; 
   bool fgzn = false;
   
@@ -316 +316 @@
 void SphereThetaPhi<T>::Limits(int_4 k,double& tetMin,double& tetMax,double& phiMin,double& phiMax) 
   {
-  int j;
+  int_4 j;
   double dphi;
   bool fgzn= false;
@@ -335 +335 @@
   
   // on recupere l'indice i de la tranche qui contient le pixel k
-  int i;
+  int_4 i;
   for( i=0; i< NTheta_; i++ ) 
     if(k < TNphi_(i+1)) break;
@@ -378 +378 @@
 int_4 SphereThetaPhi<T>::NPhi(int_4 kt) const
 {
-  int nbpix;  
+  int_4 nbpix;  
   // verification
   if((kt < 0) || (kt >= 2*NTheta_)) return(-1);
@@ -455 +455 @@
 int_4 SphereThetaPhi<T>::Index(int_4 kt,int_4 jp) const
 {
-  int k;
+  int_4 k;
   bool fgzn= false;
   
@@ -494 +494 @@
   
   // on recupere l'indice kt de la tranche qui contient le pixel k
-  int i;
+  int_4 i;
   for(i = 0; i < NTheta_; i++) 
     if(k < TNphi_(i+1)) break;
@@ -523 +523 @@
 //--
 {
-  int ntotpix,j;
+  int_4 ntotpix,j;
   
   // Decodage et controle des arguments d'appel :
@@ -634 +634 @@
     
 
-  int iring= Index(index,0);
-  int lring  = NPhi(index);
+  int_4 iring= Index(index,0);
+  int_4 lring  = NPhi(index);
 
   phi.ReSize(lring);
@@ -643 +643 @@
   PixThetaPhi(iring,Te,Fi);
   double DFi = DeuxPi/(double)NPhi(index);
-  for(int kk = 0; kk < lring; kk++)   {
+  for(int_4 kk = 0; kk < lring; kk++)   {
     value(kk)= pixels_(kk+iring);
     phi(kk)= Fi;
@@ -668 +668 @@
     throw RangeCheckError("SphereThetaPhi::GetThetaSlice(idx...)  idx out of range"); 
 
-  int iring= Index(index,0);
-  int lring  = NPhi(index);
+  int_4 iring= Index(index,0);
+  int_4 lring  = NPhi(index);
 
   pixelIndices.ReSize(lring);
@@ -675 +675 @@
   double Te= 0.;
   double Fi= 0.;
-  for(int kk = 0; kk < lring; kk++)  {
+  for(int_4 kk = 0; kk < lring; kk++)  {
     pixelIndices(kk)=kk+iring ;
     value(kk)= pixels_(kk+iring);
@@ -682 +682 @@
 }
 
-
+//! return a pointer to the specified slice pixel data
+template <class T>
+T* SphereThetaPhi<T>::GetThetaSliceDataPtr(int_4 index)
+{
+  if(index < 0 || index >= NbThetaSlices()) 
+    throw RangeCheckError("SphereThetaPhi::GetThetaSliceDataPtr(idx)  idx out of range"); 
+  return pixels_.Begin()+Index(index,0);
+}
 
 
@@ -694 +701 @@
 
   os << "... NPhi_ Values : ";
-  int i;
+  int_4 i;
   for(i=0; i < NTheta_; i++)
     {