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

[Dharmendar Reddy]

So, i need to transform the Gauss quadrature points (qp) from the reference
2D triangle to 3D using affine transformation of the form:
[x y z]^T = A x [qp_x qp_y] + [c1 c2 c3]^T; A is a 3 x 2 matrix , am i
right?

thanks
Reddy

On Sat, Nov 5, 2011 at 3:38 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:

> 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.
>



-- 
-----------------------------------------------------
Dharmendar Reddy Palle
Graduate Student
Microelectronics Research center,
University of Texas at Austin,
10100 Burnet Road, Bldg. 160
MER 2.608F, TX 78758-4445
e-mail: dharmareddy84 at gmail.com
Phone: +1-512-350-9082
United States of America.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111105/93079420/attachment.htm>