Question

这是my previous one的后续问题。我正在构建一个具有基本矢量数学功能的Array类。我想定义一个允许用户在整个数组上映射任意函数的函数。我已经可以为预定义的二进制函数（例如operator+=进行就地添加）执行此操作，但我很难看到如何对任意数量的输入执行此操作。

template<typename T, size_t ... Ns>
class Array
{
public:
    Array() { /* allocate memory of size product(Ns) * sizeof(T) */ }

    ~Array() { /* deallocate memory */ }

    inline Array<T, Ns...>& operator+=(const Array<T, Ns...> &rhs)                   
    {
        // simple case: predefined operation on exactly two Arrays
        T *p1 = data, *p2 = rhs.data;                                                
        while (p1 != data + size) { *p1 += *p2; p1++; p2++; }                                                                            
        return *this;                                                                
    }

    template<class ... Args>
    inline const Array<T, Ns...> map(T (*fn)(arg_type<Args>...), Args ... args)
    {
        // difficult case: arbitrary operations on a variable number of Arrays
    }

    private:
        T *data;
        size_t size;
}

// example usage
double f(double a, double b) { return a + b; }
Array<double,2> x, y, z;
x.map(f, y, z);

我希望这样的内容循环遍历y和z中的所有元素，将f应用于它们，并将它们存储在x中。我以为我可以在operator+=函数上对它进行建模，但是我还没有找到一个参数包扩展来让我做类似的事情（并编译）。

Answer 1

写下这个：

template<class T, class Array2>
using Array_Size_Copy = // expression that creates a Array<T,???>, where ??? is the Ns... of Array2 above

然后

template<class F, class A0, class...Arrays>
auto array_apply( F&& f, A0&& a0, Arrays&&...arrays ) {
  using R = decay_t<decltype( f(std::declval<A0>().data.front(), std::declval<Arrays>().data.front()...) )>;

  Array_Size_Copy<R, std::decay_t<A0>> retval;

  for (int i = 0; i < a0.size(); ++i) {
    retval.data[i] = f( std::forward<A0>(a0).data[i], std::forward<Arrays>(arrays).data[i]... );
  }
  return retval;
}

添加一些断言（理想情况下是静态断言），一切都是相同的大小。

上述功能可能需要成为朋友。这意味着您可能必须设置返回值而不是使用auto，但这只是将体操类型转移到声明中。

在任意数量的可变参数上应用函数

1 个答案: