[firedrake] Boundary conditions

[Lawrence Mitchell]

> On 2 Feb 2014, at 09:04, "Cotter, Colin J" <colin.cotter at imperial.ac.uk> wrote:
> 
> OK, next question is whether I can output a numpy matrix for a bilinear form assembled with boundary conditions, i.e. whether I can get the reduced matrix with the boundary rows and columns removed?

We don't give you the matrix with the bc rows and cols removed, we just zero them and slap a 1 on the diagonal, you'll have to do it yourself. Note that you can get the dof indices of the bc nodes by accessing bc.nodes where bc is a DirichletBC object.

Lawrence