feat(parallel): ensemble simulation support with Rayon#87
Conversation
This commit adds a `ParallelSolve` extension trait allowing users to batch-solve a collection of initial value problems (`IVP`) using rayon when the `rayon` feature is enabled. It ensures thread-safety without restricting core library traits for non-parallel contexts. An SDE Monte Carlo example is also included. 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 helper methods to solve many IVPs (parameter sweeps/Monte Carlo) in parallel across CPU cores using
rayon.📊 Coverage
Added a test suite
tests/parallel.rsto verify that all core solvers (ODE, DAE, SDE, and DDE) can be parallelized natively withrayon.✨ Result
A new
ParallelSolveextension trait provides thepar_solvemethod for iterators ofIVPs. It avoids modifying core traits withSendorSyncbounds by relying on Rust's auto-traits, which gracefully propagate through theIVPcomponents. An example showcasing a Monte Carlo simulation for Geometric Brownian Motion using the new parallel solver is also included.PR created automatically by Jules for task 3618549310047833983 started by @Ryan-D-Gast