Function compute_gradient_rev_with 

pub fn compute_gradient_rev_with<T, F>(
    tape: &mut Tape<T>,
    inputs: &[T],
    func: F,
) -> (T, Vec<T>)
where
    T: TapeStorage,
    F: FnOnce(&[AReal<T>]) -> AReal<T>,