start adding Jadcobian support --> not working yet

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

1 files changed, 28 insertions, 10 deletions

diff --git a/src/gr/setup_gr_gfas.maple b/src/gr/setup_gr_gfas.maple
index ede82ca..dc079cb 100644
--- a/src/gr/setup_gr_gfas.maple
+++ b/src/gr/setup_gr_gfas.maple
@@ -22,7 +22,7 @@ make_gfa('K_dd', {inert}, [1..N, 1..N], symmetric);
 # index-raised and contracted metric and extrinsic curvature
 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);
+make_gfa('K', {inert,fnd}, [], none);
 
 # xyz partial derivatives of metric
 # ... as far as Maple is concerned, these are indeed separate gridfns;
@@ -37,22 +37,40 @@ make_gfa('partial_d_ln_sqrt_g', {inert,fnd}, [1..N], none);
 make_gfa('partial_d_g_uu', {inert,fnd}, [1..N, 1..N, 1..N], symmetric3_23);
 
 # radius of horizon
-make_gfa('h', {inert, fnd}, [], none);
+make_gfa('h', {inert,fnd}, [], none);
 
 # outward-pointing *non*-unit normal (covector) 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);
 
-# 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);
+make_gfa('HA', {inert,fnd}, [], none);
+make_gfa('HB', {inert,fnd}, [], none);
+make_gfa('HC', {inert,fnd}, [], none);
+make_gfa('HD', {inert,fnd}, [], none);
+
+# LHS of apparent horizon equation
+make_gfa('H', {inert,fnd}, [], none);
+
+# Jacobian coefficients for H[ABCD]
+# these are the partial derivatives of H[ABCD]
+# ... wrt Diff(h,x_rs[x])
+make_gfa('partial_HA_wrt_partial_d_h', {inert,fnd}, [1..N_ang], none);
+make_gfa('partial_HB_wrt_partial_d_h', {inert,fnd}, [1..N_ang], none);
+make_gfa('partial_HC_wrt_partial_d_h', {inert,fnd}, [1..N_ang], none);
+make_gfa('partial_HD_wrt_partial_d_h', {inert,fnd}, [1..N_ang], none);
+# ... wrt Diff(h,x_rs[x],x_rs[y])
+make_gfa('partial_HA_wrt_partial_dd_h', {inert,fnd},
+	 [1..N_ang, 1..N_ang], symmetric);
+make_gfa('partial_HB_wrt_partial_dd_h', {inert,fnd},
+	 [1..N_ang, 1..N_ang], symmetric);
+
+# Jacobian coefficients for H itself
+make_gfa('partial_H_wrt_partial_d_h', {inert,fnd}, [1..N_ang], none);
+make_gfa('partial_H_wrt_partial_dd_h', {inert,fnd},
+	 [1..N_ang, 1..N_ang], symmetric);
 
 NULL;
 end proc;
@@ -75,7 +93,7 @@ end proc;
 # var_seq = (in) (varargs) An expression sequence of the variables to
 #			   differentiate with respect to.
 #
-`Diff/gridfn` :=
+`Diff/gridfn2` :=
 proc(operand)			# varargs
 option remember;		# performance optimization
 global