aboutsummaryrefslogtreecommitdiff
path: root/Auxiliary/Cactus/KrancNumericalTools/Ceiling/src/selectGFs.c
diff options
context:
space:
mode:
Diffstat (limited to 'Auxiliary/Cactus/KrancNumericalTools/Ceiling/src/selectGFs.c')
-rw-r--r--Auxiliary/Cactus/KrancNumericalTools/Ceiling/src/selectGFs.c94
1 files changed, 0 insertions, 94 deletions
diff --git a/Auxiliary/Cactus/KrancNumericalTools/Ceiling/src/selectGFs.c b/Auxiliary/Cactus/KrancNumericalTools/Ceiling/src/selectGFs.c
deleted file mode 100644
index 86e1196..0000000
--- a/Auxiliary/Cactus/KrancNumericalTools/Ceiling/src/selectGFs.c
+++ /dev/null
@@ -1,94 +0,0 @@
-/* $Header$ */
-
-/* this code is based on Erik Schnetter's dissipation thorn */
-
-#include <assert.h>
-#include <stdlib.h>
-
-#include "cctk.h"
-#include "cctk_Arguments.h"
-#include "cctk_Parameters.h"
-
-void CCTK_FCALL
-CCTK_FNAME(apply_check_abs) (CCTK_REAL const * const var,
- int const * const ni,
- int const * const nj,
- int const * const nk,
- CCTK_REAL const * const ceiling_value);
-
-
-void CCTK_FCALL
-CCTK_FNAME(apply_check_diff) (CCTK_REAL const * const var,
- int const * const ni,
- int const * const nj,
- int const * const nk,
- CCTK_REAL const * const ceiling_value);
-
-static void
-call_apply_check (int const varindex, char const * const optstring, void * const arg);
-
-void
-check_ceiling (CCTK_ARGUMENTS)
-{
- DECLARE_CCTK_ARGUMENTS;
- DECLARE_CCTK_PARAMETERS;
-
- CCTK_TraverseString (vars, call_apply_check, cctkGH, CCTK_GROUP_OR_VAR);
-}
-
-
-void
-call_apply_check (int const varindex, char const * const optstring, void * const arg)
-{
- cGH const * const cctkGH = (cGH const *) arg;
- DECLARE_CCTK_ARGUMENTS;
- DECLARE_CCTK_PARAMETERS;
-
- int vargroup;
- cGroup vardata;
-
- CCTK_REAL const * varptr;
- int ierr /* , terminate */ ;
-
- assert (varindex >= 0);
-
- if (verbose) {
- char * const fullvarname = CCTK_FullName (varindex);
- assert (fullvarname);
- CCTK_VInfo (CCTK_THORNSTRING,
- "Applying ceiling check to \"%s\" ",
- fullvarname);
- free (fullvarname);
- }
-
- vargroup = CCTK_GroupIndexFromVarI (varindex);
- assert (vargroup >= 0);
-
- ierr = CCTK_GroupData (vargroup, &vardata);
- assert (!ierr);
-
- assert (vardata.grouptype == CCTK_GF);
- assert (vardata.vartype == CCTK_VARIABLE_REAL);
- assert (vardata.dim == cctk_dim);
-
- varptr = CCTK_VarDataPtrI (cctkGH, 0, varindex);
- assert (varptr);
-
- if (CCTK_Equals (type, "absolute"))
- {
- CCTK_FNAME(apply_check_abs)
- (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], &ceiling_value);
- }
- else if (CCTK_Equals (type, "differential"))
- {
- CCTK_FNAME(apply_check_diff)
- (varptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2], &ceiling_value);
- }
- else
- {
- CCTK_INFO("keyword ceiling::type only allows values 'absolute' and 'differential'");
- }
-
- /* if (terminate > 0) {CCTK_TerminateNext (cctkGH);} */
-}
-