The effective use of BLAS interface for implementation of finite-element ADER-DG and finite-volume ADER-WENO methods
Fuente:
arXiv
Saved in:
| Main Author: | |
|---|---|
| Format: | Preprint |
| Published: |
2024
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866915517823975424 |
|---|---|
| author | Popov, I. S. |
| author_facet | Popov, I. S. |
| contents | Numerical methods of the ADER family, in particular finite-element ADER-DG and finite-volume ADER-WENO methods, are among the most accurate numerical methods for solving quasilinear hyperbolic PDE systems. The internal structure of ADER-DG and ADER-WENO numerical methods contains a large number of basic linear algebra operations related to matrix multiplications. The main interface of software libraries for matrix multiplications for high-performance computing is BLAS. An effective method for integration the standard functions of the BLAS interface into the implementation of these numerical methods is presented. The calculated matrices are small matrices; and this allows to use effectively JIT technologies. The proposed approach immediately operates on AoS, which allows to efficiently calculate flux, source and non-conservative terms without transposition. The obtained computational costs demonstrated that the effective implementation, based on the use of the JIT functions of the BLAS, outperformed both the implementation based on the general BLAS functions and the vanilla implementations by several orders of magnitude. The complexity of developing an implementation based on the proposed approach does not exceed the complexity of developing a vanilla implementation. Performance analysis using roofline partly explains the observed features of the decreasing of computational costs. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2409_12483 |
| institution | arXiv |
| publishDate | 2024 |
| record_format | arxiv |
| spellingShingle | The effective use of BLAS interface for implementation of finite-element ADER-DG and finite-volume ADER-WENO methods Popov, I. S. Fluid Dynamics Numerical Analysis Analysis of PDEs Applied Physics Computational Physics 65M60, 76M10, 76M12, 35L60, 35L65, 35L45, 35L67 Numerical methods of the ADER family, in particular finite-element ADER-DG and finite-volume ADER-WENO methods, are among the most accurate numerical methods for solving quasilinear hyperbolic PDE systems. The internal structure of ADER-DG and ADER-WENO numerical methods contains a large number of basic linear algebra operations related to matrix multiplications. The main interface of software libraries for matrix multiplications for high-performance computing is BLAS. An effective method for integration the standard functions of the BLAS interface into the implementation of these numerical methods is presented. The calculated matrices are small matrices; and this allows to use effectively JIT technologies. The proposed approach immediately operates on AoS, which allows to efficiently calculate flux, source and non-conservative terms without transposition. The obtained computational costs demonstrated that the effective implementation, based on the use of the JIT functions of the BLAS, outperformed both the implementation based on the general BLAS functions and the vanilla implementations by several orders of magnitude. The complexity of developing an implementation based on the proposed approach does not exceed the complexity of developing a vanilla implementation. Performance analysis using roofline partly explains the observed features of the decreasing of computational costs. |
| title | The effective use of BLAS interface for implementation of finite-element ADER-DG and finite-volume ADER-WENO methods |
| topic | Fluid Dynamics Numerical Analysis Analysis of PDEs Applied Physics Computational Physics 65M60, 76M10, 76M12, 35L60, 35L65, 35L45, 35L67 |
| url | https://arxiv.org/abs/2409.12483 |