On Thu, 12 Jun 2008, Derek Gaston wrote: > Of course putting all my thoughts together in an email made me realize > what I needed to do... all I need to do is use the usual side > assembly... and call elem->is_node_on_side(i,side)... and if it is > then that dof is on the side... then I can just set > Re(i)=soln(dof_indices[i])-bc_value like I want to... overriding > whatever values were assembled into Re before I ever add it into the > residual...
This looks good, but you might find it easier in the long run to use FE::dofs_on_side() rather than Elem::is_node_on_side; the latter might make it too easy to write code that only works on isoparametric Lagrange elements. --- 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
