conspire 0.6.0

The Rust interface to conspire.
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

macro_rules! test_explicit_fixed_step {
    ($integration: expr) => {
        use crate::math::{
            Tensor, assert_eq, assert_eq_within,
            integrate::{
                FixedStep,
                ode::explicit::test::test_explicit,
                test::{LENGTH, zero_to_one},
            },
        };
        const TIME_STEP: Scalar = 0.1;
        const TOLERANCE: Scalar = TIME_STEP;
        test_explicit!($integration);
        #[test]
        fn dxdt_eq_neg_x() -> Result<(), TestError> {
            assert_eq(&$integration.dt(), &TIME_STEP)?;
            let (time, solution, function): (Vector, Vector, _) =
                $integration.integrate(|_: Scalar, x: &Scalar| Ok(-x), &[0.0, 0.8], 1.0)?;
            time.iter()
                .zip(solution.iter().zip(function.iter()))
                .try_for_each(|(t, (y, f))| {
                    assert_eq_within(y, &(-t).exp(), TOLERANCE, TOLERANCE)?;
                    assert_eq_within(f, &-y, TOLERANCE, TOLERANCE)
                })
        }
        #[test]
        fn eval_times() -> Result<(), TestError> {
            assert_eq(&$integration.dt(), &TIME_STEP)?;
            let (time, solution, function): (Vector, Vector, _) = $integration.integrate(
                |_: Scalar, x: &Scalar| Ok(-x),
                &zero_to_one::<LENGTH>(),
                1.0,
            )?;
            time.iter()
                .zip(solution.iter().zip(function.iter()))
                .try_for_each(|(t, (y, f))| {
                    assert_eq_within(y, &(-t).exp(), TOLERANCE, TOLERANCE)?;
                    assert_eq_within(f, &-y, TOLERANCE, TOLERANCE)
                })
        }
    };
}
pub(crate) use test_explicit_fixed_step;