I was surpriced too. I used Boost.MultiArray library only and the performance was poor (may be because of bad ussage). Using standard C++ array shows better performance. Finally, fortran arrays are faster (as I heard), so I want to write function that stores a matrix in a vector.