Second pass at finding points on the level set surface. It now uses

grid arrays and works in parallel. Still needs to be generalized to work
in other modes than full mode and an algorithm has to be implemented
to figure out how many surfaces are actually present in the level set
function.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/EHFinder/trunk@78 2a26948c-0e4f-0410-aee8-f1d3e353619c

1 files changed, 27 insertions, 11 deletions

diff --git a/interface.ccl b/interface.ccl
index ed19712..d1991ec 100644
--- a/interface.ccl
+++ b/interface.ccl
@@ -9,43 +9,59 @@ USES INCLUDE: MoL.h
 
 public:
 
-CCTK_REAL level_set type=GF TimeLevels=2
+CCTK_REAL level_set TYPE=GF TIMELEVELS=2
 {
   f
 } "The scalar level set function that defines the null surface"
 
-CCTK_REAL slevel_set type=GF TimeLevels=1
+CCTK_REAL slevel_set TYPE=GF TIMELEVELS=1
 {
   sf
 } "Source for the level set function"
 
-CCTK_REAL dlevel_set type=GF TimeLevels=1
+CCTK_REAL dlevel_set TYPE=GF TIMELEVELS=1
 {
   dfx, dfy, dfz, dfsq
 } "Derivatives of the level set function"
 
-CCTK_REAL ftmp_set type=GF TimeLevels=1
+CCTK_REAL ftmp_set TYPE=GF TIMELEVELS=1
 {
   ftmp, sftmp
 } "temporary variables used in pde re-parametrization"
 
-CCTK_REAL level_bak type=GF TimeLevels=1
+CCTK_REAL level_bak TYPE=GF TIMELEVELS=1
 {
   fbak
 } "Temporary placeholder for the level set during re-parametrization"
 
-CCTK_INT eh_mask_all type=GF TimeLevels=1
+CCTK_INT eh_mask_all TYPE=GF TIMELEVELS=1
 {
   eh_mask, tm_mask
 } "Masks to define active cells"
 
-CCTK_INT mask_bak type=GF TimeLevels=1
+CCTK_INT mask_bak TYPE=GF TIMELEVELS=1
 {
   eh_mask_bak
 } "Temporary placeholder for the mask during re-parametrization"
 
-CCTK_INT rep_mask type=GF TimeLevels=1
+CCTK_REAL surface_arrays TYPE=ARRAY DIM=2 TIMELEVELS=1 SIZE=ntheta,nphi GHOSTSIZE=n_array_ghosts,n_array_ghosts DISTRIB=DEFAULT
+{
+  ctheta, cphi, rsurf, sinthetacosphi, sinthetasinphi, costheta
+} "Grid arrays for points on the surface"
+
+CCTK_REAL surface_tmp_arrays TYPE=ARRAY DIM=2 TIMELEVELS=1 SIZE=ntheta,nphi GHOSTSIZE=n_array_ghosts,n_array_ghosts DISTRIB=DEFAULT
+{
+  drsurf, interp_x, interp_y, interp_z, f_interp,
+  dfdx_interp, dfdy_interp, dfdz_interp
+} "Tmeporary grid arrays for finding points on the surface"
+
+CCTK_INT surface_int_array TYPE=ARRAY DIM=2 TIMELEVELS=1 SIZE=ntheta,nphi GHOSTSIZE=n_array_ghosts,n_array_ghosts DISTRIB=DEFAULT
+{
+  n_since_last_reduction
+} "Tmeporary integer grid array for finding points on the surface"
+
+CCTK_INT rep_mask TYPE=GF TIMELEVELS=1
 
-CCTK_INT re_param_control_pde type=SCALAR
-CCTK_INT re_param_control_approx type=SCALAR
-CCTK_REAL eh_area type=SCALAR
+CCTK_INT re_param_control_pde TYPE=SCALAR
+CCTK_INT re_param_control_approx TYPE=SCALAR
+CCTK_REAL eh_area TYPE=SCALAR