在Halide中使用split()时,我可以避免计算相同的元素吗?

时间:2015-10-21 08:39:55

标签: halide

我对Halide语言中split()的行为有疑问。

当我使用split()时,当计算区域不是分割因子的倍数时,它会在边缘处计算两次元素。例如,当计算区域为10且分割因子为4时,Halide将计算元素[0,1,2,3],[4,5,6,7]和[6,7,8,9],如跟踪trace_stores()的结果。

有没有办法在split()内循环的最后一步只计算元素[8,9]?

示例代码:

#include "Halide.h"
using namespace Halide;

#define INPUT_SIZE 10
int main(int argc, char** argv) {
    Func f("f");
    Var x("x");
    f(x) = x;

    Var xi("xi");
    f.split(x, x, xi, 4); 

    f.trace_stores();
    Image<int32_t> out = f.realize(INPUT_SIZE);
    return 0;
}

trace_stores()结果:

Store f.0(0) = 0
Store f.0(1) = 1
Store f.0(2) = 2
Store f.0(3) = 3
Store f.0(4) = 4
Store f.0(5) = 5
Store f.0(6) = 6
Store f.0(7) = 7
Store f.0(6) = 6
Store f.0(7) = 7
Store f.0(8) = 8
Store f.0(9) = 9

1 个答案:

答案 0 :(得分:1)

这可能但很难看。 Halide通常假设它可以任意重新评估Func中的点,并且输入不会与输出混淆,因此重新计算边缘附近的几个值总是安全的。

事实上,这很重要。可能还有其他方法可以实现您的目标。

无论如何,解决方法是使用显式RDom告诉Halide精确迭代的内容:

// No pure definition
f(x) = undef<int>(); 

// An update stage that does the vectorized part:
Expr w = (input.width()/4)*4;
RDom r(0, w);
f(r) = something;
f.update(0).vectorize(r, 4);

// An update stage that does the tail end:
RDom r2(input.width(), input.width() - w);
f(r2) = something;
f.update(1); // Don't vectorize the tail end