Hi there, I have a question about the finite volume formulation. In function TSComputeRHSFunction_DMPlex, when we compute the flux (line 617), why we do not multiply the face area? In the formulation of the face flux, I think we need to multiply the face area and divide the element volume.
Best, Rongliang
