ISSN 2079-6900 (Print) 
ISSN 2587-7496 (Online)

Middle Volga Mathematical Society Journal

DOI 10.15507/2079-6900.28.202602.83-98

Original article

ISSN 2079-6900 (Print)

ISSN 2587-7496 (Online)

MSC2020 65R20

Machine learning-based preconditioner to solve Poisson equation

E. A. Chekmeneva1, T. S. Khachkova2, V. V. Lisitsa2

1Novosibirsk State University (Novosibirsk, Russian Federation)

2Sobolev Institute of Mathematics (Novosibirsk, Russian Federation)

Abstract. In this paper, we present a method of constructing a preconditioner based on machine learning methods for use in the numerical solution of the Poisson equation in porous media modelling. The iterative conjugate gradient method is used to solve the problem. To precondition a system of linear algebraic equations, we propose to approximate the inverse Laplace operator using a convolutional neural network of the U-net architecture. We consider two alternative approaches to the formation of a training dataset for a neural network. The first method is based on the use of pairs of vectors and on the results of applying the Laplace operator to them. In the second method, the training pairs are residual vectors obtained by implementing the conjugate gradient method and results of applying the Laplace operator to them. The neural network learning process is based on minimizing the $L^2$-norm relative error. We show that while using each of the presented learning methods, the U-net neural network with five convolutional layers approximating the inverse Laplace operator provides insufficient accuracy to reduce the number of iterations in the conjugate gradient method. Therefore, the modified conjugate gradient method is stabilized and has an irreducible residual.

Key Words: Poisson equation, method of conjugate gradients, preconditioner, machine learning

For citation: E. A. Chekmeneva, T. S. Khachkova, V. V. Lisitsa. Machine learning-based preconditioner to solve Poisson equation. Zhurnal Srednevolzhskogo matematicheskogo obshchestva. 28:2(2026), 83–98. DOI: https://doi.org/10.15507/2079-6900.28.202602.83-98

Submitted: 01.12.2025; Revised: 24.02.2026; Accepted: 01.03.2026

Information about the authors:

Ekaterina A. Chekmeneva, student at Novosibirsk State University (2 Pirogova Street, Novosibirsk, 630090, Russia), ORCID: https://orcid.org/0009-0003-8509-0534, e.chekmeneva@g.nsu.ru

Tatyana S. Khachkova, Ph.D. (Phys. and Math.), senior researcher at the Institute of Mathematics, Siberian Branch of Russian Academy of Science (4 Koptug Avenue, Novosibirsk, 630090, Russia), ORCID: https://orcid.org/0000-0002-1595-7142, hachtanya@mail.ru

Vadim V. Lisitsa, D.Sc. (Phys. and Math.), head of laboratory at the Institute of Mathematics, Siberian Branch of Russian Academy of Science (4 Koptug Avenue, Novosibirsk, 630090, Russia), ORCID: https://orcid.org/0000-0003-3544-4878, v.v.lisitsa@math.nsc.ru

All authors have read and approved the final manuscript.

Conflict of interest: The authors declare no conflict of interest.

Creative Commons Attribution 4.0 International License This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International License.