Right... But if you think about a 2D element with local dofs as numbered here
3 ------------ 2 | | | | 0 ------------- 1 If the side between dof 2 and 3 is on the boundary, in that loop I don't think you should be adding to Ke(0,0), right? But you will be. Or have I completely over thought this? On Tue, Mar 10, 2009 at 12:54 PM, John Peterson <jwpeter...@gmail.com>wrote: > On Tue, Mar 10, 2009 at 11:51 AM, Andrea Hawkins <andjhawk...@gmail.com> > wrote: > > Hello- > > > > In looking at the examples that include boundary integration, I just had > a > > question of how something is working. > > > > For discussion I will consider example 5. The loops for test and trial > > functions run over i,j=0 to phi_face.size() which makes sense, but I'm > > wondering about the value being added to Ke(i,j) and Fe(i) as I wouldn't > > think the element dof corresponds to the side-trial function number. Is > > there something that happens in fe_face->reinit() or something? > > All the same basis functions and numbering are used, the quadrature > rule is changed (I think!) > > -- > John > ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
