From 74ef815aab12c0f9602b22b4f0542ca7c3a962be Mon Sep 17 00:00:00 2001 From: jthorn Date: Mon, 22 Apr 2002 15:10:32 +0000 Subject: change names alpha --> xcos, beta --> ycos, gamma --> zcos because maple predefines gamma = Euler's constant git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@575 f88db872-0e4f-0410-b76b-b9085cfa78c5 --- src/gr/ellipsoid.maple | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) (limited to 'src/gr/ellipsoid.maple') diff --git a/src/gr/ellipsoid.maple b/src/gr/ellipsoid.maple index c9b0ae4..0c485fa 100644 --- a/src/gr/ellipsoid.maple +++ b/src/gr/ellipsoid.maple @@ -6,10 +6,11 @@ # angular coordinate system has center (U,V,W) # # direction cosines wrt angular coordinate center are (alpha,beta,gamma) -# i.e. a point has coordinates (U+alpha*r, V+beta*r, W+gamma*r) +# but Maple predefines gamma = Euler's constant, so we use (xcos,ycos,zcos) +# instead, i.e. a point has coordinates (U+xcos*r, V+ycos*r, W+zcos*r) # # then the equation of the ellipsoid is -# (U+alpha*r - A)^2 (V+beta*r - B)^2 (W+gamma*r - C)^2 +# (U+xcos*r - A)^2 (V+ycos*r - B)^2 (W+zcos*r - C)^2 # ----------------- + ---------------- + ----------------- = 1 # a^2 b^2 c^2 # @@ -18,7 +19,7 @@ # BV = B - V # CW = C - W # -eqn := (alpha*r - AU)^2/a^2 + (beta*r - BV)^2/b^2 + (gamma*r - CW)^2/c^2 = 1; +eqn := (xcos*r - AU)^2/a^2 + (ycos*r - BV)^2/b^2 + (zcos*r - CW)^2/c^2 = 1; read "../maple/util.mm"; read "../maple/codegen2.mm"; -- cgit v1.2.3