From: Roy S. <ro...@st...> - 2008-06-12 20:08:34
|
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 |