peroxide 0.41.2

Rust comprehensive scientific computation library contains linear algebra, numerical analysis, statistics and machine learning tools with familiar syntax
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

#[macro_use]
extern crate peroxide;
use peroxide::fuga::*;

#[test]
fn test_jacobian() {
    let x = c!(1, 0);
    let j = jacobian(f, &x);
    assert_eq!(j, ml_matrix("0 1; 5 1"));
}

/// Test function
///
/// # Function
/// $$\begin{pmatrix} x^2 y \\\ 5x+ \sin y \end{pmatrix}$$
///
/// # Jacobian
/// $$\begin{pmatrix} 2xy & x^2 \\\ 5 & \cos y \end{pmatrix}$$
// Signature must stay `&Vec<AD>` to satisfy `jacobian`'s
// `Fn(&Vec<AD>) -> Vec<AD>` trait bound.
#[allow(clippy::ptr_arg)]
fn f(xs: &Vec<AD>) -> Vec<AD> {
    let x = xs[0];
    let y = xs[1];
    vec![x.powi(2) * y, 5f64 * x + y.sin()]
}