Ah. I see. So the extra guys are just set to 0. Thank you very much! Andrea
On Tue, Mar 10, 2009 at 1:02 PM, Roy Stogner <royst...@ices.utexas.edu>wrote: > > On Tue, 10 Mar 2009, Andrea Hawkins wrote: > > 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. >> > > They do correspond. > > Is there something that happens in fe_face->reinit() or something? >> > > What happens is that phi_face includes all element shape functions, > even when they're 0 on the entire face boundary. I guess the choice > was between adding a bunch of zeros vs. indexing through an extra > indirection layer, and libMesh went with simpler code over > negligably better assembly performance. > --- > Roy > ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
