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
|
#include <assert.h>
#include <stdlib.h>
#include "cctk.h"
#include "cctk_Arguments.h"
void
carpettest_slabtest (CCTK_ARGUMENTS)
{
DECLARE_CCTK_ARGUMENTS;
int vi;
CCTK_INT directions[3][3];
CCTK_INT origin[3];
CCTK_INT extent[3];
CCTK_INT hsize[3];
int mapping;
int d, dd;
int ierr;
assert (CCTK_nProcs(cctkGH) == 1);
vi = CCTK_VarIndex ("grid::y");
assert (vi >= 0);
for (d=0; d<3; ++d) {
for (dd=0; dd<3; ++dd) {
directions[d][dd] = d == dd;
}
origin[d] = cctk_lbnd[d];
extent[d] = cctk_lsh[d];
}
mapping = Hyperslab_GlobalMappingByIndex
(cctkGH, vi, 3, &directions[0][0], origin, extent, NULL, -1, NULL, hsize);
assert (mapping >= 0);
for (d=0; d<3; ++d) {
assert (hsize[d] == cctk_lsh[d]);
}
ierr = Hyperslab_Get (cctkGH, mapping, -1, vi, 0, CCTK_VARIABLE_REAL, yy);
assert (! ierr);
ierr = Hyperslab_FreeMapping (mapping);
assert (! ierr);
}
|