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
| macro VARIABLES (Point point = point, int _ig = ig, int _jg = jg, int _kg = kg)
{
double Delta = L0*_DELTA; /* cell size */
double Delta_x = Delta; /* cell size (with mapping) */
#if dimension > 1
double Delta_y = Delta; /* cell size (with mapping) */
#endif
#if dimension > 2
double Delta_z = Delta; /* cell size (with mapping) */
#endif
/* cell/face center coordinates */
double x = ((_ig + 1)/2. + _I)*Delta + X0; NOT_UNUSED(x);
#if dimension > 1
double y = ((_jg + 1)/2. + _J)*Delta + Y0;
#else
double y = 0.;
#endif
NOT_UNUSED(y);
#if dimension > 2
double z = ((_kg + 1)/2. + _K)*Delta + Z0;
#else
double z = 0.;
#endif
NOT_UNUSED(z);
/* we need this to avoid compiler warnings */
NOT_UNUSED(Delta);
NOT_UNUSED(Delta_x);
#if dimension > 1
NOT_UNUSED(Delta_y);
#endif
#if dimension > 2
NOT_UNUSED(Delta_z);
#endif
/* and this when catching FPEs */
_CATCH;
}
|
