1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33
|
#include "Halide.h"
#include <stdio.h>
using namespace Halide;
int main(int argc, char **argv) {
Func f, g, h;
Var x, y;
h(x, y) = x + y;
g(x, y) = (h(x - 1, y - 1) + h(x + 1, y + 1)) / 2;
f(x, y) = (g(x - 1, y - 1) + g(x + 1, y + 1)) / 2;
h.compute_root();
g.compute_at(f, y);
// f.trace();
Buffer<int> out = f.realize({32, 32});
for (int y = 0; y < 32; y++) {
for (int x = 0; x < 32; x++) {
if (out(x, y) != x + y) {
printf("out(%d, %d) = %d instead of %d\n", x, y, out(x, y), x + y);
return 1;
}
}
}
printf("Success!\n");
return 0;
}
|