ENH: working matrix multiply in cuda