有没有办法在犰狳中将对称矩阵乘以密集的矩阵(并使用我们有对称矩阵的事实)?我知道BLAS中的DSYMM
例程,但我处理的矩阵是mat
类型(它们不是数组),所以我想知道犰狳中是否有一些东西使用{{1来自BLAS的常规。
我无法在犰狳文档中找到任何东西。
答案 0 :(得分:2)
源代码分发中的文件include/armadillo_bits/def_blas.hpp
(BLAS接口)和include/armadillo_bits/glue_times_meat.hpp
(实现乘法)根本不会引用dsymm
(与dgemm
形成对比)所以即使你的矩阵可能是对称的,乘法也不会利用这个。
但是,如果您需要使用dsymm
,您可以使用矩阵的memptr()
方法来获取对原始连续内存的访问权并直接调用dsymm
... < / p>