num-dual 0.13.6

Generalized (hyper) dual numbers for the calculation of exact (partial) derivatives
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

#![allow(non_snake_case)]
use pyo3::prelude::*;
use pyo3::wrap_pyfunction;

mod dual;
mod dual2;
mod dual3;
mod hyperdual;
mod hyperhyperdual;

use dual::{first_derivative, gradient, jacobian};
use dual2::{hessian, second_derivative};
use dual3::third_derivative;
use hyperdual::{partial_hessian, second_partial_derivative};
use hyperhyperdual::{third_partial_derivative, third_partial_derivative_vec};

pub use dual::{
    PyDual64, PyDual64_2, PyDual64_3, PyDual64_4, PyDual64_5, PyDual64_6, PyDual64_7, PyDual64_8,
    PyDual64_9, PyDual64_10, PyDual64_11, PyDual64_12, PyDual64_13, PyDual64_14, PyDual64_15,
    PyDual64_16, PyDual64Dyn,
};
pub use dual2::{
    PyDual2_64, PyDual2_64_1, PyDual2_64_2, PyDual2_64_3, PyDual2_64_4, PyDual2_64_5, PyDual2_64_6,
    PyDual2_64_7, PyDual2_64_8, PyDual2_64_9, PyDual2_64_10, PyDual2_64_11, PyDual2_64_12,
    PyDual2_64_13, PyDual2_64_14, PyDual2_64_15, PyDual2_64_16, PyDual2_64Dyn, PyDual2Dual64,
};
pub use dual3::{PyDual3_64, PyDual3Dual64};
pub use hyperdual::{
    PyHyperDual64, PyHyperDual64_1_2, PyHyperDual64_1_3, PyHyperDual64_1_4, PyHyperDual64_1_5,
    PyHyperDual64_2_1, PyHyperDual64_2_2, PyHyperDual64_2_3, PyHyperDual64_2_4, PyHyperDual64_2_5,
    PyHyperDual64_3_1, PyHyperDual64_3_2, PyHyperDual64_3_3, PyHyperDual64_3_4, PyHyperDual64_3_5,
    PyHyperDual64_4_1, PyHyperDual64_4_2, PyHyperDual64_4_3, PyHyperDual64_4_4, PyHyperDual64_4_5,
    PyHyperDual64_5_1, PyHyperDual64_5_2, PyHyperDual64_5_3, PyHyperDual64_5_4, PyHyperDual64_5_5,
    PyHyperDual64Dyn, PyHyperDualDual64,
};
pub use hyperhyperdual::PyHyperHyperDual64;

#[pymodule]
pub fn num_dual(_py: Python, m: &Bound<'_, PyModule>) -> PyResult<()> {
    m.add("__version__", env!("CARGO_PKG_VERSION"))?;
    m.add_class::<PyDual64>()?;
    m.add_class::<PyHyperDual64>()?;
    m.add_class::<PyDual2_64>()?;
    m.add_class::<PyDual3_64>()?;
    m.add_class::<PyHyperDualDual64>()?;
    m.add_class::<PyDual2Dual64>()?;
    m.add_class::<PyDual3Dual64>()?;
    m.add_class::<PyHyperHyperDual64>()?;
    m.add_function(wrap_pyfunction!(first_derivative, m)?)?;
    m.add_function(wrap_pyfunction!(gradient, m)?)?;
    m.add_function(wrap_pyfunction!(jacobian, m)?)?;
    m.add_function(wrap_pyfunction!(second_derivative, m)?)?;
    m.add_function(wrap_pyfunction!(hessian, m)?)?;
    m.add_function(wrap_pyfunction!(third_derivative, m)?)?;
    m.add_function(wrap_pyfunction!(second_partial_derivative, m)?)?;
    m.add_function(wrap_pyfunction!(partial_hessian, m)?)?;
    m.add_function(wrap_pyfunction!(third_partial_derivative, m)?)?;
    m.add_function(wrap_pyfunction!(third_partial_derivative_vec, m)?)?;
    Ok(())
}