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
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
|
/* This file is part of the OpenLB library
*
* Copyright (C) 2014 Mathias J. Krause
* E-mail contact: info@openlb.net
* The most recent release of OpenLB can be downloaded at
* <http://www.openlb.net/>
*
* 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
* Representation of a 3d block geometry structure -- header file.
*/
#ifndef BLOCK_GEOMETRY_STRUCTURE_3D_H
#define BLOCK_GEOMETRY_STRUCTURE_3D_H
#include <vector>
#include "core/vector.h"
#include "geometry/blockGeometryStatistics3D.h"
#include "core/blockStructure3D.h"
#include "io/ostreamManager.h"
/// All OpenLB code is contained in this namespace.
namespace olb {
template<typename T> class IndicatorF3D;
/// Representation of a block geometry structure
/**
* This pure virtual class provides an interface for the classes
* block geometry and block geometry view. It presents a volume
* of voxels. Different types are given my material numbers
* which is imporant e.g. to work with different boundaries
* (like for inflow/output regions). Several methods are
* provided to rename material numbers of a type at several
* nodes at once. Further clearing methods are provided to
* rename material in order to obtain a robust mesh as well
* as a method which checks the mesh for possible errors. The
* member variable statistics provides integral information
* about the mesh and their materials.
*
* This class is intended to be derived from.
*/
template<typename T>
class BlockGeometryStructure3D {
protected:
/// Number of the cuboid, default=-1
int _iCglob;
/// Statistic class
BlockGeometryStatistics3D<T> _statistics;
/// class specific output stream
mutable OstreamManager clout;
public:
/// Constructor
BlockGeometryStructure3D(int iCglob=-1);
/// dtor
virtual ~BlockGeometryStructure3D() {};
/// Returns the underlying block structure
virtual BlockStructure3D& getBlockStructure() = 0;
/// Read only access to the global iC number which is given !=-1 if the block geometries are part of a super geometry
virtual int const& getIcGlob() const;
/// Read and write access to the statistic object
virtual BlockGeometryStatistics3D<T>& getStatistics(bool verbose=true) = 0;
/// Read only access to the statistic object
virtual BlockGeometryStatistics3D<T> const& getStatistics(bool verbose=true) const = 0;
/// Returns the position of the block origin which is the node (iX=0/iY=0/iZ=0) in physical units (meter)
virtual Vector<T,3> getOrigin() const = 0;
/// Returns the extend in x direction of the block in lattice units
virtual int getNx() const = 0;
/// Returns the extend in y direction of the block in lattice units
virtual int getNy() const = 0;
/// Returns the extend in z direction of the block in lattice units
virtual int getNz() const = 0;
/// Returns the extend of the block in lattice units
virtual Vector<int,3> const getExtend() const;
/// Returns the spacing in physical units (meter)
virtual const T getDeltaR() const = 0;
/// Transforms lattice to physical coordinates (wrapped from cuboid geometry)
virtual void getPhysR(T physR[3], const int& iX, const int& iY, const int& iZ) const = 0;
/// Transforms lattice to physical coordinates (wrapped from cuboid geometry)
virtual void getPhysR(T physR[3], const int latticeR[3]) const;
// TODO to be removed old once
/// returns the (iX,iY,iZ) entry in the 3D scalar field
virtual int getMaterial(int iX, int iY, int iZ) const = 0;
/// Returns true iff lattice location is inside block geometry
virtual bool isInside(int iX, int iY, int iZ) const;
/// Write access to a material number
virtual int& get(int iX, int iY, int iZ) = 0;
/// Read only access to a material number
virtual int const& get(int iX, int iY, int iZ) const = 0;
/// Write access to a material number
virtual int& get(std::vector<int> latticeR);
/// Read only access to a material number
virtual int const& get(std::vector<int> latticeR) const;
/// Changes all materials which are not 0 or 1 to 0 if there is no neighbour with material 1
virtual int clean(bool verbose=true);
/// Changes all materials which are not 0,1,material to 0 if there is no neighbour with material 1
virtual int clean(int material, bool verbos=true);
/// Changes all materials with material 1 to 0 if there is a neighbour with material 0
virtual int outerClean(bool verbose=true);
/// Changes all materials which are not 0 or 1 to 1 if there is a non robust constellation
virtual int innerClean(bool verbose=true);
/// Changes all materials with material fromM to 1 if there is a non robust constellation
virtual int innerClean(int fromM, bool verbose=true);
/// Returns the coordinates (iX,iY,iZ) of a voxel with a given material number (material) if there exists an neighbourhood of size (offsetX,offsetY,offsetZ) only with voxels of the given material number
virtual bool find(int material, unsigned offsetX, unsigned offsetY, unsigned offsetZ, int& iX, int& iY, int& iZ);
/// Returns true if at position (iX,iY,iZ) and in a neighbourhood of size (offsetX,offsetY,offsetZ) only voxels with a given material number (material) are there
virtual bool check(int material, int iX, int iY, int iZ, unsigned offsetX, unsigned offsetY, unsigned offsetZ);
/// Checks for errors (searches for all outer voxels (=0) with an inner voxel (=1) as a direct neighbour)
virtual bool checkForErrors(bool verbose=true) const;
/// Replaces all material numbers (fromM) to another (toM)
virtual void rename(int fromM, int toM);
/// Replaces all material numbers (fromM) to another (toM) if an indicator functor condition is fulfilled
virtual void rename(int fromM, int toM, IndicatorF3D<T>& condition);
/// Replaces all material numbers (fromM) to another (toM) if all materials in the neighbourhood (iX-offsetX,..,iX,..,ix+offsetX), .. are of the original material number (fromM)
virtual void rename(int fromM, int toM, unsigned offsetX, unsigned offsetY, unsigned offsetZ);
/// Replaces all material numbers (fromM) to another (toM) if all materials in the neighbourhood (iX+1,iX+2,..,ix+testDirection[0]), .. are of another material number (testM)
virtual void rename(int fromM, int toM, int testM, std::vector<int> testDirection);
/// Replaces all material numbers (fromM) to another (toM) if all materials in the neighbourhood (iX+discreteNormal[0],iX+2*discreteNormal[0]), .. are of another material number (testM) and if an indicator functor condition is fulfilled
virtual void rename(int fromM, int toM, int fluidM, IndicatorF3D<T>& condition, std::vector<int> discreteNormal);
/// Replaces all material numbers (fromM) to another (toM) if all materials in the neighbourhood (iX+discreteNormal[0],iX+2*discreteNormal[0]), .. are of another material number (fluidM) and if an indicator functor condition is fulfilled, the discreteNormal is computed from all fromM which fulfill the indicator functor condition
virtual void rename(int fromM, int toM, int fluidM, IndicatorF3D<T>& condition);
/// Copy a layer of material numbers inside an indicator in a discrete normal direction
virtual void copyMaterialLayer(IndicatorF3D<T>& condition, int discreteNormal[3], int numberOfLayers);
/// Replaces all material numbers (fromM) to another (toM) using a seed point and max. directions indicated by offsetX,Y,Z != 0
virtual void regionGrowing(int fromM, int toM, int seedX, int seedY, int seedZ, int offsetX, int offsetY, int offsetZ, std::map<std::vector<int>, int >* tmp=nullptr);
public:
/// Adds a pointer to the flag statistic status which is depended on the data of an block geometry
virtual void addToStatisticsList(bool* statisticStatus) = 0;
/// Removes a pointer if it exists from the flag statistic status which is depended on the data of an block geometry
virtual void removeFromStatisticsList(bool* statisticStatus) = 0;
};
} // namespace olb
#endif
|
