Dear LibMesh Developers,
Regarding the third example (http://libmesh.sourceforge.net/ex3.php), my
question is as follows:
in assemble_poission() function,
for each element, we have
"
for (unsigned int i=0; i<phi.size(); i++)
Fe(i) += JxW[qp]*fxy*phi[i][qp];
"
for each boundary side, we have correspondingly
"
for (unsigned int i=0; i<phi_face.size(); i++)
Fe(i) += JxW_face[qp]*penalty*value*phi_face[i][qp];
"
Assuming that we use tetrahedral element and linear Lagrange shape function,
in this case, for an element, there are 4 vertices, that is i=0, 1, 2, 3.
for a boundary side, how to guarantee the vertices (i=0, 1, 2) on the
boundary? Thanks a lot.
Regards,
Yujie
------------------------------------------------------------------------------
Come build with us! The BlackBerry(R) Developer Conference in SF, CA
is the only developer event you need to attend this year. Jumpstart your
developing skills, take BlackBerry mobile applications to market and stay
ahead of the curve. Join us from November 9 - 12, 2009. Register now!
http://p.sf.net/sfu/devconference
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
