在Python中实现Strassen算法的难度

时间:2014-09-16 04:14:44

标签: python algorithm recursion runtime-error strassen

我不明白如何递归调用我的代码。到目前为止,这是我的代码:

import numpy

B = [[5,5,5,5,5,5,5,5],[6,6,6,6,6,6,6,6],[7,7,7,7,7,7,7,7],[8,8,8,8,8,8,8,8],
 [9,9,9,9,9,9,9,9], [10,10,10,10,10,10,10,10],[11,11,11,11,11,11,11,11],       [12,12,12,12,12,12,12,12]]

A = [[5,5,5,5,5,5,5,5],[6,6,6,6,6,6,6,6],[7,7,7,7,7,7,7,7],[8,8,8,8,8,8,8,8],
 [1,1,1,1,1,1,1,1], [2,2,2,2,2,2,2,2],[3,3,3,3,3,3,3,3],[4,4,4,4,4,4,4,4]]

def main():
   strassen(A,B)

def strassen(A, B):
    A = numpy.asarray(A)
    B = numpy.asarray(B)
    lengthA = len(A)
    lengthB = len(B)
    if lengthA == 2:
        print "will calculate"
    else:       
        a, b = strassen(A[:lengthA//2, :lengthA//2], B[:lengthB//2, :lengthB//2])

        lengthA = lengthA//2
        lengthB = lengthB//2
        print a
        print b
        return a, b

我试图将a缩减为[[5,5],[6,6]]而将b缩减为[[5,5],[6,6]],但我收到了错误消息:

a, b = strassen(A[:lengthA//2, :lengthA//2], B[:lengthB//2, :lengthB//2])
TypeError: 'NoneType' object is not iterable. 

ab是在A和B的第二个整数矩阵除法之后形成的第一个2x2矩阵。请有人向我解释。感谢

1 个答案:

答案 0 :(得分:1)

您的递归终止条件中没有返回值。当我运行你的代码时,它打印"将计算"在给出错误之前。之后发生错误,因为最后一次调用时strassen函数没有返回值(lengthA == 2时)。