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, 33 insertions, 31 deletions

diff --git a/src/gr/maple.log b/src/gr/maple.log
index 7732b11..2e5ddd0 100644
--- a/src/gr/maple.log
+++ b/src/gr/maple.log
@@ -634,7 +634,8 @@ HA, HA__fnd, HB, HB__fnd, HC, HC__fnd, HD, HD__fnd;
     assert_fnd_exists(g_dd);
     assert_fnd_exists(g_uu, fnd);
     g_uu__fnd := linalg[inverse](g_dd);
-    if cg_flag then codegen2(g_uu__fnd, 'g_uu', "cg/inverse_metric.c")
+    if cg_flag then
+        codegen2(g_uu__fnd, 'g_uu', "../gr.cg/inverse_metric.c")
     end if;
     NULL
 end proc
@@ -661,7 +662,7 @@ HA, HA__fnd, HB, HB__fnd, HC, HC__fnd, HD, HD__fnd;
         end do
     end do;
     if cg_flag then codegen2([K__fnd, K_uu__fnd], ['K', 'K_uu'],
-        "cg/extrinsic_curvature_trace_raise.c")
+        "../gr.cg/extrinsic_curvature_trace_raise.c")
     end if;
     NULL
 end proc
@@ -692,7 +693,7 @@ HA, HA__fnd, HB, HB__fnd, HC, HC__fnd, HD, HD__fnd;
         end do
     end do;
     if cg_flag then codegen2(partial_d_g_uu__fnd, 'partial_d_g_uu',
-        "cg/inverse_metric_gradient.c")
+        "../gr.cg/inverse_metric_gradient.c")
     end if;
     NULL
 end proc
@@ -714,7 +715,7 @@ HA, HA__fnd, HB, HB__fnd, HC, HC__fnd, HD, HD__fnd;
         'g_uu[i, j]*partial_d_g_dd[k, i, j]', 'i' = 1 .. N, 'j' = 1 .. N))
     end do;
     if cg_flag then codegen2(partial_d_ln_sqrt_g__fnd,
-        'partial_d_ln_sqrt_g', "cg/metric_det_gradient.c")
+        'partial_d_ln_sqrt_g', "../gr.cg/metric_det_gradient.c")
     end if;
     NULL
 end proc
@@ -807,7 +808,8 @@ HA, HA__fnd, HB, HB__fnd, HC, HC__fnd, HD, HD__fnd;
     H__fnd :=
         HA__fnd/HD__fnd^(3/2) + HB__fnd/HD__fnd^(1/2) + HC__fnd/HD__fnd - K
         ;
-    if cg_flag then codegen2(H__fnd, 'H', "cg/horizon_function.c") end if;
+    if cg_flag then codegen2(H__fnd, 'H', "../gr.cg/horizon_function.c")
+    end if;
     NULL
 end proc
 
@@ -817,7 +819,7 @@ end proc
 
 > auxiliary(true);
 inverse_metric...
-codegen2(g_uu) --> "cg/inverse_metric.c"
+codegen2(g_uu) --> "../gr.cg/inverse_metric.c"
       --> `codegen2/input`
    convert --> equation list
       --> `codegen2/eqnlist`
@@ -826,7 +828,7 @@ codegen2(g_uu) --> "cg/inverse_metric.c"
    find temporary variables
       --> `codegen2/temps`
    convert Diff(expr,rho,sigma) --> PARTIAL_RHO_SIGMA(expr) etc
-bytes used=1000208, alloc=917336, time=0.10
+bytes used=1000148, alloc=917336, time=0.12
       --> `codegen2/fix_Diff`
    convert R_dd[2,3] --> R_dd_23 etc
       --> `codegen2/unindex`
@@ -834,12 +836,12 @@ bytes used=1000208, alloc=917336, time=0.10
       --> `codegen2/fix_rationals`
    writing C code
 extrinsic_curvature_trace_raise...
-codegen2([K, K_uu]) --> "cg/extrinsic_curvature_trace_raise.c"
+codegen2([K, K_uu]) --> "../gr.cg/extrinsic_curvature_trace_raise.c"
       --> `codegen2/input`
    convert --> equation list
       --> `codegen2/eqnlist`
    optimizing computation sequence
-bytes used=2000708, alloc=1441528, time=0.15
+bytes used=2000796, alloc=1441528, time=0.22
       --> `codegen2/optimize`
    find temporary variables
       --> `codegen2/temps`
@@ -847,39 +849,39 @@ bytes used=2000708, alloc=1441528, time=0.15
       --> `codegen2/fix_Diff`
    convert R_dd[2,3] --> R_dd_23 etc
       --> `codegen2/unindex`
-bytes used=3000940, alloc=1572576, time=0.21
+bytes used=3001776, alloc=1572576, time=0.28
    convert p/q --> RATIONAL(p/q)
       --> `codegen2/fix_rationals`
    writing C code
 > metric_derivs(true);
 inverse_metric_gradient...
-bytes used=4001284, alloc=1638100, time=0.26
-codegen2(partial_d_g_uu) --> "cg/inverse_metric_gradient.c"
+bytes used=4002196, alloc=1638100, time=0.36
+codegen2(partial_d_g_uu) --> "../gr.cg/inverse_metric_gradient.c"
       --> `codegen2/input`
    convert --> equation list
       --> `codegen2/eqnlist`
    optimizing computation sequence
-bytes used=5002720, alloc=1638100, time=0.35
+bytes used=5002380, alloc=1638100, time=0.43
       --> `codegen2/optimize`
    find temporary variables
       --> `codegen2/temps`
    convert Diff(expr,rho,sigma) --> PARTIAL_RHO_SIGMA(expr) etc
       --> `codegen2/fix_Diff`
-bytes used=6002996, alloc=1703624, time=0.41
+bytes used=6004052, alloc=1703624, time=0.50
    convert R_dd[2,3] --> R_dd_23 etc
       --> `codegen2/unindex`
-bytes used=7003240, alloc=1703624, time=0.46
+bytes used=7004228, alloc=1703624, time=0.57
    convert p/q --> RATIONAL(p/q)
       --> `codegen2/fix_rationals`
    writing C code
-bytes used=8003492, alloc=1703624, time=0.49
+bytes used=8004548, alloc=1703624, time=0.63
 metric_det_gradient...
-codegen2(partial_d_ln_sqrt_g) --> "cg/metric_det_gradient.c"
+codegen2(partial_d_ln_sqrt_g) --> "../gr.cg/metric_det_gradient.c"
       --> `codegen2/input`
    convert --> equation list
       --> `codegen2/eqnlist`
    optimizing computation sequence
-bytes used=9003660, alloc=1703624, time=0.60
+bytes used=9004804, alloc=1703624, time=0.74
       --> `codegen2/optimize`
    find temporary variables
       --> `codegen2/temps`
@@ -894,40 +896,40 @@ codegen/C/expression:   Unknown function:   RATIONAL   will be left as is.
 > horizon(true);
 non_unit_normal...
 non_unit_normal_deriv...
-bytes used=10003824, alloc=1769148, time=0.66
+bytes used=10004956, alloc=1769148, time=0.85
 horizon_function...
-bytes used=11004072, alloc=1834672, time=0.73
-codegen2(H) --> "cg/horizon_function.c"
+bytes used=11005112, alloc=1834672, time=0.93
+codegen2(H) --> "../gr.cg/horizon_function.c"
       --> `codegen2/input`
    convert --> equation list
       --> `codegen2/eqnlist`
    optimizing computation sequence
-bytes used=12004376, alloc=1834672, time=0.79
-bytes used=13004892, alloc=2031244, time=0.87
-bytes used=14005160, alloc=2031244, time=0.96
+bytes used=12005688, alloc=1834672, time=1.01
+bytes used=13005864, alloc=2031244, time=1.12
+bytes used=14006148, alloc=2031244, time=1.21
       --> `codegen2/optimize`
    find temporary variables
       --> `codegen2/temps`
    convert Diff(expr,rho,sigma) --> PARTIAL_RHO_SIGMA(expr) etc
-bytes used=15005740, alloc=2096768, time=1.01
+bytes used=15006312, alloc=2096768, time=1.28
       --> `codegen2/fix_Diff`
    convert R_dd[2,3] --> R_dd_23 etc
-bytes used=16006176, alloc=2096768, time=1.07
+bytes used=16006476, alloc=2096768, time=1.35
       --> `codegen2/unindex`
-bytes used=17006396, alloc=2096768, time=1.13
+bytes used=17006632, alloc=2096768, time=1.41
    convert p/q --> RATIONAL(p/q)
-bytes used=18006808, alloc=2096768, time=1.18
+bytes used=18006900, alloc=2096768, time=1.48
       --> `codegen2/fix_rationals`
    writing C code
 codegen/C/expression:   Unknown function:   PARTIAL_RHO   will be left as is.
 codegen/C/expression:   Unknown function:   PARTIAL_SIGMA   will be left as is.
-bytes used=19007204, alloc=2096768, time=1.22
+bytes used=19007116, alloc=2096768, time=1.54
 codegen/C/expression:   Unknown function:   PARTIAL_RHO_RHO   
 will be left as is.
 codegen/C/expression:   Unknown function:   PARTIAL_RHO_SIGMA   
 will be left as is.
 codegen/C/expression:   Unknown function:   PARTIAL_SIGMA_SIGMA   
 will be left as is.
-bytes used=20007392, alloc=2096768, time=1.38
+bytes used=20007280, alloc=2096768, time=1.66
 > quit
-bytes used=20239428, alloc=2096768, time=1.40
+bytes used=20238456, alloc=2096768, time=1.70