Add Jacobian gridfunctions to GenericFD

Storage is disabled by default.  These are here so that benchmarks and tests can be run using multipatch thorns without having to have another thorn to provide the Jacobians.  We should probably have a function to initialize them as well.

Auxiliary/Cactus/KrancNumericalTools/GenericFD/interface.ccl

@@ -45,3 +45,15 @@ USES FUNCTION MultiPatch_GetBoundarySpecification
 CCTK_INT FUNCTION \
     SymmetryTableHandleForGrid (CCTK_POINTER_TO_CONST IN cctkGH)
 REQUIRES FUNCTION SymmetryTableHandleForGrid
+
+CCTK_REAL jacobian type=GF timelevels=1 tags='checkpoint="no" prolongation="none"'
+{
+  J11, J12, J13, J21, J22, J23, J31, J32, J33
+} "Jacobian J_(ik) = da^i/dx^k of the coordinate transformation (a: local, x: global)"
+
+CCTK_REAL jacobian2 type=GF timelevels=1 tags='checkpoint="no" prolongation="none"'
+{
+  dJ111, dJ112, dJ113, dJ122, dJ123, dJ133,
+  dJ211, dJ212, dJ213, dJ222, dJ223, dJ233,
+  dJ311, dJ312, dJ313, dJ322, dJ323, dJ333,
+} "Derivative of the Jacobian  d^2[local]/d[global]^2"

Auxiliary/Cactus/KrancNumericalTools/GenericFD/param.ccl

@@ -47,3 +47,7 @@ CCTK_INT  jacobian_identity_map "Map number on which the Jacobian should not be
 {
   -1:* :: "Any integer"
 } -1
+
+CCTK_BOOLEAN allocate_genericfd_jacobian "Allocate memory for the GenericFD Jacobian (used for benchmarks)"
+{
+} "no"

Auxiliary/Cactus/KrancNumericalTools/GenericFD/schedule.ccl

@@ -3,3 +3,8 @@
 # author: S. Husa
 
 # $Header$ 
+
+if (allocate_genericfd_jacobian)
+{
+  STORAGE: jacobian, jacobian2
+}