ctx = create_some_context()
filter_buf = Buffer(ctx,
MEM.READ_ONLY | MEM.COPY_HOST_PTR,
hostbuf=GAUSSIAN_BLUR)
in_img = lena()
h, w = map(int32, in_img.shape[:2])
# in pyopencl 2018.2.2 channel orders other than RGBA cause segmentation fault
i4 = zeros((h, w, 4), dtype=uint8)
i4[:, :, 0] = in_img
in_img_buf = image_from_array(ctx, i4, 4)
fmt = ImageFormat(CHO.RGBA, CHANNEL.UNSIGNED_INT8)
out_img_buf = Image(ctx, MEM.WRITE_ONLY, fmt, shape=(w, h))
pt = perf_counter()
prg = Program(ctx, load_cl_text("convolution.cl")).build()
TIMES["Compilation"] = perf_counter() - pt
pt = perf_counter()
with CommandQueue(ctx) as queue:
prg.convolution(queue, (w, h), None, in_img_buf, out_img_buf, filter_buf,
int32(5),
Sampler(ctx, False, ADDRESS.CLAMP_TO_EDGE, FILTER.NEAREST))
TIMES["Execution"] = perf_counter() - pt
pt = perf_counter()
dest = zeros(i4.shape, dtype=uint8)
enqueue_copy(queue, dest, out_img_buf, origin=(0, 0), region=(w, h))
TIMES["Copying"] = perf_counter() - pt
in_img_buf.release()
out_img_buf.release()
filter_buf.release()
print("\n".join("%s:\t%g" % i for i in TIMES.items()))
show_img(dest[:, :, 0])
