main.m

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

clc
clear all
nnd=6;
nnc=10;
aa=[1.1+i, 2.5+3i, 2.2, 4.0, 7.6, 1.3 ,9.2,3.4, 5.3 ,1.6];
 t=[11+i, 12.3+2i, 20.7,8.9,12.6,6.9];
 nd=[1, 2 ,4, 6,8, 10];
 nc=[1, 2, 2 ,3 , 2 ,4 ,3 ,5 ,3 ,6];
for ii=1:nnc
    ss(ii)=aa(ii);
end
sillt(nnd,aa,nd,nc);
bb=iccg(aa,ss,t,nd,nc);