1 files changed, 49 insertions, 0 deletions
diff --git a/src/gr/setup_gfas.maple b/src/gr/setup_gfas.maple
new file mode 100644
index 0000000..5da990b
--- /dev/null
+++ b/src/gr/setup_gfas.maple
@@ -0,0 +1,49 @@
+# Maple code to set up all gridfn arrays
+# $Id$
+
+#
+# setup_gr_gfas - setup all the GR gfas
+#
+
+################################################################################
+
+setup_gr_gfas :=
+proc()
+global
+       @include "../gr/gr_gfas.minc";
+
+# basic metric & extrinsic curvature fields
+make_gfa('g_dd', {inert}, [1..N, 1..N], symmetric);
+make_gfa('K_dd', {inert}, [1..N, 1..N], symmetric);
+
+make_gfa('g_uu', {inert,fnd}, [1..N, 1..N], symmetric);
+make_gfa('K_uu', {inert,fnd}, [1..N, 1..N], symmetric);
+
+make_gfa('K', {inert, fnd}, [], none);
+
+# partial derivatives of metric determinant
+make_gfa('partial_d_ln_sqrt_g', {inert,fnd}, [1..N], none);
+
+# radius of horizon
+make_gfa('h', {inert, fnd}, [], none);
+
+# outward-pointing *non*-unit normal to horizon
+# and it's xyz-coordinate partial derivatives
+make_gfa('s_d', {inert,fnd}, [1..N], none);
+make_gfa('partial_d_s_d', {inert,fnd}, [1..N, 1..N], none);
+
+# outward-pointing unit normal to horizon
+make_gfa('n_u', {inert,fnd}, [1..N], none);
+
+# LHS of apparent horizon equation
+make_gfa('H', {inert, fnd}, [], none);
+
+# subexpressions for computing LHS of apparent horizon equation
+# ... these are defined by (15) in my 1996 apparent horizon finding paper
+make_gfa('HA', {inert, fnd}, [], none);
+make_gfa('HB', {inert, fnd}, [], none);
+make_gfa('HC', {inert, fnd}, [], none);
+make_gfa('HD', {inert, fnd}, [], none);
+
+NULL;
+end proc;