ida 0.1.1

A pure Rust port of the Implicit Differential-Algebraic solver from the Sundials suite.
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
56
57
58
59
60
61

use super::*;

#[test]
fn test_get_solution() {
    // --- IDAGetSolution Before:
    let t = 3623118336.24244;
    let hh = 857870592.1885694;
    let tn = 3623118336.24244;
    let kused = 4;
    let hused = 428935296.0942847;

    #[rustfmt::skip]
        let ida_phi = array![
            [ 5.716499633245077e-07, 2.286601144610028e-12, 0.9999994283477499, ],
            [ -7.779233860067279e-08, -3.111697299545603e-13, 7.779264957586927e-08, ],
            [ 2.339417551980491e-08, 9.35768837422748e-14, -2.33942692332846e-08, ],
            [ -9.503346432581604e-09, -3.801349575270522e-14, 9.503383895634436e-09, ],
            [ 7.768373161310588e-09, 3.107357755532867e-14, -7.768407422476745e-09, ],
            [ -2.242367216194777e-10, -8.970915966733762e-16, 2.242247401239887e-10, ],
        ];
    let ida_psi = array![
        428935296.0942847,
        857870592.1885694,
        1072338240.235712,
        1286805888.282854,
        1501273536.329997,
        26020582.4876316,
    ];

    //--- IDAGetSolution After:
    let yret_expect = array![
        5.716499633245077e-07,
        2.286601144610028e-12,
        0.9999994283477499,
    ];
    let ypret_expect = array![
        -1.569167478317552e-16,
        -6.276676917262037e-22,
        1.569173718962504e-16,
    ];

    let problem = Dummy {};
    let mut ida: Ida<_, linear::Dense<_>, nonlinear::Newton<_>, _> = Ida::new(
        problem,
        array![0., 0., 0.],
        array![0., 0., 0.],
        TolControlSS::new(1e-4, 1e-4),
    );

    ida.ida_hh = hh;
    ida.nlp.ida_tn = tn;
    ida.ida_kused = kused;
    ida.ida_hused = hused;
    ida.ida_phi.assign(&ida_phi);
    ida.ida_psi.assign(&ida_psi);

    ida.get_solution(t).unwrap();

    assert_nearly_eq!(ida.get_yy(), yret_expect);
    assert_nearly_eq!(ida.get_yp(), ypret_expect);
}