  We propose a flexible approach for computing the resolvent of the sum of
weakly monotone operators in real Hilbert spaces. This relies on splitting
methods where strong convergence is guaranteed. We also prove linear
convergence under Lipschitz continuity assumption. The approach is then applied
to computing the proximity operator of the sum of weakly convex functions, and
particularly to finding the best approximation to the intersection of convex
sets.
