Preconditioned proximal point algorithm (source code)

= Preconditioned proximal point algorithm

Given a <positive-definite matrix> $M$, the preconditioned proximal point map is
$$
T=(M+F)^{-1}M.
$$
It is firmly nonexpansive in the weighted <inner product> $\langle v,w\rangle_M=\langle v,Mw\rangle$ whenever $F$ is <monotone operator>[monotone].