From 94d3e79a8617f88dc0219cfdeedfa3147833719d Mon Sep 17 00:00:00 2001
From: Adrian Kummerlaender
Date: Mon, 24 Jun 2019 14:43:36 +0200
Subject: Initialize at openlb-1-3
---
src/communication/cuboidNeighbourhood3D.h | 190 ++++++++++++++++++++++++++++++
1 file changed, 190 insertions(+)
create mode 100644 src/communication/cuboidNeighbourhood3D.h
(limited to 'src/communication/cuboidNeighbourhood3D.h')
diff --git a/src/communication/cuboidNeighbourhood3D.h b/src/communication/cuboidNeighbourhood3D.h
new file mode 100644
index 0000000..639102c
--- /dev/null
+++ b/src/communication/cuboidNeighbourhood3D.h
@@ -0,0 +1,190 @@
+/* This file is part of the OpenLB library
+ *
+ * Copyright (C) 2007, 2014 Mathias J. Krause
+ * E-mail contact: info@openlb.net
+ * The most recent release of OpenLB can be downloaded at
+ *
+ *
+ * This program is free software; you can redistribute it and/or
+ * modify it under the terms of the GNU General Public License
+ * as published by the Free Software Foundation; either version 2
+ * of the License, or (at your option) any later version.
+ *
+ * This program is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ * GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public
+ * License along with this program; if not, write to the Free
+ * Software Foundation, Inc., 51 Franklin Street, Fifth Floor,
+ * Boston, MA 02110-1301, USA.
+*/
+
+/** \file
+ * The description of a 3D cuboid neighbourhood -- header file.
+ */
+
+
+#ifndef CUBOID_NEIGHBOURHOOD_3D_H
+#define CUBOID_NEIGHBOURHOOD_3D_H
+
+#include "communication/mpiManager.h"
+#include
+#include "communication/superStructure3D.h"
+
+
+/// All OpenLB code is contained in this namespace.
+namespace olb {
+
+template class SuperStructure3D;
+
+
+/// Single 3D cuboid neighbourhoods are the basic component of a
+/// 3D communicator
+/** For each cuboid a cuboid neighbourhood is defined. It stores the
+ * needed cell coordinates (Cell3D) of other cuboids. Futhermore this
+ * class provides basic initialization and communication methods
+ * for the class Communicator3D.
+ *
+ * WARNING: For unstructured grids there is an interpolation needed
+ * for the method buffer_outData which is not yet implemented!
+ *
+ * This class is not intended to be derived from.
+ */
+template
+struct Cell3D {
+ Cell3D() : latticeR {0,0,0,0} {};
+ // local position latticeR: iC, iX, iY, iZ;
+ int latticeR[4];
+ // global position physR: x, y, z;
+ T physR[3];
+
+ bool operator==(Cell3D const& rhs) const
+ {
+ return latticeR == rhs.latticeR;
+ };
+
+ /// Copy constructor
+ Cell3D(Cell3D const& rhs) = default;
+
+ /// Copy assignment
+ Cell3D& operator=(Cell3D const& rhs) = default;
+
+};
+
+
+template
+class CuboidNeighbourhood3D {
+private:
+ /// Cuboid ID
+ int _iCglob;
+ /// Number of cubboids in the structure
+ int _nC;
+ /// Delta of the cuboid
+ T _deltaC;
+ /// Reference to the super structure
+ SuperStructure3D& _superStructure;
+
+ /// Number of data to be transfered
+ int _nData;
+ /// Size of underlying data type
+ int _nDataType;
+
+ /// Internal needed cells
+ std::vector > _inCells;
+ /// External needed cells
+ std::vector > _outCells;
+ /// Internal needed neighbour cuboid and number of cells
+ std::vector _inC;
+ std::vector _inN;
+ /// External needed neighbour cuboid and number of cells
+ std::vector _outC;
+ std::vector _outN;
+ /// Buffer for the internal needed data
+ bool **_inData;
+ /// Buffer for the external needed data
+ bool **_outData;
+ /// Buffer for the internal needed data
+ T **_inDataCoordinates;
+ /// Buffer for the external needed data
+ T **_outDataCoordinates;
+ /// Indecates if there was an initialization done
+ bool _initInCNdone;
+ int* _tempInCN;
+ /// Indecates if there was an initialization done
+ bool _initOutCNdone;
+ /// Handels the MPI communication
+#ifdef PARALLEL_MODE_MPI
+ singleton::MpiNonBlockingHelper _mpiNbHelper;
+#endif
+public:
+ /// Constructor
+ //CuboidNeighbourhood2D() {};
+ /// Constructor
+ CuboidNeighbourhood3D(SuperStructure3D& superStructure, int iC);
+ /// Copy construction
+ CuboidNeighbourhood3D(CuboidNeighbourhood3D const& rhs);
+ /// Copy assignment
+ CuboidNeighbourhood3D operator=(CuboidNeighbourhood3D rhs);
+ /// Destructor
+ ~CuboidNeighbourhood3D();
+
+ /// Read only access to _inCells
+ Cell3D const& get_inCell(int i) const;
+ /// Returns the number of cells in _inCells
+ int get_inCellsSize() const;
+ /// Read only access to _inC
+ int const& get_inC(int i) const;
+ /// Returns the number of cells in _inC
+ int get_inCsize() const;
+ /// Read and write access to **_inData
+ bool** get_inData();
+ /// Read and write access to **_outData
+ bool** get_outData();
+
+ /// Adds a cell to the vector _inCells
+ void add_inCell(Cell3D cell);
+ /// Adds a cell to the vector _outCells
+ void add_outCell(Cell3D cell);
+ /// Adds a cell to the vector _inCells
+ /// if the cell is not already there and
+ /// if there is another cuboid which can deliver the information
+ void add_inCell(int iX, int iY, int iZ);
+ /// Adds all cells with the distance overlap*_delta to
+ /// the vector _inCells
+ void add_inCells(int overlap);
+ /// Initializes _inC and _inN
+ void init_inCN();
+ /// Initializes _outC and _outN
+ void init_outCN();
+ /// Initialization Helper
+#ifdef PARALLEL_MODE_MPI
+ void bufSend_inCells(singleton::MpiNonBlockingHelper& helper);
+#else
+ void bufSend_inCells() {}
+#endif
+ /// Initialization Helper
+ void recWrite_outCells();
+ /// Finishes a communication step
+ void finish_comm();
+
+ /// Buffers data to be send
+ // WARNING: Here is interpolation needed if globX, globY
+ // are not integers. This needs to be fixed if one will
+ // use unstructured grids.
+ void buffer_outData();
+ /// Sends data (pure mpi method)
+ void send_outData();
+ /// Receives data (pure mpi method)
+ void receive_inData();
+ /// Writes all data to the corresponding lattice cells
+ void write_inData();
+ /// Resets neighbourhood after initialization (init_inCN
+ /// and init_outCN)
+ void reset();
+};
+
+} // namespace olb
+
+#endif
--
cgit v1.2.3