On Mon, Oct 26, 2009 at 8:09 AM, Yujie <[email protected]> wrote: > 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?
As I understand it, this guarantee isn't required as the triangular boundary elements phi_face[] are different to the tetrahedral domain elements phi[]. ------------------------------------------------------------------------------ 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
