diff options
Diffstat (limited to 'src/gr/maple.log')
-rw-r--r-- | src/gr/maple.log | 265 |
1 files changed, 151 insertions, 114 deletions
diff --git a/src/gr/maple.log b/src/gr/maple.log index 54c762f..1a91a6c 100644 --- a/src/gr/maple.log +++ b/src/gr/maple.log @@ -1,10 +1,11 @@ +RedHat v <9 or other Linux present, starting standard mode... |\^/| Maple 7 (IBM INTEL LINUX) ._|\| |/|_. Copyright (c) 2001 by Waterloo Maple Inc. \ MAPLE / All rights reserved. Maple is a registered trademark of <____ ____> Waterloo Maple Inc. | Type ? for help. # top-level Maple file to read/run all code in this directory -# $Header: /numrelcvs/AEIDevelopment/AHFinderDirect/src/gr/doit.maple,v 1.5 2002/09/13 14:12:18 jthorn Exp $ +# $Header: /numrelcvs/AEIThorns/AHFinderDirect/src/gr/doit.maple,v 1.5 2002/09/13 14:12:18 jthorn Exp $ > > read "../maple/setup.mm"; msum := proc(fn::algebraic) @@ -205,19 +206,22 @@ rho, sigma, y_rs, y_rs_list, y_rs_set, xy_all_list, xy_all_set, inert, none, fnd, symmetric3_23, X_ud, X_ud__fnd, X_udd, X_udd__fnd, g_dd, K_dd, g_uu, g_uu__fnd, K_uu, K_uu__fnd, K, K__fnd, partial_d_g_dd, partial_d_ln_sqrt_g, partial_d_ln_sqrt_g__fnd, partial_d_g_uu, partial_d_g_uu__fnd, h, h__fnd, -s_d, s_d__fnd, partial_d_s_d, partial_d_s_d__fnd, n_u, n_u__fnd, HA, -HA__fnd, HB, HB__fnd, HC, HC__fnd, HD, HD__fnd, H, H__fnd, partial_d_HA, -partial_d_HA__fnd, partial_d_HB, partial_d_HB__fnd, partial_d_HC, -partial_d_HC__fnd, partial_d_HD, partial_d_HD__fnd, partial_d_H, -partial_d_H__fnd, partial_HA_wrt_partial_d_h, -partial_HA_wrt_partial_d_h__fnd, partial_HB_wrt_partial_d_h, -partial_HB_wrt_partial_d_h__fnd, partial_HC_wrt_partial_d_h, -partial_HC_wrt_partial_d_h__fnd, partial_HD_wrt_partial_d_h, -partial_HD_wrt_partial_d_h__fnd, partial_HA_wrt_partial_dd_h, -partial_HA_wrt_partial_dd_h__fnd, partial_HB_wrt_partial_dd_h, -partial_HB_wrt_partial_dd_h__fnd, partial_H_wrt_partial_d_h, -partial_H_wrt_partial_d_h__fnd, partial_H_wrt_partial_dd_h, -partial_H_wrt_partial_dd_h__fnd; +s_d, s_d__fnd, partial_d_s_d, partial_d_s_d__fnd, n_u, n_u__fnd, Theta_A, +Theta_A__fnd, Theta_B, Theta_B__fnd, Theta_C, Theta_C__fnd, Theta_D, +Theta_D__fnd, Theta, Theta__fnd, partial_d_Theta_A, partial_d_Theta_A__fnd, +partial_d_Theta_B, partial_d_Theta_B__fnd, partial_d_Theta_C, +partial_d_Theta_C__fnd, partial_d_Theta_D, partial_d_Theta_D__fnd, +partial_d_Theta_, partial_d_Theta__fnd, partial_Theta_A_wrt_partial_d_h, +partial_Theta_A_wrt_partial_d_h__fnd, partial_Theta_B_wrt_partial_d_h, +partial_Theta_B_wrt_partial_d_h__fnd, partial_Theta_C_wrt_partial_d_h, +partial_Theta_C_wrt_partial_d_h__fnd, partial_Theta_D_wrt_partial_d_h, +partial_Theta_D_wrt_partial_d_h__fnd, partial_Theta_A_wrt_partial_dd_h, +partial_Theta_A_wrt_partial_dd_h__fnd, partial_Theta_B_wrt_partial_dd_h, +partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_X_wrt_partial_d_h, +partial_Theta_X_wrt_partial_d_h__fnd, partial_Theta_Y_wrt_partial_d_h, +partial_Theta_Y_wrt_partial_d_h__fnd, partial_Theta_X_wrt_partial_dd_h, +partial_Theta_X_wrt_partial_dd_h__fnd, partial_Theta_Y_wrt_partial_dd_h, +partial_Theta_Y_wrt_partial_dd_h__fnd; option remember; var_list := [args[2 .. nargs]]; if type(operand, indexed) and op(0, operand) = 'X_ud' and @@ -597,9 +601,11 @@ partial_Theta_B_wrt_partial_d_h__fnd, partial_Theta_C_wrt_partial_d_h, partial_Theta_C_wrt_partial_d_h__fnd, partial_Theta_D_wrt_partial_d_h, partial_Theta_D_wrt_partial_d_h__fnd, partial_Theta_A_wrt_partial_dd_h, partial_Theta_A_wrt_partial_dd_h__fnd, partial_Theta_B_wrt_partial_dd_h, -partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_wrt_partial_d_h, -partial_Theta_wrt_partial_d_h__fnd, partial_Theta_wrt_partial_dd_h, -partial_Theta_wrt_partial_dd_h__fnd; +partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_X_wrt_partial_d_h, +partial_Theta_X_wrt_partial_d_h__fnd, partial_Theta_Y_wrt_partial_d_h, +partial_Theta_Y_wrt_partial_d_h__fnd, partial_Theta_X_wrt_partial_dd_h, +partial_Theta_X_wrt_partial_dd_h__fnd, partial_Theta_Y_wrt_partial_dd_h, +partial_Theta_Y_wrt_partial_dd_h__fnd; 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); @@ -635,9 +641,13 @@ partial_Theta_wrt_partial_dd_h__fnd; [1 .. N_ang, 1 .. N_ang], symmetric); make_gfa('partial_Theta_B_wrt_partial_dd_h', {inert, fnd}, [1 .. N_ang, 1 .. N_ang], symmetric); - make_gfa('partial_Theta_wrt_partial_d_h', {inert, fnd}, [1 .. N_ang], + make_gfa('partial_Theta_X_wrt_partial_d_h', {inert, fnd}, [1 .. N_ang], + none); + make_gfa('partial_Theta_Y_wrt_partial_d_h', {inert, fnd}, [1 .. N_ang], none); - make_gfa('partial_Theta_wrt_partial_dd_h', {inert, fnd}, + make_gfa('partial_Theta_X_wrt_partial_dd_h', {inert, fnd}, + [1 .. N_ang, 1 .. N_ang], symmetric); + make_gfa('partial_Theta_Y_wrt_partial_dd_h', {inert, fnd}, [1 .. N_ang, 1 .. N_ang], symmetric); NULL end proc @@ -660,9 +670,11 @@ partial_Theta_B_wrt_partial_d_h__fnd, partial_Theta_C_wrt_partial_d_h, partial_Theta_C_wrt_partial_d_h__fnd, partial_Theta_D_wrt_partial_d_h, partial_Theta_D_wrt_partial_d_h__fnd, partial_Theta_A_wrt_partial_dd_h, partial_Theta_A_wrt_partial_dd_h__fnd, partial_Theta_B_wrt_partial_dd_h, -partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_wrt_partial_d_h, -partial_Theta_wrt_partial_d_h__fnd, partial_Theta_wrt_partial_dd_h, -partial_Theta_wrt_partial_dd_h__fnd; +partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_X_wrt_partial_d_h, +partial_Theta_X_wrt_partial_d_h__fnd, partial_Theta_Y_wrt_partial_d_h, +partial_Theta_Y_wrt_partial_d_h__fnd, partial_Theta_X_wrt_partial_dd_h, +partial_Theta_X_wrt_partial_dd_h__fnd, partial_Theta_Y_wrt_partial_dd_h, +partial_Theta_Y_wrt_partial_dd_h__fnd; option remember; var_list := [args[2 .. nargs]]; if type(operand, indexed) and op(0, operand) = 'g_dd' and @@ -696,9 +708,11 @@ partial_Theta_B_wrt_partial_d_h__fnd, partial_Theta_C_wrt_partial_d_h, partial_Theta_C_wrt_partial_d_h__fnd, partial_Theta_D_wrt_partial_d_h, partial_Theta_D_wrt_partial_d_h__fnd, partial_Theta_A_wrt_partial_dd_h, partial_Theta_A_wrt_partial_dd_h__fnd, partial_Theta_B_wrt_partial_dd_h, -partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_wrt_partial_d_h, -partial_Theta_wrt_partial_d_h__fnd, partial_Theta_wrt_partial_dd_h, -partial_Theta_wrt_partial_dd_h__fnd; +partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_X_wrt_partial_d_h, +partial_Theta_X_wrt_partial_d_h__fnd, partial_Theta_Y_wrt_partial_d_h, +partial_Theta_Y_wrt_partial_d_h__fnd, partial_Theta_X_wrt_partial_dd_h, +partial_Theta_X_wrt_partial_dd_h__fnd, partial_Theta_Y_wrt_partial_dd_h, +partial_Theta_Y_wrt_partial_dd_h__fnd; printf("%a...\n", procname); assert_fnd_exists(g_dd); assert_fnd_exists(g_uu, fnd); @@ -727,9 +741,11 @@ partial_Theta_B_wrt_partial_d_h__fnd, partial_Theta_C_wrt_partial_d_h, partial_Theta_C_wrt_partial_d_h__fnd, partial_Theta_D_wrt_partial_d_h, partial_Theta_D_wrt_partial_d_h__fnd, partial_Theta_A_wrt_partial_dd_h, partial_Theta_A_wrt_partial_dd_h__fnd, partial_Theta_B_wrt_partial_dd_h, -partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_wrt_partial_d_h, -partial_Theta_wrt_partial_d_h__fnd, partial_Theta_wrt_partial_dd_h, -partial_Theta_wrt_partial_dd_h__fnd; +partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_X_wrt_partial_d_h, +partial_Theta_X_wrt_partial_d_h__fnd, partial_Theta_Y_wrt_partial_d_h, +partial_Theta_Y_wrt_partial_d_h__fnd, partial_Theta_X_wrt_partial_dd_h, +partial_Theta_X_wrt_partial_dd_h__fnd, partial_Theta_Y_wrt_partial_dd_h, +partial_Theta_Y_wrt_partial_dd_h__fnd; printf("%a...\n", procname); assert_fnd_exists(g_uu); assert_fnd_exists(K_dd); @@ -771,9 +787,11 @@ partial_Theta_B_wrt_partial_d_h__fnd, partial_Theta_C_wrt_partial_d_h, partial_Theta_C_wrt_partial_d_h__fnd, partial_Theta_D_wrt_partial_d_h, partial_Theta_D_wrt_partial_d_h__fnd, partial_Theta_A_wrt_partial_dd_h, partial_Theta_A_wrt_partial_dd_h__fnd, partial_Theta_B_wrt_partial_dd_h, -partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_wrt_partial_d_h, -partial_Theta_wrt_partial_d_h__fnd, partial_Theta_wrt_partial_dd_h, -partial_Theta_wrt_partial_dd_h__fnd; +partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_X_wrt_partial_d_h, +partial_Theta_X_wrt_partial_d_h__fnd, partial_Theta_Y_wrt_partial_d_h, +partial_Theta_Y_wrt_partial_d_h__fnd, partial_Theta_X_wrt_partial_dd_h, +partial_Theta_X_wrt_partial_dd_h__fnd, partial_Theta_Y_wrt_partial_dd_h, +partial_Theta_Y_wrt_partial_dd_h__fnd; printf("%a...\n", procname); assert_fnd_exists(g_dd); assert_fnd_exists(g_uu); @@ -809,9 +827,11 @@ partial_Theta_B_wrt_partial_d_h__fnd, partial_Theta_C_wrt_partial_d_h, partial_Theta_C_wrt_partial_d_h__fnd, partial_Theta_D_wrt_partial_d_h, partial_Theta_D_wrt_partial_d_h__fnd, partial_Theta_A_wrt_partial_dd_h, partial_Theta_A_wrt_partial_dd_h__fnd, partial_Theta_B_wrt_partial_dd_h, -partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_wrt_partial_d_h, -partial_Theta_wrt_partial_d_h__fnd, partial_Theta_wrt_partial_dd_h, -partial_Theta_wrt_partial_dd_h__fnd; +partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_X_wrt_partial_d_h, +partial_Theta_X_wrt_partial_d_h__fnd, partial_Theta_Y_wrt_partial_d_h, +partial_Theta_Y_wrt_partial_d_h__fnd, partial_Theta_X_wrt_partial_dd_h, +partial_Theta_X_wrt_partial_dd_h__fnd, partial_Theta_Y_wrt_partial_dd_h, +partial_Theta_Y_wrt_partial_dd_h__fnd; printf("%a...\n", procname); assert_fnd_exists(g_dd); assert_fnd_exists(g_uu); @@ -852,9 +872,11 @@ partial_Theta_B_wrt_partial_d_h__fnd, partial_Theta_C_wrt_partial_d_h, partial_Theta_C_wrt_partial_d_h__fnd, partial_Theta_D_wrt_partial_d_h, partial_Theta_D_wrt_partial_d_h__fnd, partial_Theta_A_wrt_partial_dd_h, partial_Theta_A_wrt_partial_dd_h__fnd, partial_Theta_B_wrt_partial_dd_h, -partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_wrt_partial_d_h, -partial_Theta_wrt_partial_d_h__fnd, partial_Theta_wrt_partial_dd_h, -partial_Theta_wrt_partial_dd_h__fnd; +partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_X_wrt_partial_d_h, +partial_Theta_X_wrt_partial_d_h__fnd, partial_Theta_Y_wrt_partial_d_h, +partial_Theta_Y_wrt_partial_d_h__fnd, partial_Theta_X_wrt_partial_dd_h, +partial_Theta_X_wrt_partial_dd_h__fnd, partial_Theta_Y_wrt_partial_dd_h, +partial_Theta_Y_wrt_partial_dd_h__fnd; printf("%a...\n", procname); assert_fnd_exists(h); assert_fnd_exists(X_ud); @@ -883,9 +905,11 @@ partial_Theta_B_wrt_partial_d_h__fnd, partial_Theta_C_wrt_partial_d_h, partial_Theta_C_wrt_partial_d_h__fnd, partial_Theta_D_wrt_partial_d_h, partial_Theta_D_wrt_partial_d_h__fnd, partial_Theta_A_wrt_partial_dd_h, partial_Theta_A_wrt_partial_dd_h__fnd, partial_Theta_B_wrt_partial_dd_h, -partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_wrt_partial_d_h, -partial_Theta_wrt_partial_d_h__fnd, partial_Theta_wrt_partial_dd_h, -partial_Theta_wrt_partial_dd_h__fnd; +partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_X_wrt_partial_d_h, +partial_Theta_X_wrt_partial_d_h__fnd, partial_Theta_Y_wrt_partial_d_h, +partial_Theta_Y_wrt_partial_d_h__fnd, partial_Theta_X_wrt_partial_dd_h, +partial_Theta_X_wrt_partial_dd_h__fnd, partial_Theta_Y_wrt_partial_dd_h, +partial_Theta_Y_wrt_partial_dd_h__fnd; printf("%a...\n", procname); assert_fnd_exists(h); assert_fnd_exists(X_ud); @@ -923,9 +947,11 @@ partial_Theta_B_wrt_partial_d_h__fnd, partial_Theta_C_wrt_partial_d_h, partial_Theta_C_wrt_partial_d_h__fnd, partial_Theta_D_wrt_partial_d_h, partial_Theta_D_wrt_partial_d_h__fnd, partial_Theta_A_wrt_partial_dd_h, partial_Theta_A_wrt_partial_dd_h__fnd, partial_Theta_B_wrt_partial_dd_h, -partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_wrt_partial_d_h, -partial_Theta_wrt_partial_d_h__fnd, partial_Theta_wrt_partial_dd_h, -partial_Theta_wrt_partial_dd_h__fnd; +partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_X_wrt_partial_d_h, +partial_Theta_X_wrt_partial_d_h__fnd, partial_Theta_Y_wrt_partial_d_h, +partial_Theta_Y_wrt_partial_d_h__fnd, partial_Theta_X_wrt_partial_dd_h, +partial_Theta_X_wrt_partial_dd_h__fnd, partial_Theta_Y_wrt_partial_dd_h, +partial_Theta_Y_wrt_partial_dd_h__fnd; printf("%a...\n", procname); assert_fnd_exists(g_uu); assert_fnd_exists(K_uu); @@ -958,7 +984,7 @@ partial_Theta_wrt_partial_dd_h__fnd; end proc expansion_Jacobian := proc(cg_flag::boolean) -local u, v, temp; +local u, v, temp1, temp2; global delta, N, N_ang, xx, yy, zz, x_xyz, x_xyz_list, x_xyz_set, r, r__fnd, rho, sigma, y_rs, y_rs_list, y_rs_set, xy_all_list, xy_all_set, inert, none, fnd, symmetric3_23, X_ud, X_ud__fnd, X_udd, X_udd__fnd, g_dd, K_dd, g_uu, @@ -975,9 +1001,11 @@ partial_Theta_B_wrt_partial_d_h__fnd, partial_Theta_C_wrt_partial_d_h, partial_Theta_C_wrt_partial_d_h__fnd, partial_Theta_D_wrt_partial_d_h, partial_Theta_D_wrt_partial_d_h__fnd, partial_Theta_A_wrt_partial_dd_h, partial_Theta_A_wrt_partial_dd_h__fnd, partial_Theta_B_wrt_partial_dd_h, -partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_wrt_partial_d_h, -partial_Theta_wrt_partial_d_h__fnd, partial_Theta_wrt_partial_dd_h, -partial_Theta_wrt_partial_dd_h__fnd; +partial_Theta_B_wrt_partial_dd_h__fnd, partial_Theta_X_wrt_partial_d_h, +partial_Theta_X_wrt_partial_d_h__fnd, partial_Theta_Y_wrt_partial_d_h, +partial_Theta_Y_wrt_partial_d_h__fnd, partial_Theta_X_wrt_partial_dd_h, +partial_Theta_X_wrt_partial_dd_h__fnd, partial_Theta_Y_wrt_partial_dd_h, +partial_Theta_Y_wrt_partial_dd_h__fnd; printf("%a...\n", procname); assert_fnd_exists(g_uu); assert_fnd_exists(K_uu); @@ -996,13 +1024,15 @@ partial_Theta_wrt_partial_dd_h__fnd; frontend('diff', [Theta_C__fnd, Diff(h, y_rs[u])]); partial_Theta_D_wrt_partial_d_h__fnd[u] := frontend('diff', [Theta_D__fnd, Diff(h, y_rs[u])]); - temp := 3/2*Theta_A/Theta_D^(5/2) + 1/2*Theta_B/Theta_D^(3/2) - + Theta_C/Theta_D^2; - partial_Theta_wrt_partial_d_h__fnd[u] := + temp1 := 3/2*Theta_A/Theta_D^(5/2) + 1/2*Theta_B/Theta_D^(3/2); + partial_Theta_X_wrt_partial_d_h__fnd[u] := partial_Theta_A_wrt_partial_d_h__fnd[u]/Theta_D^(3/2) + partial_Theta_B_wrt_partial_d_h__fnd[u]/Theta_D^(1/2) - + partial_Theta_C_wrt_partial_d_h__fnd[u]/Theta_D - - partial_Theta_D_wrt_partial_d_h__fnd[u]*temp + - partial_Theta_D_wrt_partial_d_h__fnd[u]*temp1; + temp2 := Theta_C/Theta_D^2; + partial_Theta_Y_wrt_partial_d_h__fnd[u] := + partial_Theta_C_wrt_partial_d_h__fnd[u]/Theta_D + - partial_Theta_D_wrt_partial_d_h__fnd[u]*temp2 end do; for u to N_ang do for v from u to N_ang do partial_Theta_A_wrt_partial_dd_h__fnd[u, v] := @@ -1011,16 +1041,22 @@ partial_Theta_wrt_partial_dd_h__fnd; partial_Theta_B_wrt_partial_dd_h__fnd[u, v] := frontend('diff', [Theta_B__fnd, Diff(h, y_rs[u], y_rs[v])]) ; - partial_Theta_wrt_partial_dd_h__fnd[u, v] := + partial_Theta_X_wrt_partial_dd_h__fnd[u, v] := partial_Theta_A_wrt_partial_dd_h__fnd[u, v]/Theta_D^(3/2) + - partial_Theta_B_wrt_partial_dd_h__fnd[u, v]/Theta_D^(1/2) + partial_Theta_B_wrt_partial_dd_h__fnd[u, v]/Theta_D^(1/2); + partial_Theta_Y_wrt_partial_dd_h__fnd[u, v] := 0 end do end do; - if cg_flag then codegen2([partial_Theta_wrt_partial_d_h__fnd, - partial_Theta_wrt_partial_dd_h__fnd], - ['partial_Theta_wrt_partial_d_h', 'partial_Theta_wrt_partial_dd_h'] - , "../gr.cg/expansion_Jacobian.c") + if cg_flag then codegen2([partial_Theta_X_wrt_partial_d_h__fnd, + partial_Theta_Y_wrt_partial_d_h__fnd, + partial_Theta_X_wrt_partial_dd_h__fnd, + partial_Theta_Y_wrt_partial_dd_h__fnd], [ + 'partial_Theta_X_wrt_partial_d_h', + 'partial_Theta_Y_wrt_partial_d_h', + 'partial_Theta_X_wrt_partial_dd_h', + 'partial_Theta_Y_wrt_partial_dd_h'], + "../gr.cg/expansion_Jacobian.c") end if; NULL end proc @@ -1040,7 +1076,7 @@ codegen2(g_uu) --> "../gr.cg/inverse_metric.c" find temporary variables --> `codegen2/temps` convert Diff(expr,rho,sigma) --> PARTIAL_RHO_SIGMA(expr) etc -bytes used=1000088, alloc=917336, time=0.12 +bytes used=1000776, alloc=917336, time=0.13 --> `codegen2/fix_Diff` convert R_dd[2,3] --> R_dd_23 etc --> `codegen2/unindex` @@ -1053,7 +1089,7 @@ codegen2([K, K_uu]) --> "../gr.cg/extrinsic_curvature_trace_raise.c" convert --> equation list --> `codegen2/eqnlist` optimizing computation sequence -bytes used=2000352, alloc=1441528, time=0.17 +bytes used=2001072, alloc=1441528, time=0.20 --> `codegen2/optimize` find temporary variables --> `codegen2/temps` @@ -1061,39 +1097,39 @@ bytes used=2000352, alloc=1441528, time=0.17 --> `codegen2/fix_Diff` convert R_dd[2,3] --> R_dd_23 etc --> `codegen2/unindex` -bytes used=3000584, alloc=1638100, time=0.22 +bytes used=3001364, alloc=1638100, time=0.26 convert p/q --> RATIONAL(p/q) --> `codegen2/fix_rationals` writing C code > curvature(true); inverse_metric_gradient... -bytes used=4000860, alloc=1703624, time=0.28 +bytes used=4001680, alloc=1703624, time=0.34 codegen2(partial_d_g_uu) --> "../gr.cg/inverse_metric_gradient.c" --> `codegen2/input` convert --> equation list --> `codegen2/eqnlist` optimizing computation sequence -bytes used=5001020, alloc=1703624, time=0.38 +bytes used=5001836, alloc=1703624, time=0.42 --> `codegen2/optimize` find temporary variables --> `codegen2/temps` convert Diff(expr,rho,sigma) --> PARTIAL_RHO_SIGMA(expr) etc -bytes used=6001324, alloc=1769148, time=0.44 +bytes used=6002212, alloc=1769148, time=0.49 --> `codegen2/fix_Diff` convert R_dd[2,3] --> R_dd_23 etc --> `codegen2/unindex` -bytes used=7001528, alloc=1769148, time=0.51 +bytes used=7002736, alloc=1769148, time=0.55 convert p/q --> RATIONAL(p/q) --> `codegen2/fix_rationals` writing C code -bytes used=8001896, alloc=1769148, time=0.57 +bytes used=8003012, alloc=1769148, time=0.62 metric_det_gradient... codegen2(partial_d_ln_sqrt_g) --> "../gr.cg/metric_det_gradient.c" --> `codegen2/input` convert --> equation list --> `codegen2/eqnlist` optimizing computation sequence -bytes used=9002080, alloc=1769148, time=0.64 +bytes used=9003256, alloc=1769148, time=0.75 --> `codegen2/optimize` find temporary variables --> `codegen2/temps` @@ -1108,29 +1144,29 @@ codegen/C/expression: Unknown function: RATIONAL will be left as is. > horizon(true); non_unit_normal... non_unit_normal_deriv... -bytes used=10002316, alloc=1769148, time=0.69 +bytes used=10005164, alloc=1769148, time=0.82 expansion... -bytes used=11002748, alloc=1834672, time=0.77 +bytes used=11005388, alloc=1834672, time=0.89 codegen2([Theta_A, Theta_B, Theta_C, Theta_D]) --> "../gr.cg/expansion.c" --> `codegen2/input` convert --> equation list --> `codegen2/eqnlist` optimizing computation sequence -bytes used=12003104, alloc=1834672, time=0.83 -bytes used=13003356, alloc=2031244, time=0.95 -bytes used=14003520, alloc=2031244, time=1.05 +bytes used=12005892, alloc=1834672, time=0.96 +bytes used=13006112, alloc=2031244, time=1.09 +bytes used=14006296, alloc=2031244, time=1.20 --> `codegen2/optimize` find temporary variables --> `codegen2/temps` convert Diff(expr,rho,sigma) --> PARTIAL_RHO_SIGMA(expr) etc -bytes used=15003896, alloc=2096768, time=1.11 +bytes used=15006844, alloc=2096768, time=1.27 --> `codegen2/fix_Diff` convert R_dd[2,3] --> R_dd_23 etc -bytes used=16004136, alloc=2096768, time=1.18 +bytes used=16007112, alloc=2096768, time=1.34 --> `codegen2/unindex` -bytes used=17004388, alloc=2096768, time=1.25 +bytes used=17007616, alloc=2096768, time=1.40 convert p/q --> RATIONAL(p/q) -bytes used=18004580, alloc=2096768, time=1.31 +bytes used=18007784, alloc=2096768, time=1.46 --> `codegen2/fix_rationals` writing C code codegen/C/expression: Unknown function: PARTIAL_RHO will be left as is. @@ -1141,57 +1177,58 @@ 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=19004732, alloc=2096768, time=1.37 -bytes used=20004972, alloc=2096768, time=1.47 +bytes used=19008004, alloc=2096768, time=1.51 +bytes used=20008500, alloc=2096768, time=1.62 expansion_Jacobian... -bytes used=21005184, alloc=2096768, time=1.56 -codegen2([partial_Theta_wrt_partial_d_h, partial_Theta_wrt_partial_dd_h]) --> "../gr.cg/expansion_Jacobian.c" +bytes used=21009180, alloc=2096768, time=1.70 +codegen2([partial_Theta_X_wrt_partial_d_h, partial_Theta_Y_wrt_partial_d_h, partial_Theta_X_wrt_partial_dd_h, partial_Theta_Y_wrt_partial_dd_h]) --> "../gr.cg/expansion_Jacobian.c" --> `codegen2/input` convert --> equation list -bytes used=22005440, alloc=2096768, time=1.63 +bytes used=22009512, alloc=2096768, time=1.77 --> `codegen2/eqnlist` optimizing computation sequence -bytes used=23005788, alloc=2162292, time=1.71 -bytes used=24006060, alloc=2293340, time=1.78 -bytes used=25006304, alloc=2293340, time=1.84 -bytes used=26006456, alloc=2293340, time=1.94 -bytes used=27007216, alloc=2293340, time=2.10 -bytes used=28007500, alloc=2293340, time=2.26 -bytes used=29007664, alloc=2293340, time=2.39 -bytes used=30007868, alloc=2424388, time=2.50 -bytes used=31008068, alloc=2424388, time=2.60 -bytes used=32009692, alloc=2555436, time=2.71 +bytes used=23009696, alloc=2162292, time=1.84 +bytes used=24009880, alloc=2293340, time=1.91 +bytes used=25010716, alloc=2293340, time=1.97 +bytes used=26010940, alloc=2293340, time=2.07 +bytes used=27011284, alloc=2293340, time=2.21 +bytes used=28011604, alloc=2293340, time=2.37 +bytes used=29011784, alloc=2293340, time=2.47 +bytes used=30012024, alloc=2424388, time=2.60 +bytes used=31012336, alloc=2424388, time=2.69 +bytes used=32012604, alloc=2489912, time=2.80 --> `codegen2/optimize` find temporary variables +bytes used=33012892, alloc=2620960, time=2.87 --> `codegen2/temps` convert Diff(expr,rho,sigma) --> PARTIAL_RHO_SIGMA(expr) etc -bytes used=33010028, alloc=2555436, time=2.77 -bytes used=34010432, alloc=2555436, time=2.84 -bytes used=35010740, alloc=2555436, time=2.91 +bytes used=34013116, alloc=2620960, time=2.94 +bytes used=35013316, alloc=2620960, time=3.01 +bytes used=36013600, alloc=2620960, time=3.08 --> `codegen2/fix_Diff` -bytes used=36021172, alloc=2555436, time=2.99 convert R_dd[2,3] --> R_dd_23 etc -bytes used=37021516, alloc=2555436, time=3.05 -bytes used=38021792, alloc=2555436, time=3.12 -bytes used=39022208, alloc=2555436, time=3.20 +bytes used=37014080, alloc=2620960, time=3.15 +bytes used=38014560, alloc=2620960, time=3.22 +bytes used=39014928, alloc=2620960, time=3.29 --> `codegen2/unindex` -bytes used=40022676, alloc=2555436, time=3.27 -bytes used=41022924, alloc=2555436, time=3.33 -bytes used=42023096, alloc=2555436, time=3.40 -bytes used=43023256, alloc=2555436, time=3.47 -bytes used=44023688, alloc=2555436, time=3.53 +bytes used=40015272, alloc=2620960, time=3.36 +bytes used=41015476, alloc=2620960, time=3.43 +bytes used=42015660, alloc=2620960, time=3.49 +bytes used=43016244, alloc=2620960, time=3.56 +bytes used=44016600, alloc=2620960, time=3.63 convert p/q --> RATIONAL(p/q) -bytes used=45023996, alloc=2555436, time=3.60 -bytes used=46024268, alloc=2555436, time=3.67 -bytes used=47024688, alloc=2555436, time=3.73 -bytes used=48025108, alloc=2555436, time=3.80 +bytes used=45016764, alloc=2620960, time=3.68 +bytes used=46016940, alloc=2620960, time=3.74 +bytes used=47017228, alloc=2620960, time=3.81 +bytes used=48017676, alloc=2620960, time=3.87 --> `codegen2/fix_rationals` writing C code -bytes used=49025260, alloc=2555436, time=3.86 -bytes used=50025744, alloc=2555436, time=3.91 -bytes used=51026008, alloc=2555436, time=4.03 -bytes used=52026168, alloc=2555436, time=4.18 -bytes used=53026328, alloc=2555436, time=4.32 -bytes used=54026584, alloc=2555436, time=4.44 +bytes used=49017944, alloc=2620960, time=3.94 +bytes used=50018104, alloc=2620960, time=3.99 +bytes used=51018308, alloc=2620960, time=4.08 +bytes used=52018568, alloc=2620960, time=4.20 +bytes used=53018900, alloc=2620960, time=4.35 +bytes used=54019192, alloc=2620960, time=4.53 +bytes used=55019432, alloc=2620960, time=4.68 > quit -bytes used=54895808, alloc=2555436, time=4.56 +bytes used=55315640, alloc=2620960, time=4.73 |