矩阵乘法的Strassen算法

标签：

伪代码

1. Strassen (A, B)
2. n = A.rows
3. let C be a n*n matrix
4. if　n == 1
5.  　　C₁₁ = A₁₁ * B₁₁
6. else
7.        S1 = B12 - B22
8.        S2=A11+A12;
9.        S2=A11+A12;
10.  　　S3=A21+A22;
11.  　　S4=B21-B11;
12.  　　S5=A11+A22;
13.  　　S6=B11+B22;
14.  　　S7=A12-A22;
15.        S8=B21+B22;
16.    　 S9=A11-A21;
17.  　　S10=B11+B12;
19.        P1 = Strassen(A11 * S1)
20.        P2 = Strassen(S2 * B22)
21.        P3 = Strassen(S3 * B11)
22.        P4 = Strassen(A22 * S4)
23.        P5 = Strassen(S5 * S6)
24.        P6 = Strassen(S7 * S8)
25.        P7 = Strassen(S9 * S10)
26.       
27.        C11 = P5 + P4 - P2 + P6
28.        C12 = P1 + P2
29.        C21 = P3 + P4
30.        C22 = P5 + P1 - P3 - P4
31.  
32. return C

实现

矩阵乘法的Strassen算法

标签：

原文地址：http://www.cnblogs.com/fanjinhua/p/4819778.html