10

1

Algorithms for fast multiplication of polynomials and integers have well-known conceptual explanations. A good survey paper is Daniel J. Bernstein's Fast Multidigit Multiplication for Mathematicians.

For example, Karatsuba's trick and FFT-based multiplication both fit into the evaluate-and-interpolate scheme. These rely on the fact that evaluation at a point is a ring homomorphism. Karatsuba's trick cheaply evaluates the product of two linear polynomials at 0, 1 and $\infty$ without multiplying them out explicitly and then interpolates the values to get the polynomial product. Similarly, the n-point FFT efficiently evaluates a polynomial at the nth roots of unity, and the inverse FFT efficiently interpolates a polynomial from its values at those same points.

You can also go beyond evaluation and use any convenient ring homomorphisms.

My question is whether Strassen's algorithm for fast matrix multiplication can be explained in conceptual terms along these lines. The core of the algorithm is Strassen's trick for computing products of 2x2 matrices over noncommutative rings. Looking at Strassen's formula, it's hard to shake the feeling that some related approach should work.

flag
2 
You might be interested in Cohn, Kleinberg, Szegedy, and Umans' Group-theoretic Algorithms for Matrix Multiplication, although I don't think it answers this specific question. – Qiaochu Yuan Sep 23 at 21:41
1 
For a diametrically opposite way (to the one desired in the question) of "explaining" things, there's this brute-force paper: ieeexplore.ieee.org/iel5/4235/5438925/… – S. Sra Sep 23 at 22:19
1 
And for an explanation / guessing of how Strassen might have obtained his equations, see the textbook material: integrator-crimea.com/ddu0177.html – S. Sra Sep 23 at 22:21

1 Answer

3

To my poor knowledge, fast matrix multiplication is not based on rings homomorphisms, but on the notion of tensorial rank. More precisely, the trilinear form $(A,B,C)\mapsto{\rm Tr}(ABC)$ extends as a linear form $L$ over $$M_{n\times m}\otimes M_{m\times p}\otimes M_{p\times n}.$$

The tensorial rank of $L$ is the minimal number $r$ in which you can write $$L=\sum_{j=1}^re_j\otimes f_j\otimes g_j,$$ with $e,f,g$ linear forms.

Notice that if you had only a pair of space, then the tensorial rank over $E\otimes F$ is just the rank of linear algebra.

link|flag

Your Answer

Get an OpenID
or

Not the answer you're looking for? Browse other questions tagged or ask your own question.