Галид: график многоступенчатого трубопровода без встраивания

Question

Я пытаюсь написать модульный многоступенчатый конвейер обработки, но у меня возникают проблемы с планированием.

Структура кода выглядит следующим образом:

#include <halide/Halide.h>

Halide::Var x, y, c;

Halide::Func producer(Halide::Func in) {
    Halide::Func producer("producer");
    producer(x, y, c) = in(x, y, c);
    return producer;
}

Halide::Func rectification(Halide::Func in, const Halide::Image<float>& rectificationMapBuffer)
{
    // Fractional pixel positions according to rectification map
    Halide::Expr x_in_frac("x_in_frac");
    Halide::Expr y_in_frac("y_in_frac");
    x_in_frac = rectificationMapBuffer(x * 2 + 0, y);
    y_in_frac = rectificationMapBuffer(x * 2 + 1, y);

    // Cast fractions down to integers. This allows to address the pixels
    // surrounding the fractional position
    Halide::Expr x_in("x_in");
    Halide::Expr y_in("y_in");
    x_in = Halide::cast(Halide::Int(32), x_in_frac);
    y_in = Halide::cast(Halide::Int(32), y_in_frac);

    // Linearly interpolate pixel values
    Halide::Func interpolate("interpolate");
    interpolate(x, y, c) =
        Halide::lerp(Halide::lerp(in(x_in + 0, y_in + 0, c), in(x_in + 1, y_in + 0, c),
                                  Halide::fract(x_in_frac)),
                     Halide::lerp(in(x_in + 0, y_in + 1, c), in(x_in + 1, y_in + 1, c),
                                  Halide::fract(x_in_frac)),
                     Halide::fract(y_in_frac));

    Halide::Func rectification("rectification");
    rectification(x, y, c) = Halide::cast(Halide::UInt(8), interpolate(x, y, c));

    return rectification;
}

Halide::Func schedule(const Halide::Image<uint8_t>& image, const Halide::Image<float>& rectificationMap) {
    Halide::Func clamped;
    clamped = Halide::BoundaryConditions::repeat_edge(image);

    Halide::Func producerFunc;
    producerFunc = producer(clamped);

    Halide::Func consumerFunc;
    consumerFunc = rectification(producerFunc, rectificationMap);

    producerFunc.compute_root();

    consumerFunc.compile_jit();

    return consumerFunc;
}

int main(int argc, char *argv[])
{
    int width = 100;
    int height = 100;
    int channels = 3;
    Halide::Image<uint8_t> input(width, height, channels);
    Halide::Image<float> rectificationMap(width * 2, height);
    Halide::Buffer output(Halide::UInt(8), width, height, channels, 0);

    Halide::Func f = schedule(input, rectificationMap);
    f.realize(output);
    return 0;
}

Теперь проблема в том, что выражение compute_root. Я получаю ошибку:

The pure definition of Function rectification calls function producer in an 
unbounded way in dimension 0

Когда я удаляю часть compute_root, функция продюсера встраивается, и я не вижу никаких проблем.

Я попытался добавить ограничения.bound в функцию расписания, но это не помогло. Может кто-нибудь помочь мне понять, что означает эта ошибка?