This is another adventure into the world of LuaJIT and OpenCL.
This has been the driving force behind my
LuaJIT + OpenCL library of bindings, wrappers, and classes.
I'm also writing CL options into my
Krylov solvers lua library.
This is a port of an original version done in
C++.
My original plan was to use JFNK to solve the Einstein field equations. This has finally made some progress to converging, albeit very slowly, and so I embarked on a quest to find more power.
For the LuaJIT version, in absence of a JFNK OpenCL port, I started my first attempt with Newton descent. Unfortunately I implemented the descent with respect to the continuous equations (see
here).
There's a heavy amount of
templated OpenCL code generation too.