title: non-square C-order matrices in cuBLAS ( numba )
body: <p>I'm trying to use the cuBLAS functions in Anaconda's Numba package and having an issue. I need the input matrices to be in C-order. The output can be in Fortran order.</p>
  
  <p>I can run the example script provided with the package, <a href="https://github.com/ContinuumIO/numbapro-examples/blob/master/cublas/gemm.py" rel="noreferrer">here</a>. The script has two functions, <code>gemm_v1</code> and <code>gemm_v2</code>. In <code>gemm_v1</code>, the user has to create the input matrices in Fortran order. In <code>gemm_v2</code>, they can be passed to the cuda implementation of GEMM and transposed on the device. I can get these examples to work with square matrices. However, I can't figure out how to get <code>gemm_v2</code> to work with non-square input matrices. Is there a way to work with C-order input matrices that are non-square?</p>
  
  <p><strong>Note:</strong><br>
  Ideally, both the input and output matrices would stay on the device after the call to GEMM to be used in other calculations ( this is part of an iterative method ).</p>