olb::BlockGeometryStatistics2D< T > Class Template Reference

#include <blockGeometryStatistics2D.h>

Collaboration diagram for olb::BlockGeometryStatistics2D< T >:

Public Member Functions
	BlockGeometryStatistics2D (BlockGeometry< T, 2 > *blockGeometry)
bool &	getStatisticsStatus ()
	Read and write access to a flag, which indicates if an uptate is needed (=true)
bool const &	getStatisticsStatus () const
	Read only access to a flag, which indicates if an uptate is needed (=true)
std::map< int, int >	getMaterial2n ()
	Returns the map with the numbers of voxels for each material.
std::map< int, int >	getMaterial2n () const
void	update (bool verbose=true)
	Updates the statistics if it is really needed.
int	getNmaterials ()
	Returns the number of different materials.
int	getNmaterials () const
int	getNvoxel (int material)
	Returns the number of voxels for a given material number.
int	getNvoxel (int material) const
int	getNvoxel ()
	Returns the number of voxels with material!=0.
int	getNvoxel () const
std::vector< int >	getMinLatticeR (int material)
	Returns the min. lattice position in each direction.
std::vector< int >	getMinLatticeR (int material) const
std::vector< int >	getMaxLatticeR (int material)
	Returns the max. lattice position in each direction.
std::vector< int >	getMaxLatticeR (int material) const
std::vector< T >	getMinPhysR (int material) const
	Returns the min. phys position in each direction.
std::vector< T >	getMaxPhysR (int material) const
	Returns the max. phys position in each direction.
std::vector< T >	getLatticeExtend (int material)
	Returns the lattice extend as length in each direction.
std::vector< T >	getLatticeExtend (int material) const
std::vector< T >	getPhysExtend (int material)
	Returns the phys extend as length in each direction.
std::vector< T >	getPhysExtend (int material) const
std::vector< T >	getPhysRadius (int material)
	Returns the phys radius as length in each direction.
std::vector< T >	getPhysRadius (int material) const
std::vector< T >	getCenterPhysR (int material)
	Returns the center position.
std::vector< T >	getCenterPhysR (int material) const
std::vector< int >	getType (const int *input) const
	Returns the boundary type which is characterized by a discrete normal (c.f. Zimny)
std::vector< int >	getType (int iX, int iY) const
std::vector< int >	computeNormal (int iX, int iY)
	Returns normal that points into the fluid for paraxial surfaces.
std::vector< int >	computeNormal (int iX, int iY) const
std::vector< T >	computeNormal (int material)
	Returns normal that points into the fluid for paraxial surfaces.
std::vector< T >	computeNormal (int material) const
std::vector< int >	computeDiscreteNormal (int material, T maxNorm=1.1)
	Returns discrete normal with norm maxNorm that points into the fluid for paraxial surfaces maxNorm=1.1 implies only normals parallel to the axises.
std::vector< int >	computeDiscreteNormal (int material, T maxNorm=1.1) const
bool	check (int material, int iX, int iY, unsigned offsetX, unsigned offsetY)
bool	check (int material, int iX, int iY, unsigned offsetX, unsigned offsetY) const
bool	find (int material, unsigned offsetX, unsigned offsetY, int &iX, int &iY)
bool	find (int material, unsigned offsetX, unsigned offsetY, int &iX, int &iY) const
void	print ()
	Prints some statistic information, i.e. the number of voxels and min. max. physical position for each different material.
void	print () const

Detailed Description

template<typename T>
class olb::BlockGeometryStatistics2D< T >

Definition at line 53 of file blockGeometryStatistics2D.h.

Constructor & Destructor Documentation

BlockGeometryStatistics2D()

template<typename T >

olb::BlockGeometryStatistics2D< T >::BlockGeometryStatistics2D

(

BlockGeometry< T, 2 > *

blockGeometry

)

Definition at line 43 of file blockGeometryStatistics2D.hh.

  : _blockGeometry(blockGeometry),
    clout(std::cout,"BlockGeometryStatistics2D")
{
  _statisticsUpdateNeeded = true;
}

Member Function Documentation

check() [1/2]

template<typename T >

bool olb::BlockGeometryStatistics2D< T >::check	(	int	material,
		int	iX,
		int	iY,
		unsigned	offsetX,
		unsigned	offsetY )

Definition at line 509 of file blockGeometryStatistics2D.hh.

{
  return const_this->check(material, iX, iY, offsetX, offsetY);
}

check() [2/2]

template<typename T >

bool olb::BlockGeometryStatistics2D< T >::check	(	int	material,
		int	iX,
		int	iY,
		unsigned	offsetX,
		unsigned	offsetY ) const

Definition at line 516 of file blockGeometryStatistics2D.hh.

{
  bool found = true;
  for (int iOffsetX = -offsetX; iOffsetX <= (int) offsetX; ++iOffsetX) {
    for (int iOffsetY = -offsetY; iOffsetY <= (int) offsetY; ++iOffsetY) {
      if (_blockGeometry->getMaterial({iX + iOffsetX, iY + iOffsetY}) != material) {
        found = false;
      }
    }
  }
  return found;
}

computeDiscreteNormal() [1/2]

template<typename T >

std::vector< int > olb::BlockGeometryStatistics2D< T >::computeDiscreteNormal	(	int	material,
		T	maxNorm = 1.1 )

Returns discrete normal with norm maxNorm that points into the fluid for paraxial surfaces maxNorm=1.1 implies only normals parallel to the axises.

Definition at line 480 of file blockGeometryStatistics2D.hh.

{
  return const_this->computeDiscreteNormal(material, maxNorm);
}

computeDiscreteNormal() [2/2]

template<typename T >

std::vector< int > olb::BlockGeometryStatistics2D< T >::computeDiscreteNormal	(	int	material,
		T	maxNorm = 1.1 ) const

Definition at line 486 of file blockGeometryStatistics2D.hh.

{
  std::vector<T> normal = computeNormal(material);
  std::vector<int> discreteNormal(2,int(0));
 
  T smallestAngle = T(0);
  for (int iX = -1; iX<=1; iX++) {
    for (int iY = -1; iY<=1; iY++) {
      T norm = util::sqrt(iX*iX+iY*iY);
      if (norm>0.&& norm<maxNorm) {
        T angle = (iX*normal[0] + iY*normal[1])/norm;
        if (angle>=smallestAngle) {
          smallestAngle=angle;
          discreteNormal[0] = iX;
          discreteNormal[1] = iY;
        }
      }
    }
  }
  return discreteNormal;
}

References olb::norm(), and olb::util::sqrt().

Here is the call graph for this function:

computeNormal() [1/4]

template<typename T >

std::vector< int > olb::BlockGeometryStatistics2D< T >::computeNormal	(	int	iX,
		int	iY )

Returns normal that points into the fluid for paraxial surfaces.

Definition at line 418 of file blockGeometryStatistics2D.hh.

{
  return const_this->computeNormal(iX, iY);
}

computeNormal() [2/4]

template<typename T >

std::vector< int > olb::BlockGeometryStatistics2D< T >::computeNormal	(	int	iX,
		int	iY ) const

Definition at line 424 of file blockGeometryStatistics2D.hh.

{
  std::vector<int> normal (2,int(0));
 
  if (iX != 0) {
    if (_blockGeometry->getMaterial({iX - 1, iY}) == 1) {
      normal[0] = -1;
    }
  }
  if (iX != _nX - 1) {
    if (_blockGeometry->getMaterial({iX + 1, iY}) == 1) {
      normal[0] = 1;
    }
  }
  if (iY != 0) {
    if (_blockGeometry->getMaterial({iX, iY - 1}) == 1) {
      normal[1] = -1;
    }
  }
  if (iY != _nY - 1) {
    if (_blockGeometry->getMaterial({iX, iY + 1}) == 1) {
      normal[1] = 1;
    }
  }
  return normal;
}

computeNormal() [3/4]

template<typename T >

std::vector< T > olb::BlockGeometryStatistics2D< T >::computeNormal

(

int

material

)

Returns normal that points into the fluid for paraxial surfaces.

Definition at line 452 of file blockGeometryStatistics2D.hh.

{
  return const_this->computeNormal(material);
}

computeNormal() [4/4]

template<typename T >

std::vector< T > olb::BlockGeometryStatistics2D< T >::computeNormal

(

int

material

)

const

Definition at line 458 of file blockGeometryStatistics2D.hh.

{
  std::vector<T> normal (2,int(0));
  std::vector<int> minC = getMinLatticeR(material);
  std::vector<int> maxC = getMaxLatticeR(material);
  for (int iX = minC[0]; iX<=maxC[0]; iX++) {
    for (int iY = minC[1]; iY<=maxC[1]; iY++) {
      if (_blockGeometry->getMaterial({iX,iY}) == material) {
        normal[0]+=computeNormal(iX,iY)[0];
        normal[1]+=computeNormal(iX,iY)[1];
      }
    }
  }
  T norm = util::sqrt(normal[0]*normal[0]+normal[1]*normal[1]);
  if (norm>0.) {
    normal[0]/=norm;
    normal[1]/=norm;
  }
  return normal;
}

References olb::norm(), and olb::util::sqrt().

Here is the call graph for this function:

find() [1/2]

template<typename T >

bool olb::BlockGeometryStatistics2D< T >::find	(	int	material,
		unsigned	offsetX,
		unsigned	offsetY,
		int &	iX,
		int &	iY )

Definition at line 531 of file blockGeometryStatistics2D.hh.

{
  return const_this->find(material, offsetX, offsetY, foundX, foundY);
}

References olb::BlockGeometryStatistics2D< T >::find().

Here is the call graph for this function:

Here is the caller graph for this function:

find() [2/2]

template<typename T >

bool olb::BlockGeometryStatistics2D< T >::find	(	int	material,
		unsigned	offsetX,
		unsigned	offsetY,
		int &	iX,
		int &	iY ) const

Definition at line 538 of file blockGeometryStatistics2D.hh.

{
  bool found = false;
  for (foundX = 0; foundX < _nX; foundX++) {
    for (foundY = 0; foundY < _nY; foundY++) {
      found = check(material, foundX, foundY, offsetX, offsetY);
      if (found) {
        return found;
      }
    }
  }
  return found;
}

getCenterPhysR() [1/2]

template<typename T >

std::vector< T > olb::BlockGeometryStatistics2D< T >::getCenterPhysR

(

int

material

)

Returns the center position.

Definition at line 262 of file blockGeometryStatistics2D.hh.

{
  update();
  return const_this->getCenterPhysR(material);
}

getCenterPhysR() [2/2]

template<typename T >

std::vector< T > olb::BlockGeometryStatistics2D< T >::getCenterPhysR

(

int

material

)

const

Definition at line 269 of file blockGeometryStatistics2D.hh.

{
  std::vector<T> center;
  for (int iDim=0; iDim<2; iDim++) {
    center.push_back(getMinPhysR(material)[iDim] + getPhysRadius(material)[iDim]);
  }
  return center;
}

getLatticeExtend() [1/2]

template<typename T >

std::vector< T > olb::BlockGeometryStatistics2D< T >::getLatticeExtend

(

int

material

)

Returns the lattice extend as length in each direction.

Definition at line 205 of file blockGeometryStatistics2D.hh.

{
  update();
  return const_this->getLatticeExtend(material);
}

getLatticeExtend() [2/2]

template<typename T >

std::vector< T > olb::BlockGeometryStatistics2D< T >::getLatticeExtend

(

int

material

)

const

Definition at line 212 of file blockGeometryStatistics2D.hh.

{
  try {
    std::vector<T> extend;
    for (int iDim = 0; iDim < 2; iDim++) {
      extend.push_back(_material2max.at(material)[iDim] - _material2min.at(material)[iDim]);
    }
    return extend;
  }
  catch (std::out_of_range& ex) {
    std::vector<T> null;
    return null;
  }
}

getMaterial2n() [1/2]

template<typename T >

std::map< int, int > olb::BlockGeometryStatistics2D< T >::getMaterial2n

(

)

Returns the map with the numbers of voxels for each material.

Definition at line 121 of file blockGeometryStatistics2D.hh.

{
  update();
  return const_this->getMaterial2n();
}

getMaterial2n() [2/2]

template<typename T >

std::map< int, int > olb::BlockGeometryStatistics2D< T >::getMaterial2n

(

)

const

Definition at line 128 of file blockGeometryStatistics2D.hh.

{
  return _material2n;
}

getMaxLatticeR() [1/2]

template<typename T >

std::vector< int > olb::BlockGeometryStatistics2D< T >::getMaxLatticeR

(

int

material

)

Returns the max. lattice position in each direction.

Definition at line 170 of file blockGeometryStatistics2D.hh.

{
  update();
  return const_this->getMaxLatticeR(material);
}

getMaxLatticeR() [2/2]

template<typename T >

std::vector< int > olb::BlockGeometryStatistics2D< T >::getMaxLatticeR

(

int

material

)

const

Definition at line 177 of file blockGeometryStatistics2D.hh.

{
  try {
    return _material2max.at(material);
  }
  catch (std::out_of_range& ex) {
    std::vector<int> null;
    return null;
  }
}

getMaxPhysR()

template<typename T >

std::vector< T > olb::BlockGeometryStatistics2D< T >::getMaxPhysR

(

int

material

)

const

Returns the max. phys position in each direction.

Definition at line 197 of file blockGeometryStatistics2D.hh.

{
  std::vector<T> tmp(2,T());
  _blockGeometry->getPhysR(&(tmp[0]), &(getMaxLatticeR(material)[0]));
  return tmp;
}

getMinLatticeR() [1/2]

template<typename T >

std::vector< int > olb::BlockGeometryStatistics2D< T >::getMinLatticeR

(

int

material

)

Returns the min. lattice position in each direction.

Definition at line 151 of file blockGeometryStatistics2D.hh.

{
  update();
  return const_this->getMinLatticeR(material);
}

getMinLatticeR() [2/2]

template<typename T >

std::vector< int > olb::BlockGeometryStatistics2D< T >::getMinLatticeR

(

int

material

)

const

Definition at line 158 of file blockGeometryStatistics2D.hh.

{
  try {
    return _material2min.at(material);
  }
  catch (std::out_of_range& ex) {
    std::vector<int> null;
    return null;
  }
}

getMinPhysR()

template<typename T >

std::vector< T > olb::BlockGeometryStatistics2D< T >::getMinPhysR

(

int

material

)

const

Returns the min. phys position in each direction.

Definition at line 189 of file blockGeometryStatistics2D.hh.

{
  std::vector<T> tmp(2,T());
  _blockGeometry->getPhysR(&(tmp[0]), &(getMinLatticeR(material)[0]));
  return tmp;
}

getNmaterials() [1/2]

template<typename T >

int olb::BlockGeometryStatistics2D< T >::getNmaterials

(

)

Returns the number of different materials.

Definition at line 90 of file blockGeometryStatistics2D.hh.

{
  update();
  return const_this->getNmaterials();
}

getNmaterials() [2/2]

template<typename T >

int olb::BlockGeometryStatistics2D< T >::getNmaterials

(

)

const

Definition at line 97 of file blockGeometryStatistics2D.hh.

{
  return _nMaterials;
}

getNvoxel() [1/4]

template<typename T >

int olb::BlockGeometryStatistics2D< T >::getNvoxel

(

)

Returns the number of voxels with material!=0.

Definition at line 134 of file blockGeometryStatistics2D.hh.

{
  update();
  return const_this->getNvoxel();
}

getNvoxel() [2/4]

template<typename T >

int olb::BlockGeometryStatistics2D< T >::getNvoxel

(

)

const

Definition at line 141 of file blockGeometryStatistics2D.hh.

{
  int total = 0;
  for (const auto& material : _material2n ) {
    total += material.second;
  }
  return total;
}

getNvoxel() [3/4]

template<typename T >

int olb::BlockGeometryStatistics2D< T >::getNvoxel

(

int

material

)

Returns the number of voxels for a given material number.

Definition at line 103 of file blockGeometryStatistics2D.hh.

{
  update();
  return const_this->getNvoxel(material);
}

getNvoxel() [4/4]

template<typename T >

int olb::BlockGeometryStatistics2D< T >::getNvoxel

(

int

material

)

const

Definition at line 110 of file blockGeometryStatistics2D.hh.

{
  try {
    return _material2n.at(material);
  }
  catch (std::out_of_range& ex) {
    return 0;
  }
}

getPhysExtend() [1/2]

template<typename T >

std::vector< T > olb::BlockGeometryStatistics2D< T >::getPhysExtend

(

int

material

)

Returns the phys extend as length in each direction.

Definition at line 228 of file blockGeometryStatistics2D.hh.

{
  update();
  return const_this->getPhysExtend(material);
}

getPhysExtend() [2/2]

template<typename T >

std::vector< T > olb::BlockGeometryStatistics2D< T >::getPhysExtend

(

int

material

)

const

Definition at line 235 of file blockGeometryStatistics2D.hh.

{
  std::vector<T> extend;
  for (int iDim = 0; iDim < 2; iDim++) {
    extend.push_back(getMaxPhysR(material)[iDim] - getMinPhysR(material)[iDim]);
  }
  return extend;
}

getPhysRadius() [1/2]

template<typename T >

std::vector< T > olb::BlockGeometryStatistics2D< T >::getPhysRadius

(

int

material

)

Returns the phys radius as length in each direction.

Definition at line 245 of file blockGeometryStatistics2D.hh.

{
  update();
  return const_this->getPhysRadius(material);
}

getPhysRadius() [2/2]

template<typename T >

std::vector< T > olb::BlockGeometryStatistics2D< T >::getPhysRadius

(

int

material

)

const

Definition at line 252 of file blockGeometryStatistics2D.hh.

{
  std::vector<T> radius;
  for (int iDim=0; iDim<2; iDim++) {
    radius.push_back((getMaxPhysR(material)[iDim] - getMinPhysR(material)[iDim])/2.);
  }
  return radius;
}

getStatisticsStatus() [1/2]

template<typename T >

bool & olb::BlockGeometryStatistics2D< T >::getStatisticsStatus

(

)

Read and write access to a flag, which indicates if an uptate is needed (=true)

Definition at line 52 of file blockGeometryStatistics2D.hh.

{
  return _statisticsUpdateNeeded;
}

getStatisticsStatus() [2/2]

template<typename T >

bool const & olb::BlockGeometryStatistics2D< T >::getStatisticsStatus

(

)

const

Read only access to a flag, which indicates if an uptate is needed (=true)

Definition at line 58 of file blockGeometryStatistics2D.hh.

{
  return _statisticsUpdateNeeded;
}

getType() [1/2]

template<typename T >

std::vector< int > olb::BlockGeometryStatistics2D< T >::getType

(

const int *

input

)

const

Returns the boundary type which is characterized by a discrete normal (c.f. Zimny)

Definition at line 279 of file blockGeometryStatistics2D.hh.

{
  return const_this->getType(input[0], input[1]);
}

getType() [2/2]

template<typename T >

std::vector< int > olb::BlockGeometryStatistics2D< T >::getType	(	int	iX,
		int	iY ) const

boundary0N and boundary 0P

boundary1N and boundary1P

externalCornerNN and externalCornerNP

externalCornerPN and externalCornerPP

internalCornerNN and internalCornerNP

internalCornerPN and internalCornerPP

Definition at line 285 of file blockGeometryStatistics2D.hh.

{
  std::vector<int> discreteNormal(3, 0);
 
  if (_blockGeometry->getMaterial({iX, iY}) != 1
      && _blockGeometry->getMaterial({iX, iY}) != 0) {
 
    if (_blockGeometry->getMaterial({iX, iY + 1}) != 1
        && _blockGeometry->getMaterial({iX, iY + 1}) != 0
        && _blockGeometry->getMaterial({iX, iY - 1}) != 1
        && _blockGeometry->getMaterial({iX, iY - 1}) != 0) {
 
      if (_blockGeometry->getMaterial({iX + 1, iY}) == 1) {
        discreteNormal[0] = 0;
        discreteNormal[1] = -1;
        discreteNormal[2] = 0;
      }
 
      if (_blockGeometry->getMaterial({iX - 1, iY}) == 1) {
        discreteNormal[0] = 0;
        discreteNormal[1] = 1;
        discreteNormal[2] = 0;
      }
    }
 
    if (_blockGeometry->getMaterial({iX + 1, iY}) != 1
        && _blockGeometry->getMaterial({iX + 1, iY}) != 0
        && _blockGeometry->getMaterial({iX - 1, iY}) != 1
        && _blockGeometry->getMaterial({iX - 1, iY}) != 0) {
 
      if (_blockGeometry->getMaterial({iX, iY + 1}) == 1) {
        discreteNormal[0] = 0;
        discreteNormal[1] = 0;
        discreteNormal[2] = -1;
      }
 
      if (_blockGeometry->getMaterial({iX, iY - 1}) == 1) {
        discreteNormal[0] = 0;
        discreteNormal[1] = 0;
        discreteNormal[2] = 1;
      }
    }
 
    if (   _blockGeometry->getMaterial({iX + 1, iY}) != 1
        && _blockGeometry->getMaterial({iX + 1, iY}) != 0) {
 
      if (_blockGeometry->getMaterial({iX, iY + 1}) != 1
          && _blockGeometry->getMaterial({iX, iY + 1}) != 0
          && _blockGeometry->getMaterial({iX + 1, iY + 1}) == 1) {
        discreteNormal[0] = 1;
        discreteNormal[1] = -1;
        discreteNormal[2] = -1;
      }
 
      if (_blockGeometry->getMaterial({iX, iY - 1}) != 1
          && _blockGeometry->getMaterial({iX, iY - 1}) != 0
          && _blockGeometry->getMaterial({iX + 1, iY - 1}) == 1) {
        discreteNormal[0] = 1;
        discreteNormal[1] = -1;
        discreteNormal[2] = 1;
      }
    }
 
    if (_blockGeometry->getMaterial({iX - 1, iY}) != 1
        && _blockGeometry->getMaterial({iX - 1, iY}) != 0) {
 
      if (_blockGeometry->getMaterial({iX, iY + 1}) != 1
          && _blockGeometry->getMaterial({iX, iY + 1}) != 0
          && _blockGeometry->getMaterial({iX - 1, iY + 1}) == 1) {
        discreteNormal[0] = 1;
        discreteNormal[1] = 1;
        discreteNormal[2] = -1;
      }
 
      if (_blockGeometry->getMaterial({iX, iY - 1}) != 1
          && _blockGeometry->getMaterial({iX, iY - 1}) != 0
          && _blockGeometry->getMaterial({iX - 1, iY - 1}) == 1) {
        discreteNormal[0] = 1;
        discreteNormal[1] = 1;
        discreteNormal[2] = 1;
      }
    }
 
    if (_blockGeometry->getMaterial({iX - 1, iY}) != 1
        && _blockGeometry->getMaterial({iX - 1, iY}) != 0) {
 
      if (_blockGeometry->getMaterial({iX, iY - 1}) != 1
          && _blockGeometry->getMaterial({iX, iY - 1}) != 0
          && _blockGeometry->getMaterial({iX - 1, iY - 1}) == 0) {
        discreteNormal[0] = 2;
        discreteNormal[1] = -1;
        discreteNormal[2] = -1;
      }
 
      if (_blockGeometry->getMaterial({iX, iY + 1}) != 1
          && _blockGeometry->getMaterial({iX, iY + 1}) != 0
          && _blockGeometry->getMaterial({iX - 1, iY + 1}) == 0) {
        discreteNormal[0] = 2;
        discreteNormal[1] = -1;
        discreteNormal[2] = 1;
      }
    }
 
    if (_blockGeometry->getMaterial({iX + 1, iY}) != 1
        && _blockGeometry->getMaterial({iX + 1, iY}) != 0) {
 
      if (_blockGeometry->getMaterial({iX, iY - 1}) != 1
          && _blockGeometry->getMaterial({iX, iY - 1}) != 0
          && _blockGeometry->getMaterial({iX + 1, iY - 1}) == 0) {
        discreteNormal[0] = 2;
        discreteNormal[1] = 1;
        discreteNormal[2] = -1;
      }
 
      if (_blockGeometry->getMaterial({iX, iY + 1}) != 1
          && _blockGeometry->getMaterial({iX, iY + 1}) != 0
          && _blockGeometry->getMaterial({iX + 1, iY + 1}) == 0) {
        discreteNormal[0] = 2;
        discreteNormal[1] = 1;
        discreteNormal[2] = 1;
      }
    }
  }
  return discreteNormal;
}

print() [1/2]

template<typename T >

void olb::BlockGeometryStatistics2D< T >::print

(

)

Prints some statistic information, i.e. the number of voxels and min. max. physical position for each different material.

Definition at line 554 of file blockGeometryStatistics2D.hh.

{
  update();
  return const_this->print();
}

print() [2/2]

template<typename T >

void olb::BlockGeometryStatistics2D< T >::print

(

)

const

Definition at line 561 of file blockGeometryStatistics2D.hh.

{
  try {
    std::map<int, int>::iterator iter;
    for (const auto& material : _material2n) {
      clout << "materialNumber=" << material.first
            << "; count=" << material.second
            << "; minLatticeR=(" << _material2min.at(material.first)[0] <<","<< _material2min.at(material.first)[1] <<")"
            << "; maxLatticeR=(" << _material2max.at(material.first)[0] <<","<< _material2max.at(material.first)[1] <<")"
            << std::endl;
    }
  }
  catch (std::out_of_range& ex)
  { }
}

update()

template<typename T >

void olb::BlockGeometryStatistics2D< T >::update

(

bool

verbose = true

)

Updates the statistics if it is really needed.

Definition at line 65 of file blockGeometryStatistics2D.hh.

{
  const_this = const_cast<const BlockGeometryStatistics2D<T>*>(this);
 
  if (getStatisticsStatus() ) {
    _material2n.clear();
    _blockGeometry->forCoreSpatialLocations([&](auto iX, auto iY) {
      takeStatistics(iX, iY);
    });
 
    _nMaterials=int();
    std::map<int, int>::iterator iter;
    for (iter = _material2n.begin(); iter != _material2n.end(); iter++) {
      _nMaterials++;
    }
 
    if (verbose) {
      clout << "updated" << std::endl;
    }
    getStatisticsStatus() = false;
  }
}

---

The documentation for this class was generated from the following files:

• src/geometry/blockGeometryStatistics2D.h
• src/geometry/blockGeometryStatistics2D.hh