https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Which-algorithm-is-implemented-in-DGEMM/m-p/1131382#M25613 <P>Interestingly, I've been unable to find an answer to this simple question. What is the algorithm that is used for matrix-matrix multiplications (e.g., DGEMM) in MKL? Is is classical (O(N^3)), Strassen (O(N^2.7)), or something else? Thanks.</P> Tue, 26 Sep 2017 18:52:34 GMT Laasner__Raul 2017-09-26T18:52:34Z https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Which-algorithm-is-implemented-in-DGEMM/m-p/1131382#M25613 <P>Interestingly, I've been unable to find an answer to this simple question. What is the algorithm that is used for matrix-matrix multiplications (e.g., DGEMM) in MKL? Is is classical (O(N^3)), Strassen (O(N^2.7)), or something else? Thanks.</P> Tue, 26 Sep 2017 18:52:34 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Which-algorithm-is-implemented-in-DGEMM/m-p/1131382#M25613 Laasner__Raul 2017-09-26T18:52:34Z https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Which-algorithm-is-implemented-in-DGEMM/m-p/1131383#M25614 <P>Hi Raul,</P> <P>I am afraid BLAS standard gemm uses classical O(N<SUP>3</SUP>), for algorithm design, you could follow Netlib gemm source code. Intel MKL optimized BLAS routines with SIMD instruction sets, do some work to fit data into the caches enabling contiguous, aligned accesses.</P> <P>Here's another algorithm for matrix matrix multiplication, call 3M. It split a complex matrix into two matrices, performs 3 GEMM and 4 matrix additions. For other algorithm, like Winograd which implemented for NN convolution kernel in MKL-DNN.</P> Wed, 27 Sep 2017 02:07:00 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Which-algorithm-is-implemented-in-DGEMM/m-p/1131383#M25614 Zhen_Z_Intel 2017-09-27T02:07:00Z