Strassen's algorithm: Difference between revisions

m
Line 331:
===Recursive===
Output is the same as the dynamically padded version.
<lang Julia>function Strassen(A::Matrix, B::Matrix)
n = size(A, 1)
if n == 1
Line 348:
P2 = Strassen(A11 + A22, B11 + B22)
P3 = Strassen(A11 - A21, B11 + B12)
P4 = Strassen(A11 + A12, Matrix(B22))
P5 = Strassen(Matrix(A11), B12 - B22)
P6 = Strassen(Matrix(A22), B21 - B11)
P7 = Strassen(A21 + A22, Matrix(B11))
 
C11 = P1 + P2 - P4 + P6
4,104

edits