Finite-budget Newton iteration for the nonsingular polar factor

Input a stored square complex matrix $A\in {\mathbb{C}}^{n\times n}$, $n\ge 1$, which is assumed nonsingular, a finite iteration budget $N\ge 0$, and finite Householder-QR controls including a nonnegative integer candidate-operation budget $K_H$. Set the candidate $X_0=A$. For $k=0,\dots ,N-1$:

1. Call Finite Householder QR solve for a Padé denominator with explicit statuses on $X_k^{\ast }Y=I$, passing $K_H$ to its instrumented kernel dispatcher. If it does not return verified-solution, return its failure status with the last candidate.
2. Compute $X_{k+1}=\mathrm{fl}((X_k+Y)\mathrm{/}2)$ componentwise in the operationwise model of Operationwise real and complex floating point arithmetic with rounded square roots and fused multiply-add. A range failure returns range-failure with the last candidate.

After $N$ completed steps, return $X_N$ with status budget-exhausted and the number of completed steps. The exact Newton sequence is analyzed separately; this procedure makes no forward-error or floating-point convergence claim.