From: Harald H. Soleng (harald.soleng_at_[hidden])
Date: 2007-02-15 08:49:45
Gunter Winkler wrote:
> On Thursday 15 February 2007 10:04, Harald H. Soleng wrote:
>> 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?
> I think it is important to check the orientation of the symmetric matrix (row
> major, column major) and the stored part (lower/upper triangle). Did you use
> this? (The default can be found in fwd.hpp)
The packed storage of symmetric_matrix is in fact incompatible with BLAS
dsymm. dsymm assumes full N x N matrix
filled either in the upper or lower part.