From bd0fc470a7c3c64f5a14dc021b2d39fbd431a45d Mon Sep 17 00:00:00 2001
From: schnetter <schnetter@54511f98-0e4f-0410-826e-eb8b393f5a1e>
Date: Wed, 5 Oct 2005 19:38:55 +0000
Subject: #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
---
 src/volume_integral.c | 9 ++++++---
 1 file 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++)
-- 
cgit v1.2.3