An Augmented Lagrangian Uzawa Iterative Method for Solving Double Saddle-Point Systems With Semidefinite (2,2) Block and Its Application to DLM FD Method for Elliptic Interface Problems. Commun. Comput. Phys. 2021, 30 (1), 124-143. https://doi.org/10.4208/cicp.OA-2020-0084.