diff options
| -rw-r--r-- | interface.ccl | 28 | ||||
| -rw-r--r-- | src/README | 2 | ||||
| -rw-r--r-- | src/driver/make.code.defn | 2 |
3 files changed, 30 insertions, 2 deletions
diff --git a/interface.ccl b/interface.ccl index 0ac8643..83cc539 100644 --- a/interface.ccl +++ b/interface.ccl @@ -7,7 +7,33 @@ inherits: Grid ADMBase StaticConformal SpaceMask IO always # include file provided by SpaceMask USES INCLUDE: SpaceMask.h -# aliased function maybe provided by DriftCorrect +# provides aliased functions for various other thorns +# FIXME: this really belongs in Erik Schnetter's "generic surface" thorn +# ... returns local coordinate origin for a given horizon +SUBROUTINE HorizonLocalCoordinateOrigin \ + (CCTK_INT IN horizon_number, \ + CCTK_REAL OUT origin_x, CCTK_REAL OUT origin_y, CCTK_REAL OUT origin_z) +# ... computes horizon radius in the direction of each (x,y,z) point, +# or -1.0 if this horizon wasn't found +# ... more precisely: +# For each (x,y,z), consider the ray from the local coordinate +# origin through (x,y,z). This function computes the Euclidean +# distance between the local coordinate origin and this ray's +# intersection with the horizon, or -1.0 if this horizon wasn't +# found the last time we searched for it. +SUBROUTINE HorizonRadiusInDirection \ + (CCTK_INT IN horizon_number, \ + CCTK_INT IN N_points, \ + CCTK_REAL IN ARRAY x, CCTK_REAL IN ARRAY y, CCTK_REAL IN ARRAY z, \ + CCTK_REAL OUT ARRAY radius) +PROVIDES FUNCTION HorizonLocalCoordinateOrigin \ + WITH AHFinderDirect_local_coordinate_origin \ + LANGUAGE C +PROVIDES FUNCTION HorizonRadiusInDirection \ + WITH AHFinderDirect_radius_in_direction \ + LANGUAGE C + +# uses aliased function maybe provided by DriftCorrect SUBROUTINE SetDriftCorrectPosition(CCTK_POINTER IN GH, CCTK_REAL IN x, \ CCTK_REAL IN y, \ CCTK_REAL IN z) @@ -22,6 +22,8 @@ within Cactus; these are listed in top-down order (no code depends on, or includes header files from, higher-level code): driver/ contains high-level driver routines to solve the Theta(h) = 0 equations and interface to the rest of Cactus + (all the aliased-function stuff to interface with other + thorns also lives here) gr/ contains relativity code; all knowledge of the actual apparent horizon equation lives in the code in this directory gr.cg/ contains Maple-generated C code (generated from Maple diff --git a/src/driver/make.code.defn b/src/driver/make.code.defn index 352a0cb..d26a076 100644 --- a/src/driver/make.code.defn +++ b/src/driver/make.code.defn @@ -6,7 +6,7 @@ SRCS = state.cc \ setup.cc find_horizons.cc \ initial_guess.cc Newton.cc \ io.cc misc-driver.cc \ - announce.cc mask.cc \ + announce.cc horizon_radius.cc mask.cc \ BH_diagnostics.cc horizon_sequence.cc # Subdirectories containing source files |