SpGEMM: Added CUDA first implementation. Needs further tweaks.

Only works for matrices with up to 256 nonzeros per row.
Extension for general sizes on the way.