结构矩阵分析家庭作业：编程项目与LDLT分解优化.pdfVIP

家庭作业9

1.根据给定的模板完成编程项目。2.第119页上的LDLT分解代码可以改进。

你能找到如何改进吗？（提示：内层do循环中的一些操作可以移出循环。）

diag(1:ncol)=(/(Kcol(j)%row(j),j=1,ncol)/)

doj=2,ncol

row1j=lbound(Kcol(j)%row,1)

doi=row1j,j-1

row_1=max(row1j,lbound(Kcol(i)%row,1))

s=sum(diag(row_1:i-1)*Kcol(i)%row(row_1:i-1)*Kcol(j)%row(row_1:i-1))

Kcol(j)%row(i)=Kcol(j)%row(i)-s

enddo

Kcol(j)%row(row1j:j-1)=Kcol(j)%row(row1j:j-1)/diag(row1j:j-1)

ssum(diag(row1j:j-1)*Kcol(j)%row(row1j:j-1)**2)

diag(j)=diag(j)-s

enddo

把除法提取出来，再进行逐元操作，不知道是否可以使算法加速？