From 2d323c5f7b3725edd86bb0e4bb8c0ab01878a249 Mon Sep 17 00:00:00 2001 From: pollney Date: Wed, 18 Oct 2006 12:04:19 +0000 Subject: Added center_offset parameter to allow the initial data BHs to be shifted around the grid. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@70 b2a53a04-0f4f-0410-87ed-f9f25ced00cf --- src/TwoPunctures.c | 20 ++++++++++++++++++++ 1 file changed, 20 insertions(+) (limited to 'src') diff --git a/src/TwoPunctures.c b/src/TwoPunctures.c index cc48a82..a5164f6 100644 --- a/src/TwoPunctures.c +++ b/src/TwoPunctures.c @@ -306,6 +306,16 @@ TwoPunctures (CCTK_ARGUMENTS) imax[d] = cctk_lsh[d]; } + 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, -- cgit v1.2.3