Add test for Householder::append

Author: termoshtt

src/krylov/householder.rs

@@ -80,3 +80,26 @@ impl<A: Scalar + Lapack> Orthogonalizer for Householder<A> {
         a
     }
 }
+
+#[cfg(test)]
+mod tests {
+    use super::*;
+    use crate::assert::*;
+
+    #[test]
+    fn householder_append() {
+        let mut householder = Householder::new(3);
+        let coef = householder.append(array![0.0, 1.0, 0.0], 1e-9).unwrap();
+        close_l2(&coef, &array![1.0], 1e-9).unwrap();
+
+        let coef = householder.append(array![1.0, 1.0, 0.0], 1e-9).unwrap();
+        close_l2(&coef, &array![1.0, 1.0], 1e-9).unwrap();
+
+        assert!(householder.append(array![1.0, 2.0, 0.0], 1e-9).is_err());
+
+        if let Err(coef) = householder.append(array![1.0, 2.0, 0.0], 1e-9) {
+            close_l2(&coef, &array![2.0, 1.0, 0.0], 1e-9).unwrap();
+        }
+    }
+
+}