diff options
author | hawke <hawke@eec4d7dc-71c2-46d6-addf-10296150bf52> | 2005-11-10 14:05:35 +0000 |
---|---|---|
committer | hawke <hawke@eec4d7dc-71c2-46d6-addf-10296150bf52> | 2005-11-10 14:05:35 +0000 |
commit | 77424202a473eae863478cbd7c901a0c494d652d (patch) | |
tree | 1e678e64edfa876af9e77dae37c6703950171930 /src/SymInterp.c | |
parent | 2b42ff6020d0ab36b11b65466ef6850f7ec68242 (diff) |
Zeroth attempt at doing excision with Cartoon. Now requires the
Spacemask thorn from CactusEinstein.
Initial tests suggest perfect second order convergence but not quite
perfect fourth order convergence, at least at reasonable
resolutions. Convergence is still better than 3rd order, though.
Excision with the new spacemask has been added but not tested.
git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/Cartoon2D/trunk@97 eec4d7dc-71c2-46d6-addf-10296150bf52
Diffstat (limited to 'src/SymInterp.c')
-rw-r--r-- | src/SymInterp.c | 7 |
1 files changed, 7 insertions, 0 deletions
diff --git a/src/SymInterp.c b/src/SymInterp.c index 0efc291..24668ab 100644 --- a/src/SymInterp.c +++ b/src/SymInterp.c @@ -360,6 +360,13 @@ Cartoon2D_SymmetryInterpolate (CCTK_POINTER_TO_CONST const cctkGH_, for (m=0; m<N_output_arrays; ++m) { new_output_arrays[m] = malloc (N_interp_points * sizeof(CCTK_REAL)); assert (new_output_arrays[m]); + { + CCTK_INT mm; + for (mm = 0; mm < N_interp_points; ++mm) + { + ((CCTK_REAL*)new_output_arrays[m])[mm] = 0; + } + } } |