On Thu, 12 Jun 2008, Roy Stogner wrote: > On Thu, 12 Jun 2008, Derek Gaston wrote: > >> It's kind of hard to see how to go from dofs_on_side back to the >> mapping to Re()... you have to know what "node number" those dofs map >> to for the current element. Or am I missing something? > > IIRC dofs_on_side gives you local dof numbers; you still need the > result of a dof_indices call to look up the corresponding global dof > indices.
Oh, and to really be general this implies that your residual terms aren't just the pointwise difference between current value and desired value at node i, they're the sum of the side integral of the difference between current and desired values weighted by test function i. So basically you still evaluate those integrals, but you use the dofs_on_side information to decide which test functions to skip. --- Roy ------------------------------------------------------------------------- Check out the new SourceForge.net Marketplace. It's the best place to buy or sell services for just about anything Open Source. http://sourceforge.net/services/buy/index.php _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
