循环通过圣人的矩阵行

时间:2011-12-12 14:04:25

标签: python matrix numpy sage

我正在尝试在圣人中编写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;

1 个答案:

答案 0 :(得分:4)

你让事情变得太复杂了。如果W不是稀疏矩阵,则可以执行

for row in W:

由于您还需要行索引,因此可以使用Python的内置enumerate

for i, row in enumerate(W):

或(uglier)

for i in xrange(len(W.shape[0])):