We develop a mixed-precision iterative refinement framework for solving low-rank Lyapunov matrix equations AX + XAT + W =0, where W=LLT or W=LSLT. Via rounding error analysis of the algorithms we derive sufficient conditions for the attainable normwise residuals in different precision settings and show how the algorithmic parameters should be chosen. Using the sign function Newton iteration as the solver, we show that reduced precisions, such as the half precision, can be used as the solver precision (with unit roundoff uₛ) to accelerate the solution of Lyapunov equations of condition number up to 1/uₛ without compromising its quality.
Benner et al. (Thu,) studied this question.