aboutsummaryrefslogtreecommitdiff
path: root/src/ParamChecker.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/ParamChecker.c')
-rw-r--r--src/ParamChecker.c13
1 files changed, 13 insertions, 0 deletions
diff --git a/src/ParamChecker.c b/src/ParamChecker.c
index c18fcbe..6b24083 100644
--- a/src/ParamChecker.c
+++ b/src/ParamChecker.c
@@ -10,6 +10,7 @@
#include <stdio.h>
#include <stdlib.h>
+#include <math.h>
#include "cctk.h"
#include "cctk_Arguments.h"
@@ -47,6 +48,18 @@ void ParamChecker(CCTK_ARGUMENTS)
CCTK_INFO(message);
free(message);
}
+ if (CCTK_Equals(initial_data,"kerr") == 1)
+ {
+ CCTK_REAL tmp;
+ CCTK_INFO("Kerr black hole");
+ message = (char *)malloc(200*sizeof(char));
+ sprintf(message," m = %f, a = %f",mass,a_Kerr);
+ CCTK_INFO(message);
+ tmp=pow((mass*mass-a_Kerr*a_Kerr)/4.0,.5);
+ sprintf(message," throat at %f",tmp);
+ CCTK_INFO(message);
+ free(message);
+ }
else if (CCTK_Equals(initial_data,"bl_bh") == 1)
{
CCTK_INFO("Brill Lindquist black holes");