A metaprogramming approach to generating optimized code for algorithms in linear algebra