Re: [ublas] symmetric_matrix and BLAS dsymm

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.

Harald