c++ - 如何使用 std::valarray 来存储/操作连续的二维数组？

Question

如何使用 std::valarray 来存储/操作二维数组？

我想看一个二维数组的例子，其中的元素由行/列索引访问。像这样的伪代码:

matrix(i,j) = 42;

如何初始化这样一个数组的例子也很好。

我已经知道 Boost.MultiArray、Boost.uBlas 和 Blitz++。

请随意回答为什么我不应该在我的用例中使用 valarray。但是，我希望多维数组的内存是一个连续的(列 x 行) block 。没有 Java 风格的嵌套数组。

Answer 1

在我的脑海中:

template <class element_type>
class matrix
{
public:
    matrix(size_t width, size_t height): m_stride(width), m_height(height), m_storage(width*height) {  }

    element_type &operator()(size_t row, size_t column)
    {
        // column major
        return m_storage[std::slice(column, m_height, m_stride)][row];

        // row major
        return m_storage[std::slice(row, m_stride, m_height)][column];
    }

private:
    std::valarray<element_type> m_storage;
    size_t m_stride;
    size_t m_height;
};

std::valarray 提供了许多有趣的方式来访问元素，通过切片、掩码、多维切片或间接表。参见 std::slice_array、std::gslice_array、std::mask_array 和 std::indirect_array更多详情。