Assemble form(s) and apply any given boundary conditions in a symmetric fashion and return tensor(s).
The standard application of boundary conditions does not necessarily preserve the symmetry of the assembled matrix. In order to perserve symmetry in a system of equations with boundary conditions, one may use the alternative assemble_system instead of multiple calls to assemble.
Examples of usage
For instance, the statements
A = assemble(a) b = assemble(L) bc.apply(A, b)can alternatively be carried out by
A, b = assemble_system(a, L, bc)The statement above is valid even if bc is a list of DirichletBC instances. For more info and options, see assemble.