From 3923e58ef3b314bea51f1ec3ed40652a040ad2ad Mon Sep 17 00:00:00 2001 From: pollney Date: Mon, 28 Aug 2006 17:06:49 +0000 Subject: Added initial lapse profile "brownsville" accordint to the specification in Phys. Rev. D 74, 041501 (2006). git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@67 b2a53a04-0f4f-0410-87ed-f9f25ced00cf --- src/TwoPunctures.c | 13 ++++++++++--- 1 file changed, 10 insertions(+), 3 deletions(-) (limited to 'src') diff --git a/src/TwoPunctures.c b/src/TwoPunctures.c index 8352890..f1fcf4f 100644 --- a/src/TwoPunctures.c +++ b/src/TwoPunctures.c @@ -200,7 +200,7 @@ TwoPunctures (CCTK_ARGUMENTS) enum GRID_SETUP_METHOD { GSM_Taylor_expansion, GSM_evaluation }; enum GRID_SETUP_METHOD gsm; - int antisymmetric_lapse, averaged_lapse, pmn_lapse; + int antisymmetric_lapse, averaged_lapse, pmn_lapse, brownsville_lapse; int nvar = 1, n1 = npoints_A, n2 = npoints_B, n3 = npoints_phi; @@ -273,9 +273,14 @@ TwoPunctures (CCTK_ARGUMENTS) antisymmetric_lapse = CCTK_EQUALS(initial_lapse, "twopunctures-antisymmetric"); averaged_lapse = CCTK_EQUALS(initial_lapse, "twopunctures-averaged"); pmn_lapse = CCTK_EQUALS(initial_lapse, "psi^n"); - if (pmn_lapse) + if (pmn_lapse) CCTK_VInfo(CCTK_THORNSTRING, "Setting initial lapse to psi^%f profile.", (double)initial_lapse_psi_exponent); + brownsville_lapse = CCTK_EQUALS(initial_lapse, "brownsville"); + if (brownsville_lapse) + CCTK_VInfo(CCTK_THORNSTRING, + "Setting initial lapse to a Brownsville-style profile.", + (double)initial_lapse_psi_exponent); CCTK_INFO ("Interpolating result"); if (CCTK_EQUALS(metric_type, "static conformal")) { @@ -382,7 +387,7 @@ TwoPunctures (CCTK_ARGUMENTS) if (multiply_old_lapse) old_alp = alp[ind]; - if ((*conformal_state > 0) || (pmn_lapse)) { + if ((*conformal_state > 0) || (pmn_lapse) || (brownsville_lapse)) { CCTK_REAL xp, yp, zp, rp, ir; CCTK_REAL s1, s3, s5; @@ -474,6 +479,8 @@ TwoPunctures (CCTK_ARGUMENTS) if (pmn_lapse) alp[ind] = pow(p, initial_lapse_psi_exponent); + if (brownsville_lapse) + alp[ind] = 2.0/(1.0+pow(p, initial_lapse_psi_exponent)); } /* if conformal-state > 0 */ -- cgit v1.2.3