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

This class serializes the cells which are marked by indicator. More...

#include <serialization.h>

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

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

Public Member Functions
	SparseGeometrySerializer (SuperGeometry< S, dim > &superGeometry, FunctorPtr< SuperIndicatorF< S, dim > > &&indicator)
	SparseGeometrySerializer (SuperGeometry< S, dim > &superGeometry, FunctorPtr< IndicatorF< S, dim > > &&indicator)
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.

Detailed Description

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

This class serializes the cells which are marked by indicator.

It produces a smaller control vector, but access is more expensive compared to SimpleGeometrySerializer.

Template Parameters

S	Floating point type
dim	Geometrical dimension

Definition at line 209 of file serialization.h.

Constructor & Destructor Documentation

SparseGeometrySerializer() [1/2]

template<typename S , unsigned dim>

olb::opti::SparseGeometrySerializer< S, dim >::SparseGeometrySerializer ( SuperGeometry< S, dim > & superGeometry, FunctorPtr< SuperIndicatorF< S, dim > > && indicator )

inline

Definition at line 220 of file serialization.h.

   : _indicator(std::move(indicator))
  {
    const auto cGeometry = superGeometry.getCuboidDecomposition();
    L latticeR;
    for (unsigned iC = 0; iC < cGeometry.size(); ++iC) {
      latticeR[0] = iC;
      const int nX = cGeometry.get(iC).getNx();
      const int nY = cGeometry.get(iC).getNy();
      for (int iX=0; iX<nX; iX++) {
        latticeR[1] = iX;
        for (int iY=0; iY<nY; iY++) {
          latticeR[2] = iY;
          if constexpr (dim == 2) {
            if (_indicator(latticeR.data())) {
              _coords.push_back(latticeR);
            }
          } else {
            const int nZ = cGeometry.get(iC).getNz();
            for (int iZ=0; iZ<nZ; iZ++) {
              latticeR[3] = iZ;
 
              if (_indicator(latticeR.data())) {
                _coords.push_back(latticeR);
              }
            }
          }
        }
      }
    }
  }

References olb::Vector< T, Size >::data(), and olb::SuperStructure< T, D >::getCuboidDecomposition().

Here is the call graph for this function:

SparseGeometrySerializer() [2/2]

template<typename S , unsigned dim>

olb::opti::SparseGeometrySerializer< S, dim >::SparseGeometrySerializer ( SuperGeometry< S, dim > & superGeometry, FunctorPtr< IndicatorF< S, dim > > && indicator )

inline

Definition at line 253 of file serialization.h.

   : SparseGeometrySerializer(superGeometry,
     new SuperIndicatorFfromIndicatorF<S,dim>(std::move(indicator), superGeometry))
  { }

Member Function Documentation

getLatticeR()

template<typename S , unsigned dim>

LatticeR< dim+1 > olb::opti::SparseGeometrySerializer< 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 283 of file serialization.h.

                                                              {
#ifdef OLB_DEBUG
    // check if index lies in accessible range
    if (index >= _coords.size()) {
      OstreamManager clout (std::cout, "SparseGeometrySerializer");
      clout << "Warning: the passed index is higher than the number of cells "
        << "which are marked by indicator." << std::endl;
      std::exit(1);
    }
#endif
    return _coords[index];
  }

getNoCells()

template<typename S , unsigned dim>

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

inlineoverridevirtual

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

Definition at line 296 of file serialization.h.

                                       {
    return _coords.size();
  }

getSerializedCellIndex()

template<typename S , unsigned dim>

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

inlineoverridevirtual

Compute serialized cell index from lattice coordinates.

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

Definition at line 260 of file serialization.h.

                                                                        {
#ifdef OLB_DEBUG
    // check if cell lies in indicated domain
    bool isInside;
    _indicator(&isInside, latticeR);
    if (! isInside) {
      OstreamManager clout (std::cout, "SparseGeometrySerializer");
      clout << "Warning: the passed cell does not lie in the indicated domain "
        << "and can therefore not be accessed." << std::endl;
      std::exit(1);
    }
#endif
    const L point (latticeR);
    const auto upper = std::upper_bound(
      _coords.begin(), _coords.end(), point, [](const L& a, const L& b){
      return lex_smaller(a, b);
    });
    return std::distance(_coords.begin(), upper - 1);
  }

---

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

• src/optimization/solver/serialization.h