my_qmr_cmplxsym.m

来自「求解线性系统的Krylov方法的工具箱」· M 代码 · 共 9 行

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 =