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

Jed Brown jedbrown at mcs.anl.gov
Sat Nov 5 15:38:02 CDT 2011


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>


More information about the petsc-users mailing list