diff options
Diffstat (limited to 'src/Misner_multiple.c')
-rw-r--r-- | src/Misner_multiple.c | 19 |
1 files changed, 10 insertions, 9 deletions
diff --git a/src/Misner_multiple.c b/src/Misner_multiple.c index aa50249..d53a3b9 100644 --- a/src/Misner_multiple.c +++ b/src/Misner_multiple.c @@ -3,11 +3,12 @@ @date @author Carsten Gundlach @desc - Set up initial data for multiple Misner black holes + Set up initial data for multiple Misner black holes @enddesc + @version $Id$ @@*/ -#include <math.h> +#include <string.h> #include "cctk.h" #include "cctk_Arguments.h" @@ -15,6 +16,10 @@ #include "CactusEinstein/Einstein/src/Einstein.h" +static char *rcsid = "$Header$"; +CCTK_FILEVERSION(CactusEinstein_IDAnalyticBH_Misner_multiple_c) + + #define SQR(a) ((a)*(a)) void Misner_init(int n, CCTK_REAL mu, int terms); @@ -26,13 +31,9 @@ void Misner_multiple(CCTK_ARGUMENTS); @date @author Carsten Gundlach @desc - Set up initial data for multiple Misner black holes + Set up initial data for multiple Misner black holes @enddesc - @calls - @history - - @endhistory - + @calls MisnerEvalPsi @@*/ void Misner_multiple(CCTK_ARGUMENTS) { @@ -148,7 +149,7 @@ void Misner_multiple(CCTK_ARGUMENTS) { for(i = 0; i < npoints; i++) { - gxx[i] = pow(psi[i], 4); + gxx[i] = psi[i] * psi[i] * psi[i] * psi[i]; gyy[i] = gxx[i]; gzz[i] = gxx[i]; } |