for i in range(n):
    sum += A[row, i] * b[i]
c[row] = sum
A = np.array(np.random.random((m, n)), dtype=np.float32)
B = np.array(np.random.random(m), dtype=np.float32)
C = np.empty_like(B)
s = time()
dA = cuda.to_device(A)
dB = cuda.to_device(B)
dC = cuda.to_device(C)
cu_matrix_vector[(m + 511) / 512, 512](dA, dB, dC)
dC.to_host()
print(C)
e = time()
tcuda = e - s
Numba Matrix Vector multiplication