From 84e6c44587215f6ca45e2b2829e7b3eccd73981b Mon Sep 17 00:00:00 2001 From: rhaas Date: Tue, 6 Aug 2013 01:34:57 +0000 Subject: correct vector component computation during interpolation 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 --- src/interpolate.c | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) (limited to 'src') 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")) { -- cgit v1.2.3