Detailed Description

Classes to be used together with the PDFReconstruction class to reconstruct PDFs.

Each reconstructor must exactly implement the member function void operator()( const cell_idx_t & x, const cell_idx_t & y, const cell_idx_t & z, IBlock * const block, PdfField_T * const pdfField ); that reconstructs all PDFs in a specific cell with indices x,y,z on a given block.

Different variants are available:

EquilibriumReconstructor: Determines a local average density and sets the PDFs to their equilibrium values based on this density and the body's velocity
EquilibriumAndNonEquilibriumReconstructor: First reconstructs the equilibrium values with the help of the EquilibriumReconstructor. Then extrapolates the non-equilibrium part from a neighboring cell and adds it to the reconstructed PDFs for better accuracy. Defaults to EquilibriumReconstructor if no suitable cell for extrapolation is available.
ExtrapolationReconstructor: Extrapolates the PDFs of three or two neighboring cells that lie in extrapolation direction. Optionally, a special treatment after the extrapolation can be used that sets certain moments directly (only D3Q19). Defaults to EquilibriumAndNonEquilibriumReconstructor if not enough cells in this direction are available.

EquilibriumAndNonEquilibriumReconstructor and ExtrapolationReconstructor need an extrapolation direction which is provided by the ExtrapolationDirectionFinder that is to be chosen (see ExtrapolationDirectionFinder.h).

#include <Reconstructor.h>

Public Types
using	PdfField_T = lbm::PdfField< LatticeModel_T >

using	BodyField_T = Field< pe::BodyID, 1 >

Public Member Functions
	EquilibriumReconstructor (const shared_ptr< StructuredBlockStorage > &blockStorage, const BlockDataID &boundaryHandlingID, const BlockDataID &bodyFieldID)

void	operator() (const cell_idx_t &x, const cell_idx_t &y, const cell_idx_t &z, IBlock const block, PdfField_T const pdfField)

Private Member Functions
void	setLocalEquilibrium (const cell_idx_t &x, const cell_idx_t &y, const cell_idx_t &z, IBlock const block, PdfField_T const pdfField)

real_t	getLocalAverageDensity (const cell_idx_t &x, const cell_idx_t &y, const cell_idx_t &z, IBlock const block, PdfField_T const pdfField) const

Private Attributes
shared_ptr< StructuredBlockStorage >	blockStorage_

const BlockDataID	boundaryHandlingID_

const BlockDataID	bodyFieldID_