Reply To: Compabilitities of superLattticeIntegralF2D for DynOmegaD2Q9Descriptor
OpenLB – Open Source Lattice Boltzmann Code › Forums › on OpenLB › General Topics › Compabilitities of superLattticeIntegralF2D for DynOmegaD2Q9Descriptor › Reply To: Compabilitities of superLattticeIntegralF2D for DynOmegaD2Q9Descriptor
September 13, 2019 at 10:52 am
#4541
jonathan
Participant
An even easier way of achieving a static solution would be to just specify the field type inside the physFactor calculation.
_physFactor = -this->_blockLattice.get(input[0]+_overlap,input[1]+_overlap).template getFieldPointer<descriptors::OMEGA>()[0] * descriptors::invCs2<T,DESCRIPTOR>() / dt * this->_converter.getPhysDensity() * this->_converter.getPhysViscosity();
The method described in my first answer is more generic and necessary if you want to have easily swappable descriptors.
Best,
Jonathan