On Wed, 31 Oct 2012, Kirk, Benjamin (JSC-EG311) wrote: > dof_map.heterogenously_constrain_element_matrix_and_vector (Ke, Fe, > dof_indices_disp); > > But for the life of me, they are *not* in the output of the linear solver. > The output of the following code
Sure they are, you're just erasing them. When I originally wrote the heterogenous constraints code, I didn't realize that DofMap::extract_local_vector was recalculating constrained values, so I didn't update it to account for the fact that constraints can now have right hand sides. Sorry about that; this probably counts as a bug worth fixing ASAP. --- Roy ------------------------------------------------------------------------------ Everyone hates slow websites. So do we. Make your web apps faster with AppDynamics Download AppDynamics Lite for free today: http://p.sf.net/sfu/appdyn_sfd2d_oct _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
