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
//! 
//! This crate provides a library for performing automatic differentiation.
//!
//! # Examples
//!
//! The following example differentiates a 1D function defined by a closure.
//!
//! ```rust
//! # extern crate autodiff;
//! # use autodiff::*;
//! # fn main() {
//!     // Define a function `f(x) = e^{-0.5*x^2}`
//!     let f = |x: Num| (-x * x / Num::cst(2.0)).exp();
//!
//!     // Differentiate `f` at zero.
//!     println!("{}", diff(f, 0.0)); // prints `0`
//! #   assert_eq!(diff(f, 0.0), 0.0);
//! # }
//! ```
//!
//! To compute the gradient of a function, use the function `grad` as follows:
//!
//! ```rust
//! # extern crate autodiff;
//! # use autodiff::*;
//! # fn main() {
//!     // Define a function `f(x,y) = x*y^2`
//!     let f = |x: &[Num]| x[0] * x[1] * x[1];
//!
//!     // Differentiate `f` at `(1,2)`.
//!     let g = grad(f, &vec![1.0, 2.0]);
//!     println!("({}, {})", g[0], g[1]); // prints `(4, 4)`
//! #   assert_eq!(g, vec![4.0, 4.0]);
//! # }
//! ```
//!
//! Compute a specific derivative of a multi-variable function:
//! 
//! ```rust
//! # extern crate autodiff;
//! # use autodiff::*;
//! # fn main() {
//!     // Define a function `f(x,y) = x*y^2`.
//!     let f = |v: &[Num]| v[0] * v[1] * v[1];
//! 
//!     // Differentiate `f` at `(1,2)` with respect to `x` (the first unknown) only.
//!     let v = vec![
//!         Num::var(1.0), // Create a variable.
//!         Num::cst(2.0), // Create a constant.
//!     ];
//!     println!("{}", f(&v).deriv()); // prints `4`
//! #   assert_eq!(f(&v).deriv(), 4.0);
//! # }
//! ```
extern crate num_traits;

pub mod forward_autodiff;

pub use crate::forward_autodiff::*;

// Re-export useful traits for performing computations.
pub use num_traits::{Float, FloatConst, NumCast, One, ToPrimitive, Zero};