Efficient primal--dual splitting methods for a Poisson-constrained JKO scheme for Poisson-Nernst-Planck models
To efficiently solve the resulting large-scale constrained problems, preconditioned and transformed primal--dual algorithms equipped with tailored fast dual solvers that exploit the coupled block structure of the PDE constraints are developed.