这是我第一次提出问题,如果我做错了什么,请原谅我。
我正在尝试在 matlab 中编写 Strassen 的算法,它似乎可以工作,但速度很慢(取决于截止时间,对于 64x64 矩阵,它已经可能需要一秒钟以上的时间)。
它比简单的实现(有 3 个循环)要慢。我做错了什么吗?下面是函数,输入的大小已经正确(2^n)
function [ D ] = ActualStrassen( PaddedA, PaddedB )
[m n] = size(PaddedA);
if n < 5
D = PaddedA*PaddedB;
else
A11 = PaddedA( 1:n/2 , 1:n/2 );
A12 = PaddedA( 1:n/2 , n/2+1:end);
A21 = PaddedA( n/2+1:end , 1:n/2 );
A22 = PaddedA( n/2+1:end , n/2+1:end);
B11 = PaddedB( 1:n/2 , 1:n/2 );
B12 = PaddedB( 1:n/2 , n/2+1:end);
B21 = PaddedB( n/2+1:end , 1:n/2 );
B22 = PaddedB( n/2+1:end , n/2+1:end);
M1=ActualStrassen( A11 + A22 , B11 + B22);
M2=ActualStrassen( A21 + A22 , B11 );
M3=ActualStrassen( A11 , B12 - B22);
M4=ActualStrassen( A22 , B21 - B11);
M5=ActualStrassen( A11 + A12 , B22 );
M6=ActualStrassen( A21 - A11 , B11 + B12);
M7=ActualStrassen( A12 - A22 , B21 + B22);
D11= M1 + M4 - M5 + M7;
D12= M3 + M5;
D21= M2 + M4;
D22= M1 - M2 + M3 + M6;
D = [D11, D12; D21, D22];
end
end