spsolve 0.1.0

Solve sparse systems of linear equations
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

use std::time::Instant;

#[cfg(feature = "cpuprofiler")]
use cpuprofiler::PROFILER;

use spsolve::Solver;

fn main() {
    let trans = false;
    // let (n, a_p, a_i, a_x) = spsolve::matrix::activsg2000_bbus(!trans);
    let (n, a_p, a_i, a_x) = spsolve::matrix::activsg10k_bbus(!trans);
    // let (n, a_p, a_i, a_x) = spsolve::matrix::activsg25k_bbus(!trans);
    let mut b: Vec<f64> = (0..n).map(|i| 1.0 + i as f64 / n as f64).collect();

    let solver = spsolve::rlu::RLU::default();

    let t0 = Instant::now();

    #[cfg(feature = "cpuprofiler")]
    PROFILER.lock().unwrap().start("./spsolve.profile").unwrap();

    solver.solve(n, &a_i, &a_p, &a_x, &mut b, trans).unwrap();

    #[cfg(feature = "cpuprofiler")]
    PROFILER.lock().unwrap().stop().unwrap();

    println!("t = {:#?}", t0.elapsed());
}