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

5 files changed, 33 insertions, 29 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;
diff --git a/src/maple/coeffs.minc b/src/maple/coeffs.minc
new file mode 100644
index 0000000..c7ff962
--- /dev/null
+++ b/src/maple/coeffs.minc
@@ -0,0 +1,7 @@
+# coeffs.minc -- global variables for coordinate-transformation coefficients
+# $Id$
+
+# coordinate-transformation coefficient gfas
+X_ud, X_ud__fnd,
+X_udd, X_udd_fnd
+		# no comma
diff --git a/src/maple/coords.maple b/src/maple/coords.maple
index 4e8953d..24404bf 100644
--- a/src/maple/coords.maple
+++ b/src/maple/coords.maple
@@ -3,7 +3,6 @@
 
 #
 # setup_coords - set up coordinates
-# setup_coord_gfas - setup up coordinate-transformation gfas
 #
 
 ###############################################################################
@@ -48,21 +47,3 @@ x_all_list := [op(x_xyz_list), op(x_rrs_list)];
 
 NULL;
 end proc;
-
-###############################################################################
-
-#
-# This function sets up the coordinate-transformation gfas
-#
-setup_coord_gfas :=
-proc()
-global
-@include "coords.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,fnd}, [A_min..N, 1..N], none);
-make_gfa('X_udd', {inert,fnd}, [A_min..N, 1..N, 1..N], symmetric3_23);
-
-NULL;
-end proc;
diff --git a/src/maple/coords.minc b/src/maple/coords.minc
index 532bfa0..6ab0c9d 100644
--- a/src/maple/coords.minc
+++ b/src/maple/coords.minc
@@ -10,9 +10,5 @@ x_xyz, x_xyz_list,	# array and list of (xx,yy,zz) coordinates
 x_rrs, x_rrs_list,	# array and list of (r,rho,sigma) coordinates
 r, rho, sigma,		# individual patch coordinates
 xx, yy, zz,		# Cartesian coordinates
-x_all_list,		# list of all coordinates
-
-# coordinate-transformation gfas
-X_ud, X_ud__fnd,
-X_udd, X_udd_fnd
-		# no comma
+x_all_list		# list of all coordinates
+	# no comma
diff --git a/src/maple/setup.maple b/src/maple/setup.maple
index 8602d4e..ef5b7a3 100644
--- a/src/maple/setup.maple
+++ b/src/maple/setup.maple
@@ -21,4 +21,7 @@ setup_coords();
 read "../maple/Diff.mm";
 
 read "../maple/gfa.mm";
-setup_coord_gfas();
+
+read "../maple/coeffs.mm";
+setup_coeff_gfas();
+setup_coeffs();