我正在进行doolittle分解。我的功能很好(没有错误),但我没有收到正确的结果。我确信函数中的一切都很好,除了2-4个变量,我可能会放错位置。所以我怀疑这里是否有一个函数或命令来打印这个函数正在使用的所有计算,一步一步,所以我可以理解出了什么问题。
我正在使用Python 3.6.2,Geany,Windows 10。
以下是代码:
def doolittle(A):
n = len(A)
b = 0
c = 0
L = [[0.0] * n for s in range(n)]
U = [[0.0] * n for s in range(n)]
for k in range(n):
L[k][k] = 1.0
for j in range(n):
b += L[k][j] * U[j][k]
U[k][j] = A[k][j] - b
for i in range(k+1,n):
c += L[k][i] * U[i][k]
L[k][j] = (A[k][j] - c)/U[k][k]
return(L,U)
A = [[1, -3, 2], [-2,8,-1], [4, -6, 5]]
print(doolittle(A))
答案 0 :(得分:2)
import ipdb
ipdb.set_trace()
尝试在想要输入功能的地方使用。