← Back to team overview

dolfin team mailing list archive

Preserving Symmetry in applying boundaryconditions

 

Hello,

I have a simple question. I want to apply Dirichlet boundary
conditions and mantain the symmetry of the matrix. Previously I was
using "assemble_system" which does this. Now I have a situation, in
which I need to assemble a symmetric matrix A separately from vector
b. If I use DirichletBC::apply(A,b), I loose the symmetry in A. Can
someone tell me what to use in this case?

Thanks,
-Jehanzeb


Follow ups