Gauss-Seidel is done ( I think - need testing )!  Solver's constructor needs to construct the B-Vector to get the whole eq.