TY - JOUR

T1 - A new formulation for the numerical proof of the existence of solutions to elliptic problems

AU - Sekine, Kouta

AU - Nakao, Mitsuhiro T.

AU - Oishi, Shin'ichi

N1 - Publisher Copyright:
Copyright © 2019, The Authors. All rights reserved.
Copyright:
Copyright 2020 Elsevier B.V., All rights reserved.

PY - 2019/10/1

Y1 - 2019/10/1

N2 - Infinite-dimensional Newton methods can be effectively used to derive numerical proofs of the existence of solutions to partial differential equations (PDEs). In computer-Assisted proofs of PDEs, the original problem is transformed into the infinite Newton-Type fixed point equation w = L-1F(u) + L-1G(w), where L is a linearized operator, F(u) is a residual, and G(w) is a local Lipschitz term. Therefore, the estimations of kL-1F(u)k and kL-1G(w)k play major roles in the verification procedures. In this paper, using a similar concept as the `Schur complement' for matrix problems, we represent the inverse operator L-1as an infinite-dimensional operator matrix that can be decomposed into two parts, one finite dimensional and one infinite dimensional. This operator matrix yields a new effective realization of the infinite-dimensional Newton method, enabling a more efficient verification procedure compared with existing methods for the solution of elliptic PDEs. We present some numerical examples that confirm the usefulness of the proposed method. Related results obtained from the representation of the operator matrix as L-1are presented in the appendix.MSC Codes 65G20, 65N30, 35J25

AB - Infinite-dimensional Newton methods can be effectively used to derive numerical proofs of the existence of solutions to partial differential equations (PDEs). In computer-Assisted proofs of PDEs, the original problem is transformed into the infinite Newton-Type fixed point equation w = L-1F(u) + L-1G(w), where L is a linearized operator, F(u) is a residual, and G(w) is a local Lipschitz term. Therefore, the estimations of kL-1F(u)k and kL-1G(w)k play major roles in the verification procedures. In this paper, using a similar concept as the `Schur complement' for matrix problems, we represent the inverse operator L-1as an infinite-dimensional operator matrix that can be decomposed into two parts, one finite dimensional and one infinite dimensional. This operator matrix yields a new effective realization of the infinite-dimensional Newton method, enabling a more efficient verification procedure compared with existing methods for the solution of elliptic PDEs. We present some numerical examples that confirm the usefulness of the proposed method. Related results obtained from the representation of the operator matrix as L-1are presented in the appendix.MSC Codes 65G20, 65N30, 35J25

UR - http://www.scopus.com/inward/record.url?scp=85095145708&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85095145708&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:85095145708

JO - Nuclear Physics A

JF - Nuclear Physics A

SN - 0375-9474

ER -