split out X_ud and X_udd coeffs into new coeffs.* files

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

1 files changed, 20 insertions, 3 deletions

diff --git a/src/maple/coeffs.maple b/src/maple/coeffs.maple
index b45114b..de07b8e 100644
--- a/src/maple/coeffs.maple
+++ b/src/maple/coeffs.maple
@@ -1,10 +1,27 @@
 # coeffs.maple - set up coordinate coefficients
 # $Id$
 #
-# setup_X_ud
-# setup_X_udd
+# setup_coeff_gfas - create the coordinate-transformation coefficient gfas
+# setup_coeffs - compute the coordinate-transformation coefficients themselves
 #
 
-################################################################################
+###############################################################################
 
+#
+# This function sets up the coordinate-transformation gfas.
+# At present these are only inert -- Maple doesn't worry about
+# their fnd expressions (ugly trig functions).
+#
+setup_coeff_gfas :=
+proc()
+global
+@include "coords.minc",
+@include "coeffs.minc";
+
+# 1st ($X^u_i$) and 2nd ($X^u_{ij}$) partial derivs
+# of (r,rho,sigma) wrt (xx,yy,zz)
+make_gfa('X_ud', {inert}, [A_min..N, 1..N], none);
+make_gfa('X_udd', {inert}, [A_min..N, 1..N, 1..N], symmetric3_23);
 
+NULL;
+end proc;