diff options
| -rw-r--r-- | src/interpolate.c | 6 |
1 files changed, 5 insertions, 1 deletions
diff --git a/src/interpolate.c b/src/interpolate.c index b1d618f..9c6a392 100644 --- a/src/interpolate.c +++ b/src/interpolate.c @@ -701,7 +701,11 @@ Rot90_SymmetryInterpolate (CCTK_POINTER_TO_CONST const cctkGH_, } assert (thetensor[num_derivs][time_level][basevar] == tensortype); assert (thevars[num_derivs][time_level][basevar]); - assert (thevars[num_derivs][time_level][basevar][var] == -1); + /* This does not hold if the caller requests the same + interpolation to be done into different output arrays. This + may happen e.g. when CarpetInterp needs to differentiate in + time. This is arguably a performance bug in CarpetInterp. */ + /* assert (thevars[num_derivs][time_level][basevar][var] == -1); */ thevars[num_derivs][time_level][basevar][var] = m; } /* for m */ |