Bug fix

Author: termoshtt

src/triangular.rs

@@ -47,7 +47,7 @@ where
     So: DataMut<Elem = A> + DataOwned,
 {
     fn solve_triangular_into(&self, uplo: UPLO, diag: Diag, mut b: ArrayBase<So, Ix2>) -> Result<ArrayBase<So, Ix2>> {
-        self.solve_triangular(uplo, diag, &mut b)?;
+        self.solve_triangular_mut(uplo, diag, &mut b)?;
         Ok(b)
     }
 }

tests/triangular.rs

@@ -29,12 +29,11 @@ where
 {
     println!("a = {:?}", &a);
     println!("b = {:?}", &b);
-    let ans = b.clone();
-    let x = a.solve_triangular(uplo, Diag::NonUnit, b).unwrap();
+    let x = a.solve_triangular(uplo, Diag::NonUnit, &b).unwrap();
     println!("x = {:?}", &x);
     let b_ = a.dot(&x);
     println!("Ax = {:?}", &b_);
-    assert_close_l2!(&b_, &ans, tol);
+    assert_close_l2!(&b_, &b, tol);
 }
 
 #[test]