smi.m

来自「iccg不完全乔尔斯基分解共轭梯度求解大规模稀疏方程组」· M 代码 · 共 20 行

function v=smi(ss,x,nd,nc)
nnd=6;
nnc=10;
for ii=1:nnd
    if(ii==1)
        mi=1;
    end
    if(ii>1)
        mi=nd(ii-1)+1;
    end
    ni=nd(ii);
    m1=ni-1;
    y(ii)=x(ii)*ss(ni);
    for k=mi:m1
        ns=nc(k);
        y(ii)=y(ii)+x(ns)*ss(k);
        y(ns)=y(ns)+x(ii)*ss(k);
    end
end
v=y;