<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, May 20, 2015 at 10:27 PM, Rongliang Chen <span dir="ltr"><<a href="mailto:rongliang.chan@gmail.com" target="_blank">rongliang.chan@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div><div><div>Hi there,<br><br></div>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. '<br></div></div></div></blockquote><div><br></div><div>The normals are area weighted. I will indicate this in the documentation:</div><div><br></div><div><a href="https://bitbucket.org/petsc/petsc/src/5a737092961693f09ed5214c2133c4a7b76f72d7/src/dm/impls/plex/plexgeometry.c?at=master#cl-1208">https://bitbucket.org/petsc/petsc/src/5a737092961693f09ed5214c2133c4a7b76f72d7/src/dm/impls/plex/plexgeometry.c?at=master#cl-1208</a><br></div><div> </div><div>  Thanks,</div><div><br></div><div>     Matt</div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div><div></div>Best,<br></div>Rongliang<br></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>