Solver-in-the-Loop joint operator learning: fractional Laplace-Beltrami features for interface reconstruction

Abstract

In this work, we propose a joint operator learning method for reconstructing images of conductivity coefficients from boundary data. Inspired by the idea of employing partial differential equation (PDE) solvers as preconditioners for this inverse problem, we investigate a ``solver-in-the-loop'' training mechanism. It allows the interaction of learnable parameters integrated in a PDE solver module and those in neural networks for reconstructing images. Specifically, we employ a fractional Laplace-Beltrami operator with a learnable fractional order, which transforms boundary data into high-dimensional features. These features then serve as input to a neural network, significantly improving reconstruction accuracy. For this purpose, a Learning-Automated FEM (LA-FEM) package, facilitating this ``solver-in-the-loop'' property, is developed with PyTorch as a backend. The new LA-FEM module conveniently allows the auto-differentiation regarding an objective function to freely propagate through the PDE solver from the forward problem and the coupled neural networks for the inverse problem.

References in corpus (15)

• B-PINNs: Bayesian Physics-Informed Neural Networks for Forward and Inverse PDE Problems with Noisy Data
• Ten equivalent definitions of the fractional Laplace operator
• Deep D-bar: Real time Electrical Impedance Tomography Imaging with Deep Neural Networks
• Physics-informed graph neural Galerkin networks: A unified framework for solving PDE-governed forward and inverse problems
• Automated derivation of the adjoint of high-level transient finite element programs
• Solving Electrical Impedance Tomography with Deep Learning
• Solving inverse-PDE problems with physics-aware neural networks
• Numerical Solution of Inverse Problems by Weak Adversarial Networks
• BCR-Net: a neural network based on the nonstandard wavelet form
• Relaxed Gauss-Newton methods with applications to electrical impedance tomography
• Residual-based error correction for neural operator accelerated infinite-dimensional Bayesian inverse problems
• A Fixed Mesh Method With Immersed Finite Elements for Solving Interface Inverse Problems
• Recovering the Potential and Order in One-Dimensional Time-Fractional Diffusion with Unknown Initial Condition and Source
• Adaptive Reconstruction for Electrical Impedance Tomography with a Piecewise Constant Conductivity
• A data-assisted two-stage method for the inverse random source problem