Fix Inverse for LUFactorized

This fixes the result of `a.factorize_into()?.inv()?` when `a` is
column-major.

Author: jturner314

ndarray-linalg/src/solve.rs

@@ -323,8 +323,15 @@ where
     type Output = Array2<A>;
 
     fn inv(&self) -> Result<Array2<A>> {
+        // Preserve the existing layout. This is required to obtain the correct
+        // result, because the result of `A::inv` is layout-dependent.
+        let a = if self.a.is_standard_layout() {
+            replicate(&self.a)
+        } else {
+            replicate(&self.a.t()).reversed_axes()
+        };
         let f = LUFactorized {
-            a: replicate(&self.a),
+            a,
             ipiv: self.ipiv.clone(),
         };
         f.inv_into()