olb::SuperGeometryStatistics3D< T > Class Template Reference

#include <superGeometryStatistics3D.h>

Collaboration diagram for olb::SuperGeometryStatistics3D< T >:

Public Member Functions
	SuperGeometryStatistics3D (SuperGeometry< T, 3 > *superGeometry)
	Constructor.
	SuperGeometryStatistics3D (SuperGeometryStatistics3D const &rhs)
	Copy constructor.
SuperGeometryStatistics3D< T > &	operator= (SuperGeometryStatistics3D const &rhs)
	Copy assignment.
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)
void	update (bool verbose=false)
	Updates the statistics if it is really needed.
int	getNmaterials ()
	Returns the number of different materials.
int	getNmaterials () const
std::size_t	getNvoxel (int material)
	Returns the number of voxels for a given material number.
std::size_t	getNvoxel (int material) const
std::size_t	getNvoxel ()
	Returns the number of voxels with material!=0.
std::size_t	getNvoxel () const
olb::Vector< T, 3 >	getMinPhysR (int material)
	Returns the min. phys position in each direction corresponding to material number.
olb::Vector< T, 3 >	getMinPhysR (int material) const
olb::Vector< T, 3 >	getMinPhysR ()
	Returns the min. phys position in each direction corresponding to all non-zero material numbers.
olb::Vector< T, 3 >	getMinPhysR () const
olb::Vector< T, 3 >	getMaxPhysR (int material)
	Returns the max. phys position in each direction corresponding to material number.
olb::Vector< T, 3 >	getMaxPhysR (int material) const
olb::Vector< T, 3 >	getMaxPhysR ()
	Returns the max. phys position in each direction corresponding to all non-zero material numbers.
olb::Vector< T, 3 >	getMaxPhysR () const
olb::Vector< T, 3 >	getPhysExtend (int material)
	Returns the phys extend as length in each direction.
std::vector< T >	getPhysExtend (int material) const
olb::Vector< T, 3 >	getPhysRadius (int material)
	Returns the phys radius as length in each direction.
olb::Vector< T, 3 >	getPhysRadius (int material) const
olb::Vector< T, 3 >	getCenterPhysR (int material)
	Returns the center position.
olb::Vector< T, 3 >	getCenterPhysR (int material) const
olb::Vector< int, 3 >	getType (int iC, int iX, int iY, int iZ)
	Returns the boundary type which is characterized by a discrte normal (c.f. Zimny)
olb::Vector< int, 3 >	getType (int iC, int iX, int iY, int iZ) const
olb::Vector< T, 3 >	computeNormal (int material)
	Returns normal that points into the fluid for paraxial surfaces.
olb::Vector< T, 3 >	computeNormal (int material) const
olb::Vector< int, 3 >	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.
olb::Vector< int, 3 >	computeDiscreteNormal (int material, T maxNorm=1.1) const
T	computeMaxPhysDistance (int material) const
	Returns util::sqrt( maxX^2 + maxY^2 + maxZ^2 ) max over a certain material number.
T	computeMaxPhysDistance () const
	Returns util::sqrt( maxX^2 + maxY^2 + maxZ^2 ) max over all material numbers.
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::SuperGeometryStatistics3D< T >

Definition at line 57 of file superGeometryStatistics3D.h.

Constructor & Destructor Documentation

SuperGeometryStatistics3D() [1/2]

template<typename T >

olb::SuperGeometryStatistics3D< T >::SuperGeometryStatistics3D

(

SuperGeometry< T, 3 > *

superGeometry

)

Constructor.

Definition at line 42 of file superGeometryStatistics3D.hh.

  : _superGeometry(superGeometry), _statisticsUpdateNeeded(true),
    _nMaterials(0), _minOverMaterial(0), _maxOverMaterial(0),
    clout(std::cout,"SuperGeometryStatistics3D")
{
}

SuperGeometryStatistics3D() [2/2]

template<typename T >

olb::SuperGeometryStatistics3D< T >::SuperGeometryStatistics3D

(

SuperGeometryStatistics3D< T > const &

rhs

)

Copy constructor.

Definition at line 50 of file superGeometryStatistics3D.hh.

  : _superGeometry(rhs._superGeometry), _statisticsUpdateNeeded(true),
    _nMaterials(rhs._nMaterials),
    _minOverMaterial(rhs._minOverMaterial), _maxOverMaterial(rhs._maxOverMaterial),
    clout(std::cout,"SuperGeometryStatistics3D")
{
}

Member Function Documentation

computeDiscreteNormal() [1/2]

template<typename T >

olb::Vector< int, 3 > olb::SuperGeometryStatistics3D< 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 509 of file superGeometryStatistics3D.hh.

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

computeDiscreteNormal() [2/2]

template<typename T >

olb::Vector< int, 3 > olb::SuperGeometryStatistics3D< T >::computeDiscreteNormal	(	int	material,
		T	maxNorm = 1.1 ) const

Definition at line 516 of file superGeometryStatistics3D.hh.

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

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

Here is the call graph for this function:

computeMaxPhysDistance() [1/2]

template<typename T >

T olb::SuperGeometryStatistics3D< T >::computeMaxPhysDistance

(

)

const

Returns util::sqrt( maxX^2 + maxY^2 + maxZ^2 ) max over all material numbers.

Definition at line 549 of file superGeometryStatistics3D.hh.

{
  Vector<T,3> vec(getMaxPhysR()[0] -getMinPhysR()[0], getMaxPhysR()[1] -getMinPhysR()[1], getMaxPhysR()[2] -getMinPhysR()[2]);
  return norm(vec);
}

References olb::norm().

Here is the call graph for this function:

computeMaxPhysDistance() [2/2]

template<typename T >

T olb::SuperGeometryStatistics3D< T >::computeMaxPhysDistance

(

int

material

)

const

Returns util::sqrt( maxX^2 + maxY^2 + maxZ^2 ) max over a certain material number.

Definition at line 542 of file superGeometryStatistics3D.hh.

{
  Vector<T,3> vec(getMaxPhysR(material)[0] -getMinPhysR(material)[0], getMaxPhysR(material)[1] -getMinPhysR(material)[1], getMaxPhysR(material)[2] -getMinPhysR(material)[2]);
  return norm(vec);
}

References olb::norm().

Here is the call graph for this function:

computeNormal() [1/2]

template<typename T >

olb::Vector< T, 3 > olb::SuperGeometryStatistics3D< T >::computeNormal

(

int

material

)

Returns normal that points into the fluid for paraxial surfaces.

Definition at line 466 of file superGeometryStatistics3D.hh.

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

computeNormal() [2/2]

template<typename T >

olb::Vector< T, 3 > olb::SuperGeometryStatistics3D< T >::computeNormal

(

int

material

)

const

Definition at line 473 of file superGeometryStatistics3D.hh.

{
  std::vector<T> normal (3,int());
  for (int iCloc=0; iCloc<_superGeometry->getLoadBalancer().size(); iCloc++) {
    for (int iDim=0; iDim<3; iDim++) {
      if (_superGeometry->getBlockGeometry(iCloc).getStatistics().getNvoxel(material)!=0) {
        normal[iDim] += _superGeometry->getBlockGeometry(iCloc).getStatistics().computeNormal(material)[iDim]*_superGeometry->getBlockGeometry(iCloc).getStatistics().getNvoxel(material);
      }
    }
  }
 
#ifdef PARALLEL_MODE_MPI
  for (int iDim=0; iDim<3; iDim++) {
    singleton::mpi().reduceAndBcast((normal[iDim]), MPI_SUM);
  }
#endif
 
  if (getNvoxel(material) == 0) {
    std::cerr << "Unkown material number: " << material << std::endl;
    std::exit(-1);
  }
 
  for (int iDim=0; iDim<3; iDim++) {
    normal[iDim] /= getNvoxel(material);
  }
 
  T norm = util::sqrt(normal[0]*normal[0]+normal[1]*normal[1]+normal[2]*normal[2]);
  if (norm>0.) {
    normal[0]/=norm;
    normal[1]/=norm;
    normal[2]/=norm;
  }
  return normal;
}

References olb::singleton::mpi(), olb::norm(), olb::singleton::MpiManager::reduceAndBcast(), and olb::util::sqrt().

Here is the call graph for this function:

getCenterPhysR() [1/2]

template<typename T >

olb::Vector< T, 3 > olb::SuperGeometryStatistics3D< T >::getCenterPhysR

(

int

material

)

Returns the center position.

Definition at line 434 of file superGeometryStatistics3D.hh.

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

getCenterPhysR() [2/2]

template<typename T >

olb::Vector< T, 3 > olb::SuperGeometryStatistics3D< T >::getCenterPhysR

(

int

material

)

const

Definition at line 441 of file superGeometryStatistics3D.hh.

{
  olb::Vector<T, 3> center;
  for (int iDim=0; iDim<3; iDim++) {
    center[iDim] = getMinPhysR(material)[iDim] + getPhysRadius(material)[iDim];
  }
  return center;
}

getMaxPhysR() [1/4]

template<typename T >

olb::Vector< T, 3 > olb::SuperGeometryStatistics3D< T >::getMaxPhysR

(

)

Returns the max. phys position in each direction corresponding to all non-zero material numbers.

Definition at line 381 of file superGeometryStatistics3D.hh.

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

getMaxPhysR() [2/4]

template<typename T >

olb::Vector< T, 3 > olb::SuperGeometryStatistics3D< T >::getMaxPhysR

(

)

const

Definition at line 388 of file superGeometryStatistics3D.hh.

{
  return _maxOverMaterial;
}

getMaxPhysR() [3/4]

template<typename T >

olb::Vector< T, 3 > olb::SuperGeometryStatistics3D< T >::getMaxPhysR

(

int

material

)

Returns the max. phys position in each direction corresponding to material number.

Definition at line 363 of file superGeometryStatistics3D.hh.

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

getMaxPhysR() [4/4]

template<typename T >

olb::Vector< T, 3 > olb::SuperGeometryStatistics3D< T >::getMaxPhysR

(

int

material

)

const

Definition at line 370 of file superGeometryStatistics3D.hh.

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

getMinPhysR() [1/4]

template<typename T >

olb::Vector< T, 3 > olb::SuperGeometryStatistics3D< T >::getMinPhysR

(

)

Returns the min. phys position in each direction corresponding to all non-zero material numbers.

Definition at line 350 of file superGeometryStatistics3D.hh.

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

getMinPhysR() [2/4]

template<typename T >

olb::Vector< T, 3 > olb::SuperGeometryStatistics3D< T >::getMinPhysR

(

)

const

Definition at line 357 of file superGeometryStatistics3D.hh.

{
  return _minOverMaterial;
}

getMinPhysR() [3/4]

template<typename T >

olb::Vector< T, 3 > olb::SuperGeometryStatistics3D< T >::getMinPhysR

(

int

material

)

Returns the min. phys position in each direction corresponding to material number.

Definition at line 331 of file superGeometryStatistics3D.hh.

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

getMinPhysR() [4/4]

template<typename T >

olb::Vector< T, 3 > olb::SuperGeometryStatistics3D< T >::getMinPhysR

(

int

material

)

const

Definition at line 338 of file superGeometryStatistics3D.hh.

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

getNmaterials() [1/2]

template<typename T >

int olb::SuperGeometryStatistics3D< T >::getNmaterials

(

)

Returns the number of different materials.

Definition at line 281 of file superGeometryStatistics3D.hh.

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

getNmaterials() [2/2]

template<typename T >

int olb::SuperGeometryStatistics3D< T >::getNmaterials

(

)

const

Definition at line 288 of file superGeometryStatistics3D.hh.

{
  return _nMaterials;
}

getNvoxel() [1/4]

template<typename T >

std::size_t olb::SuperGeometryStatistics3D< T >::getNvoxel

(

)

Returns the number of voxels with material!=0.

Definition at line 312 of file superGeometryStatistics3D.hh.

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

getNvoxel() [2/4]

template<typename T >

std::size_t olb::SuperGeometryStatistics3D< T >::getNvoxel

(

)

const

Definition at line 319 of file superGeometryStatistics3D.hh.

{
  std::size_t total = 0;
  for (const auto& material : _material2n) {
    if (material.first!=0) {
      total += material.second;
    }
  }
  return total;
}

getNvoxel() [3/4]

template<typename T >

std::size_t olb::SuperGeometryStatistics3D< T >::getNvoxel

(

int

material

)

Returns the number of voxels for a given material number.

Definition at line 294 of file superGeometryStatistics3D.hh.

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

getNvoxel() [4/4]

template<typename T >

std::size_t olb::SuperGeometryStatistics3D< T >::getNvoxel

(

int

material

)

const

Definition at line 301 of file superGeometryStatistics3D.hh.

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

getPhysExtend() [1/2]

template<typename T >

olb::Vector< T, 3 > olb::SuperGeometryStatistics3D< T >::getPhysExtend

(

int

material

)

Returns the phys extend as length in each direction.

Definition at line 394 of file superGeometryStatistics3D.hh.

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

getPhysExtend() [2/2]

template<typename T >

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

(

int

material

)

const

Definition at line 401 of file superGeometryStatistics3D.hh.

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

getPhysRadius() [1/2]

template<typename T >

olb::Vector< T, 3 > olb::SuperGeometryStatistics3D< T >::getPhysRadius

(

int

material

)

Returns the phys radius as length in each direction.

Definition at line 416 of file superGeometryStatistics3D.hh.

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

getPhysRadius() [2/2]

template<typename T >

olb::Vector< T, 3 > olb::SuperGeometryStatistics3D< T >::getPhysRadius

(

int

material

)

const

Definition at line 423 of file superGeometryStatistics3D.hh.

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

getStatisticsStatus() [1/2]

template<typename T >

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

(

)

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

Definition at line 71 of file superGeometryStatistics3D.hh.

{
  return _statisticsUpdateNeeded;
}

getStatisticsStatus() [2/2]

template<typename T >

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

(

)

const

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

Definition at line 77 of file superGeometryStatistics3D.hh.

{
  return _statisticsUpdateNeeded;
}

getType() [1/2]

template<typename T >

olb::Vector< int, 3 > olb::SuperGeometryStatistics3D< T >::getType	(	int	iC,
		int	iX,
		int	iY,
		int	iZ )

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

Definition at line 451 of file superGeometryStatistics3D.hh.

{
  update();
  return const_this->getType(iC, iX, iY, iZ);
}

getType() [2/2]

template<typename T >

olb::Vector< int, 3 > olb::SuperGeometryStatistics3D< T >::getType	(	int	iC,
		int	iX,
		int	iY,
		int	iZ ) const

Definition at line 458 of file superGeometryStatistics3D.hh.

{
  int iCloc=_superGeometry->getLoadBalancer().loc(iC);
  olb::Vector<int, 3> discreteNormal = _superGeometry->getBlockGeometry(iCloc).getStatistics().getType(iX, iY, iZ);
  return discreteNormal;
}

operator=()

template<typename T >

SuperGeometryStatistics3D< T > & olb::SuperGeometryStatistics3D< T >::operator=

(

SuperGeometryStatistics3D< T > const &

rhs

)

Copy assignment.

Definition at line 59 of file superGeometryStatistics3D.hh.

{
  _superGeometry = rhs._superGeometry;
  _statisticsUpdateNeeded = true;
  _nMaterials = rhs._nMaterials;
  _minOverMaterial = rhs._minOverMaterial;
  _maxOverMaterial = rhs._maxOverMaterial;
  return *this;
}

print() [1/2]

template<typename T >

void olb::SuperGeometryStatistics3D< T >::print

(

)

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

Definition at line 556 of file superGeometryStatistics3D.hh.

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

print() [2/2]

template<typename T >

void olb::SuperGeometryStatistics3D< T >::print

(

)

const

Definition at line 563 of file superGeometryStatistics3D.hh.

{
  try {
    std::size_t nCells = 0;
    for (const auto& material : _material2n) {
      clout << "materialNumber=" << material.first
            << "; count=" << material.second
            << "; minPhysR=(" << _material2min.at(material.first)[0] <<","<< _material2min.at(material.first)[1] <<","<< _material2min.at(material.first)[2] <<")"
            << "; maxPhysR=(" << _material2max.at(material.first)[0] <<","<< _material2max.at(material.first)[1] <<","<< _material2max.at(material.first)[2] <<")"
            << std::endl;
      nCells += material.second;
    }
    clout << "countTotal[1e6]=" << nCells / 1.e6 << std::endl;
  }
  catch (std::out_of_range& ex) {
  }
}

update()

template<typename T >

void olb::SuperGeometryStatistics3D< T >::update

(

bool

verbose = false

)

Updates the statistics if it is really needed.

Definition at line 84 of file superGeometryStatistics3D.hh.

{
  const_this = const_cast<const SuperGeometryStatistics3D<T>*>(this);
 
#ifdef PARALLEL_MODE_MPI
  // This needs to be done with integers due to undesired behaviour with bool and LOR implicated by the MPI standard
  int updateReallyNeededGlobal = 0;
  if (_statisticsUpdateNeeded) {
    updateReallyNeededGlobal = 1;
  }
  singleton::mpi().reduceAndBcast(updateReallyNeededGlobal, MPI_SUM);
 
  if (updateReallyNeededGlobal>0) {
    _statisticsUpdateNeeded = true;
  }
  //singleton::mpi().reduceAndBcast(_statisticsUpdateNeeded, MPI_LOR);
#endif
 
  // check if update is really needed
  if (_statisticsUpdateNeeded ) {
    int updateReallyNeeded = 0;
    for (int iCloc=0; iCloc<_superGeometry->getLoadBalancer().size(); iCloc++) {
      if (_superGeometry->getBlockGeometry(iCloc).getStatistics().getStatisticsStatus() ) {
        auto& blockGeometry = const_cast<BlockGeometry<T,3>&>(_superGeometry->getBlockGeometry(iCloc));
        blockGeometry.getStatistics().update(false);
        updateReallyNeeded++;
      }
    }
 
 
#ifdef PARALLEL_MODE_MPI
    singleton::mpi().reduceAndBcast(updateReallyNeeded, MPI_SUM);
#endif
 
    if (updateReallyNeeded==0) {
      _statisticsUpdateNeeded = false;
      //      clout << "almost updated" << std::endl;
      return;
    }
 
 
    // get total number of different materials right
    _nMaterials = int();
    {
      std::set<int> tmpMaterials{};
      for (int iCloc=0; iCloc<_superGeometry->getLoadBalancer().size(); iCloc++) {
        const auto& blockMaterial2n = _superGeometry->getBlockGeometry(iCloc).getStatistics().getMaterial2n();
        for (auto [material, _] : blockMaterial2n) {
          tmpMaterials.insert(material);
        }
      }
      _nMaterials = tmpMaterials.size();
    }
 
    _material2n = std::map<int, std::size_t>();
 
#ifdef PARALLEL_MODE_MPI
    singleton::mpi().reduceAndBcast(_nMaterials, MPI_SUM);
#endif
 
    // store the number and min., max. possition for each rank
    for (int iCloc=0; iCloc<_superGeometry->getLoadBalancer().size(); iCloc++) {
      std::map<int, std::size_t> material2n = _superGeometry->getBlockGeometry(iCloc).getStatistics().getMaterial2n();
      std::map<int, std::size_t>::iterator iter;
 
      for (iter = material2n.begin(); iter != material2n.end(); iter++) {
        if (iter->second!=0) {
          std::vector<T> minPhysR = _superGeometry->getBlockGeometry(iCloc).getStatistics().getMinPhysR(iter->first);
          std::vector<T> maxPhysR = _superGeometry->getBlockGeometry(iCloc).getStatistics().getMaxPhysR(iter->first);
          if (_material2n.count(iter->first) == 0) {
            _material2n[iter->first] = iter->second;
            _material2min[iter->first] = minPhysR;
            _material2max[iter->first] = maxPhysR;
          }
          else {
            _material2n[iter->first] += iter->second;
            for (int iDim=0; iDim<3; iDim++) {
              if (_material2min[iter->first][iDim] > minPhysR[iDim]) {
                _material2min[iter->first][iDim] = minPhysR[iDim];
              }
              if (_material2max[iter->first][iDim] < maxPhysR[iDim]) {
                _material2max[iter->first][iDim] = maxPhysR[iDim];
              }
            }
          }
        }
      }
    }
 
    // store the number and min., max. possition for all ranks
#ifdef PARALLEL_MODE_MPI
    std::ptrdiff_t materials[_nMaterials];
    std::ptrdiff_t materialsInBuf[_nMaterials];
    std::ptrdiff_t materialCount[_nMaterials];
    std::ptrdiff_t materialCountInBuf[_nMaterials];
    T materialMinR[3*_nMaterials];
    T materialMaxR[3*_nMaterials];
    T materialMinRinBuf[3*_nMaterials];
    T materialMaxRinBuf[3*_nMaterials];
 
    for (int iM=0; iM<_nMaterials; iM++) {
      materials[iM]=-1;
      materialCount[iM]=0;
      for (int iDim=0; iDim<3; iDim++) {
        materialMinRinBuf[3*iM + iDim] = T();
        materialMaxRinBuf[3*iM + iDim] = T();
      }
    }
    std::size_t counter = 0;
    std::map<int, std::size_t>::iterator iMaterial;
    for (iMaterial = _material2n.begin(); iMaterial != _material2n.end(); iMaterial++) {
      materials[counter] = iMaterial->first;
      materialCount[counter] = iMaterial->second;
      for (int iDim=0; iDim<3; iDim++) {
        materialMinR[3*counter + iDim] = _material2min[iMaterial->first][iDim];
        materialMaxR[3*counter + iDim] = _material2max[iMaterial->first][iDim];
      }
      counter++;
    }
 
    for (int iRank=1; iRank<singleton::mpi().getSize(); iRank++) {
      int myRank = singleton::mpi().getRank();
      singleton::mpi().sendRecv(materials, materialsInBuf, _nMaterials,
                                (myRank+iRank)%singleton::mpi().getSize(),
                                (myRank-iRank+singleton::mpi().getSize())%singleton::mpi().getSize(), 0);
      singleton::mpi().sendRecv(materialCount, materialCountInBuf, _nMaterials,
                                (myRank+iRank)%singleton::mpi().getSize(),
                                (myRank-iRank+singleton::mpi().getSize())%singleton::mpi().getSize(), 1);
      singleton::mpi().sendRecv(materialMinR, materialMinRinBuf, 3*_nMaterials,
                                (myRank+iRank)%singleton::mpi().getSize(),
                                (myRank-iRank+singleton::mpi().getSize())%singleton::mpi().getSize(), 2);
      singleton::mpi().sendRecv(materialMaxR, materialMaxRinBuf, 3*_nMaterials,
                                (myRank+iRank)%singleton::mpi().getSize(),
                                (myRank-iRank+singleton::mpi().getSize())%singleton::mpi().getSize(), 3);
      for (int iM=0; iM<_nMaterials; iM++) {
        if (materialsInBuf[iM]!=-1) {
          std::vector<T> minPhysR(3,T());
          std::vector<T> maxPhysR(3,T());
          for (int iDim=0; iDim<3; iDim++) {
            minPhysR[iDim] = materialMinRinBuf[3*iM + iDim];
            maxPhysR[iDim] = materialMaxRinBuf[3*iM + iDim];
          }
          if (_material2n.count(materialsInBuf[iM]) == 0) {
            _material2n[materialsInBuf[iM]] = materialCountInBuf[iM];
            _material2min[materialsInBuf[iM]] = minPhysR;
            _material2max[materialsInBuf[iM]] = maxPhysR;
          }
          else {
            _material2n[materialsInBuf[iM]] += materialCountInBuf[iM];
            for (int iDim=0; iDim<3; iDim++) {
              if (_material2min[materialsInBuf[iM]][iDim] > minPhysR[iDim]) {
                _material2min[materialsInBuf[iM]][iDim] = minPhysR[iDim];
              }
              if (_material2max[materialsInBuf[iM]][iDim] < maxPhysR[iDim]) {
                _material2max[materialsInBuf[iM]][iDim] = maxPhysR[iDim];
              }
            }
          }
        }
      }
    }
#endif
 
    // update _minOverMaterial and _maxOverMaterial
    typename std::map< int, olb::Vector<T, 3> >::const_iterator iter;
    // find componentwise minimal extension over all material numbers
    for ( int iDim = 0; iDim < 3; iDim++ ) {
      // minimum
      for ( iter = _material2min.begin(); iter != _material2min.end(); iter++ ) {
        if ( iter->first != 0 ) { // only relevant material number are considered
          if ( _minOverMaterial[iDim] > iter->second[iDim] ) {
            _minOverMaterial[iDim] = iter->second[iDim];
          }
        }
      }
      // maximum
      for ( iter = _material2max.begin(); iter != _material2max.end(); iter++ ) {
        if ( iter->first != 0 ) { // only relevant material number are considered
          if ( _maxOverMaterial[iDim] < iter->second[iDim] ) {
            _maxOverMaterial[iDim] = iter->second[iDim];
          }
        }
      }
    }
 
    //clout.setMultiOutput(true);
//    print();
    //clout.setMultiOutput(false);
 
    if (verbose) {
      clout << "updated" << std::endl;
    }
    _statisticsUpdateNeeded = false;
  }
}

References olb::singleton::MpiManager::getRank(), olb::singleton::MpiManager::getSize(), olb::BlockGeometry< T, D >::getStatistics(), olb::singleton::mpi(), olb::singleton::MpiManager::reduceAndBcast(), and olb::singleton::MpiManager::sendRecv().

Here is the call graph for this function:

---

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

• src/geometry/superGeometryStatistics3D.h
• src/geometry/superGeometryStatistics3D.hh