Recursive vectorized computation of the vector $p$-norm
Fuente:
arXiv
Saved in:
| Main Author: | |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866914556072165376 |
|---|---|
| author | Novaković, Vedran |
| author_facet | Novaković, Vedran |
| contents | Recursive algorithms for computing the Frobenius norm of a real array are proposed, based on hypot, a hypotenuse function. Comparing their relative accuracy bounds with those of the BLAS routine DNRM2 it is shown that the proposed algorithms could in many cases be significantly more accurate. The scalar recursive algorithms are vectorized with the Intel's vector instructions to achieve performance comparable to DNRM2, and are further parallelized with OpenCilk. Some scalar algorithms are unconditionally bitwise reproducible, while the reproducibility of the vector ones depends on the vector width. A modification of the proposed algorithms to compute the vector $p$-norm is also presented. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2509_06220 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Recursive vectorized computation of the vector $p$-norm Novaković, Vedran Numerical Analysis 65F35 (Primary) 65Y05, 65G50 (Secondary) G.1.3; G.4 Recursive algorithms for computing the Frobenius norm of a real array are proposed, based on hypot, a hypotenuse function. Comparing their relative accuracy bounds with those of the BLAS routine DNRM2 it is shown that the proposed algorithms could in many cases be significantly more accurate. The scalar recursive algorithms are vectorized with the Intel's vector instructions to achieve performance comparable to DNRM2, and are further parallelized with OpenCilk. Some scalar algorithms are unconditionally bitwise reproducible, while the reproducibility of the vector ones depends on the vector width. A modification of the proposed algorithms to compute the vector $p$-norm is also presented. |
| title | Recursive vectorized computation of the vector $p$-norm |
| topic | Numerical Analysis 65F35 (Primary) 65Y05, 65G50 (Secondary) G.1.3; G.4 |
| url | https://arxiv.org/abs/2509.06220 |