my_qmr_cmplxsym.m

求解线性系统的Krylov方法的工具箱

function [x,iter,flg,res] = my_qmr_cmplxsym(A,b,tol,maxit)
n = length(A);
x = zeros(n,1);
v0 = zeros(n,1);
p0 = zeros(n,1);
p_1 = zeros(n,1);
v1 = b-A*x;
beta = sqrt(v1.'*v1);
tau =