olb::AnalyticalPorosityVolumeF< T > Class Template Reference

#include <analyticalPorosityVolumeF.h>

Inheritance diagram for olb::AnalyticalPorosityVolumeF< T >:

Collaboration diagram for olb::AnalyticalPorosityVolumeF< T >:

Public Member Functions
	AnalyticalPorosityVolumeF (std::string fileName, T spacing)
Vector< T, 3 >	getPhysShape () const
bool	operator() (T output[], const T physR[]) override
Public Member Functions inherited from olb::AnalyticalF< 3, T, T >
AnalyticalF< D, T, T > &	operator- (AnalyticalF< D, T, T > &rhs)
AnalyticalF< D, T, T > &	operator+ (AnalyticalF< D, T, T > &rhs)
AnalyticalF< D, T, T > &	operator* (AnalyticalF< D, T, T > &rhs)
AnalyticalF< D, T, T > &	operator/ (AnalyticalF< D, T, T > &rhs)
Public Member Functions inherited from olb::GenericF< T, S >
virtual	~GenericF ()=default
int	getSourceDim () const
	read only access to member variable _m
int	getTargetDim () const
	read only access to member variable _n
std::string &	getName ()
	read and write access to name
std::string const &	getName () const
	read only access to name
virtual bool	operator() (T output[], const S input[])=0
	has to be implemented for 'every' derived class
bool	operator() (T output[])
	wrapper that call the pure virtual operator() (T output[], const S input[]) from above
bool	operator() (T output[], S input0)
bool	operator() (T output[], S input0, S input1)
bool	operator() (T output[], S input0, S input1, S input2)
bool	operator() (T output[], S input0, S input1, S input2, S input3)

Additional Inherited Members
Public Types inherited from olb::AnalyticalF< 3, T, T >
using	identity_functor_type
Public Types inherited from olb::GenericF< T, S >
using	targetType = T
using	sourceType = S
Public Attributes inherited from olb::GenericF< T, S >
std::shared_ptr< GenericF< T, S > >	_ptrCalcC
	memory management, frees resouces (calcClass)
Static Public Attributes inherited from olb::AnalyticalF< 3, T, T >
static constexpr unsigned	dim
Protected Member Functions inherited from olb::AnalyticalF< 3, T, T >
	AnalyticalF (int n)
Protected Member Functions inherited from olb::GenericF< T, S >
	GenericF (int targetDim, int sourceDim)

Detailed Description

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

Definition at line 42 of file analyticalPorosityVolumeF.h.

Constructor & Destructor Documentation

AnalyticalPorosityVolumeF()

template<typename T >

olb::AnalyticalPorosityVolumeF< T >::AnalyticalPorosityVolumeF ( std::string fileName, T spacing )

inline

Definition at line 56 of file analyticalPorosityVolumeF.h.

      : AnalyticalF<3,T,T>(1), _spacing{spacing}
  {
    OstreamManager clout(std::cout, "PorosityVolumeImporter");
    _fileName = fileName;
    #ifdef FEATURE_VDB
    if (isVDBFile(_fileName)) {
      openvdb::initialize();
      openvdb::io::File _vdbFile(_fileName);
      // Read the OpenVDB file
      _vdbFile.open();
      grids = _vdbFile.getGrids();
      auto grid = openvdb::gridPtrCast<openvdb::FloatGrid>(*grids->begin());
      _vdbFile.close();
      openvdb::CoordBBox bounds = grid->evalActiveVoxelBoundingBox();
      _shape[0] = bounds.max().x() - bounds.min().x();
      _shape[1] = bounds.max().y() - bounds.min().y();
      _shape[2] = bounds.max().z() - bounds.min().z();
    }
    #endif
    #if not defined(FEATURE_VTK)
    if (isVTKFile(_fileName))
    {
      std::cerr << "To use the VTK format, add VTK to the FEATURES list in config.mk";
      exit(1);
    }
    #endif
    #if defined(FEATURE_VTK)
    if (isVTKFile(_fileName)) {
      _reader->SetFileName(fileName.c_str());
      _reader->Update();
      _reader->SetScalarsName(_reader->GetScalarsNameInFile(0));
      auto* _data = _reader->GetOutput();
      _data->GetDimensions(_shape.data());
    }
    #endif
    #ifndef FEATURE_VDB
    if (isVDBFile(_fileName))
    {
      std::cerr << "To use the VDB format, add VDB to the FEATURES list in config.mk";
      exit(1);
    }
    #endif
  }

References olb::Vector< T, D >::data().

Here is the call graph for this function:

Member Function Documentation

getPhysShape()

template<typename T >

Vector< T, 3 > olb::AnalyticalPorosityVolumeF< T >::getPhysShape ( ) const

inline

Definition at line 101 of file analyticalPorosityVolumeF.h.

                                   {
    return _shape * _spacing;
  }

operator()()

template<typename T >

bool olb::AnalyticalPorosityVolumeF< T >::operator() ( T output[], const T physR[] )

inlineoverride

Definition at line 105 of file analyticalPorosityVolumeF.h.

                                                       {
    #ifdef FEATURE_VDB
    if (isVDBFile(_fileName)) {
      int iX;
      int iY;
      int iZ;
      iX = util::floor(physR[0] / _spacing);
      iY = util::floor(physR[1] / _spacing);
      iZ = util::floor(physR[2] / _spacing);
      if (iX >= 0 && iY >= 0 && iZ >= 0 && iX < _shape[0] && iY < _shape[1] && iZ < _shape[2]) {
        auto grid = openvdb::gridPtrCast<openvdb::FloatGrid>(*grids->begin());
        openvdb::Coord location(iX, iY, iZ);
        openvdb::FloatGrid::Accessor accessor = grid->getAccessor();
        output[0] =accessor.getValue(location);
      } else {
        output[0]=0;
      }
      return true;
    }
    #endif
    #if defined(FEATURE_VTK)
    if(isVTKFile(_fileName)) {
      int iX;
      int iY;
      int iZ;
      iX = util::floor(physR[0] / _spacing);
      iY = util::floor(physR[1] / _spacing);
      iZ = util::floor(physR[2] / _spacing);
      if (iX >= 0 && iY >= 0 && iZ >= 0 && iX < _shape[0] && iY < _shape[1] && iZ < _shape[2]) {
        auto* _data = _reader->GetOutput();
        output[0] = *static_cast<T*>(_data->GetScalarPointer(iX, iY, iZ));
      } else {
        output[0] = 0;
      }
      return true;
    }
    #endif
    return false;
  }

References olb::util::floor().

Here is the call graph for this function:

---

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

• src/io/analyticalPorosityVolumeF.h