2025/08/23 by Leplat, Valentin, Ahmadi-Asl, Salman, Pan, JunJun +1
#FOS: Mathematics #Numerical Analysis (math.NA)
paper · doi:10.48550/arxiv.2508.16979
We develop quaternion--native iterative methods for computing the Moore--Penrose (MP) pseudoinverse of quaternion matrices and analyze their convergence. Our starting point is a damped Newton--Schulz (NS) iteration tailored to noncommutativity: we enforce the appropriate left/right identities for rectangular inputs and prove convergence directly in ℍ under a simple spectral scaling. We then derive higher--order (hyperpower) NS schemes with exact residual recurrences that yield order-p local convergence, together with factorizations that reduce the number of s× s quaternion products per iteration. Beyond NS, we introduce a randomized sketch--and--project method (RSP--Q), a hybrid RSP+NS scheme that interleaves inexpensive randomized projections with an exact hyperpower step, and a matrix--form conjugate gradient on the normal equations (CGNE--Q). All algorithms operate directly in ℍ (no real or complex embeddings) and are matrix--free.