我正在尝试在圣人中编写Graham-Schmidt进程,但无法弄清楚如何循环遍历数组的行。
def graham_schmidt(W):
a=0
U=W
for i in W.dims()[0]:# this is the not working part
print w
a=a+1
for j in xrange(0,-2):
a=a+1
U[i]=U[i]-(transpose(U[j])*w)/(transpose(U[j])*U[j])*U[j]
return a;
答案 0 :(得分:4)
你让事情变得太复杂了。如果W
不是稀疏矩阵,则可以执行
for row in W:
由于您还需要行索引,因此可以使用Python的内置enumerate
:
for i, row in enumerate(W):
或(uglier)
for i in xrange(len(W.shape[0])):