我在做杜利特尔分解。我的函数正常(没有错误),但是我没有收到正确的结果。我确信函数中的一切都很好,除了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))发布于 2017-10-03 14:50:47
import ipdb
ipdb.set_trace()试着在你想要输入函数的地方使用它。
https://stackoverflow.com/questions/46538596
复制相似问题