dorsal/arxiv
View SchemaAnalysis and Efficient Sylvester-Based Implementation of a Dimension-Split ETD2RK Scheme for Multidimensional Reaction-Diffusion Equations
| Authors | Ibrahim O. Sarumi |
|---|---|
| Categories | |
| ArXiv ID | 2601.06849vv1 |
| URL | https://arxiv.org/abs/2601.06849 |
| License | http://arxiv.org/licenses/nonexclusive-distrib/1.0/ |
Abstract
We propose and analyze a second-order, dimension-split exponential time differencing Runge--Kutta scheme (ETD2RK-DS) for multidimensional reaction--diffusion equations in two and three spatial dimensions. Under mild assumptions on the nonlinear source term, we establish uniform stability bounds and prove second-order temporal convergence for the underlying dimension-split scheme. To enable efficient implementation, we employ Pad\'e approximations of the matrix exponential, converting each required matrix-exponential--vector product into the solution of a shifted linear system. A convergence analysis of the resulting Pad\'e-based ETD2RK-DS formulation is provided. We derive explicit and reproducible tensor-slicing and reshaping algorithms that realize the dimension-splitting strategy, decomposing multidimensional systems into collections of independent one-dimensional problems. This leads to a reduction of the dominant per-time-step computational cost from $\mathcal{O}(m^3)$ to $\mathcal{O}(m^2)$ in two dimensions and from $\mathcal{O}(m^5)$ to $\mathcal{O}(m^3)$ in three dimensions when compared with banded LU solvers for the unsplit problem, where $m$ denotes the number of grid points per spatial direction. Furthermore, we develop a Sylvester-equation reformulation of the resulting one-dimensional systems, enabling a highly efficient spectral implementation based on reusable eigendecompositions, matrix--vector multiplications, and Hadamard divisions. Numerical experiments in two and three dimensions, including a coupled FitzHugh--Nagumo system, confirm the second-order temporal accuracy, stability of the underlying scheme, and scalability of the proposed ETD2RK-DS framework, as well as the substantial computational advantages of the Sylvester-based implementation over classical LU-based solvers.
{
"annotation_id": "e832ed71-ee0a-4c22-bc49-86fee08bfc96",
"date_created": "2026-02-17T05:53:08.801000Z",
"date_modified": "2026-02-17T05:53:08.801000Z",
"file_hash": "9f40fdf5a63eb5e69eaeb89b00caa645804c599b833a106387e5cd302cb49471",
"private": false,
"record": {
"abstract": "We propose and analyze a second-order, dimension-split exponential time differencing Runge--Kutta scheme (ETD2RK-DS) for multidimensional reaction--diffusion equations in two and three spatial dimensions. Under mild assumptions on the nonlinear source term, we establish uniform stability bounds and prove second-order temporal convergence for the underlying dimension-split scheme.\n To enable efficient implementation, we employ Pad\\\u0027e approximations of the matrix exponential, converting each required matrix-exponential--vector product into the solution of a shifted linear system. A convergence analysis of the resulting Pad\\\u0027e-based ETD2RK-DS formulation is provided. We derive explicit and reproducible tensor-slicing and reshaping algorithms that realize the dimension-splitting strategy, decomposing multidimensional systems into collections of independent one-dimensional problems. This leads to a reduction of the dominant per-time-step computational cost from $\\mathcal{O}(m^3)$ to $\\mathcal{O}(m^2)$ in two dimensions and from $\\mathcal{O}(m^5)$ to $\\mathcal{O}(m^3)$ in three dimensions when compared with banded LU solvers for the unsplit problem, where $m$ denotes the number of grid points per spatial direction.\n Furthermore, we develop a Sylvester-equation reformulation of the resulting one-dimensional systems, enabling a highly efficient spectral implementation based on reusable eigendecompositions, matrix--vector multiplications, and Hadamard divisions. Numerical experiments in two and three dimensions, including a coupled FitzHugh--Nagumo system, confirm the second-order temporal accuracy, stability of the underlying scheme, and scalability of the proposed ETD2RK-DS framework, as well as the substantial computational advantages of the Sylvester-based implementation over classical LU-based solvers.",
"arxiv_id": "2601.06849",
"authors": [
"Ibrahim O. Sarumi"
],
"categories": [
"math.NA",
"cs.NA"
],
"license": "http://arxiv.org/licenses/nonexclusive-distrib/1.0/",
"title": "Analysis and Efficient Sylvester-Based Implementation of a Dimension-Split ETD2RK Scheme for Multidimensional Reaction-Diffusion Equations",
"url": "https://arxiv.org/abs/2601.06849",
"version": "v1"
},
"schema_id": "dorsal/arxiv",
"source": {
"execution_id": "143567dc-0368-4ea9-9bc7-8dfff3ca0bfb",
"id": "arXiv Dataset",
"type": "Model",
"variant": "snapshot-2026-01-17",
"version": "0.1.0"
},
"user_id": 1000002
}