Question

我正在尝试构建一个代码，我必须在堆中声明一个大型数组。同时我将使用boost库来执行一些矩阵计算（如Fill a symmetric matrix using an array 中所示）。

我的局限是两个：我将处理大型数组和矩阵，所以我必须在堆上声明所有内容，我必须使用数组而不是向量。

然而，对于许多人来说，我面临着相当微不足道的问题......当填充矩阵时，最后一个元素无法正确填充。所以虽然我希望得到

[3,3]((0,1,3),(1,2,4),(3,4,5))

代码的输出是

[3,3]((0,1,3),(1,2,4),(3,4,2.6681e-315))

我正在ROOT6编译此代码。我认为这与此无关，我只是提到完成。

代码的一小部分示例如下

#include <iterator>
#include <iostream>
#include <fstream>

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

using namespace std;

int test_boost () {

    using namespace boost::numeric::ublas;
    symmetric_matrix<double, upper> m_sym1 (3, 3);

    float* filler = new float[6];
    for (int i = 0; i<6; ++i) filler[i] = i;
    float const* in1 = filler;

    for (size_t i = 0; i < m_sym1.size1(); ++ i)
        for (size_t j = 0; j <= i && in1 != &filler[5]; ++ j)
            m_sym1 (i, j) = *in1++;
    delete[] filler;
    std::cout << m_sym1 << std::endl;
    return 0;
}

如何解决这个问题？

Answer 1

数组和指针不是类类型的对象，它们没有成员。您已经有float *，它是filler。

float const* in1 = filler; // adding const is always allowed

Answer 2

我已经通过将&filler[5]更改为&filler[6]来最终解决问题。所以有效的版本见下文

#include <iterator>
#include <iostream>
#include <fstream>

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

using namespace std;

int test_boost () {

    using namespace boost::numeric::ublas;
    symmetric_matrix<double, upper> m_sym1 (3, 3);

    float* filler = new float[6];
    for (int i = 0; i<6; ++i) filler[i] = i;
    float const* in1 = filler;

    for (size_t i = 0; i < m_sym1.size1(); ++ i)
        for (size_t j = 0; j <= i && in1 != &filler[6]; ++ j)
            m_sym1 (i, j) = *in1++;
    delete[] filler;
    std::cout << m_sym1 << std::endl;
    return 0;
}

运行此代码会产生以下输出

[3,3]((0,1,3),(1,2,4),(3,4,5))

c ++ - 使用存储在堆上的数组填充对称矩阵

2 个答案: