Add the gpu implementation of the solver. It is not parallelized yet (and very slow).