Skip to content

Commit

Permalink
WIP: add second derivative test (#4)
Browse files Browse the repository at this point in the history
  • Loading branch information
jedbrown committed Apr 14, 2024
1 parent 7df4710 commit 1625aa8
Show file tree
Hide file tree
Showing 2 changed files with 31 additions and 1 deletion.
2 changes: 1 addition & 1 deletion samples/tests/mod.rs
Original file line number Diff line number Diff line change
Expand Up @@ -2,7 +2,7 @@

mod forward;
mod reverse;

mod second;
mod higher;
mod neohookean;
mod traits;
30 changes: 30 additions & 0 deletions samples/tests/second/mod.rs
Original file line number Diff line number Diff line change
@@ -0,0 +1,30 @@
samples::test! {
forward_of_reverse;
// ANCHOR: forward_of_reverse
#[autodiff(df, Reverse, Duplicated, Active)]
fn f(x: &[f32; 2]) -> f32 {
x[0] * x[0] + x[1] * x[0]
}

#[autodiff(dg, Forward, Dual, Dual)]
fn g(x: &[f32; 2], gx: &mut [f32; 2]) {
let mut bx = [0.0, 0.0];
df(x, &mut bx, 1.0);
*gx = bx;
}

fn main() {
let x = [2.0, 3.0];
let fx = f(&x);
assert_eq!(10.0, fx);

let mut gx = [0.0, 0.0];
g(&x, &mut gx);
assert_eq!([7.0, 2.0], gx);

let dx = [1.0, 0.0];
let mut dgx = [0.0, 0.0];
dg(&x, &dx, &mut gx, &mut dgx);
}
// ANCHOR_END: forward_of_reverse
}

0 comments on commit 1625aa8

Please sign in to comment.