diff options
author | schnetter <schnetter@54511f98-0e4f-0410-826e-eb8b393f5a1e> | 2005-10-05 19:38:55 +0000 |
---|---|---|
committer | schnetter <schnetter@54511f98-0e4f-0410-826e-eb8b393f5a1e> | 2005-10-05 19:38:55 +0000 |
commit | bd0fc470a7c3c64f5a14dc021b2d39fbd431a45d (patch) | |
tree | 795c11c5249f503b9135c61e9e897d522cf6092c | |
parent | b56eefb92972385433caae60b97b6957ec6aa2f8 (diff) |
#include <math.h> for sqrt()
#define PI to be more accurate (but then, who cares...)
git-svn-id: http://svn.aei.mpg.de/numrel/AEIThorns/ADMMass/trunk@3 54511f98-0e4f-0410-826e-eb8b393f5a1e
-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++) |