Yes, that's necessary; see the last bit of the doc-string for assemble(): If ``bcs`` is supplied and ``f`` is a 2-form, the rows and columns of the resulting :class:`.Matrix` corresponding to boundary nodes will be set to 0 and the diagonal entries to 1. If ``f`` is a 1-form, the vector entries at boundary nodes are set to the boundary condition values. Cheers Stephan On 31/12/15 20:01, Justin Chang wrote:
Hi all,
Say I have this:
a = dot(grad(u),grad(v))*dx L = f*v*dx bcs = DirichletBC(...)
A = assemble(a) b = assemble(L) solver = LinearSolver(A, ...) solver.solve(x,b)
If I want to apply the DirichletBC's, I would use "A = assemble(a,bcs=bcs)." Is it necessary to also do this for b? Would it have any affect?
Thanks, Justin
_______________________________________________ firedrake mailing list firedrake@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/firedrake