diff options
-rw-r--r-- | src/volume_integral.c | 9 |
1 files changed, 6 insertions, 3 deletions
diff --git a/src/volume_integral.c b/src/volume_integral.c index 0a9a1dd..7db9d6d 100644 --- a/src/volume_integral.c +++ b/src/volume_integral.c @@ -1,4 +1,5 @@ +#include <math.h> #include <stdio.h> #include <stdlib.h> @@ -8,6 +9,8 @@ #include "SpaceMask.h" +#define PI 3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342117068 + void ADMMass_Volume(CCTK_ARGUMENTS) { DECLARE_CCTK_ARGUMENTS @@ -109,11 +112,11 @@ void ADMMass_Volume(CCTK_ARGUMENTS) ( dg[ti][tk][tj] - dg[ti][tj][tk] ); } ADMMass_VolumeMass_pot_x[ijk] *= alp[ijk] * sqrt(DETG_DETG) / - 16 / 3.14159265; + 16 / PI; ADMMass_VolumeMass_pot_y[ijk] *= alp[ijk] * sqrt(DETG_DETG) / - 16 / 3.14159265; + 16 / PI; ADMMass_VolumeMass_pot_z[ijk] *= alp[ijk] * sqrt(DETG_DETG) / - 16 / 3.14159265; + 16 / PI; } for(i=4; i<cctk_lsh[0]-4; i++) for(j=4; j<cctk_lsh[1]-4; j++) |