English

Iterative Refinement with Low-Precision Posits

Numerical Analysis 2024-08-28 v2 Numerical Analysis

Abstract

This research investigates using a mixed-precision iterative refinement method using posit numbers instead of the standard IEEE floating-point format. The method is applied to solve a general linear system represented by the equation Ax=bAx = b, where AA is a large sparse matrix. Various scaling techniques, such as row and column equilibration, map the matrix entries to higher-density regions of machine numbers before performing the O(n3)O(n^3) factorization operation. Low-precision LU factorization followed by forward/backward substitution provides an initial estimate. The results demonstrate that a 16-bit posit configuration combined with equilibration produces accuracy comparable to IEEE half-precision (fp16), indicating a potential for achieving a balance between efficiency and accuracy.

Keywords

Cite

@article{arxiv.2408.13400,
  title  = {Iterative Refinement with Low-Precision Posits},
  author = {James Quinlan and E. Theodore L. Omtzigt},
  journal= {arXiv preprint arXiv:2408.13400},
  year   = {2024}
}

Comments

preprint CoNGA'24

R2 v1 2026-06-28T18:22:40.028Z