|
|
Ublas : |
From: Harald H. Soleng (harald.soleng_at_[hidden])
Date: 2007-02-15 04:04:33
In an application I would like to replace a few computationally
intensive functions with BLAS functions.
For example, when replacing prod(const matrix<double>&,const
matrix<double>&) with dgemm and linking with gotoblas
I get a speed increase of a factor of between 10 and 30 depending on the
platform. I tried to replace
prod(const symmetric_matrix<double>&,const matrix<double>&) with dsymm
but with no success. I just don't get the correct
results. Any suggestions?
Harald
-- Harald H. Soleng, Dr. Philos., MSCC 12398 Chief Research Scientist Norwegian Computing Center Phone (direct): +47 2285 2671 Phone (central): +47 2285 2500 Fax (central): +47 2269 7660 http://nr.no/~soleng/