Inconsistency in writing the residual norms between solvers

[codecircuit]

• Set the solver verbosity to 1 (solver_par->verbose = 1)
• Set number of iterations to N (solver_par->maxiter = N)
• execute solver PBICG, which writes N + 1 residual norms to solver_par->res_vec from 0 to N inclusively
• execute solver PGMRES, which writes only N residual norms to solver_par->res_vec from 1 to N inclusively, where solver_par->res_vec[1] = 0 always.

‌