aboutsummaryrefslogtreecommitdiff
path: root/src/adjust.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/adjust.c')
-rw-r--r--src/adjust.c90
1 files changed, 0 insertions, 90 deletions
diff --git a/src/adjust.c b/src/adjust.c
deleted file mode 100644
index 43a2df0..0000000
--- a/src/adjust.c
+++ /dev/null
@@ -1,90 +0,0 @@
-/* $Header$ */
-
-#include <assert.h>
-#include <stdio.h>
-
-#include "cctk.h"
-#include "cctk_Arguments.h"
-#include "cctk_Parameters.h"
-
-
-
-void
-BHTracker_Adjust (CCTK_ARGUMENTS);
-
-static void
-setpar (char const * name, int index, char const * thorn, CCTK_REAL value);
-
-
-
-void
-BHTracker_Adjust (CCTK_ARGUMENTS)
-{
- DECLARE_CCTK_PARAMETERS;
- DECLARE_CCTK_ARGUMENTS;
-
- int n;
- int si, ri;
- char name[100], val[100];
- int icnt, ierr;
-
- CCTK_REAL xpos,ypos,zpos;
-
- if ((cctk_iteration-1) < track_first) return;
- if ((cctk_iteration-1) % track_every != 0) return;
-
-
-
- CCTK_INFO ("Tracking...");
-
- for (n=0; n<num_tracked; ++n) {
-
- si = surface_index[n];
- assert (si>=0 && si<nsurfaces);
-
- ri = region_index[n];
- assert (ri>=0 && ri<num_offsets);
-
- if (sf_valid[si] > 0) {
-
- CCTK_VInfo (CCTK_THORNSTRING,
- " Object #%d has moved to [%g,%g,%g]", n,
- (double)sf_centroid_x[si], (double)sf_centroid_y[si],
- (double)sf_centroid_z[si]);
-
- /* calculate new positions */
- xpos = (sf_centroid_x[si] - initial_x[si]);
- ypos = (sf_centroid_y[si] - initial_y[si]);
- zpos = (sf_centroid_z[si] - initial_z[si]);
-
- setpar ("offsetx", ri, "CarpetRegrid", xpos);
- setpar ("offsety", ri, "CarpetRegrid", ypos);
- setpar ("offsetz", ri, "CarpetRegrid", zpos);
-
- } else {
-
- CCTK_VInfo (CCTK_THORNSTRING,
- " Object #%d cannot be detected at this time", n);
-
- }
- }
-}
-
-
-
-static void
-setpar (char const * const name, int const index, char const * const thorn,
- CCTK_REAL const value)
-{
- char nambuf[100], valbuf[100];
- int icnt, ierr;
-
- icnt = snprintf (nambuf, sizeof nambuf, "%s[%d]", name, index);
- assert (icnt>=0 && icnt<sizeof nambuf);
-
- icnt = snprintf (valbuf, sizeof valbuf, "%g", value);
- assert (icnt>=0 && icnt<sizeof valbuf);
-
- ierr = CCTK_ParameterSet (nambuf, thorn, valbuf);
- assert (! ierr);
-}