Hierarchical interpolative factorization for elliptic operators: integral equations (Q2812291)

The authors study the methods to approximate solving the integral equation NEWLINE\[NEWLINE a(x) u(x)+ b(x)\int_{\Omega}K(\|x-y\|u(y)d\Omega (y)=f(x),\qquad x\in \Omega \subset \mathbb R^{d},NEWLINE\]NEWLINE derived from the boundary value problem NEWLINE\[NEWLINE \nabla ( a(x) \nabla u(x)) = f(x) \qquad x\in \Omega\subset \mathbb R^{d}, \,d= 2,3 NEWLINE\]NEWLINE Discretization of the integral equation, using the different approximate schemes (collocation, Galerkin method and other), leads to a system of linear equations NEWLINE\[NEWLINE Ax = f, \qquad A \in C^{N\times N}.NEWLINE\]NEWLINE The main idea of this paper is to construct algorithms which permit to obtain the matrix \(A\) in factorized form \(A= LV\) or \(A=LDV\). Also an analysis of the complexity of the constructed algorithms and the numerical results for three concrete examples are given.