convolution.chpl

proc convShape() {}

iter regions(dom: domain(2,int), k: int, stride: int) {
    const (h,w) = convShape(dom.shape,k,stride,padding=0);
    for (i,j) in {0..#h,0..#w} {
        yield dom[i * k * stride..#k,j * k * stride..#k];
    }
}