diff options
| author | diener <diener@2a26948c-0e4f-0410-aee8-f1d3e353619c> | 2003-02-07 23:08:15 +0000 |
|---|---|---|
| committer | diener <diener@2a26948c-0e4f-0410-aee8-f1d3e353619c> | 2003-02-07 23:08:15 +0000 |
| commit | 6f5c1c3a019cf4f116874944ed30a5098d5f50d5 (patch) | |
| tree | ebb8ac56dc13fe43285cd3e5fec382ec15c5b07b /interface.ccl | |
| parent | 4727e99931b092d07fae4e3188b45ae31bea927e (diff) | |
Added support for finding out how many surfaces are present and then do the
integrals over all surfaces. Still only supports full mode. Only weakly tested.
There is a problem with the outputs in multiprocessor mode. The results of
the integration are stored in 1D grid arrays. However if I use DISTRIB=DEFAULT
(taking care only to store results on the correct processor), the
output fails if the number of processors are bigger than the number of
elements in the grid array (essentially some processors contains a chunk
of length zero). If I on the other hand uses DISTRIB=CONSTANT, it works
but on multiple processors i get multiple copies of the data in the output.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/EHFinder/trunk@90 2a26948c-0e4f-0410-aee8-f1d3e353619c
Diffstat (limited to 'interface.ccl')
| -rw-r--r-- | interface.ccl | 33 |
1 files changed, 28 insertions, 5 deletions
diff --git a/interface.ccl b/interface.ccl index f4fe515..9af0506 100644 --- a/interface.ccl +++ b/interface.ccl @@ -53,7 +53,8 @@ CCTK_INT surface_counter TYPE=SCALAR CCTK_INT points_counter TYPE=SCALAR CCTK_INT more_surfaces TYPE=SCALAR CCTK_INT more_points TYPE=SCALAR - +CCTK_INT integrate_counter TYPE=SCALAR + 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, sintheta, costheta, sinphi, cosphi, @@ -90,9 +91,31 @@ 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_REAL eh_area2 TYPE=SCALAR -CCTK_REAL eh_centroid TYPE=SCALAR + +#CCTK_REAL eh_area TYPE=ARRAY DIM=1 TIMELEVELS=1 SIZE=maximum_surface_number GHOSTSIZE=0 DISTRIB=DEFAULT + +CCTK_REAL eh_area TYPE=ARRAY DIM=1 TIMELEVELS=1 SIZE=maximum_surface_number DISTRIB=CONSTANT + +#CCTK_REAL eh_area[maximum_surface_number] TYPE=SCALAR + +#CCTK_REAL eh_area2 TYPE=ARRAY DIM=1 TIMELEVELS=1 SIZE=maximum_surface_number GHOSTSIZE=0 DISTRIB=DEFAULT #"The area of the event horizon(s)" +CCTK_REAL eh_area2 TYPE=ARRAY DIM=1 TIMELEVELS=1 SIZE=maximum_surface_number DISTRIB=CONSTANT #"The area of the event horizon(s)" + +#CCTK_REAL eh_area2[maximum_surface_number] TYPE=SCALAR +#CCTK_REAL eh_area3 TYPE=ARRAY DIM=1 TIMELEVELS=1 SIZE=maximum_surface_number DISTRIB=CONSTANT + +CCTK_REAL eh_centroid TYPE=ARRAY DIM=1 TIMELEVELS=1 SIZE=maximum_surface_number DISTRIB=CONSTANT { eh_centroid_x, eh_centroid_y, eh_centroid_z -} "The centroid of the event horizon" +} "The centroid of the event horizon(s)" + +CCTK_REAL eh_centroid2 TYPE=ARRAY DIM=1 TIMELEVELS=1 SIZE=maximum_surface_number DISTRIB=CONSTANT +{ + eh_centroid2_x, eh_centroid2_y, eh_centroid2_z +} "The centroid of the event horizon(s)" + +#CCTK_REAL eh_centroid3 TYPE=ARRAY DIM=1 TIMELEVELS=1 SIZE=maximum_surface_number DISTRIB=CONSTANT +#{ +# eh_centroid3_x, eh_centroid3_y, eh_centroid3_z +#} "The centroid of the event horizon(s)" + |