I took a old solver from about 6 years ago and finally updated it to opencl, so much faster

I didn’t bother to do a exact speed comparison, but its pretty substantial.

I then moved on to doing this in opencl on a 3D grid, not sure the implementation is correct but it generates some interesting results.

It took me a bit to change the 5 point stencil to a 7 point for the volume grid, I originally tried it with connected point grid but things got really weird trying to reuse the 2D versions code which is based on connected points