Boundaries don't work with OpenCL

Boundary conditions work fine on CPU and CUDA, but they don't have any effect on OpenCL. Here is a Poiseuille flow test case that shows that, just change the target in the third block. gpu_boundary.ipynb