rlu 0.7.0

Sparse LU Decomposition (Gilbert-Peierls)
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55

#[cfg(feature = "debug")]
use std::io::Write;

#[cfg(feature = "debug")]
const MIN_WIDTH: usize = 5;

#[cfg(feature = "debug")]
const PADDING: usize = 1;

#[cfg(feature = "debug")]
const FLOAT_CONFIG: pretty_dtoa::FmtFloatConfig = pretty_dtoa::FmtFloatConfig::default()
    .add_point_zero(false)
    .max_significant_digits(6);

#[cfg(feature = "debug")]
pub fn matrix_table<I: crate::Int, S: crate::Scalar>(m: &crate::Matrix<I, S>) -> String {
    let mut tw = tabwriter::TabWriter::new(vec![])
        .minwidth(MIN_WIDTH)
        .padding(PADDING)
        .alignment(tabwriter::Alignment::Right);

    let n = m.len();

    for r in 0..n {
        for c in 0..n {
            match m[c].iter().find(|(i, _)| i.to_index() == r) {
                None => {
                    tw.write(b"-").unwrap();
                }
                Some((_, x)) => {
                    let s = x.pretty_string(FLOAT_CONFIG);
                    tw.write(s.as_bytes()).unwrap();
                }
            }
            if c == n - 1 {
                tw.write(b"\t\n").unwrap();
            } else {
                tw.write(b"\t").unwrap();
            }
        }
    }
    String::from_utf8(tw.into_inner().unwrap()).unwrap()
}

#[cfg(feature = "debug")]
macro_rules! debug {
    ($( $args:expr ),*) => { println!( $( $args ),* ); }
}

#[cfg(not(feature = "debug"))]
macro_rules! debug {
    ($( $args:expr ),*) => {};
}

pub(crate) use debug;