diff options
| author | Anton Khirnov <anton@khirnov.net> | 2012-12-04 12:51:11 +0100 |
|---|---|---|
| committer | Anton Khirnov <anton@khirnov.net> | 2012-12-04 12:51:11 +0100 |
| commit | e0131415f228602579c6ff2423e0b84828cef094 (patch) | |
| tree | cf7919be09509d522ba107b246fbe4080adb6639 | |
| parent | 308395f8a0f0fcd28ad50a467ef6f47ff2e35510 (diff) | |
Fix extra multiplication by gamma.
| -rw-r--r-- | src/trumpet.c | 8 |
1 files 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); |