diff options
Diffstat (limited to 'src/interpolate.c')
-rw-r--r-- | src/interpolate.c | 9 |
1 files changed, 7 insertions, 2 deletions
diff --git a/src/interpolate.c b/src/interpolate.c index 2bdb183..80a11c1 100644 --- a/src/interpolate.c +++ b/src/interpolate.c @@ -577,8 +577,9 @@ Rot90_SymmetryInterpolate (CCTK_POINTER_TO_CONST const cctkGH_, for (m=0; m<N_output_arrays; ++m) { assert (operand_indices[m]>=0 && operand_indices[m]<N_input_arrays); output_array_indices[m] = input_array_indices[operand_indices[m]]; - assert (output_array_indices[m]>=0 - && output_array_indices[m]<CCTK_NumVars()); + assert (output_array_indices[m]==-1 + || (output_array_indices[m]>=0 + && output_array_indices[m]<CCTK_NumVars())); } @@ -608,6 +609,7 @@ Rot90_SymmetryInterpolate (CCTK_POINTER_TO_CONST const cctkGH_, assert (thevar); for (m=0; m<N_output_arrays; ++m) { + if (output_array_indices[m]!=-1) { int vi, gi; int numvars, firstvar; @@ -854,12 +856,14 @@ Rot90_SymmetryInterpolate (CCTK_POINTER_TO_CONST const cctkGH_, == -1); thevars[num_time_derivs][num_derivs][time_level][basevar][var] = m; + } } /* for m */ /* Loop over all output arrays */ for (m=0; m<N_output_arrays; ++m) { + if (output_array_indices[m]!=-1) { int num_time_derivs; int num_derivs; @@ -970,6 +974,7 @@ Rot90_SymmetryInterpolate (CCTK_POINTER_TO_CONST const cctkGH_, } /* for n */ + } } /* for m */ |