A comparison of mixed precision iterative refinement approaches for least-squares problems

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Carson, Erin, Daužickaitė, Ieva
Format: Preprint
Published: 2024
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866909458070765568
author Carson, Erin
Daužickaitė, Ieva
author_facet Carson, Erin
Daužickaitė, Ieva
contents Various approaches to iterative refinement (IR) for least-squares problems have been proposed in the literature and it may not be clear which approach is suitable for a given problem. We consider three approaches to IR for least-squares problems when two precisions are used and review their theoretical guarantees, known shortcomings and when the method can be expected to recognize that the correct solution has been found, and extend uniform precision analysis for an IR approach based on the semi-normal equations to the two-precision case. We focus on the situation where it is desired to refine the solution to the working precision level. It is shown that the IR methods exhibit different sensitivities to the conditioning of the problem and the size of the least-squares residual, which should be taken into account when choosing the IR approach. We also discuss a new approach that is based on solving multiple least-squares problems.
format Preprint
id arxiv_https___arxiv_org_abs_2405_18363
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle A comparison of mixed precision iterative refinement approaches for least-squares problems
Carson, Erin
Daužickaitė, Ieva
Numerical Analysis
Various approaches to iterative refinement (IR) for least-squares problems have been proposed in the literature and it may not be clear which approach is suitable for a given problem. We consider three approaches to IR for least-squares problems when two precisions are used and review their theoretical guarantees, known shortcomings and when the method can be expected to recognize that the correct solution has been found, and extend uniform precision analysis for an IR approach based on the semi-normal equations to the two-precision case. We focus on the situation where it is desired to refine the solution to the working precision level. It is shown that the IR methods exhibit different sensitivities to the conditioning of the problem and the size of the least-squares residual, which should be taken into account when choosing the IR approach. We also discuss a new approach that is based on solving multiple least-squares problems.
title A comparison of mixed precision iterative refinement approaches for least-squares problems
topic Numerical Analysis
url https://arxiv.org/abs/2405.18363