aboutsummaryrefslogtreecommitdiff
path: root/src/gr/horizon.maple
diff options
context:
space:
mode:
authorjthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5>2002-04-13 16:53:10 +0000
committerjthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5>2002-04-13 16:53:10 +0000
commit39fe29f9a5d91e149a794b22ba4d61cf07836fa8 (patch)
tree7e09140ee4ce1313154365800281fced5a5f7738 /src/gr/horizon.maple
parentfc79a77f60a3eb527a40b53b343f40f01fbd3da9 (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.maple14
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;