j9263178

> What's your setting? Inline? Yes, running with `%matplotlib inline` will reproduce the bug.

cuSolver uses column major so we always need to do some transpose to get the desired result with [cusolverDngeqrf()](https://docs.nvidia.com/cuda/cusolver/index.html#cuSolverDN-lt-t-gt-geqrf) and [cusolverDnorgqr()](https://docs.nvidia.com/cuda/cusolver/index.html#cuSolverDN-lt-t-gt-orgqr), maybe it is better we use Xgelqf and Xorglq...