numdiff 0.7.3

Numerical differentiation via forward-mode automatic differentiation and finite difference approximations.
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
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93

use crate::HyperDual;
use linalg_traits::{Scalar, Vector};

/// Trait to create a vector of hyper-dual numbers.
pub trait HyperDualVector<S, V>
where
    S: Scalar,
    V: Vector<S>,
{
    /// Convert this vector of scalars to a vector of hyper-dual numbers.
    ///
    /// # Returns
    ///
    /// A copy of this vector with each element converted to a hyper-dual number (with `ε1`, `ε2`,
    /// and `ε1ε2` coefficients equal to `0.0`).
    fn to_hyper_dual_vector(self) -> V::VectorT<HyperDual>;
}

impl<S, V> HyperDualVector<S, V> for V
where
    S: Scalar,
    V: Vector<S>,
{
    fn to_hyper_dual_vector(self) -> V::VectorT<HyperDual> {
        let mut vec_hyper_dual = V::VectorT::new_with_length(self.len());
        for i in 0..self.len() {
            vec_hyper_dual.vset(
                i,
                HyperDual::new(self.vget(i).to_f64().unwrap(), 0.0, 0.0, 0.0),
            );
        }
        vec_hyper_dual
    }
}

#[cfg(test)]
mod tests {
    use super::*;
    use nalgebra::{SVector, dvector};
    use ndarray::array;

    #[test]
    fn test_vec() {
        let vec = vec![1.0, 2.0, 3.0];
        assert_eq!(
            vec.to_hyper_dual_vector(),
            vec![
                HyperDual::new(1.0, 0.0, 0.0, 0.0),
                HyperDual::new(2.0, 0.0, 0.0, 0.0),
                HyperDual::new(3.0, 0.0, 0.0, 0.0)
            ]
        );
    }

    #[test]
    fn test_nalgebra_dvector() {
        let vec = dvector![1.0, 2.0, 3.0];
        assert_eq!(
            vec.to_hyper_dual_vector(),
            dvector![
                HyperDual::new(1.0, 0.0, 0.0, 0.0),
                HyperDual::new(2.0, 0.0, 0.0, 0.0),
                HyperDual::new(3.0, 0.0, 0.0, 0.0)
            ]
        );
    }

    #[test]
    fn test_nalgebra_svector() {
        let vec = SVector::<f64, 3>::from_row_slice(&[1.0, 2.0, 3.0]);
        assert_eq!(
            vec.to_hyper_dual_vector(),
            SVector::<HyperDual, 3>::from_row_slice(&[
                HyperDual::new(1.0, 0.0, 0.0, 0.0),
                HyperDual::new(2.0, 0.0, 0.0, 0.0),
                HyperDual::new(3.0, 0.0, 0.0, 0.0)
            ])
        );
    }

    #[test]
    fn test_ndarray_array1() {
        let vec = array![1.0, 2.0, 3.0];
        assert_eq!(
            vec.to_hyper_dual_vector(),
            vec![
                HyperDual::new(1.0, 0.0, 0.0, 0.0),
                HyperDual::new(2.0, 0.0, 0.0, 0.0),
                HyperDual::new(3.0, 0.0, 0.0, 0.0)
            ]
        );
    }
}