Handle the case of 0 interpolation points correctly. malloc can

return NULL in this case.


git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/RotatingSymmetry90/trunk@38 c3c03602-0f4f-0410-b3fa-d2c81c8a7dc5

2 files changed, 4 insertions, 2 deletions

diff --git a/src/interpolate.c b/src/interpolate.c
index 80a11c1..30a18f8 100644
--- a/src/interpolate.c
+++ b/src/interpolate.c
@@ -455,7 +455,7 @@ Rot90_SymmetryInterpolate (CCTK_POINTER_TO_CONST const cctkGH_,
   /* Copy coordinates */
   for (d=0; d<DIM; ++d) {
     new_interp_coords[d] = malloc (N_interp_points * sizeof(CCTK_REAL));
-    assert (new_interp_coords[d]);
+    assert (N_interp_points == 0 || new_interp_coords[d]);
   }
   
   /* Fold coordinates */
@@ -487,7 +487,7 @@ Rot90_SymmetryInterpolate (CCTK_POINTER_TO_CONST const cctkGH_,
   assert (new_output_arrays);
   for (m=0; m<N_output_arrays; ++m) {
     new_output_arrays[m] = malloc (N_interp_points * sizeof(CCTK_REAL));
-    assert (new_output_arrays[m]);
+    assert (N_interp_points == 0 || new_output_arrays[m]);
   }
   
   
diff --git a/src/rotatingsymmetry90.c b/src/rotatingsymmetry90.c
index 8e89689..51150c4 100644
--- a/src/rotatingsymmetry90.c
+++ b/src/rotatingsymmetry90.c
@@ -594,6 +594,8 @@ void Rot90_ApplyBC (CCTK_ARGUMENTS)
   nvars = Boundary_SelectedGVs (cctkGH, 0, 0, 0, 0, 0, 0);
   assert (nvars>=0);
   
+  if (nvars==0) return;
+  
   indices = malloc (nvars * sizeof *indices);
   assert (indices);
   faces = malloc (nvars * sizeof *faces);