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

1 files changed, 10 insertions, 10 deletions

diff --git a/src/gr/ellipsoid.c b/src/gr/ellipsoid.c
index b600b76..49a6457 100644
--- a/src/gr/ellipsoid.c
+++ b/src/gr/ellipsoid.c
@@ -4,25 +4,25 @@ fp t55;
       t1 = a*a;
       t2 = b*b;
       t3 = t1*t2;
-      t5 = t3*0.5772156649015329*CW;
+      t5 = t3*zcos*CW;
       t6 = c*c;
       t7 = t1*t6;
-      t9 = t7*beta*BV;
+      t9 = t7*ycos*BV;
       t10 = t2*t6;
-      t12 = t10*alpha*AU;
-      t28 = alpha*alpha;
+      t12 = t10*xcos*AU;
+      t28 = xcos*xcos;
       t30 = CW*CW;
       t33 = BV*BV;
       t35 = t10*t28;
-      t36 = beta*beta;
+      t36 = ycos*ycos;
       t40 = AU*AU;
       t42 = t7*t36;
-      t43 = 0.5772156649015329*0.5772156649015329;
+      t43 = zcos*zcos;
       t48 = t3*t43;
-      t49 = -2.0*t1*0.5772156649015329*CW*beta*BV-2.0*t2*0.5772156649015329*CW*
-alpha*AU-2.0*t6*beta*BV*alpha*AU+t2*t28*t30+t6*t28*t33-t35+t1*t36*t30+t6*t36*
-t40-t42+t1*t43*t33+t2*t43*t40-t48;
+      t49 = -2.0*t1*zcos*CW*ycos*BV-2.0*t2*zcos*CW*xcos*AU-2.0*t6*ycos*BV*xcos*
+AU+t2*t28*t30+t6*t28*t33-t35+t1*t36*t30+t6*t36*t40-t42+t1*t43*t33+t2*t43*t40-
+t48;
       t52 = sqrt(-t3*t6*t49);
-      t55 = 1/(t48+t42+t35);
+      t55 = 1/(t35+t42+t48);
       r_plus = (t5+t9+t12+t52)*t55;
       r_minus = (t5+t9+t12-t52)*t55;