1 // --------------------------------------------------------------------------
4 // Markus Wittmann, 2016-2017
5 // RRZE, University of Erlangen-Nuremberg, Germany
6 // markus.wittmann -at- fau.de or hpc -at- rrze.fau.de
9 // LSS, University of Erlangen-Nuremberg, Germany
11 // This file is part of the Lattice Boltzmann Benchmark Kernels (LbmBenchKernels).
13 // LbmBenchKernels is free software: you can redistribute it and/or modify
14 // it under the terms of the GNU General Public License as published by
15 // the Free Software Foundation, either version 3 of the License, or
16 // (at your option) any later version.
18 // LbmBenchKernels is distributed in the hope that it will be useful,
19 // but WITHOUT ANY WARRANTY; without even the implied warranty of
20 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
21 // GNU General Public License for more details.
23 // You should have received a copy of the GNU General Public License
24 // along with LbmBenchKernels. If not, see <http://www.gnu.org/licenses/>.
26 // --------------------------------------------------------------------------
29 // Dumps the layers [zStart, zStop] of lattice as ASCII.
30 // Specify zStart = -1 and zStop = -1 as begin/end of lattice.
32 void LatDumpAscii(LatticeDesc * ld, int zStart, int zStop)
36 const char strLatCellType[] = "X.IxO"; // X = Obstacle, . = Fluid, I = inlet, O = outlet
38 int localZStart = zStart;
39 int localZStop = zStop;
41 int * dims = ld->Dims;
42 LatticeT * lattice = ld->Lattice;
44 if (localZStart == -1) localZStart = 0;
45 if (localZStop == -1) localZStop = dims[2] - 1;
47 for (int z = localZStop; z >= localZStart; --z) {
48 printf("plane % 2d\n", z);
50 for (int y = dims[1] - 1; y >= 0; --y) {
52 for (int x = 0; x < dims[0]; ++x) {
53 printf("%c", strLatCellType[lattice[L_INDEX_4(dims, x, y, z)]]);