0

这是我第一次提出问题,如果我做错了什么,请原谅我。

我正在尝试在 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
4

0 回答 0