1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
|
/* $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);} */
}
|