I attached a simple implementation of a matrix-matrix product (gemm.hpp) which has less

then 250 lines of code.  I also attached a simple test programm (matprod.cc).  The test

Programm has a simple wrapper for the implementation in gemm.hpp.  For a real applications

you should use traits like for the BLAS bindings.  You can compile it with

g++ -O3 -Wall -std=c++11 -I $BOOST_PATH matprod.cc

where BOOST_PATH contains your boost implementation.  Let me know if this works

for you.

Cheers,

Michael