From e0131415f228602579c6ff2423e0b84828cef094 Mon Sep 17 00:00:00 2001 From: Anton Khirnov Date: Tue, 4 Dec 2012 12:51:11 +0100 Subject: Fix extra multiplication by gamma. --- src/trumpet.c | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/trumpet.c b/src/trumpet.c index 7857160..40ce181 100644 --- a/src/trumpet.c +++ b/src/trumpet.c @@ -25,7 +25,7 @@ /* * isotropic/coordinate radius */ -#define ISO_R(x, y, z, gamma) (sqrt(SQR(gamma*x) + SQR(y) + SQR(z) + EPS)) +#define ISO_R(x, y, z) (sqrt(SQR(x) + SQR(y) + SQR(z) + EPS)) #define TRUMPET_ALPHA(R) (sqrt(1 - 2*MASS/R + SQR(TRUMPET_CONST)/SQR(SQR(R)))) @@ -195,7 +195,7 @@ void trumpet_data(CCTK_ARGUMENTS) for (int i = 0; i < cctk_lsh[0]; i++) { int index = CCTK_GFINDEX3D(cctkGH, i, j, k); CCTK_REAL xx = gamma*x[index], yy = y[index], zz = z[index]; - CCTK_REAL r = ISO_R(xx, yy, zz, gamma); + CCTK_REAL r = ISO_R(xx, yy, zz); CCTK_REAL R = gsl_spline_eval(spline, r, acc); CCTK_REAL alpha = TRUMPET_ALPHA(R); @@ -247,7 +247,7 @@ void trumpet_lapse(CCTK_ARGUMENTS) for (int i = 0; i < cctk_lsh[0]; i++) { int index = CCTK_GFINDEX3D(cctkGH, i, j, k); CCTK_REAL xx = gamma*x[index], yy = y[index], zz = z[index]; - CCTK_REAL r = ISO_R(xx, yy, zz, gamma); + CCTK_REAL r = ISO_R(xx, yy, zz); CCTK_REAL R = gsl_spline_eval(spline, r, acc); CCTK_REAL alpha = TRUMPET_ALPHA(R); @@ -288,7 +288,7 @@ void trumpet_shift(CCTK_ARGUMENTS) for (int i = 0; i < cctk_lsh[0]; i++) { int index = CCTK_GFINDEX3D(cctkGH, i, j, k); CCTK_REAL xx = gamma*x[index], yy = y[index], zz = z[index]; - CCTK_REAL r = ISO_R(xx, yy, zz, gamma); + CCTK_REAL r = ISO_R(xx, yy, zz); CCTK_REAL R = gsl_spline_eval(spline, r, acc); CCTK_REAL alpha = TRUMPET_ALPHA(R); -- cgit v1.2.3