# Efficient implementations of the modified Gram-Schmidt orthogonalization with a non-standard inner product

## TL;DR

This paper introduces optimized implementations of the modified Gram-Schmidt orthogonalization for non-standard inner products, reducing computational cost and maintaining high accuracy, with practical benefits demonstrated through experiments.

## Contribution

It proposes $n$-matrix-vector multiplication implementations of MGS for non-standard inner products, improving efficiency and accuracy over naive methods.

## Key findings

- The HA-type implementation achieves high accuracy with error bounds.
- The HP-type implementation offers high performance with reduced computation.
- Numerical experiments show competitive advantages in cost and accuracy.

## Abstract

The modified Gram-Schmidt (MGS) orthogonalization is one of the most well-used algorithms for computing the thin QR factorization. MGS can be straightforwardly extended to a non-standard inner product with respect to a symmetric positive definite matrix $A$. For the thin QR factorization of an $m \times n$ matrix with the non-standard inner product, a naive implementation of MGS requires $2n$ matrix-vector multiplications (MV) with respect to $A$. In this paper, we propose $n$-MV implementations: a high accuracy (HA) type and a high performance (HP) type, of MGS. We also provide error bounds of the HA-type implementation. Numerical experiments and analysis indicate that the proposed implementations have competitive advantages over the naive implementation in terms of both computational cost and accuracy.