如何使用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 (将#修改为@）