C++ code to solve linear systems with GMRES solver