Recursive vectorized computation of the vector $p$-norm
Abstract
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 -norm is also presented.
Cite
@article{arxiv.2509.06220,
title = {Recursive vectorized computation of the vector $p$-norm},
author = {Vedran Novaković},
journal= {arXiv preprint arXiv:2509.06220},
year = {2026}
}
Comments
Preprint of an article published in Parallel Processing Letters, Vol. 36, No. 01n02, 2026, Article 2650002. Copyright World Scientific Publishing Company