2

我在 C++ 中遇到了这个错误。我正在尝试使用 multi_array 实现 Strassen 矩阵乘法。我将一个数组分配给另一个相同维度的数组。像那样A11[i][j][k] = A[i][j][k]。我认为原因是那种台词。

断言失败：(size_type(idx - index_bases[0]) < extents[0])，函数访问，文件 /usr/local/Cellar/boost/1.65.1/include/boost/multi_array/base.hpp，第 136 行。中止陷阱：6

你知道原因吗？这个错误是什么意思？

typedef boost::multi_array<int_type, 3> array_type;
array_type::extent_gen extents;

array_type A(boost::extents[size][size][noc]);
array_type B(boost::extents[size][size][noc]);
array_type C(boost::extents[size][size][noc]);

 std::fill( A.origin(), A.origin() + A.num_elements(), 0 );
 std::fill( B.origin(), B.origin() + B.num_elements(), 0 );
 std::fill( C.origin(), C.origin() + C.num_elements(), 0 );

array_type Strr(int size,int noc,array_type A,array_type B, array_type C) {

    if(size == 2) {  //2-order
       C=Matrix_Multiply(size,noc, A, B, C);

    } else {
        //
        for(int i=0; i<size/2; i++) {
            for(int j=0; j<size/2; j++) {
                 for(int k=0; k<noc; j++) {

                A11[i][j][k] = A[i][j][k] ;

                A12[i][j][k]  = A[i][j+size/2][k] ;

        }

    }

}

我的代码是这样的：我不知道问题是什么。

错误：断言失败：(size_type(idx - index_bases[0]) < extents[0])，函数访问，文件 /usr/local/Cellar/boost/1.65.1/include/boost/multi_array/base.hpp，行136.

4

1 回答 1

2

在最里面的循环中，您有：

        for (int k = 0; k < noc; j++) {

您的意思一定++k是++j：

        for (int k = 0; k < noc; ++k) {

我也会简化 main ：

int dim[] = {size,size,noc};
array_type A(dim), B(dim), C(dim);

默认情况下进行值初始化。

的想法multi_array是数组自我描述，而不是您传递单独的参数（size和noc例如）：

array_type Strr(array_type A, array_type B) {
    static_assert(array_type::dimensionality == 3, "static invariant");
    size_t size = A.shape()[0];
    size_t noc  = A.shape()[2];

    assert(A.shape()[0] ==  A.shape()[1]);
    assert(std::equal_range(A.shape(), A.shape()+3, B.shape()));
    assert(std::equal_range(A.shape(), A.shape()+3, C.shape()));
于 2017-11-07T14:19:32.067 回答