Working!! (non mpi CG)

CG solver needs to set boundary conditions within the algorithm.
Currently we pass the Funken IndexVector of boundary nodes (100%
transparent via the template), since converting to Flens vector would
be too ugly.