Eigen 自定义表达式类型（或 CWise Op 仿函数？）

Question

考虑这个简单的函数（一种伪代码）：

template <typename Derived>
void f (const MatrixBase<Derived>& input1,
        const MatrixBase<Derived>& input2,
        const MatrixBase<Derived>& input3,
        Vec6& output)
{
    part1(output) = <an expr using the three inputs>;
    part2(output) = <another expr using the three inputs>;
}

其中part1 和part2 是分别选择output 的前3 行和后3 行的块表达式。我知道output总是一个有6个系数的列向量。

我想避免传递 output 参数。但我也想避免返回明确的Vec6，以避免临时副本。

因此，我想返回其中一个神奇的特征表达式对象，它只存储对输入的引用，并且仅在需要时才实际进行计算。

我的问题：

1. 值得吗？鉴于它是一个“小”的 6 维向量（例如，文档说使用小的 3 维向量，编译器通常能够优化掉临时变量，因此没有太多避免它们的意义）

2. 我可以使用自定义函子和CWiseNullaryOp 实现我的返回表达式类型，还是真的需要自定义表达式类型？

谢谢

---

编辑：所以我们确定临时对象已被优化（参见 cmets）；这回答了问题 1，并使问题 2 不相关。

但是，为了好奇和更多地了解 Eigen，是否有人对我的函数的可行表达式类型有一些提示？ （假设值得设计一个）

Answer 1

假设 Vec6 的设计“正确”（参见 例如 Rule Of Zero），那么如果您的函数 f() 采用以下形式

Vec6 f(const MatrixBase<Derived>& input1,
       const MatrixBase<Derived>& input2,
       const MatrixBase<Derived>& input3)
{
    Vec6 answer;

    //Manipulate answer
    //...

    return answer; 
}

然后在大多数现代编译器上启用优化不会发生复制。请参阅关于返回值优化here 的 Wikipedia 部分。

（要使用您的对象和编译器测试理论，您可以将 Vec6 包装在您自己的类中，让所有构造函数在执行时打印，并检查您是否获得了预期的构造函数调用次数。）