olb::opti::SimpleGeometrySerializer< S, dim > Class Template Reference

This class serializes the cells inside the geometry. More...

#include <serialization.h>

Inheritance diagram for olb::opti::SimpleGeometrySerializer< S, dim >:

Collaboration diagram for olb::opti::SimpleGeometrySerializer< S, dim >:

Public Member Functions
	SimpleGeometrySerializer (CuboidDecomposition< S, dim > &cGeometry)
	SimpleGeometrySerializer (SuperGeometry< S, dim > &sGeometry)
std::size_t	getSerializedCellIndex (const int latticeR[]) const override
	Compute serialized cell index from lattice coordinates.
LatticeR< dim+1 >	getLatticeR (std::size_t index) const override
	Get lattice coordinates from serialized cell index.
unsigned	getNoCells () const override
Public Member Functions inherited from olb::opti::GeometrySerializer< S, dim >
std::size_t	getSerializedCellIndex (LatticeR< dim+1 > latticeR) const
	Compute serialized cell index from lattice coordinates.
std::size_t	getSerializedComponentIndex (const int latticeR[], unsigned iD, unsigned fieldDim) const
	Get index of field component from lattice coordinates and component index.
std::size_t	getSerializedComponentIndex (LatticeR< dim+1 > latticeR, unsigned iD, unsigned fieldDim) const
	Get index of field component from lattice coordinates and component index.
std::size_t	getSerializedComponentIndex (LatticeAndFieldR coords, unsigned fieldDim) const
	Get index of field component from lattice coordinates and component index.
LatticeAndFieldR	getLatticeAndFieldR (std::size_t index, unsigned fieldDim) const
	Get lattice coordinates and field component from serialized field index.

Protected Attributes
CuboidDecomposition< S, dim > &	_cGeometry
const unsigned	_noCuboids

Detailed Description

template<typename S, unsigned dim>
class olb::opti::SimpleGeometrySerializer< S, dim >

This class serializes the cells inside the geometry.

Formula for serialized cell index (3d): cuboid_offset + NX*NY*z + NX*y + x Formula for serialized field component index: fieldDim * cellIndex + fieldComponent

This class is not intended to be used in performance-relevant context. Index computations could be accellerated if desired.

Definition at line 114 of file serialization.h.

Constructor & Destructor Documentation

SimpleGeometrySerializer() [1/2]

template<typename S , unsigned dim>

olb::opti::SimpleGeometrySerializer< S, dim >::SimpleGeometrySerializer ( CuboidDecomposition< S, dim > & cGeometry )

inline

Definition at line 126 of file serialization.h.

   : _cGeometry(cGeometry),
     _noCuboids(_cGeometry.size())
  {
    _cuboidSizes.reserve(_noCuboids);
    _offsets.reserve(_noCuboids);
    for (unsigned i = 0; i < _noCuboids; ++i) {
      _cuboidSizes.push_back(_cGeometry.get(i).getLatticeVolume());
    }
    _offsets.push_back(0);
    for (unsigned i = 1; i < _noCuboids; ++i) {
      _offsets.push_back(_offsets[i-1] + _cuboidSizes[i-1]);
    }
 
    const auto mc = _cGeometry.getMotherCuboid();
    if constexpr (dim == 2) {
      _noCells = (mc.getNx() + 1) * (mc.getNy() + 1);
    } else {
      _noCells = (mc.getNx() + 1) * (mc.getNy() + 1) * (mc.getNz() + 1);
    }
  }

References olb::opti::SimpleGeometrySerializer< S, dim >::_cGeometry, olb::opti::SimpleGeometrySerializer< S, dim >::_noCuboids, olb::CuboidDecomposition< T, D >::get(), olb::Cuboid< T, D >::getLatticeVolume(), and olb::CuboidDecomposition< T, D >::getMotherCuboid().

Here is the call graph for this function:

SimpleGeometrySerializer() [2/2]

template<typename S , unsigned dim>

olb::opti::SimpleGeometrySerializer< S, dim >::SimpleGeometrySerializer ( SuperGeometry< S, dim > & sGeometry )

inline

Definition at line 148 of file serialization.h.

   : SimpleGeometrySerializer(sGeometry.getCuboidDecomposition())
  { }

Member Function Documentation

getLatticeR()

template<typename S , unsigned dim>

LatticeR< dim+1 > olb::opti::SimpleGeometrySerializer< S, dim >::getLatticeR ( std::size_t index ) const

inlineoverridevirtual

Get lattice coordinates from serialized cell index.

Implements olb::opti::GeometrySerializer< S, dim >.

Definition at line 171 of file serialization.h.

                                                              {
    const auto cuboidIt = std::upper_bound(_offsets.begin(), _offsets.end(), index) - 1;
    LatticeR<dim+1> res;
    res[0] = std::distance(_offsets.begin(), cuboidIt);
 
    const Cuboid<S,dim>& c = _cGeometry.get(res[0]);
    const std::size_t nX = c.getNx();
    if constexpr (dim == 2) {
      // index = cuboid_offset + NX*y + x
      const auto divByNx = std::ldiv((long int) index - *cuboidIt, (long int) nX);
      res[2] = divByNx.quot;
      res[1] = divByNx.rem;
    } else {
      // index = cuboid_offset + NX*NY*z + NX*y + x
      const std::size_t nY = c.getNy();
      const auto divByNxNy = std::ldiv((long int) index - *cuboidIt, (long int) nX*nY);
      res[3] = divByNxNy.quot;
      const auto divByNx = std::ldiv(divByNxNy.rem, (long int) nX);
      res[2] = divByNx.quot;
      res[1] = divByNx.rem;
    }
    return res;
  }

References olb::opti::SimpleGeometrySerializer< S, dim >::_cGeometry, and olb::CuboidDecomposition< T, D >::get().

Here is the call graph for this function:

getNoCells()

template<typename S , unsigned dim>

unsigned olb::opti::SimpleGeometrySerializer< S, dim >::getNoCells ( ) const

inlineoverridevirtual

Implements olb::opti::GeometrySerializer< S, dim >.

Definition at line 195 of file serialization.h.

                                       {
    return _noCells;
  }

getSerializedCellIndex()

template<typename S , unsigned dim>

std::size_t olb::opti::SimpleGeometrySerializer< S, dim >::getSerializedCellIndex ( const int latticeR[] ) const

inlineoverridevirtual

Compute serialized cell index from lattice coordinates.

Implements olb::opti::GeometrySerializer< S, dim >.

Definition at line 153 of file serialization.h.

                                                                        {
    const Cuboid<S,dim>& c = _cGeometry.get(latticeR[0]);
    const int nX = c.getNx();
    const int nY = c.getNy();
    std::size_t res;
    if constexpr (dim == 2) {
      // cuboid_offset + NX*y + x
      res = _offsets[latticeR[0]] + nX*latticeR[2] + latticeR[1];
    } else {
      // cuboid_offset + NX*NY*z + NX*y + x
      res = _offsets[latticeR[0]] + nX*nY*latticeR[3] + nX*latticeR[2] + latticeR[1];
    }
    return res;
  }

References olb::opti::SimpleGeometrySerializer< S, dim >::_cGeometry, and olb::CuboidDecomposition< T, D >::get().

Here is the call graph for this function:

Member Data Documentation

_cGeometry

template<typename S , unsigned dim>

CuboidDecomposition<S,dim>& olb::opti::SimpleGeometrySerializer< S, dim >::_cGeometry

protected

Definition at line 117 of file serialization.h.

_noCuboids

template<typename S , unsigned dim>

const unsigned olb::opti::SimpleGeometrySerializer< S, dim >::_noCuboids

protected

Definition at line 118 of file serialization.h.

---

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

• src/optimization/solver/serialization.h