Registered Member
|
Hi,
I'd like to multiply a symmetric matrix (using selfAdjointView<Upper>) with a matrix. I can I do please? Thanks |
Moderator
|
|
Registered Member
|
thanks but sorry I wasn't clear.
I have an array double* ap which represents a symmetric matrix (as in lapack). So I map it like this Map<SelfAdjointView<double,Upper> > apMap (t,n); And i have a double matrix B and I want to do C=apMap*B is it possible? (if I write it like this it doesn't work the operator * doesn't exist...) |
Moderator
|
is ap in a packed format? if so Eigen cannot handle it. Otherwise, you should rather do:
Map<MatrixXd> (ap,n,n).selfadjointView<Upper>() * B; |
Registered Member
|
|
Moderator
|
The packed format is not suited for fast matrix-matrix products.
|
Registered users: Baidu [Spider], Bing [Bot], Google [Bot]