diff options
author | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2002-04-13 16:53:10 +0000 |
---|---|---|
committer | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2002-04-13 16:53:10 +0000 |
commit | 39fe29f9a5d91e149a794b22ba4d61cf07836fa8 (patch) | |
tree | 7e09140ee4ce1313154365800281fced5a5f7738 /src/gr/horizon.maple | |
parent | fc79a77f60a3eb527a40b53b343f40f01fbd3da9 (diff) |
small tweaks to work with latest changes to ../maple/codegen2.maple
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@512 f88db872-0e4f-0410-b76b-b9085cfa78c5
Diffstat (limited to 'src/gr/horizon.maple')
-rw-r--r-- | src/gr/horizon.maple | 14 |
1 files changed, 8 insertions, 6 deletions
diff --git a/src/gr/horizon.maple b/src/gr/horizon.maple index 89d5537..d610844 100644 --- a/src/gr/horizon.maple +++ b/src/gr/horizon.maple @@ -18,8 +18,8 @@ horizon := proc(cg_flag::boolean) -non_unit_normal(cg_flag); -non_unit_normal_deriv(cg_flag); +non_unit_normal(); +non_unit_normal_deriv(); horizon_function(cg_flag); NULL; @@ -30,10 +30,10 @@ end proc; # # This function computes the non-unit outward-pointing normal (covector) # to the horizon, s_d = n_u(h), using the equation on p.7.1 of my apparent -# horizon finding notes. +# horizon finding notes. This function does *NOT* generate any C code. # non_unit_normal := -proc(cg_flag::boolean) +proc() global @include "../maple/coords.minc", @include "../maple/gfa.minc", @@ -60,7 +60,8 @@ end proc; # # This function computes the 1st derivatives of the non-unit outward-pointing # normal (covector) to the horizon, using the equation on p.7.2 of my -# apparent horizon finding notes. +# apparent horizon finding notes. This function does *NOT* generate any +# C code. # non_unit_normal_deriv := proc(cg_flag::boolean) @@ -106,6 +107,7 @@ end proc; # as a function of 1st and 2nd angular partial derivatives of the # apparent horizon radius $h$, in the Maple gridfn array H.fnd , # using equation (15) in my 1996 apparent horizon finding paper. +# This function also optionally generates C code for this computation. # horizon_function := proc(cg_flag::boolean) @@ -153,7 +155,7 @@ HD__fnd := msum('g_uu[i,j]*s_d__fnd[i]*s_d__fnd[j]', 'i'=1..N, 'j'=1..N); 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', H__fnd, "../gr.cg/H.c"); + then codegen2(H__fnd, 'H', "cg/horizon.c"); fi; NULL; |