dlu.m

来自「gmres.m为迭代算法GMRES的实现」· M 代码 · 共 15 行

function [D,L,U]=dlu(A) %dlu分解子函数 
D=zeros(size(A)); 
L=D; 
U=D; 

for i=1:size(A,1) 
    D(i,i)=A(i,i); 
end 
for i=1:size(A,1) 
    for j=1:size(A,2) 
        L(i,j)=-(i <j)*A(i,j); 
        U(i,j)=-(i>j)*A(i,j); 
    end 
end;