Solve for Ax=B where A is a sparse matrix

Does anyone know of a library or example of openCL code which will solve for Ax=B where A is large and sparse? I do not want to calculate the inverse of A as it would be very large and dense.

I would look at ViennaCL or ArrayFire as a starting point. I don’t know their specific capabilities, but their general page looks like they might be a good starting point.