[petsc-users] [FEM doubt] Not Related to PETSc

[Jed Brown]

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>