Skip to content

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

#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