CODEWITHSUNDEEP
c++performancevisual-c++armadilloblas
MAh2014
MAh2014

Reputation: 147

A * B computation when B is a symmetric matrix in armadillo

Is there any way to multiply a symmetric matrix by a dense one in armadillo(and use the fact that we have a symmetric matrix)? I know about DSYMM Routine in BLAS,but the matrices I'm dealing with are mat type(they are not arrays),So I wonder if there is some thing in armadillo using DSYMM Routine from BLAS. I couldn't find any thing in armadillo documentation.

Upvotes: 0

Views: 199

Answers (1)

ewcz
ewcz

Reputation: 13087

The files include/armadillo_bits/def_blas.hpp (BLAS interface) and include/armadillo_bits/glue_times_meat.hpp (implementing the multiplication) in the source distribution do not reference dsymm at all (in contrast to dgemm) so even though your matrices might be symmetric, the multiplication won't take advantage of this.

However, should you need to use dsymm, you might use the memptr() method of the matrices to obtain access to the raw contiguous memory and call dsymm directly yourself...

Upvotes: 2

Related Questions