Stage-Parallel Implicit Runge--Kutta methods via low-rank matrix equation corrections

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Durastante, Fabio, Mazza, Mariarosa
Format: Preprint
Published: 2025
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866915498706337792
author Durastante, Fabio
Mazza, Mariarosa
author_facet Durastante, Fabio
Mazza, Mariarosa
contents Implicit Runge--Kutta (IRK) methods are highly effective for solving stiff ordinary differential equations (ODEs) but can be computationally expensive for large-scale problems due to the need of solving coupled algebraic equations at each step. This study improves IRK efficiency by leveraging parallelism to decouple stage computations and reduce communication overhead, specifically we stably decouple a perturbed version of the stage system of equations and recover the exact solution by solving a Sylvester matrix equation with an explicitly known low-rank right-hand side. Two IRK families -- symmetric methods and collocation methods -- are analyzed, with extensions to nonlinear problems using a simplified Newton method. Implementation details, shared memory parallel code, and numerical examples, particularly for ODEs from spatially discretized PDEs, demonstrate the efficiency of the proposed IRK technique.
format Preprint
id arxiv_https___arxiv_org_abs_2505_17719
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Stage-Parallel Implicit Runge--Kutta methods via low-rank matrix equation corrections
Durastante, Fabio
Mazza, Mariarosa
Numerical Analysis
Implicit Runge--Kutta (IRK) methods are highly effective for solving stiff ordinary differential equations (ODEs) but can be computationally expensive for large-scale problems due to the need of solving coupled algebraic equations at each step. This study improves IRK efficiency by leveraging parallelism to decouple stage computations and reduce communication overhead, specifically we stably decouple a perturbed version of the stage system of equations and recover the exact solution by solving a Sylvester matrix equation with an explicitly known low-rank right-hand side. Two IRK families -- symmetric methods and collocation methods -- are analyzed, with extensions to nonlinear problems using a simplified Newton method. Implementation details, shared memory parallel code, and numerical examples, particularly for ODEs from spatially discretized PDEs, demonstrate the efficiency of the proposed IRK technique.
title Stage-Parallel Implicit Runge--Kutta methods via low-rank matrix equation corrections
topic Numerical Analysis
url https://arxiv.org/abs/2505.17719