pycuda path for elementwise matrix multiplication works, but is a bit slower than numpy with ATLAS and is commented out