在凹形矩阵的行中循环

Question

我试图在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;

Answer 1

你把事情搞得太复杂了。如果W不是一个稀疏矩阵，那么可以这样做

for row in W:

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

^{pr2}$

或者（更丑）

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