Initial work on BCs
Much more BC
Down to the Jacobian
Down to the linear solve
complete implementation section