Update test for non-square b

Author: termoshtt

tests/triangular.rs

@@ -34,9 +34,6 @@ fn test2d<A, Sa, Sb, Tol>(uplo: UPLO, a: ArrayBase<Sa, Ix2>, b: ArrayBase<Sb, Ix
     println!("x = {:?}", &x);
     let b_ = a.dot(&x);
     println!("Ax = {:?}", &b_);
-    println!("A^Tx = {:?}", a.t().dot(&x));
-    println!("Ax^T = {:?}", a.dot(&x.t()));
-    println!("(Ax^T)^T = {:?}", a.dot(&x.t()).t());
     assert_close_l2!(&b_, &ans, tol);
 }
 
@@ -74,64 +71,56 @@ fn triangular_1d_upper_t() {
 
 #[test]
 fn triangular_2d_upper() {
-    let n = 3;
-    let b: Array2<f64> = random_square(n);
-    let a: Array2<f64> = random_square(n).into_triangular(UPLO::Upper);
+    let b: Array2<f64> = random_matrix(3, 4);
+    let a: Array2<f64> = random_square(3).into_triangular(UPLO::Upper);
     test2d(UPLO::Upper, a, b, 1e-7);
 }
 
 #[test]
 fn triangular_2d_lower() {
-    let n = 3;
-    let b: Array2<f64> = random_square(n);
-    let a: Array2<f64> = random_square(n).into_triangular(UPLO::Lower);
+    let b: Array2<f64> = random_matrix(3, 4);
+    let a: Array2<f64> = random_square(3).into_triangular(UPLO::Lower);
     test2d(UPLO::Lower, a, b, 1e-7);
 }
 
 #[test]
 fn triangular_2d_lower_t() {
-    let n = 3;
-    let b: Array2<f64> = random_square(n);
-    let a: Array2<f64> = random_square(n).into_triangular(UPLO::Lower).reversed_axes();
+    let b: Array2<f64> = random_matrix(3, 4);
+    let a: Array2<f64> = random_square(3).into_triangular(UPLO::Lower).reversed_axes();
     test2d(UPLO::Upper, a, b, 1e-7);
 }
 
 #[test]
 fn triangular_2d_upper_t() {
-    let n = 3;
-    let b: Array2<f64> = random_square(n);
-    let a: Array2<f64> = random_square(n).into_triangular(UPLO::Upper).reversed_axes();
+    let b: Array2<f64> = random_matrix(3, 4);
+    let a: Array2<f64> = random_square(3).into_triangular(UPLO::Upper).reversed_axes();
     test2d(UPLO::Lower, a, b, 1e-7);
 }
 
 #[test]
 fn triangular_2d_upper_bt() {
-    let n = 3;
-    let b: Array2<f64> = random_square(n).reversed_axes();
-    let a: Array2<f64> = random_square(n).into_triangular(UPLO::Upper);
+    let b: Array2<f64> = random_matrix(4, 3).reversed_axes();
+    let a: Array2<f64> = random_square(3).into_triangular(UPLO::Upper);
     test2d(UPLO::Upper, a, b, 1e-7);
 }
 
 #[test]
 fn triangular_2d_lower_bt() {
-    let n = 3;
-    let b: Array2<f64> = random_square(n).reversed_axes();
-    let a: Array2<f64> = random_square(n).into_triangular(UPLO::Lower);
+    let b: Array2<f64> = random_matrix(4, 3).reversed_axes();
+    let a: Array2<f64> = random_square(3).into_triangular(UPLO::Lower);
     test2d(UPLO::Lower, a, b, 1e-7);
 }
 
 #[test]
 fn triangular_2d_lower_t_bt() {
-    let n = 3;
-    let b: Array2<f64> = random_square(n).reversed_axes();
-    let a: Array2<f64> = random_square(n).into_triangular(UPLO::Lower).reversed_axes();
+    let b: Array2<f64> = random_matrix(4, 3).reversed_axes();
+    let a: Array2<f64> = random_square(3).into_triangular(UPLO::Lower).reversed_axes();
     test2d(UPLO::Upper, a, b, 1e-7);
 }
 
 #[test]
 fn triangular_2d_upper_t_bt() {
-    let n = 3;
-    let b: Array2<f64> = random_square(n).reversed_axes();
-    let a: Array2<f64> = random_square(n).into_triangular(UPLO::Upper).reversed_axes();
+    let b: Array2<f64> = random_matrix(4, 3).reversed_axes();
+    let a: Array2<f64> = random_square(3).into_triangular(UPLO::Upper).reversed_axes();
     test2d(UPLO::Lower, a, b, 1e-7);
 }