如何解决如何使用Eigen有效地进行SYR2K操作?
我想使用Eigen对Matrix进行第2级更新,并将性能与另一个BLAS库进行比较。
通常称为syr2k
操作。 (http://www.netlib.org/lapack/explore-html/db/def/group__complex__blas__level3_gaa8320d51ded07cd3038db237fd400547.html)
我尝试了以下代码,如Eigen建议的等级2更新
M.selfadjointView<Eigen::Upper>().rankUpdate(u,v,s);
(https://eigen.tuxfamily.org/dox/group__QuickRefPage.html)
问题是,在这种情况下,u
和v
必须是向量,我想使用矩阵。我可以简单地做C = alpha*A*B.transpose() + alpha*B*A.transpose() + beta*C.template selfadjointView<Eigen::Upper>();
,但问题是:
我可以做得更好吗?或者换句话说,我应该使用一种更快的方法吗?
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。