This is a pretty common topic in many courses in software performance engineering — the methods to optimize matrix multiplication are well known, the problem is simple enough to be accessible to students, and the results are surprising to novices because the speedups are so large.
See also the discussion in this thread of some of the techniques: Julia matrix-multiplication performance