src/Lattice.c

   1 // --------------------------------------------------------------------------
   2 //
   3 // Copyright
   4 //   Markus Wittmann, 2016-2017
   5 //   RRZE, University of Erlangen-Nuremberg, Germany
   6 //   markus.wittmann -at- fau.de or hpc -at- rrze.fau.de
   7 //
   8 //   Viktor Haag, 2016
   9 //   LSS, University of Erlangen-Nuremberg, Germany
  10 //
  11 //  This file is part of the Lattice Boltzmann Benchmark Kernels (LbmBenchKernels).
  12 //
  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.
  17 //
  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.
  22 //
  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/>.
  25 //
  26 // --------------------------------------------------------------------------
  27 #include "Lattice.h"
  28
  29 // Dumps the layers [zStart, zStop] of lattice as ASCII.
  30 // Specify zStart = -1 and zStop = -1 as begin/end of lattice.
  31
  32 void LatDumpAscii(LatticeDesc * ld, int zStart, int zStop)
  33 {
  34         Assert(ld != NULL);
  35
  36         const char strLatCellType[] = "X.IxO"; // X = Obstacle, . = Fluid, I = inlet, O = outlet
  37
  38         int localZStart = zStart;
  39         int localZStop  = zStop;
  40
  41         int * dims = ld->Dims;
  42         LatticeT * lattice = ld->Lattice;
  43
  44         if (localZStart == -1) localZStart = 0;
  45         if (localZStop  == -1) localZStop  = dims[2] - 1;
  46
  47         for (int z = localZStop; z >= localZStart; --z) {
  48                 printf("plane % 2d\n", z);
  49
  50                 for (int y = dims[1] - 1; y >= 0; --y) {
  51                         printf(" %2d  ", y);
  52                         for (int x = 0; x < dims[0]; ++x) {
  53                                 printf("%c", strLatCellType[lattice[L_INDEX_4(dims, x, y, z)]]);
  54                         }
  55                         printf("\n");
  56                 }
  57         }
  58 }
  59