From 13ce8affa4d740cb44201a2da7f9a3499d71082e Mon Sep 17 00:00:00 2001
From: diener <diener@2a26948c-0e4f-0410-aee8-f1d3e353619c>
Date: Fri, 24 Jan 2003 19:55:43 +0000
Subject: 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
---
 interface.ccl | 38 +++++++++++++++++++++++++++-----------
 1 file changed, 27 insertions(+), 11 deletions(-)

(limited to 'interface.ccl')

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
-- 
cgit v1.2.3