diff options
Diffstat (limited to 'src/gr/horizon.maple')
| -rw-r--r-- | src/gr/horizon.maple | 34 |
1 files changed, 21 insertions, 13 deletions
diff --git a/src/gr/horizon.maple b/src/gr/horizon.maple index 7284253..423642f 100644 --- a/src/gr/horizon.maple +++ b/src/gr/horizon.maple @@ -258,7 +258,7 @@ global @include "../maple/gfa.minc", @include "../gr/gr_gfas.minc"; local u,v, - temp; + temp1,temp2; printf("%a...\n", procname); @@ -284,14 +284,16 @@ assert_fnd_exists(Theta_D); := frontend('diff', [Theta_D__fnd, Diff(h,y_rs[u])]); # equation (A1a) in my 1996 apparent horizon finding paper - 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; # Jacobian coefficients of Theta_[AB] and Theta wrt Diff(h,y_rs[u],y_rs[v]) @@ -305,17 +307,23 @@ assert_fnd_exists(Theta_D); := frontend('diff', [Theta_B__fnd, Diff(h,y_rs[u],y_rs[v])]); # equation (A1b) in my 1996 apparent horizon finding paper - 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'], + 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"); fi; |