Question

我正在编写一堆（〜打）算法，以下列模式迭代处理向量：

ArrayXd prev;
ArrayXd curr;

prev = some_initial_vector();
for (i = 0; i < N; ++i) {
    // do lots of stuff, many lines to compute curr
    // use lots of algorithm specific variables initialised above
    ...
    prev = curr;
}
return curr;

我希望有一种方法可以将curr的值的整个历史记录作为ArrayXXd中的行返回。

我尝试通过编写两个展示curr句柄的类来解决此问题，一个作为ArrayXd &，另一个作为Block<ArrayXXd, 1, -1>，但失败为it's not possible to reassign Blocks。

解决这个问题的好方法是什么？也许我可以将Block或ArrayXd自己存储在std::vector中，然后存储在convert that to an ArrayXXd中。

修改：添加了示例输入，输出

struct NotAccumulator {
    typedef ArrayXd return_type;
    ArrayXd curr;
    ArrayXd prev;
    void record () {}
    return_type result() {
        return prev;
    }
};

struct RowAccumulator {
    typedef ArrayXXd return_type;
    ArrayXd curr;
    ArrayXd prev;
    RowAccumulator(const uint N) {
        history.reserve(N);
    }

    void record () {
        history.push_back(curr);
    }
    return_type result () {
        uint rows = history.size();
        uint cols = history[0].size();
        ArrayXXd result_matrix (rows, cols);
        for(uint i = 0; i < rows; ++i) {
            result_matrix.row(i) = Map<ArrayXd> (history[i].data(), cols);
        }
        return result_matrix;
    }

private:
    std::vector<ArrayXd> history;
};

template <typename Accumulator>
typename Accumulator::return_type add_one(const ArrayXd & start, const uint how_many_times, Accumulator & u) {
    u.prev = start;
    for (uint i = 0; i < how_many_times; ++i) {
        u.curr = 1 + u.prev;
        u.record();
        u.prev = u.curr;
    }
    return u.result();
}

ArrayXd start (3);
start << 1, 0, -1;
NotAccumulator notAccumulator;
RowAccumulator rowAccumulator (5);
cout << add_one(start, 5, notAccumulator) << endl;
// outputs 6 5 4
cout << add_one(start, 5, rowAccumulator) << endl;
// outputs 2 1 0\n 3 2 1\n ... 6 5 4

Answer 1

如果您想尽可能避免数据副本，可以提前分配N - 行ArrayXXd并避免使用return。

prev和curr等缓冲区也是不必要的。以下代码还假定事先知道迭代次数N。

void GenerateHistory(int N, ArrayXXd* result) {
    result.row(0) = gen_row0_without_prev();
    for(int curr=1; curr<N; curr++) {
        result.row(curr) = gen_row_with_prev_row(result.row(curr-1));
    }
}

你可以用

调用这个函数

int N = 100;
int m = 20;
ArrayXXd history(N, m);
GenerateHistory(N, &history);

修改

让我们把复制问题放在一边。这是一个应根据您的代码示例满足要求的更新代码。

ArrayXXd add_one(const ArrayXd& start, int num_iterations, bool acc) {
    if (acc) {
        ArrayXXd result = start;
        for(int i=0; i<num_iterations; i++) {
            result = result + 1;
        }
        return result;
    } else {
        ArrayXXd result(num_iterations, start.cols());
        result.row(0)=start+1;
        for(int i=1; i<num_iterations; i++) {
            result.row(i)=result.row(i-1);
        }
        return result;
    }
}

ArrayXd start (3);
start << 1, 0, -1;
bool acc = true;
cout << add_one(start, 5, !acc) << endl;
// outputs 6 5 4
cout << add_one(start, 5, acc) << endl;
// outputs 2 1 0\n 3 2 1\n ... 6 5 4

返回所有行，或仅返回最后一行

1 个答案:

修改