Reply To: Modifying field values for neighboring cells inside a postprocessor

[Danial.Khazaeipoul]

I did though about storing the cell IDs, instead of storing pointers to the cell interfaces. However, I would need to create a Cell interface to be able to modify a global data which requires a blockLattice to be passed as an argument:

Cell<T,DESCRIPTOR> cell(blockLattice, iCell);

From the user-guide, I can see that the OperatorScope can be set to PerBlock instead of PerCell or PerCellWithParameters. However, I am not sure if this works with the freeSurfacePostProcessor or if OpenLB provides a similar method for block lattice operator such as:

template <typename T, typename DESCRIPTOR>
template <typename BLOCK>
void TestPostProcessor3D<T, DESCRIPTOR>::apply(BLOCK& block)
{}