Home Science On the sharp linear convergence rate of the...
Science

On the sharp linear convergence rate of the circumcentered--reflection method on subspaces

Key Points

arXiv:2606.07888v1 Announce Type: cross Abstract: For two subspaces $U,V\subseteq\RR^n$, the circumcentered--reflection method (CRM) of Behling, Bello-Cruz, and Santos~\cite{BBS2018} computes the projection onto $U\cap V$ using only the reflections across $U$ and $V$, with known linear-convergence rate $c_F$, the cosine of the Friedrichs angle. We prove that, when CRM is initialized in $V$, it contracts at the strictly smaller rate...

arXiv:2606.07888v1 Announce Type: cross Abstract: For two subspaces $U,V\subseteq\RR^n$, the circumcentered--reflection method (CRM) of Behling, Bello-Cruz, and Santos~\cite{BBS2018} computes the projection onto $U\cap V$ using only the reflections across $U$ and $V$, with known linear-convergence rate $c_F$, the cosine of the Friedrichs angle. We prove that, when CRM is initialized in $V$, it contracts at the strictly smaller rate $\rho_V=(\sin^2\theta_p-\sin^2\theta_F)/(\sin^2\theta_p+\sin^2\theta_F)$, where $\theta_F\in(0,\pi/2]$ is the Friedrichs angle and $\theta_p\in[\theta_F,\pi/2]$ the largest principal angle between $U$ and $V$. The bound is sharp, attained on an explicit ray in $V$, and optimal among parameter-free single-step iterations. The constant itself is not new: Bauschke, Bello-Cruz, Nghia, Phan, and Wang~\cite{BBNPW2016} identified it as the optimal rate of the relaxed alternating-projection family and of their adaptive linesearch map $B_T$. Our contribution is that the parameter-free geometric circumcenter attains it as well, via Kantorovich's inequality applied to a single self-adjoint operator on $V$. Restricted to $V$, CRM coincides pointwise with the linesearch maps $A_T$ and $B_T$ from the Gubin--Polyak--Raik framework~\cite{GPR1967}. We further prove $\rho_V<\pi/2$, with one-step convergence exactly when $\theta_F=\theta_p$. Over-reflecting either or both of $R_U$, $R_V$ inside the circumcenter does not help. Going faster than $\rho_V$ universally requires memory: Chebyshev semi-iteration applied to $P_VP_U$ attains a strictly smaller rate, beating $\rho_V$ by a factor at most $2$, attained in the limit $\theta_F\to\theta_p$.
Behling (LOCATION) Bello-Cruz (PERSON) Friedrichs (PERSON) Nghia, Phan (ORG) Kantorovich (ORG) Restricted (ORG) \theta_F=\theta_p$. (ORG)
Originally published by arXiv CS Read original →