Add documentation to boundary conditions Correct DirichletBC::fillBC to use column basis for solution