Changeset 2973 in Sophya for trunk/SophyaLib/SkyMap/HEALPixUtils.cc

Timestamp:

Jun 20, 2006, 6:01:49 PM (19 years ago)

Author:

ansari

Message:

1/ Nettoyage+commentaires/doxygen ds Vector3d, UnitVector, LongLat,
SphereCoordSys ...
2/ Ajout de la classe angle pour faciliter les conversions rad<>deg<>arcmin
dans le fichier vector3d.h .cc
3/ nettoyage/uniformisation methodes print pour pixelmap, ajout de la
methode PixelMap<T>::Show()
4/ Ajout methodes SphericalMap<T>::HasSymThetaSlice() et
GetSymThetaSliceIndex(int_4 idx) et leurs implementations pour
SphereHEALPix et SphereThetaPhi en vue de l'optimisation du calcul
transforme Ylm
5/ Ajout methode ResolToSizeIndex ds SphereThetaPhi , SphereHEALPix et
SphereECP

Reza , 20 Juin 2006

File:

• trunk/SophyaLib/SkyMap/HEALPixUtils.cc (modified) (1 diff)

trunk/SophyaLib/SkyMap/HEALPixUtils.cc

@@ -687 +687 @@
 
  
+/*! 
+  \brief return the size index value corresponding to resolution \b res  
+  Computes the size index m providing a resolution better or equal the 
+  specified value \b res (in radian)
+
+  - m : size-index 
+  - number of pixels for 4Pi : 12*m*m
+
+  \sa SOPHYA::Angle
+*/
+int_4 HEALPix::ResolToSizeIndex(double res)
+{
+  int_4 m = 1;
+  double cres = sqrt(M_PI/3.); // 4PI/12
+  while (cres > 1.001*res) {
+    m *= 2;
+    cres = sqrt(4.*M_PI/(12*m*m));
+  }
+  return m;
+}