作为一个小练习,我开始在
python中使用数字代码,我正在尝试制作LDLT算法.只是为了“弄湿脚”.
但是我似乎缺乏对numpy数组的基本理解.请参阅以下示例:
def ldlt(Matrix):
import numpy
(NRow,NCol) = Matrix.shape
for col in range(NCol):
Tmp = 1/Matrix[col,col]
for D in range(col+1,NCol):
Matrix[col,D] = Matrix[D,col]*Tmp
if __name__ == '__main__':
import numpy
A = numpy.array([[2,-1,0],[-1,2,-1],[0,2]])
ldlt(A)
这个例子不是我正在处理的完整代码.但是,尝试运行它,并在Matrix [col,D] = …设置断点.
我对第一次评估的期望是第0行第1列(起始值-1)被设置为等于= -1 *(1/2)= -0.5.
但是在运行代码时,它似乎设置为0.为什么?
必须有一些基本的东西,我还没有真正理解?
在此先感谢所有帮助我的人.
编辑1:
Python Ver.:3.3
Tmp.:变为0.5(由我的调试器报告).