Function check_matrix

examples/clippy_verify.rs:169–205

()

}

fn check_matrix() {
    // py_matrix (api.rs:54-55)
    let m = py_matrix(vec![
        vec![1.0_f64, 2.0, 3.0, 4.0],
        vec![5.0, 6.0, 7.0, 8.0],
        vec![9.0, 10.0, 11.0, 12.0],
    ]);
    print_kv("py_matrix 3x4 hash", hash_matrix(&m));
    // change_shape (matrix.rs:897, 905)
    let m_col = m.change_shape();
    print_kv("change_shape Row->Col hash", hash_matrix(&m_col));
    let m_back = m_col.change_shape();
    print_kv("change_shape Col->Row hash", hash_matrix(&m_back));
    // diag (matrix.rs:1105)
    let sq = ml_matrix("1 2 3; 4 5 6; 7 8 9");
    let d = sq.diag();
    print_kv("diag 3x3", hash_f64s(&d));
    // to_vec (matrix.rs:1184)
    let v_of_v = sq.to_vec();
    let flat: Vec<f64> = v_of_v.into_iter().flatten().collect();
    print_kv("to_vec hash", hash_f64s(&flat));
    // col_reduce / row_reduce (matrix.rs:2863, 2874)
    let cr = sq.col_reduce(|c| c.iter().sum::<f64>());
    let rr = sq.row_reduce(|c| c.iter().sum::<f64>());
    print_kv("col_reduce sum", hash_f64s(&cr));
    print_kv("row_reduce sum", hash_f64s(&rr));
    // is_symmetric (matrix.rs:3660, 3683)
    let sym = ml_matrix("1 2 3; 2 5 6; 3 6 9");
    let nonsym = ml_matrix("1 2 3; 4 5 6; 7 8 9");
    print_kv("is_symmetric sym", sym.is_symmetric());
    print_kv("is_symmetric nonsym", nonsym.is_symmetric());
    // matrix.rs:764-765 - inspect; this site is inside `py_matrix`
    // creation path already exercised above. Cover additionally with a
    // non-square shape so both row/col loops fire.
    let rect = py_matrix(vec![vec![1.0_f64, 2.0], vec![3.0, 4.0], vec![5.0, 6.0]]);
    print_kv("py_matrix 3x2 hash", hash_matrix(&rect));
}

fn check_polynomial() {
    let p = poly(vec![1.0, -2.0, 0.5, 3.0]); // 1 - 2x + 0.5x^2 + 3x^3