Boost logo

Boost :

From: Dimitar Dimitrov (mail_mitko_at_[hidden])
Date: 2007-10-11 03:20:25

 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,
 Check out the hottest 2008 models today at Yahoo! Autos.

Boost list run by bdawes at, gregod at, cpdaniel at, john at