dolfin team mailing list archive
-
dolfin team
-
Mailing list archive
-
Message #15758
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