On Wed, May 20, 2015 at 10:27 PM, Rongliang Chen <[email protected]> wrote:
> 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. ' > The normals are area weighted. I will indicate this in the documentation: https://bitbucket.org/petsc/petsc/src/5a737092961693f09ed5214c2133c4a7b76f72d7/src/dm/impls/plex/plexgeometry.c?at=master#cl-1208 Thanks, Matt Best, > Rongliang > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener
