[boost] [ublas] multiplication of triangular matrices

Hello, I have to efficiently multiply A'*A, where A is an upper triangular matrix with dimension 150. Defining A as a ublas::triangular_matrix<double, upper> and using the function "prod" is very slow (I used -DNDEBUG flag and noalias). If I directly define ublas::matrix<double> A, and use atlas::gemm I get better results. My question is: Is there a way to use the fact that A is triangular in order to speed up the latter computation. (I tried using triangular_matrix as an input for atlas::gemm but I got many compilation errors :)) Thank you, Dimitar --------------------------------- Check out the hottest 2008 models today at Yahoo! Autos.