Я пытаюсь реализовать функцию без матрицы для использования с сопряженным градиентом в Eigen
.Моя матричная функция - это в основном приложение нескольких матриц последовательно к вектору, поэтому я предпочитаю реализовывать свой собственный матричный объект и делать это самостоятельно.
На веб-странице Eigen
есть хороший учебник о том, как этого добиться: https://eigen.tuxfamily.org/dox/group__MatrixfreeSolverExample.html
Этот подход требует открытия пространства имен Eigen::internal
, чтобы добавить некоторые специализации для моего собственного объекта матрицы,AFAIK, писать вещи в пространствах имен других библиотек - плохая практика, и она должна быть реализована с использованием точек настройки.
Возможно ли реализовать это с помощью точек настройки?Eigen
в этом вопросе не так?