Changeset 2973 in Sophya for trunk/SophyaLib/SkyMap/vector3d.h

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/vector3d.h (modified) (6 diffs)

trunk/SophyaLib/SkyMap/vector3d.h

@@ -1 +1 @@
 //   3-D Geometry 
 //        B. Revenu, G. Le Meur   2000
+//        R. Ansari 2006
 // DAPNIA/SPP (Saclay) / CEA    LAL - IN2P3/CNRS  (Orsay)
 
@@ -15 +16 @@
 #include "longlat.h"
 
+
+namespace SOPHYA { 
+
+//! Class to ease angle conversions (radian <> degree <> arcmin <> arcsec)
+class Angle {
+public:
+  enum AngleUnit { Radian, Degree, ArcMin, ArcSec };
+  //! Constructor with specification of angle value in radian 
+  Angle(double val=0.) { _angrad = val; }
+  //! Constructor with specification of angle value and unit
+  Angle(double val, Angle::AngleUnit un);
+  //! Copy constructor
+  Angle(Angle const& a) { _angrad = a._angrad; }
+
+  //! Conversion to radian
+  inline double ToRadian() const { return _angrad; }
+  //! Conversion to degree
+  inline double ToDegree() const { return _angrad*_rad2deg; }
+  //! Conversion to arcmin
+  inline double ToArcMin() const { return _angrad*_rad2min; }
+  //! Conversion to arcsec
+  inline double ToArcSec() const { return _angrad*_rad2sec; }
+
+  //! return the angle value in radian
+  inline operator double () const { return _angrad; }
+
+protected:  
+  double _angrad;  // angle in radians 
+
+  static double _deg2rad;  // deg -> radian conversion factor
+  static double _rad2deg;  // rad -> degree conversion factor
+  static double _rad2min;  // rad -> arcmin conversion factor
+  static double _rad2sec;  // rad -> arcmin conversion factor
+  
+};
+
 /*
   Geometrie en dimension 3. 
@@ -20 +57 @@
   et en coordonnees spheriques theta,phi
   pour les rotations (angles d'Euler) ma source est 
+  B. Revenu / G. Le Meur
   "Classical Mechanics" 2nd edition, H. Goldstein, Addison Wesley
 */
-/*!    3-D geometry. 
-
-    All computations are made with angles in radians and with spherical
-    coordinates theta, phi.
-
-    Concerning Euler's angles, the reference is :
- 
-    "Classical Mechanics" 2nd edition, H. Goldstein, Addison Wesley
-*/
-
-namespace SOPHYA { 
  
 class Vector3d 
@@ -65 +92 @@
 
   // ecart angulaire entre 2 vecteurs dans [0,Pi]
-  /*!   angular gap between 2 vectors in [0,Pi] */
+  //!   angular gap between 2 vectors in [0,Pi] 
   virtual double SepAngle(const Vector3d&) const;
 
   // produit vectoriel
-  /*!    vector product */
-  virtual Vector3d Vect(const Vector3d&) const;
+  //! return the vector product (*this)^v2
+  virtual Vector3d Vect(const Vector3d& v2) const;
 
   // vecteur perpendiculaire de meme phi
-  /*!    perpendicular vector, with equal phi */
+  //! return the perpendicular vector, with equal phi 
   virtual Vector3d VperpPhi() const;
 
   // vecteur perpendiculaire de meme theta
-  /*!    perpendicular vector with equal theta */
+  //! return the perpendicular vector, with equal theta 
   virtual Vector3d VperpTheta() const;
 
@@ -84 +111 @@
 
   // rotations d'Euler
-  /*!    Euler's rotations */
-  // rotations d Euler
+  //! Perform   Euler's rotations 
   virtual Vector3d Euler(double, double, double) const;
 
   // rotation inverse
-  /*!    inverse rotation */
+  //! perform   inverse Euler rotation 
   Vector3d InvEuler(double, double, double) const;
 
   // rotation d'angle phi autour d'un axe omega (regle du tire-bouchon)
-  /*!    rotation of angle phi around an axis omega (Maxwell's rule) */
+  //! perform rotation of angle phi around an axis omega (Maxwell's rule) 
   Vector3d Rotate(const Vector3d& omega,double phi) const;
 
@@ -148 +174 @@
 }
 
+
 } // namespace SOPHYA