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 _______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
