An ADMM-Newton-CNN Numerical Approach to a TV Model for Identifying Discontinuous Diffusion Coefficients in Elliptic Equations: Convex Case with Gradient Observations

by   Wenyi Tian, et al.

Identifying the discontinuous diffusion coefficient in an elliptic equation with observation data of the gradient of the solution is an important nonlinear and ill-posed inverse problem. Models with total variational (TV) regularization have been widely studied for this problem, while the theoretically required nonsmoothness property of the TV regularization and the hidden convexity of the models are usually sacrificed when numerical schemes are considered in the literature. In this paper, we show that the favorable nonsmoothness and convexity properties can be entirely kept if the well-known alternating direction method of multipliers (ADMM) is applied to the TV-regularized models, hence it is meaningful to consider designing numerical schemes based on the ADMM. Moreover, we show that one of the ADMM subproblems can be well solved by the active-set Newton method along with the Schur complement reduction method, and the other one can be efficiently solved by the deep convolutional neural network (CNN). The resulting ADMM-Newton-CNN approach is demonstrated to be easily implementable and very efficient even for higher-dimensional spaces with fine mesh discretization.



There are no comments yet.


page 16

page 20

page 23


Identifying source term in the subdiffusion equation with L^2-TV regularization

In this paper, we consider the inverse source problem for the time-fract...

Inertial nonconvex alternating minimizations for the image deblurring

In image processing, Total Variation (TV) regularization models are comm...

A Dual Symmetric Gauss-Seidel Alternating Direction Method of Multipliers for Hyperspectral Sparse Unmixing

Since sparse unmixing has emerged as a promising approach to hyperspectr...

An ADMM-LAP method for total variation blind deconvolution of adaptive optics retinal images

Adaptive optics (AO) corrected flood imaging of the retina is a popular ...

Bilinear Constraint based ADMM for Mixed Poisson-Gaussian Noise Removal

In this paper, we propose new operator-splitting algorithms for the tota...

Graph- and finite element-based total variation models for the inverse problem in diffuse optical tomography

Total variation (TV) is a powerful regularization method that has been w...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.