sparse_cmp.m

来自「改进平方根法求解对称正定线性方程组」· M 代码 · 共 22 行

function [Acmp,B]=sparse_cmp(A,n)
B=1;          %初始化半带宽
for i=1:n
    for j=n:-1:1
        if A(i,j)~=0
            temp=j;      % 每行从尾往头遍历查找非零元素确定每行带宽
            break;
        end
    end
    if temp-i+1>B
        B=temp-i+1;      % 判断是否为最大半带宽
    end
end
Acmp=zeros(n,B);         % 根据既得带宽定义压缩存储矩阵的维数
for i=1:n
    for j=1:B
        if i+j-1<=n
            Acmp(i,j)=A(i,i+j-1);
        end
    end
end
return