Saved in:
Bibliographic Details
Main Authors: Al-Sayed, Mouhamad, Bernard, Samuel, Marzorati, Arsène, Rouzaud-Cornabas, Jonathan
Format: Preprint
Published: 2026
Subjects:
Online Access:https://arxiv.org/abs/2605.23727
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866918518666035200
author Al-Sayed, Mouhamad
Bernard, Samuel
Marzorati, Arsène
Rouzaud-Cornabas, Jonathan
author_facet Al-Sayed, Mouhamad
Bernard, Samuel
Marzorati, Arsène
Rouzaud-Cornabas, Jonathan
contents Mixed-precision methods combine low and high precision arithmetics to exploit low precision computational speed and high precision accuracy. Large ODE systems that contain many heterogeneous interactions lead to a high computational cost that could be tackled with mixed-precision solvers. We tested mixedprecision versions of the Bogacki-Shampine 3(2) Runge-Kutta pair over three benchmark systems: coupled linear oscillators, the Kuramoto model and a circadian clock model. Our study is performed in a way that can be adapted to any finite-precision format, software architecture and numerical scheme. We found that mixed-precision solvers can preserve most of the high-precision solver accuracy under a wide range of solver tolerances. Moreover, mixed-precision solver accuracy improves with system size, reaching levels equivalent to high-precision solvers in small system size. We also observed that mixed-precision arithmetic does not impact the number of evaluation in a way that balances the benefit of fast operations in low precision. Taken together, these results show that mixed-precision methods can offer significant computational speed-up at little or no loss of accuracy in large coupled ODE systems.
format Preprint
id arxiv_https___arxiv_org_abs_2605_23727
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle Mixed-Precision in adaptive Runge-Kutta method for large ODE systems
Al-Sayed, Mouhamad
Bernard, Samuel
Marzorati, Arsène
Rouzaud-Cornabas, Jonathan
Numerical Analysis
Mixed-precision methods combine low and high precision arithmetics to exploit low precision computational speed and high precision accuracy. Large ODE systems that contain many heterogeneous interactions lead to a high computational cost that could be tackled with mixed-precision solvers. We tested mixedprecision versions of the Bogacki-Shampine 3(2) Runge-Kutta pair over three benchmark systems: coupled linear oscillators, the Kuramoto model and a circadian clock model. Our study is performed in a way that can be adapted to any finite-precision format, software architecture and numerical scheme. We found that mixed-precision solvers can preserve most of the high-precision solver accuracy under a wide range of solver tolerances. Moreover, mixed-precision solver accuracy improves with system size, reaching levels equivalent to high-precision solvers in small system size. We also observed that mixed-precision arithmetic does not impact the number of evaluation in a way that balances the benefit of fast operations in low precision. Taken together, these results show that mixed-precision methods can offer significant computational speed-up at little or no loss of accuracy in large coupled ODE systems.
title Mixed-Precision in adaptive Runge-Kutta method for large ODE systems
topic Numerical Analysis
url https://arxiv.org/abs/2605.23727