Hi all,
    In case of no hanging nodes contraint, and the boundary condition is
only Dirichlet, should the asymmetric_constraint_rows always be set to true
in the function of
DofMap::heterogenously_constraint_element_matrix_and_vector?  Thanks.

I am looking at the libmesh examples introduction_ex4. Changing
asymmetric_constraint_rows from true to false will give wrong solution.
 dof_map.heterogenously_constrain_element_matrix_and_vector (Ke, Fe,
dof_indices,true);

void DofMap::heterogenously_constrain_element_matrix_and_vector
  (DenseMatrix<Number>& matrix,
   DenseVector<Number>& rhs,
   std::vector<dof_id_type>& elem_dofs,
   bool asymmetric_constraint_rows)
------------------------------------------------------------------------------
Minimize network downtime and maximize team effectiveness.
Reduce network management and security costs.Learn how to hire 
the most talented Cisco Certified professionals. Visit the 
Employer Resources Portal
http://www.cisco.com/web/learning/employer_resources/index.html
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to