def qmrs(*args, **kwargs): return krylov.qmrs(*args, **kwargs)
x = np.zeros(n, 'd') info, iter, relres = minres(A, b, x, 1e-6, 1000) print 'minres, K_none: ', info, iter, relres, resid(A, b, x) x = np.zeros(n, 'd') info, iter, relres = minres(A, b, x, 1e-6, 1000, K_diag) print 'minres, K_diag: ', info, iter, relres, resid(A, b, x) x = np.zeros(n, 'd') info, iter, relres = minres(A, b, x, 1e-6, 1000, K_jac) print 'minres, K_jac: ', info, iter, relres, resid(A, b, x) x = np.zeros(n, 'd') info, iter, relres = minres(A, b, x, 1e-6, 1000, K_ssor) print 'minres, K_ssor: ', info, iter, relres, resid(A, b, x) x = np.zeros(n, 'd') info, iter, relres = qmrs(A, b, x, 1e-6, 1000) print 'qmrs, K_none: ', info, iter, relres, resid(A, b, x) x = np.zeros(n, 'd') info, iter, relres = qmrs(A, b, x, 1e-6, 1000, K_diag) print 'qmrs, K_diag: ', info, iter, relres, resid(A, b, x) x = np.zeros(n, 'd') info, iter, relres = qmrs(A, b, x, 1e-6, 1000, K_jac) print 'qmrs, K_jac: ', info, iter, relres, resid(A, b, x) x = np.zeros(n, 'd') info, iter, relres = qmrs(A, b, x, 1e-6, 1000, K_ssor) print 'qmrs, K_ssor: ', info, iter, relres, resid(A, b, x) x = np.zeros(n, 'd') info, iter, relres = cgs(A, b, x, 1e-6, 1000) print 'cgs, K_none: ', info, iter, relres, resid(A, b, x) x = np.zeros(n, 'd')