An official website of the United States government
We studied the least-squares ReLU neural network (LSNN) method for solving a linear advection-reaction equation with discontinuous solution in [Z. Cai et al., J. Comput. Phys., 443 (2021), 110514]. The method is based on a least-squares formulation and uses a new class of approximating functions: ReLU neural network (NN) functions. A critical and additional component of the LSNN method, differing from other NN-based methods, is the introduction of a properly designed and physics preserved discrete differential operator. In this paper, we study the LSNN method for problems with discontinuity interfaces. First, we show that ReLU NN functions with depth \(\lceil \log\_2(d+1)\rceil+1\) can approximate any \(d\)-dimensional step function on a discontinuity interface generated by a vector field as streamlines with any prescribed accuracy. By decomposing the solution into continuous and discontinuous parts, we prove theoretically that the discretization error of the LSNN method using ReLU NN functions with depth \(\lceil \log\_2(d+1)\rceil+1\) is mainly determined by the continuous part of the solution provided that the solution jump is constant. Numerical results for both two- and three-dimensional test problems with various discontinuity interfaces show that the LSNN method with enough layers is accurate and does not exhibit the common Gibbs phenomena along discontinuity interfaces.
more »
« less
On a class of splines free of Gibbs phenomenon
When interpolating data with certain regularity, spline functions are useful. They are defined as piecewise polynomials that satisfy certain regularity conditions at the joints. In the literature about splines it is possible to find several references that study the apparition of Gibbs phenomenon close to jump discontinuities in the results obtained by spline interpolation. This work is devoted to the construction and analysis of a new nonlinear technique that allows to improve the accuracy of splines near jump discontinuities eliminating the Gibbs phenomenon. The adaption is easily attained through a nonlinear modification of the right hand side of the system of equations of the spline, that contains divided differences. The modification is based on the use of a new limiter specifically designed to attain adaption close to jumps in the function. The new limiter can be seen as a nonlinear weighted mean that has better adaption properties than the linear weighted mean. We will prove that the nonlinear modification introduced in the spline keeps the maximum theoretical accuracy in all the domain except at the intervals that contain a jump discontinuity, where Gibbs oscillations are eliminated. Diffusion is introduced, but this is fine if the discontinuity appears due to a discretization of a high gradient with not enough accuracy. The new technique is introduced for cubic splines, but the theory presented allows to generalize the results very easily to splines of any order. The experiments presented satisfy the theoretical aspects analyzed in the paper.
more »
« less
- Award ID(s):
- 1719410
- PAR ID:
- 10226110
- Date Published:
- Journal Name:
- ESAIM: Mathematical Modelling and Numerical Analysis
- Volume:
- 55
- Issue:
- Supplement
- ISSN:
- 0764-583X
- Page Range / eLocation ID:
- S29 to S64
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Abstract The M waves introduced by Burridge and Willis (1969) are emitted by the surface of a self-similarly expanding elliptical crack, and they give Rayleigh waves at the corresponding crack speed. In the analysis for the self-similarly expanding spherical inclusion with phase change (dynamic Eshelby problem) the M waves are related to the waves obtained on the basis of the dynamic Green’s function containing the contribution from the latest wavelets emitted by the expanding boundary of phase discontinuity, and they satisfy the Hadamard jump conditions for compatibility and linear momentum across the moving phase boundary of discontinuity. In the interior of the expanding inclusion they create a “lacuna” with zero particle velocity by canceling the effect of the P and S. It is shown that the “lacuna” and Eshelby properties are also valid for a Newtonian fluid undergoing phase change in a self-similarly expanding ellipsoidal region of a fluid with different viscosity.more » « less
-
Abstract We consider a surveillance-evasion game in an environment with obstacles. In such an environment, a mobile pursuer seeks to maintain visibility with a mobile evader, who tries to hide from the pursuer in the shortest time possible. In this two-player zero-sum game setting, we study the discontinuities of the value of the game near the boundary of the target set, where the players cannot see each other (the non-visibility region). In particular, we describe the transition between the usable part of the boundary of the target (where the value vanishes) and the non-usable part (where the value is positive). We show that the value exhibits different behaviour depending on the regularity of the obstacles. Namely, we prove that the boundary profile is continuous in the case of smooth obstacles and that it exhibits a jump discontinuity when the obstacle contains corners. Moreover, we prove that, in the latter case, there is a semi-permeable barrier emanating from the interface between the usable and the non-usable part of the boundary of the target set.more » « less
-
Abstract A mesh refinement method is described for solving optimal control problems using Legendre‐Gauss‐Radau collocation. The method detects discontinuities in the control solution by employing an edge detection scheme based on jump function approximations. When discontinuities are identified, the mesh is refined with a targetedh‐refinement approach whereby the discontinuity locations are bracketed with mesh points. The remaining smooth portions of the mesh are refined using previously developed techniques. The method is demonstrated on two examples, and results indicate that the method solves optimal control problems with discontinuous control solutions using fewer mesh refinement iterations and less computation time when compared with previously developed methods.more » « less
-
Abstract Lindblad dynamics and other open-system dynamics provide a promising path towards efficient Gibbs sampling on quantum computers. In these proposals, the Lindbladian is obtained via an algorithmic construction akin to designing an artificial thermostat in classical Monte Carlo or molecular dynamics methods, rather than being treated as an approximation to weakly coupled system-bath unitary dynamics. Recently, Chen, Kastoryano, and Gilyén (arXiv:2311.09207) introduced the first efficiently implementable Lindbladian satisfying the Kubo–Martin–Schwinger (KMS) detailed balance condition, which ensures that the Gibbs state is a fixed point of the dynamics and is applicable to non-commuting Hamiltonians. This Gibbs sampler uses a continuously parameterized set of jump operators, and the energy resolution required for implementing each jump operator depends only logarithmically on the precision and the mixing time. In this work, we build upon the structural characterization of KMS detailed balanced Lindbladians by Fagnola and Umanità, and develop a family of efficient quantum Gibbs samplers using a finite set of jump operators (the number can be as few as one), akin to the classical Markov chain-based sampling algorithm. Compared to the existing works, our quantum Gibbs samplers have a comparable quantum simulation cost but with greater design flexibility and a much simpler implementation and error analysis. Moreover, it encompasses the construction of Chen, Kastoryano, and Gilyén as a special instance.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1051/m2an/2020021
