我正在尝试实现与Eigen
中的共轭梯度一起使用的无矩阵函数。我的矩阵函数基本上是将多个矩阵连续应用于一个向量,因此我更喜欢实现自己的矩阵对象,并自行执行此操作。
Eigen
的网页上有一个很好的教程,介绍如何实现此目的:
https://eigen.tuxfamily.org/dox/group__MatrixfreeSolverExample.html
此方法需要打开Eigen::internal
命名空间,以为我自己的矩阵对象添加一些特殊化。 AFAIK,在其他库名称空间中编写内容是一种不好的做法,应使用自定义点来实现。
是否可以使用定制点来实现? Eigen
在这件事上错了吗?