A gcc fix in LinearizerBase.
Implement ParalutionLinearMatrixSolver<Scalar>::get_residual().
Make PARALUTION compliant with other solvers in terms of initial_guess usage.