https://gcc.gnu.org/bugzilla/show_bug.cgi?id=68600
--- Comment #9 from Thomas Koenig <tkoenig at gcc dot gnu.org> --- > I took the example code found in > http://wiki.cs.utexas.edu/rvdg/HowToOptimizeGemm/ where the register based > vector computations are explicitly called via the SSE registers and > converted it to use the builtin gcc vector extensions. I had to experiment > a little to get some of the equivalent operations of the original code. Nice one, I think we can use this as a basis for a new library function in 7.1. It probably will be necessary to add some zero-padding in some places where the matix dimensions are not divisible by four.