xyst

Skip applying BCs in CG solver if not needed

This set of changes allows the conjugate gradients solver to be
initialized (before a new solve) without applying the boundary
conditions on the matrix and rhs. In ChoCG this is only necessary during
startup and in the first time step, but not afterwards if the mesh does
not change.