diff options
| author | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2002-04-22 15:10:32 +0000 |
|---|---|---|
| committer | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2002-04-22 15:10:32 +0000 |
| commit | 74ef815aab12c0f9602b22b4f0542ca7c3a962be (patch) | |
| tree | d6ab428f9b7bf9f14c19f981f929a08087bad310 /src/gr/ellipsoid.maple | |
| parent | ca337a6176d82455b04e3552079f389852f49813 (diff) | |
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
Diffstat (limited to 'src/gr/ellipsoid.maple')
| -rw-r--r-- | src/gr/ellipsoid.maple | 7 |
1 files changed, 4 insertions, 3 deletions
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"; |