David and Nasos,

Thank you for the interest.

I heard about Armadillo and Eigen3 but I have no idea of functionalities provided by them.

For what concerns my code, one possible "weak-point" is that many operation rely on LAPACK (through boost-numeric_binding).

While I'm not sure this is compliant with uBLAS philosophy (e.g., I think to the C reimplementation of BLAS subrourines), my primary target was to bring some MATLAB/Octave functionalities to uBLAS without reinventing the wheel and the use of LAPACK appeared to me a good trade-off between high-quality code, speed of development and portability.