correct vector component computation during interpolationHEADmaster

before (since rev 89) there was a bug where always the x component of a vector
vel[] would be interpolated. An assert() would trigger if multiple vector
components were interpolated at the same time.


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

1 files changed, 5 insertions, 1 deletions

diff --git a/src/interpolate.c b/src/interpolate.c
index 4087119..985c419 100644
--- a/src/interpolate.c
+++ b/src/interpolate.c
@@ -710,7 +710,11 @@ Rot90_SymmetryInterpolate (CCTK_POINTER_TO_CONST const cctkGH,
               (numvars == DIM && vectorlength == DIM));
       tensortype = &vector;
       basevar = firstvar;
-      var = (vi - basevar) / vectorlength;
+      if (numvars == vectorlength) {
+        var = (vi - basevar);
+      } else {
+        var = (vi - basevar) / vectorlength;
+      }
     } else if (CCTK_EQUALS (tensortypealias, "4u")
                || CCTK_EQUALS (tensortypealias, "4d"))
     {