gpt4 book ai didi

c++ - 调整 Eigen::Ref 大小的解决方法

转载 作者:行者123 更新时间:2023-11-30 05:09:08 26 4
gpt4 key购买 nike

我想使用 Eigen::Ref 使非模板函数使用 Eigen::Matrix 参数。我的问题是,在这些函数中,我可能必须调整 Eigen::Ref 引用的矩阵的大小。我明白,一般来说,不应调整 Eigen::Ref 的大小,因为它可以映射到表达式或矩阵 block ,但就我而言,我确信我的 Eigen::Ref 背后是一个 Eigen::Matrix。

为了说明这一点:

#include "Eigen/Dense"

void add(Eigen::Ref<Eigen::MatrixXd> M, const Eigen::Ref<const Eigen::MatrixXd> &A, const Eigen::Ref<const Eigen::MatrixXd> &B) {
M=A+B;
}

int main() {
Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic, Eigen::ColMajor, 32, 32> M(2,3);
Eigen::Matrix<double, 2, 2> A;
Eigen::Matrix<double, 2, 2> B;

add(M,A,B);
}

在运行时给出:

void Eigen::DenseBase<Derived>::resize(Eigen::Index, Eigen::Index) [with Derived = Eigen::Ref<Eigen::Matrix<double, -1, -1> >; Eigen::Index = long int]: Assertion `rows == this->rows() && cols == this->cols() && "DenseBase::resize() does not actually allow to resize."' failed.

我试图欺骗它:

void add(Eigen::Ref<Eigen::MatrixXd> M, const Eigen::Ref<const Eigen::MatrixXd> &A, const Eigen::Ref<const Eigen::MatrixXd> &B) {
Eigen::Ref<Eigen::Matrix<double,2,2>> MM(M);
MM=A+B;
}

但我在运行时得到了:

Eigen::internal::variable_if_dynamic<T, Value>::variable_if_dynamic(T) [with T = long int; int Value = 2]: Assertion `v == T(Value)' failed.

那么,我该怎么做才能处理这个问题呢?在 Eigen 文档中,针对使用 MatrixBase 作为参数的模板函数解决了调整大小的问题,但对于 Eigen::Ref?

最佳答案

这是一个使用成员函数指针和粗暴转换的 hacky 解决方案:

#include <iostream>
#include <Eigen/Core>
template<class MatrixType>
struct ResizableRef
{
typedef typename MatrixType::Scalar Scalar;
class MatrixDummy;
typedef void (MatrixDummy::*ResizeFunctor)(Eigen::Index rows, Eigen::Index Cols);
typedef Scalar* (MatrixDummy::*DataGetter)();

MatrixDummy *m;
const ResizeFunctor resizer;
const DataGetter getData;

template<class Derived>
ResizableRef(Eigen::MatrixBase<Derived>& M)
: m(reinterpret_cast<MatrixDummy*>(&M))
, resizer(reinterpret_cast<ResizeFunctor>((void (Derived::*)(Eigen::Index, Eigen::Index)) &Derived::resize))
, getData(reinterpret_cast<DataGetter>((Scalar* (Derived::*)()) &Derived::data))
{ }

template<class Derived>
ResizableRef& operator=(const Eigen::EigenBase<Derived>& other)
{
(m->*resizer)(other.rows(), other.cols());
MatrixType::Map((m->*getData)(), other.rows(), other.cols()) = other;
}
};

void foo(ResizableRef<Eigen::MatrixXd> A)
{
A = Eigen::Matrix2d::Identity();
}

int main(int argc, char *argv[])
{
using namespace Eigen;
MatrixXd A;
Matrix<double, Dynamic, Dynamic, Eigen::ColMajor, 20, 12> B;
Matrix<double, 2, Dynamic, Eigen::ColMajor, 2, 4> C;
Matrix2d D;
foo(A);
foo(B);
foo(C);
foo(D);

std::cout << A << "\n\n" << B << "\n\n" << C << "\n\n" << D << '\n';
}

这可能违反了严格的别名规则,我通常建议重新考虑您的设计。然而,它应该在没有不必要的运行时分配的情况下工作,并且它可以安全地防止一些错误的使用:

    MatrixXf fail1;
Matrix3d fail2;
foo(fail1); // fails at compile time
foo(fail2); // fails at runtime

关于c++ - 调整 Eigen::Ref 大小的解决方法,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/46339455/

26 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com