Changeset 2985 in Sophya for trunk/SophyaLib/SkyMap/spherethetaphi.cc

Timestamp:

Jun 21, 2006, 5:55:29 PM (19 years ago)

Author:

ansari

Message:

1/ Suite codage HEALPix-NEST , test avec transforme Ylm OK
2/ Petites modifs et corrections diverses (ajout SizeIndex2Resol() par exemple)

Reza 21 Juin 2006

File:

• trunk/SophyaLib/SkyMap/spherethetaphi.cc (modified) (8 diffs)

trunk/SophyaLib/SkyMap/spherethetaphi.cc

@@ -631 +631 @@
 
   if(index < 0 || index >= NbThetaSlices()) 
-    throw RangeCheckError("SphereThetaPhi::GetThetaSlice() index out of range"); 
+    throw RangeCheckError("SphereThetaPhi::GetThetaSlice(idx...) index out of range"); 
     
 
@@ -641 +641 @@
   double Te= 0.;
   double Fi= 0.;
-  for(int kk = 0; kk < lring; kk++) 
-    {
-      PixThetaPhi(kk+iring,Te,Fi);
-      phi(kk)= Fi;
-      value(kk)= PixVal(kk+iring);
-    }
+  PixThetaPhi(iring,Te,Fi);
+  double DFi = DeuxPi/(double)NPhi(index);
+  for(int kk = 0; kk < lring; kk++)   {
+    value(kk)= pixels_(kk+iring);
+    phi(kk)= Fi;
+    Fi += DFi;
+  }
   theta= Te;
 }
@@ -665 +666 @@
 
   if(index < 0 || index >= NbThetaSlices()) 
-    {
-       throw RangeCheckError("SphereThetaPhi::PIxVal Pixel index out of range"); 
- 
-    }
+    throw RangeCheckError("SphereThetaPhi::GetThetaSlice(idx...)  idx out of range"); 
 
   int iring= Index(index,0);
@@ -677 +675 @@
   double Te= 0.;
   double Fi= 0.;
-  for(int kk = 0; kk < lring; kk++) 
-    {
-      pixelIndices(kk)=kk+iring ;
-      value(kk)= PixVal(kk+iring);
-    }
-      PixThetaPhi(iring,theta,phi0);
+  for(int kk = 0; kk < lring; kk++)  {
+    pixelIndices(kk)=kk+iring ;
+    value(kk)= pixels_(kk+iring);
+  }
+  PixThetaPhi(iring,theta,phi0);
 }
 
@@ -692 +689 @@
 {
   Show(os);
+  os << "SphereThetaPhi<T> NTheta_= " << NTheta_ << " NPix_    = " << NPix_ 
+     << " Omega_  =  " << Omega_   << endl;
   if(this->mInfo_) os << "  DVList Info= " << *(this->mInfo_) << endl;
-  //
-  os << " NTheta_= " << NTheta_ << endl;
-  os << " NPix_    = " << NPix_   << endl;
-  os << " Omega_  =  " << Omega_   << endl;
-
-  os << " contenu de NPhi_ : ";
+
+  os << "... NPhi_ Values : ";
   int i;
   for(i=0; i < NTheta_; i++)
@@ -707 +702 @@
   os << endl;
 
-  os << " contenu de Theta_ : ";
+  os << "... Theta_ Values : ";
   for(i=0; i < NTheta_+1; i++)
     {
@@ -715 +710 @@
   os << endl;
 
-  os << " contenu de TNphi_ : ";
+  os << "... TNphi_ Values : ";
   for(i=0; i < NTheta_+1; i++)
     {
@@ -723 +718 @@
   os << endl;
 
-  os << " contenu de pixels : ";
+  os << "... Pixel Values : ";
   for(i=0; i < NPix_; i++)
     {