C++ Эйген: рекурсивные функции, принимающие любую матрицу класса

Question

Я хочу иметь рекурсивную функцию

template <typename Derived> 
void f(Eigen::MatrixBase<Derived>& m) 
{ 
    size_t blockRows = ... 
    size_t blockCols = ... 
    .... 
    f(m.block(0, 0, blockRows, blockCols)); 
} 

К сожалению, это приводит к бесконечному времени компиляции рекурсии.

Первый вызов будет

f<Eigen::MatrixBase<Derived> > 

Второй будет

f<Eigen::Block<Eigen::MatrixBase<Derived>, ... > > 

Третий вызов будет

f<Eigen::Block<Eigen::Block<Eigen::MatrixBase<Derived>, ... >, ... > > 

Каждый раз, когда блок блока запрашивается ,

Какова наилучшая практика для реализации рекурсивных функций в Eigen, которые все еще работают с любым типом матрицы Eigen?

Я думаю, я должен использовать какой-то тип, который все еще обертывает одну и ту же часть памяти, но не является шаблоном выражения и оценивается.

Answer 1

Вы можете отдать свой рекурсивный блок к Eigen::Ref, чтобы избежать бесконечного типа экземпляра:

Ref<MatrixXd> bl = m.block(0, 0, blockRows, blockCols); 
f(bl); 

Чтобы оставаться полностью родовым, вы можете заменить MatrixXd на typename Derived::PlainObject.