如何创建一个常量boost矩阵？

Question

如何创建一个常量boost矩阵？

以下方法不起作用：

const boost::numeric::ublas::matrix<double> arrayM(1, 3) = { {1.0, 2.0, 3.0} };

- Journey Man

---

它不起作用是因为你声明了一个常量矩阵，然后试图使用可变赋值运算符对其进行赋值。下面GMan的解决方案将可以满足你的需求。 - Dan Story

我应该对我的问题表达更清晰：有没有一种类似于我们创建const数组的方式来创建它的方法？ 例如：int billy [5] = { 16, 2, 77, 40, 12071 }; - Journey Man

Answer 1

通常类似于以下内容：

typedef boost::numeric::ublas::matrix<double> matrix_type;

const matrix_type get_matrix(void)
{
    matrix_type result(1, 3);
    result(0, 0) = 1;
    result(0, 1) = 2;
    result(0, 2) = 3;

    return result;
}

const matrix_type arrayM = get_matrix();

您也可以尝试类似以下的方法（大多数未经过测试）：

#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/io.hpp>

template <typename T, typename L = boost::numeric::ublas::row_major,
            typename A = boost::numeric::ublas::unbounded_array<T> >
class matrix_builder
{
public:
    // types
    typedef boost::numeric::ublas::matrix<T, L, A> matrix_type;
    typedef typename matrix_type::size_type size_type;

    // creation
    matrix_builder(size_type pRows, size_type pColumns) :
    mMatrix(pRows, pColumns),
    mRow(0),
    mColumn(0)
    {}

    matrix_builder& operator()(const T& pValue)
    {
        mMatrix(mRow, mColumn) = pValue;
        if (++mColumn == mMatrix.size2())
        {
            mColumn = 0;
            mRow++;
        }

        return *this;
    }

    // access
    operator const matrix_type&(void) const
    {
        return mMatrix;
    }

private:
    // non copyable
    matrix_builder(const matrix_builder&);
    matrix_builder& operator=(const matrix_builder&);

    // members
    matrix_type mMatrix;
    size_type mRow;
    size_type mColumn;
};

typedef boost::numeric::ublas::matrix<double> matrix_type;

static const matrix_type m1 = matrix_builder<double>(3, 1)
                                (1)(2)(3);

static const matrix_type m2 = matrix_builder<double>(3, 3)
                                (1)(2)(3)
                                (4)(5)(6)
                                (7)(8)(9);

int main(void)
{
    std::cout << m1 << std::endl;
    std::cout << m2 << std::endl;
}

相同的想法，更加通用。同时也更具视觉效果，这可能很不错。

- GManNickG

---

虽然这样可以工作，但我需要为每个常量矩阵类型创建一个函数。 - Journey Man

@Venkata：我现在会写一个更通用的方法。 - GManNickG

@Stacked：谢谢。我承认可能还有更好的方法，但这应该适用于我能想到的任何情况。 - GManNickG