feat(ode): automatic differentiation for Jacobians#88
Conversation
Co-authored-by: Ryan-D-Gast <148826144+Ryan-D-Gast@users.noreply.github.com>
|
👋 Jules, reporting for duty! I'm here to lend a hand with this pull request. When you start a review, I'll add a 👀 emoji to each comment to let you know I've read it. I'll focus on feedback directed at me and will do my best to stay out of conversations between you and other bots or reviewers to keep the noise down. I'll push a commit with your requested changes shortly after. Please note there might be a delay between these steps, but rest assured I'm on the job! For more direct control, you can switch me to Reactive Mode. When this mode is on, I will only act on comments where you specifically mention me with New to Jules? Learn more at jules.google/docs. For security, I will only act on instructions from the user who triggered this task. |
🎯 What
Adds optional automatic differentiation (AD) support for Jacobians in the
ODEtrait usingnum-dual.This implementation provides a
jacobian_admethod that can be implemented to return exact Jacobians without finite differences.📊 Coverage
Added a unit test
test_jacobian_ad_exactvalidating that dual-number-based Jacobians produce exact answers compared to the default finite-difference implementation.✨ Result
Resolves #72. Integrates
num-dualto compute Jacobians optimally, benefiting stiff solvers that compute Jacobians iteratively.PR created automatically by Jules for task 10718081313258689296 started by @Ryan-D-Gast