I wanted to try coding a python program that used a Strassen algorithm to multiply two extremely large matrices. I my first hunch was to use numpy and after a bit of googling I found my way to this SO question.

http://stackoverflow.com/questions/3218645/handling-large-dense-matrices-in-python

People who answer their own questions are nice, but anyway this HDF5 library looks promising. I’m going to spend some time today an maybe next weekend attempting to code this.