On Thu, 4 Sep 2014, Miguel Angel Salazar de Troya wrote: > I'm using a FEMSystem with a NewtonSolver. I realized that in the FEMSystem, > in the assembly routines, DofMap::constrain_element_matrix_and_vector > is called with asymmetric_constraint_rows = false. Also, because it's a > nonlinear solver. DofMap::enforce_constraints_exactly(system) is called as > well. Would this ensure me a symmetric matrix even for the refinement case?
Assuming your element matrices are symmetric before refinement, this should keep them symmetric after refinement, yes. --- Roy ------------------------------------------------------------------------------ Slashdot TV. Video for Nerds. Stuff that matters. http://tv.slashdot.org/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
