Maple code to compute r(angular coords) for points on offset ellipsoid

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@562 f88db872-0e4f-0410-b76b-b9085cfa78c5

1 files changed, 28 insertions, 0 deletions

diff --git a/src/gr/ellipsoid.c b/src/gr/ellipsoid.c
new file mode 100644
index 0000000..b600b76
--- /dev/null
+++ b/src/gr/ellipsoid.c
@@ -0,0 +1,28 @@
+fp t1, t2, t3, t5, t6, t7, t9, t10, t12, t28;
+fp t30, t33, t35, t36, t40, t42, t43, t48, t49, t52;
+fp t55;
+      t1 = a*a;
+      t2 = b*b;
+      t3 = t1*t2;
+      t5 = t3*0.5772156649015329*CW;
+      t6 = c*c;
+      t7 = t1*t6;
+      t9 = t7*beta*BV;
+      t10 = t2*t6;
+      t12 = t10*alpha*AU;
+      t28 = alpha*alpha;
+      t30 = CW*CW;
+      t33 = BV*BV;
+      t35 = t10*t28;
+      t36 = beta*beta;
+      t40 = AU*AU;
+      t42 = t7*t36;
+      t43 = 0.5772156649015329*0.5772156649015329;
+      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;
+      t52 = sqrt(-t3*t6*t49);
+      t55 = 1/(t48+t42+t35);
+      r_plus = (t5+t9+t12+t52)*t55;
+      r_minus = (t5+t9+t12-t52)*t55;