diff options
Diffstat (limited to 'src/ParamChecker.c')
-rw-r--r-- | src/ParamChecker.c | 13 |
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"); |