diff options
| author | pollney <pollney@b2a53a04-0f4f-0410-87ed-f9f25ced00cf> | 2006-08-28 17:06:49 +0000 |
|---|---|---|
| committer | pollney <pollney@b2a53a04-0f4f-0410-87ed-f9f25ced00cf> | 2006-08-28 17:06:49 +0000 |
| commit | 3923e58ef3b314bea51f1ec3ed40652a040ad2ad (patch) | |
| tree | 847733e5c20e7270c49bebbcb246a7b2d87a9e61 /src | |
| parent | ee8ea69303dc134d9ed7849ee6c51105778d109f (diff) | |
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
Diffstat (limited to 'src')
| -rw-r--r-- | src/TwoPunctures.c | 13 |
1 files changed, 10 insertions, 3 deletions
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 */ |