LSQR algorithm with PyOpenCL

As part of the FreeCAD-Ship project I created a LSQR algorithm to solve Linear systems of equations, accelerated with OpenCL using PyOpenCL.

LSQR is a great, and relatively recent, iterative method to solve linear systems which can work with unsymmetrical rectangular matrix of an arbitrary range. Fortunately this method has also the property to fit so quite well to the GPU computing paradigm.