Changeset 1177 in Sophya

Timestamp:

Sep 4, 2000, 12:15:57 PM (25 years ago)

Author:

ansari

Message:

bug constructeur Vector3d(theta, phi) corrige.

dominique

File:

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

trunk/SophyaLib/SkyMap/vector3d.cc

@@ -43 +43 @@
 //-- 
 {
-  _theta=mod(theta,M_PI); // dans [0;pi]
+ // _theta=mod(theta,M_PI); // dans [0;pi]
+ // Version precedente fausse: _theta=M_PI est valide. Or mod(M_PI,M_PI)=0!
+ // De plus theta>pi ou <0 n'a pas de sens. Dominique Yvon
+        if( (theta<0.) || (theta>M_PI) )
+        {       string exmsg = "Wrong initialisation of theta in Vector3d::Vector3d(double theta, double phi)";
+                throw( ParmError(exmsg) );
+        }
+  _theta=theta; // dans [0;pi]
   _phi=mod(phi,pi2); // dans [0;2pi]
   ThetaPhi2xyz();