On Sat, Nov 5, 2011 at 14:34, Dharmendar Reddy <dharmareddy84 at gmail.com>wrote:
> I am writing a 3D FEM code for learning purpose (experimenting with object > oriented concepts in Fortran 2003). Can some one tell me (pseudo code) how > to implement a non homogenous Neumann boundary condition. You can also > point me to a book. I am using tetrahedral elements for solving Poisson > equation. I am using a 4 point Gauss quadrature, and linear basis functions > for cell volume integrals. I am confused on how to do the integration on > the facet. You need to compute an integral over the face. It appears in the weak form. Any book on finite element methods should cover this. -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111105/6fcf3811/attachment.htm>
