I was going to post something similar, with one addition: One performance problem with Blitz was that the expression templates obfuscated the code enough that

I was going to post something similar, with one addition: One performance problem with Blitz was that the expression templates obfuscated the code enough that the compiler bailed on doing SIMD vectorization of the evaluation loops. That was fixed this year at least for the Intel compiler (gcc doesn't have pragmas to control vectorization), so now performance can even exceed Fortran for large-ish arrays. I made some plots while fiddling with this if you are interested: http://governator.ucsc.edu/filer/blitzbench_r1845/blitzcomp.html http://governator.ucsc.edu/filer/blitzbench_r1845/blitzcomp....