diff options
-rw-r--r-- | param.ccl | 5 | ||||
-rw-r--r-- | src/TwoPunctures.c | 20 |
2 files changed, 25 insertions, 0 deletions
@@ -122,6 +122,11 @@ REAL par_S_minus[3] "spin of the m- puncture" (*:*) :: "" } 0.0 +REAL center_offset[3] "offset b=0 to position (x,y,z)" +{ + (*:*) :: "" +} 0.0 + REAL initial_lapse_psi_exponent "Exponent n for psi^-n initial lapse profile" { (*:*) :: "Should be negative" diff --git a/src/TwoPunctures.c b/src/TwoPunctures.c index cc48a82..a5164f6 100644 --- a/src/TwoPunctures.c +++ b/src/TwoPunctures.c @@ -307,6 +307,16 @@ TwoPunctures (CCTK_ARGUMENTS) } for (k = imin[2]; k < imax[2]; ++k) + for (j = imin[1]; j < imax[1]; ++j) + for (i = imin[0]; i < imax[0]; ++i) + { + int ijk = CCTK_GFINDEX3D(cctkGH, i, j, k); + x[ijk] -= center_offset[0]; + y[ijk] -= center_offset[1]; + z[ijk] -= center_offset[2]; + } + + for (k = imin[2]; k < imax[2]; ++k) { for (j = imin[1]; j < imax[1]; ++j) { @@ -543,6 +553,16 @@ TwoPunctures (CCTK_ARGUMENTS) } /* for j */ } /* for k */ + for (k = imin[2]; k < imax[2]; ++k) + for (j = imin[1]; j < imax[1]; ++j) + for (i = imin[0]; i < imax[0]; ++i) + { + int ijk = CCTK_GFINDEX3D(cctkGH, i, j, k); + x[ijk] += center_offset[0]; + y[ijk] += center_offset[1]; + z[ijk] += center_offset[2]; + } + if (use_sources && rescale_sources) { Rescale_Sources(cctkGH, |