B为犰狳中对称矩阵时的A * B计算

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

有什么方法可以在犰狳中将对称矩阵乘以稠密矩阵(并利用我们有一个对称矩阵这一事实)?我知道 BLAS 中的 DSYMM 例程,但我正在处理的矩阵是 mat 类型(它们不是数组),所以我想知道犰狳中是否有一些东西使用 DSYMM来自 BLAS 的例程。 我在犰狳文档中找不到任何内容。

源分发中的文件 include/armadillo_bits/def_blas.hpp(BLAS 接口)和 include/armadillo_bits/glue_times_meat.hpp(实现乘法)根本不引用 dsymm(与 dgemm 形成对比) ) 因此即使您的矩阵可能是对称的,乘法也不会利用这一点。

但是,如果您需要使用 dsymm,您可以使用矩阵的 memptr() 方法来获取对原始连续内存的访问权,然后自己直接调用 dsymm。 .