Re: [Libmesh-users] Stiffness matrix in Adaptive Mesh Refinement

[Roy S. <roy...@ic...>]

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