正在学习Cholesky分解
需要return一个上三角矩阵
因为是第一次接触python,感觉自己的代码写的很复杂,想问问大家该怎么变得更简洁
希望得到类似于这个帖子的第一个代码
https://blog.csdn.net/ReDreamme/article/details/108745255?spm=1001.2014.3001.5501
以下是我自己的代码
def Cholesky(A):
len=A.shape[0]
R=np.zeros((len,len))
for i in range(0,len):
for j in range(i, len):
temp=A[i,j]
k=i-1
while(k>=0):
temp-=R[k,j]*R[k,i]
k=k-1
if(i==j):
R[i,j]=np.sqrt(temp)
else:
R[i,j]=temp/R[i,i]
return R