## WELL-CONDITIONED DIPOLE-TYPE METHOD OF FUNDAMENTAL SOLUTIONS: DERIVATION AND ITS MATHEMATICAL ANALYSIS

## KOYA SAKAKIBARA

Abstract. In this paper, we examine the dipole-type method of fundamental solutions, which can be conceptualized as a discretization of the 'singularity-removed' double-layer potential. We present a method for removing the ill-conditionality, which was previously considered a significant challenge, and provide a mathematical analysis in the context of disk regions. Moreover, we extend the proposed method to the general Jordan region using conformal mapping, demonstrating the efficacy of the proposed method through numerical experiments.

## Contents

| 1.         | Introduction                                     |   1 |
|------------|--------------------------------------------------|-----|
| 2.         | Dipole simulation method                         |   5 |
| 3.         | DSM-QR                                           |   5 |
| 4.         | Mathematical analysis of DSM-QR                  |   7 |
| 4.1.       | Concrete expression of basis functions of DSM-QR |   7 |
| 4.2.       | Unique existence and condition number            |   9 |
| 4.3.       | Error estimate                                   |  12 |
| 5.         | Numerical experiments                            |  15 |
| 5.1.       | Disk                                             |  15 |
| 5.2.       | Jordan region                                    |  16 |
| 6.         | Concluding remarks                               |  17 |
| References | References                                       |  19 |

## 1. Introduction

The method of fundamental solutions (MFS, for short) is a mesh-free numerical solver for partial differential equations, and it has been actively used especially in the field of engineering. Its idea is very simple. Let us consider the boundary-value problem for the linear partial differential operator L ; more precisely,

$$\begin{cases} \mathcal { L } u = 0 & \text{in } \Omega, \\ \mathcal { B } g = f & \text{on } \partial \Omega, \end{cases}$$

̸

where Ω ⊂ R 2 denotes a region such that R 2 \ Ω = ∅ . The second equation represents the linear boundary condition such as Dirichlet ( B = id), Neumann ( B = ∂ ν ), or Robin ( B = id + γ∂ ν ). Suppose that the partial differential operator L has a fundamental solution E ; that is, E satisfies

L

E

=

δ

in the distributional sense. The MFS offers an approximate solution to the problem (1.1) as in the following procedure.

2020 Mathematics Subject Classification. 65N80, 65N35, 65N12, 35J05, 35J08.

Key words and phrases. Method of fundamental solutions, Double-layer potential, QR decomposition, Illconditionality.

- (i) Choose N ∈ N , and arrange the singular points { y k } N k =1 ⊂ R 2 \ Ω 'suitably'.
- (ii) Seek an approximate solution u ( N ) of the form

$$u ^ { ( N ) } ( x ) = \sum _ { k = 1 } ^ { N } Q _ { k } E ( x - y _ { k } ), \ \ x \in \mathbb { R } ^ { 2 } \, \sum \{ y _ { 1 }, y _ { 2 }, \dots, y _ { N } \}.$$

- (iii) Determine the coefficients { Q k } N k =1 by the collocation method; that is, choose P ∈ N so that P ≥ N , arrange the collocation points { x j } P j =1 ⊂ ∂ Ω 'suitably', and impose the collocation equations:

$$u ^ { ( N ) } ( x _ { j } ) = f ( x _ { j } ), \ \ j = 1, 2, \dots, P,$$

which is equivalent to the linear system

$$\mathbf G Q = f,$$

$$\begin{matrix} \end{matrix}$$

where

$$\mathbf G \coloneqq ( \mathcal { B } E ( x _ { j } - y _ { k } ) ) _ { \stackrel { j = 1, 2, \dots, P } { k = 1, 2, \dots, N } } \in \mathbb { R } ^ { P \times N }, \ \mathbf Q \coloneqq ( Q _ { k } ) _ { k = 1 } ^ { N } \in \mathbb { R } ^ { N }, \ \ f \coloneqq ( f ( x _ { j } ) ) _ { j = 1 } ^ { P } \in \mathbb { R } ^ { P }.$$

As previously stated, in contrast to the finite element method and the finite difference method, the MFS does not necessitate the partitioning of the region into a mesh. An approximate solution can be derived by simply positioning points in a suitable manner on the boundary and exterior of the region. Consequently, the MFS can be regarded as a mesh-free numerical method. It has been demonstrated in numerous prior studies that when the singular and collocation points are correctly positioned and the boundary data is real-analytic, the approximation error exhibits a noteworthy property of decaying exponentially with respect to the number of points [13, 10, 11, 12, 19, 14, 16, 23, 4, 17, 18, 15, 21, 22, 7]. Nevertheless, the outcomes of these mathematical investigations remain inadequate. For instance, there is a paucity of findings pertaining to the existence of approximate solutions in multiply-connected regions in the plane, error estimation, and analysis in high-dimensional spaces. Furthermore, the actual numerical computation with the MFS is subject to the issue of ill-conditionality. This refers to the fact that the condition number of the coefficient matrix G of a linear system (1.2) increases exponentially with respect to the number of points, indicating that the linear system (1.2) is exceedingly challenging to solve. Numerous prior studies have devised MFS-specific pre-processing techniques based on numerical linear algebra (see, for instance, [6, 5]). Nevertheless, the issue of ill-conditionality persists, and it has been identified as a significant impediment in actual numerical computation.

In light of these considerations, Antunes put forth an alternative approach in his paper [2], namely, modifying the basis function. This concept was previously proposed in [9] within the context of the radial basis function method.

The infinite series expansion obtained by the Taylor series expansion of the basis function is truncated by a finite number of terms (e.g., up to the M -th order terms) to obtain

$$\Theta ( r, \theta ) = \text{BDF} ( r, \theta ),$$

where the singular points are arranged as y k = ε -1 (cos ϕ , k sin ϕ k ) ⊤ and

$$\mathbf B \coloneqq \begin{pmatrix} 1 & \sin \phi _ { 1 } & \cos \phi _ { 1 } & \cdots & \sin M \phi _ { 1 } & \cos M \phi _ { 1 } \\ 1 & \sin \phi _ { 2 } & \cos \phi _ { 2 } & \cdots & \sin M \phi _ { 2 } & \cos M \phi _ { 2 } \\ \vdots & \vdots & \vdots & \vdots & \vdots \\ 1 & \sin \phi _ { N } & \cos \phi _ { N } & \cdots & \sin M \phi _ { N } & \cos M \phi _ { N } \\ / & \quad \bar { \ } M \ \bar { \ } M \ \end{pmatrix} \in \mathbb { R } ^ { N \times ( 2 M + 1 ) },$$

$$\mathbf D \coloneqq \text{diag} \left ( \log \varepsilon, \varepsilon, \varepsilon, \dots, \frac { \varepsilon ^ { M } } { M }, \frac { \varepsilon ^ { M } } { M } \right ) \in \mathbb { R } ^ { ( 2 M + 1 ) \times ( 2 M + 1 ) },$$

$$F ( r, \theta ) \coloneqq \begin{pmatrix} 1 \\ r \sin \theta \\ r \cos \theta \\ \vdots \\ r ^ { M } \cos M \theta \\ r ^ { M } \sin M \theta \end{pmatrix} \in \mathbb { R } ^ { 2 M + 1 }, \\ \sl i s \, \text{chosen to satisfy} \, 2 M + 1 > \, N. \, \text{In this mat}$$

and M is chosen to satisfy 2 M + 1 &gt; N . In this matrix representation, matrix B represents a well condition, whereas the ill-conditionality mainly stems from matrix D . Accordingly, the matrix B is decomposed via the QR algorithm to yield the representation B = QR . The N -th leading principal submatrix of D is subsequently designated as D N . The new basis function is then expressed as

$$\Psi ( r, \theta ) = ( \psi _ { k } ( r, \theta ) ) _ { k = 1 } ^ { N } \coloneqq D _ { N } ^ { - 1 } Q ^ { \top } \Theta ( r, \theta ) = D _ { N } ^ { - 1 } R D F ( r, \theta ).$$

Subsequently, he put forth the proposition of MFS-QR as a MFS-based system, founded upon this novel basis function. The approximate solution is expressed in the form of

$$u ^ { ( N ) } ( x ) = u ^ { ( N ) } ( r, \theta ) = \sum _ { k = 1 } ^ { N } Q _ { k } \psi _ { k } ( r, \theta ), \quad x = \begin{pmatrix} r \cos \theta \\ r \sin \theta \end{pmatrix}.$$

The coefficients of the linear combination are determined by the collocation method.

In paper [2], it was numerically shown that the ill-conditionality of the coefficient matrix can be completely eliminated for the disk domain, and that for other domains, the ill-conditionality cannot be completely eliminated, but can be reduced to some extent. Later, in paper [3], a method based on singular value decomposition (MFS-SVD) instead of QR decomposition was developed, and it was numerically confirmed that the ill-conditionality can be removed when the domain is star-shaped. It should be noted that in both cases, the order of convergence is the same as that of the original MFS. In other words, MFS-QR and MFS-SVD achieve the exponential decay of errors, which is one of the major advantages of MFS, while eliminating one of the major problems of MFS, i.e., the problem of illconditionality. Therefore, it is expected that MFS-QR and MFS-SVD, rather than ordinary MFS, will become more important numerical methods in the future, but unfortunately, there are no theoretical guarantees to support the numerical experimental results obtained in previous studies.

It is also important to consider whether the construction of the approximate solution of the MFS is a sensible approach in the first place. In a sense, the approximate solution obtained through the MFS can be regarded as a 'singularity-removed' discretization of the single-layer potential. However, as is well established in the theory of potential problems [8], the single-layer potential should be employed for problems with Neumann boundary conditions, whereas the double-layer potential is more appropriate for problems with Dirichlet boundary conditions. In light of this awareness of the problem, the dipoletype method of fundamental solutions (also known as the dipole simulation method) was proposed. For further details, please refer to § 2. Henceforth, the dipole-type method of fundamental solutions will be referred to as the DSM. In the DSM, the approximate solution is constructed using the normal derivative of the fundamental solution as the basis function in lieu of the fundamental solution. In the original theoretical analysis of the MFS, it was necessary to make some mathematically unnatural assumptions. However, with the DSM, a similar theoretical analysis can be performed without the need to make unnatural assumptions. Therefore, it is more natural to use the DSM when dealing with Dirichlet boundary conditions. Nevertheless, even in the DSM, the problem of ill-conditionality persists, and to the best of the author's knowledge, there is no known solution to this problem.

The objective of this paper is to put forth a novel DSM, the DSM-QR, which is founded upon the tenets of the MFS-QR and eliminates the problematic ill-conditionality. Additionally, this paper endeavors to conduct a fundamental mathematical analysis of the DSM-QR. With the intention of establishing a robust foundation for future mathematical analysis, this paper concentrates on the Dirichlet boundary value problem

$$\begin{cases} \triangle u = 0 & \text{in } \mathbb { B } _ { \rho }, \\ u = f & \text{on } \partial \mathbb { B } _ { \rho }, \end{cases}$$

for the Laplace equation in the disk domain, where B ρ denotes the disk in the complex plane with radius ρ having the origin as its center; that is,

$$\mathbb { B } _ { \rho } \coloneqq \{ z \in \mathbb { C } \, | \, | z | < \rho \}.$$

The collocation points { x j } P j =1 and singular points { y k } N k =1 are to be placed as

$$\Lambda ^ { \prime }, \Lambda ^ { - 1 } \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \ \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \cdot \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \quad \Lambda \end{cases}$$

where R &gt; ρ and R 2 and C are identified here and hereafter. Moreover, set

$$\kappa \coloneqq \frac { \rho } { R }.$$

The following three theorems represent the core of this paper's theoretical framework.

Theorem 1.1. The DSM-QR gives a unique approximate solution.

Denote by cond ( 2 A ) the condition number of A induced by the ℓ 2 -norm.

Theorem 1.2. The condition number cond ( 2 G ) of the coefficient matrix G is of O(1) as N → ∞ . More precisely, it holds that

$$\text{cond} _ { 2 } ( \mathbf G ) = \frac { 2 } { 1 + \kappa ^ { N - 2 } }, \ \ N \in \mathbb { N }.$$

↦

Theorem 1.3. (i) Suppose that the Fourier series ∑ n ∈ Z f n e i nθ of θ → f ( ρ e i θ ) is absolutely convergent. Then, u ( N ) uniformly converges to u on B ρ ; that is,

$$\| u - u ^ { ( N ) } \| _ { L ^ { \infty } ( \mathbb { B } ) } \longrightarrow 0 \ \ ( N \to \infty ).$$

- (ii) Suppose there exists an α &gt; 1 such that | f n | = O( | n | -α ) ( | n | →∞ ) holds. Then, there exists a positive constant C , independent of N , such that

$$\| u - u ^ { ( N ) } \| _ { L ^ { \infty } ( \mathbb { B } _ { \rho } ) } \leq C N ^ { - \alpha + 1 }$$

holds.

- (iii) Suppose that there exists an a ∈ (0 1) , such that | f n | = O( a | n | ) holds as | n | →∞ . Then, there exists a positive constant C , independent of N , such that

$$\| u - u ^ { ( N ) } \| _ { L ^ { \infty } ( \mathbb { B } _ { \rho } ) } \leq C \times \begin{cases} a ^ { N / 2 } & \text{if $a>\kappa^{2}$}, \\ N \kappa ^ { N } & \text{if $a=\kappa^{2}$}, \\ \kappa ^ { N } & \text{if $a<\kappa^{2}$}. \end{cases}$$

The following is a description of the organization of this paper. In Section 2, we undertake a review of the MFS from the perspective of potential theory, and subsequently derive a tabular expression for the approximate solution by the DSM. In Section 3, we put forth the DSM-QR, which is based on the concept of the MFS-QR. Section 4 presents a mathematical analysis of DSM-QR. First, we provide a detailed account of the basis functions utilized in the DSM-QR through the presentation of concrete calculations pertaining to QR decomposition. Subsequently, the eigenvalues of the coefficient matrix are obtained in Lemma 4.2, which also demonstrates that the condition number of the coefficient matrix is of O(1). Then, we provide a proof of the error estimate presented in Theorem 1.3. In Section 5, we present the results of numerical experiments. First, we consider the case of the disk region and demonstrate that the problematic ill-conditionality is entirely eliminated and that the error decays exponentially, as predicted by the main theorem of this paper. Next, we extend the DSM-QR to the case of Jordan regions, although no mathematical analysis is available. Specifically, DSM-QR is extended to the general Jordan region by mapping point configurations in the disk region to the Jordan region by a conformal mapping, and by viewing basis functions as perturbations of basis functions in

the disk region. The utility of this approach is demonstrated through numerical experiments. Section 6 provides a summary of the contents of this paper and suggests avenues for future research.

## 2. Dipole simulation method

In this section, we briefly explain the idea of the dipole simulation method (DSM).

Let Γ be a Jordan curve surrounding Ω; more precisely, Ω ⊂ Int Ω[Γ], where Ω[Γ] denotes a Jordan region with boundary Γ and Int Ω[Γ] denotes the interior of Ω[Γ]. Suppose that there exists a function Q defined on Γ such that the solution u to the Dirichlet boundary value problem (1.3) can be expressed as

$$u ( z ) = \int _ { \Gamma } E ( z - \zeta ) Q ( \zeta ) \, \mathrm d \sigma ( \zeta ) = \colon \text{SL} [ Q ] ( z ),$$

which can regarded as a 'singularity-removed' single-layer potential. Here, E = E z ( ) = -(2 π ) -1 log | z | is the fundamental solution to the Laplace operator -△ . If we consider a linear combination of Dirac delta distributions Q ( N ) ( z ) := ∑ N k =1 Q δ ζ k ( -ζ k ) with ζ k ∈ Γ ( k = 1 2 , , . . . , N ), we formally obtain that

$$\text{SL} [ Q ^ { ( N ) } ] ( z ) = \sum _ { k = 1 } ^ { N } Q _ { k } E ( z - \zeta _ { k } ),$$

which is an approximate solution by the MFS. Namely, the MFS can be regarded as a discretization of the 'singularity-removed' single-layer potential representation.

However, it is natural to consider the double-layer potential, not single-layer potential, when we deal with the Dirichlet boundary condition (see, for instance, [8]). Recall that the dipole-layer potential is given by replacing the integral kernel E z ( -ζ ) with its normal derivative as

$$\text{DL} [ Q ] ( z ) \coloneqq \int _ { \Gamma } \frac { \partial E } { \partial \nu ( \zeta ) } ( z - \zeta ) \, \text{d} \zeta,$$

where ν ζ ( ) denotes the unit outward normal vector to Γ. Then, we formally obtain

$$\text{dL} [ Q ^ { ( N ) } ] ( z ) = \sum _ { k = 1 } ^ { N } Q _ { k } \frac { \partial E } { \partial \nu ( \zeta _ { k } ) } ( z - \zeta _ { k } ).$$

Therefore, it is natural to consider that a function of the form (2.1) can also approximate the solution to the Dirichlet boundary value problem (1.3). This is the idea of the DSM. More precisely, the algorithm of the DSM can be described as follows:

- (i) Choose N ∈ N , and arrange the singular points { ζ k } N k =1 ⊂ C \ Ω and the dipole moments { ν k } N k =1 ⊂ S 1 suitably, where S 1 := { z ∈ C | | z | = 1 } denotes the unit circle.
- (ii) Seek an approximate solution u ( N ) of the form

$$u ^ { ( N ) } ( z ) = \sum _ { k = 1 } ^ { N } Q _ { k } \frac { \partial E } { \partial \nu _ { k } } ( z - \zeta _ { k } ).$$

- (iii) Determine the coefficients { Q k } N k =1 by the collocation method; that is, arrange the collocation points { z j } P j =1 ⊂ ∂ Ω suitably, and impose the collocation equations:

$$u ^ { ( N ) } ( z _ { j } ) = f ( z _ { j }$$

## 3. DSM-QR

In this section, we develop DSM-QR algorithm. In particular, we replace the basis functions { ( ∂E/∂ν k )( · -ζ k ) } N k =1 with new ones by a similar idea to the MFS-QR.

Let us arrange the singular points ζ k on the circle with radius R uniformly as

$$\zeta _ { k } \coloneqq R \omega ^ { k } \ ( k$$

Regarding the circle | z | = R as Γ, it is natural to define dipole moments ν k as

$$\nu _ { k } \coloneqq \frac { \zeta _ { k } } {$$

Then, the basis function of the DSM can be computed as

$$\frac { \partial E } { \partial \nu _ { k } } ( z - \z$$

For z = e r i θ and ζ = R e i ϕ with 0 ≤ r ≤ ρ &lt; R , we formally have

$$\begin{array} {$$

where κ := ρ/R ∈ (0 1) and , s := r/ρ ∈ [0 , 1]. Basis functions (3.1) of the DSM can be obtained by taking ρ = ϕ k := 2 πk/N ( k = 1 2 , , . . . , N ). Namely,

$$\text{taking } \rho = \phi _ { k } \coloneqq 2 \pi$$

$$\begin{matrix} \\ \\ \end{matrix}$$

Hereafter, we omit the constant factor 1 / (2 πR ).

Let M ∈ N satisfying 2 M +1 &gt; N , and truncate the above infinite series expansion at the M -th order. More precisely, set

$$\begin{array} {$$

$$\mathbf D \coloneqq \text{diag} \left ( 1, \kappa$$

$$\mathbf D \coloneqq \text{diag} \left ( 1, \kappa$$

We perform the QR decomposition of B as B = QR , where Q ∈ R N × N is an orthogonal matrix and R ∈ R N × (2 M +1) is an upper triangular matrix. Furthermore, let D N be the N -th leading principal submatrix of D . Then, we define new basis functions as

$$( 3. 2 ) & & \Psi ( s, \theta ) = ( \psi _ { k } ( s, \theta ) ) _ { k = 1 } ^ { N } \coloneqq D _ { N } ^ { - 1 } Q ^ { \top } Q R D F ( s, \theta ) = D _ { N } ^ { - 1 } R D F ( s, \theta ).$$

Hence, the DSM-QR gives an approximate solution to problem (1.3) of the form

$$u ^ { ( N ) } ( z ) = u ^ { ( N ) } ( r, \theta ) = \sum _ { k = 1 } ^ { N } Q _ { k } \psi _ { k } ( r / \rho, \theta ), \ \ z = r \text{e} ^ { i \theta }.$$

Coefficients { Q k } N k =1 are determined by the collocation method; that is,

$$u ^ { ( N ) } ( \rho e ^ { i \theta _ { j } } ) = f ( \rho e ^ { i \theta _ { j } } ), \ \theta _ { j } \coloneqq \frac { 2 \pi j } { P }, \ \ j = 1, 2, \dots, P,$$

which is equivalent to the linear system

$$\text{G} Q = f,$$

$$\psi$$

where

$$\mathbf G \coloneqq ( \psi _ { k } ( 1, \theta _ { j } ) ) _ { \substack { j = 1, 2, \dots, P \\ k = 1, 2, \dots, N } } \in \mathbb { R } ^ { P \times N }, \ \ Q \coloneqq ( Q _ { k } ) _ { k = 1 } ^ { N } \in \mathbb { R } ^ { N }, \ \ f \coloneqq ( f ( \rho e ^ { i \theta _ { j } } ) ) _ { j = 1 } ^ { P } \in \mathbb { R } ^ { P }.$$

4. Mathematical analysis of DSM-QR

- 4.1. Concrete expression of basis functions of DSM-QR. Hereafter, we choose M = N -1 and suppose that N is odd and greater than 1 for the sake of simplicity. We explicitly compute the QR decomposition of B . To this end, define vectors c j , s j ∈ R N by

$$c _ { j } \coloneqq \begin{pmatrix} \cos j \phi _ { 1 } \\ \cos j \phi _ { 2 } \\ \vdots \\ \cos j \phi _ { N } \end{pmatrix}, \quad s _ { j } \coloneqq \begin{pmatrix} \sin j \phi _ { 1 } \\ \sin j \phi _ { 2 } \\ \vdots \\ \sin j \phi _ { N } \end{pmatrix}$$

for j ∈ Z . Then, the matrix B can be represented as

$$\mathbf B = \left ( c _ { 0 } \ s _ { 1 } \ c _ { 1 } \ \cdots \ s _ { N } \ c _ { N } \right ).$$

Lemma 4.1. It holds that

̸

$$\Phi _ { 1 } \colon \, \iota \ n v u s \ u n u \\ \langle c _ { j }, c _ { j ^ { \prime } } \rangle _ { \ell ^ { 2 } } & = \langle s _ { j }, s _ { j ^ { \prime } } \rangle _ { \ell ^ { 2 } } = 0 \quad \text{for $0\leq j,j^{\prime} \leq \frac{N-1}{2}$ with $j\neq j^{\prime}$,} \\ \langle c _ { j }, s _ { j ^ { \prime } } \rangle _ { \ell ^ { 2 } } & = 0 \quad \text{for $0\leq j,j^{\prime} \leq \frac{N-1}{2}$,}$$

$$& \quad \cdots \cdots \cdots \cdots \\ \| c _ { 0 } \| _ { \ell ^ { 2 } } & = \sqrt { N }, \\ \| c _ { j } \| _ { \ell ^ { 2 } } & = \| s _ { j } \| _ { \ell ^ { 2 } } = \sqrt { \frac { N } { 2 } } \quad \text{ for } 1 \leq j \leq \frac { N - 1 } { 2 }.$$

Furthermore, it holds that

$$c _ { N - j } = c _ { j }, \ \ s _ { N - j$$

Since the above lemma immediately follows from the elementary relation

̸

$$\sum _ { l = 1 } ^ { N } \omega ^ { j l } = \$$

we omit its proof. Here and hereafter, the congruence ≡ is taken in modulo N . Define

$$q _ { j } & \coloneqq \begin{cases} \frac { c _ { l$$

Then, Lemma 4.1 yields that

$$\begin{array} {$$

Therefore, setting

$$Q \coloneqq \left ( q _ { 0 } \ q _ { 1 } \ q _$$

$$\begin{array} {$$

we have B = QR . Since Q is an orthogonal matrix by Lemma 4.1 and R is an upper triangular matrix, this is the desired QR decomposition of B . A direct computation shows that

$$| ^ { - 1 }$$

$$\text{matrix}, this is the desired QR decomposition of } B. \ A \text$$

Therefore, the basis functions Ψ = ( ψ k ) N k =1 of the DSM-QR, defined as (3.2), can be explicitly given by

$$\psi _ { k } ( s, \theta ) & = \begin{cases} \$$

4.2. Unique existence and condition number. In this subsection, we prove that the linear system (3.3) is uniquely solvable. We will establish mathematical theory for P = αN with α ∈ N . When α = 1, the linear system (3.3) is understood in a standard manner. When α ≥ 2, we seek a solution to the linear sytem (3.3) in the least-squares sense; that is,

$$Q \in \underset { \tilde { Q } \in \mathbb { R } ^$$

which is equivalent to

$$G ^ { \top } G Q = G ^ { \top } f.$$

Hereafter, we focus on the case where α = 1. Other cases α ≥ 2 can be shown in a similar manner. We examine properties of the coefficient matrix G G ⊤ .

Lemma 4.2. The matrix G G ⊤ ∈ R N × N is diagonal and its diagonal elements are given by

$$[ G ^ { \top } \mathbf G ] _ { k k } = \begin{cases} N ^ { 2 } & \text{if $k=1$,} \\ \frac { N ^ { 2 } } { 4 } ( 1 + \kappa ^ { N - 2 m } ) ^ { 2 } & \text{if $k=2m$ or $k=2m+1$ with $m=1,2,\dots,\frac{N-1}{2}$.} \\ - \cdots - \cdots - \cdots - \cdots - \cdots \end{cases}.$$

Proof. By the definition of G , we have

$$\mathbf G ^ { \top } \mathbf G = \left ( \sum _ { j = 1 } ^ { N } \psi _ { k } ( 1, \theta _ { j } ) \psi _ { l } ( 1, \theta _ { j } ) \right ) _ { k, l } \in \mathbb { R } ^ { N \times N },$$

and

$$& \text{and} \\ & ( 4. 3 ) \quad \psi _ { k } ( 1, \theta _ { j } ) = \begin{cases} \sqrt { N } & \text{if $k=1$,} \\ \sqrt { \frac { N } { 2 } } \left [ \sin m \theta _ { j } - \kappa ^ { N - 2 m } \sin ( N - m ) \theta _ { j } \right ] & \text{if $k=2m$, $m=1,2,\dots,\frac{N-1}{2}$,} \\ \sqrt { \frac { N } { 2 } } \left [ \cos m \theta _ { j } + \kappa ^ { N - 2 m } \cos ( N - m ) \theta _ { j } \right ] & \text{if $k=2m+1$, $m=1,2,\dots,\frac{N-1}{2}$.} \\ \end{cases}$$

Using an elementary relation (4.1), it can be easily shown that the following hold:

$$\sum _ { j = 1 } ^ { N } \sin m \theta _ { j } = 0, \ \ m \in \mathbb { Z },$$

$$\sum _ { j = 1 } ^ { N } \cos m \theta _ { j } = \begin{cases} N & \text{if $m\equiv 0$,} \\ 0 & \text{otherwise}, \end{cases}$$

$$( 4. 6 ) & \sum _ { j = 1 } ^ { j = 1 } \cos m \theta _ { j } \cos n \theta _ { j } = \begin{cases} \frac { N } { 2 } & \text{if $(m-n\equiv 0)\nmid(m+n\equiv 0)$,} \\ N & \text{if $(m-n\equiv 0)\wedge(m+n\equiv 0)$,} \\ 0 & \text{otherwise,} \end{cases}$$

̸

̸

$$( 4. 7 ) & \sum _ { j = 1 } ^ { N } \sin m \theta _ { j } \sin n \theta _ { j } = \begin{cases} \frac { N } { 2 } & \text{if $(m-n\equiv 0)\wedge(m+n\not\equiv 0$)}, \\ - \frac { N } { 2 } & \text{if $(m+n\equiv 0)\wedge(m-n\not\equiv 0$)}, \\ 0 & \text{otherwise}, \end{cases}$$

$$\sum _ { j = 1 } ^ { N } \sin m \theta _ { j } \cos n \theta _ { j } = 0, \ \ m, n \in \mathbb { Z },$$

where the symbol ⊻ denotes the exclusive or.

Off-diagonal elements can be computed as in the following manner. (1 2 , n ) and (1 , 2 n +1) elements for n = 1 2 , , . . . , ( N -1) / 2 are

$$[ { \mathbf G } ^ { \top } { \mathbf G } ] _ { 1, 2 n } = \frac { N } { \sqrt { 2 } } \sum _ { j = 1 } ^ { N } \left [ \sin n \theta _ { j } - \kappa ^ { N - 2 n } \sin ( N - n ) \theta _ { j } \right ] = 0,$$

$$[ G ^ { \top } G ] _ { 1, 2 n + 1 } = \frac { N } { \sqrt { 2 } } \sum _ { j = 1 } ^ { N } \left [ \cos n \theta _ { j } + \kappa ^ { N - 2 n } \cos ( N - n ) \theta _ { j } \right ] = 0$$

$$[ { \mathbf G } ^ { \top } { \mathbf G } ] _ { 1, 2 n } & = \frac { N } { \sqrt { 2 } } \sum _ { j = 1 } ^ { \prime \prime } \left [ \sin n \theta _ { j } - \kappa ^ { N - 2 n } \sin ( N - n ) \theta _ { j } \right ] = 0, \\ [ { \mathbf G } ^ { \top } { \mathbf G } ] _ { 1, 2 n + 1 } & = \frac { N } { \sqrt { 2 } } \sum _ { i = 1 } ^ { N } \left [ \cos n \theta _ { j } + \kappa ^ { N - 2 n } \cos ( N - n ) \theta _ { j } \right ] = 0$$

by (4.4) and (4.5). (2 m, 2 n ) elements for m = 1 2 , , . . . , ( N -1) / 2 and n = m +1 , . . . , ( N -1) / 2 are

$$\mathfrak { v y } \, \left ( \mathfrak { d } { - 1 } \right ) \, \text{and} \, \left ( \mathfrak { d } { \mathfrak { e } } \right ). \, \left ( \mathfrak { e } { m }, \, \mathfrak { e } { n } \right ) \, \text{enums to} \, m \, - 1, \, \mathfrak { e }, \, \dots, \, \left ( \mathfrak { v } \, - 1 \right ) / 2 \, \text{and} \, n \, - m \, + 1, \, \dots, \, \left ( \mathfrak { v } \, - 1 \right ) / 2 \, \text{ are} \\ [ G ^ { \top } G ] _ { 2 m, 2 n } & = \frac { N } { 2 } \sum _ { j = 1 } ^ { N } \left [ \sin m \theta _ { j } \sin n \theta _ { j } - \kappa ^ { N - 2 n } \sin m \theta _ { j } \sin ( N - n ) \theta _ { j } \\ & \quad - \kappa ^ { N - 2 m } \sin ( N - m ) \theta _ { j } \sin n \theta _ { j } + \kappa ^ { 2 ( N - ( m + n ) ) } \sin ( N - m ) \theta _ { j } \sin ( N - n ) \theta _ { j } \right ] \\ & = 0 \\ \mathfrak { d } { - 1 } \right. \, / \, \mathfrak { d } { \mathfrak { e } } \, \mathfrak { m } \, + 1 \, \mathfrak { m } \, \mathfrak { d } { \mathfrak { e } } \, \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m n } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } \, + \mathfrak { m m } = 0$$

by (4.7). (2 m, 2 n +1) elements for m = 1 2 , , . . . , ( N -1) / 2 and n = m,... , ( N -1) / 2 are

$$\L ^ { \diamondsuit \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon \\ [ G ^ { \top } G ] _ { 2 m, 2 n + 1 } & = \frac { N } { 2 } \sum _ { j = 1 } ^ { N } \left [ \sin m \theta _ { j } \cos n \theta _ { j } + \kappa ^ { N - 2 n } \sin m \theta _ { j } \cos ( N - n ) \theta _ { j } \\ & \quad - \kappa ^ { N - 2 m } \sin ( N - m ) \theta _ { j } \cos n \theta _ { j } - \kappa ^ { 2 ( N - ( m + n ) ) } \sin ( N - m ) \theta _ { j } \cos ( N - n ) \theta _ { j } \right ]$$

= 0

by (4.8). (2 m +1 2 , n +1) elements for m = 1 2 , , . . . , ( N -1) / 2 and n = m +1 , . . . , ( N -1) / 2 are

$$[ G ^ { \top } G ] _ { 2 m + 1, 2 n + 1 } = \frac { N } { 2 } \sum _ { j = 1 } ^ { N } [ \cos m \theta _ { j } \cos n \theta _ { j } + \kappa ^ { N - 2 n } \cos m \theta _ { j } \cos ( N - n ) \theta _ { j }$$

$$+ \kappa ^ { N - 2 m } \cos ( N - m ) \theta _ { j } \cos n \theta _ { j } + \kappa ^ { 2 ( N - ( m + n ) ) } \cos ( N - m ) \theta _ { j } \cos ( N - n ) \theta _ { j } \Big |$$

$$[ G ^ { \top } \text{G} ] _ { 2 m + 1, 2 n + 1 } = \frac { N } { 2 } \sum _ { j = 1 } \left [ \cos m \theta _ { j } \cos n \theta _ { j } + \kappa ^ { N - 2 n } \cos m \theta _ { j } \cos ( N - n ) \theta _ { j } \\ + \kappa ^ { N - 2 m } \cos ( N - m ) \theta _ { i } \cos n \theta _ { i } + \kappa ^ { 2 ( N - ( m + n ) ) } \cos ( N - m ) \theta _ { i } \cos ( N - n ) \theta _ { i } \right ]$$

by (4.6). In the same way in computing (2 m, 2 n +1) elements, we can show that (2 m +1 2 , n ) elements are equal to 0 for m = 1 2 , , . . . , ( N -1) / 2 and n = m + 1 , . . . , ( N -1) / 2. Summarizing the above, we have shown that upper triangular elements of G G ⊤ are equal to 0. Since G G ⊤ is symmetric, it follows that G G ⊤ is diagonal.

Concerning diagonal elements, we have by (4.6) and (4.7) that

$$\text{follows that $G^{T}G$ is diagonal.} \\ [ G ^ { T } G ] _ { 1 1 } = \sum _ { j = 1 } ^ { N } N = N ^ { 2 }, \\ [ G ^ { T } G ] _ { 2 m, 2 m } = \frac { N } { 2 } \sum _ { j = 1 } ^ { N } \left [ \sin ^ { 2 } m \theta _ { j } - 2 \kappa ^ { N - 2 m } \sin m \theta _ { j } \sin ( N - m ) \theta _ { j } + \kappa ^ { 2 ( N - 2 m ) } \sin ^ { 2 } ( N - m ) \theta _ { j } \right ] \\ = \frac { N ^ { 2 } } { 4 } ( 1 + \kappa ^ { N - 2 m } ) ^ { 2 }, \\ [ G ^ { T } G ] _ { 2 m + 1, 2 m + 1 } = \frac { N } { 2 } \sum _ { j = 1 } ^ { N } \left [ \cos ^ { 2 } m \theta _ { j } + 2 \kappa ^ { N - 2 m } \cos m \theta _ { j } \cos ( N - m ) \theta _ { j } + \kappa ^ { 2 ( N - 2 m ) } \cos ^ { 2 } ( N - m ) \theta _ { j } \right ] \\ = \frac { N ^ { 2 } } { 4 } ( 1 + \kappa ^ { N - 2 m } ) ^ { 2 } \\ \text{for $m=1,2,\dots,(N-1)/2$.} \quad \Box \\ \text{The above lemma shows that $G^{T}G$ is non-singular since all the eigenvalues of $G^{T}G$ are positive.}$$

for m = 1 2 , , . . . , ( N -1) / 2.

The above lemma shows that G G ⊤ is non-singular since all the eigenvalues of G G ⊤ are positive. Namely, Theorem 1.1 has been proved.

Proof of Theorem 1.2. It is well known that the relation

$$\text{cond} _ { 2 } ( \mathbf A ) = \sqrt { \frac { \lambda _ { \max } ( \mathbf A ^ { \top } \mathbf A ) } { \lambda _ { \min } ( \mathbf A ^ { \top } \mathbf A ) } }$$

holds, where λ max ( A A ⊤ ) and λ min ( A A ⊤ ) are the maximal and minimal (by moduli) eigenvalues of A A ⊤ , respectively. Lemma 4.2 shows that

$$\lambda _ { \max } ( G ^ { \top } G ) = N ^ { 2 }, \ \lambda _ { \min } ( G ^ { \top } G ) = \frac { N ^ { 2 } } { 4 } ( 1 + \kappa ^ { N - 2 } ) ^ { 2 }$$

hold. Therefore, the relation (4.9) yields that

$$\text{cond} _ { 2 } ( \mathbf G ) = \sqrt { \frac { \lambda _ { \max } ( \mathbf G ^ { \top } \mathbf G ) } { \lambda _ { \min } ( \mathbf G ^ { \top } \mathbf G ) } } } = \frac { 2 } { 1 + \kappa ^ { N - 2 } } \longrightarrow 2 \ ( N \to \infty ).$$

$$<_Python_> 
\infty
                                                                                                                                                                                                        
                                                                                                                                                                                                      
\$$

The above corollary implies that the coefficient matrix G in the collocation equation (3.3) is wellconditioned, which is a quite difference compared with the original DSM.

4.3. Error estimate. Since G G ⊤ is diagonal by Lemma 4.2, the solution to (4.2) is explicitly given by

$$Q _ { k } & = \frac { 1 } { [ G ^ { \top } G ] _ { k k } } \sum _ { j = 1 } ^ { N } \psi _ { k } ( 1, \theta _ { j } ) f ( \rho e ^ { i \theta _ { j } } ), \quad k = 1, 2, \dots, N.$$

To derive the error estimate, we write down the exact solution u to the problem (1.3) using the Fourier series expansion as

$$u ( r \mathrm e ^ { \mathrm i \theta } ) = \sum _ { n \in \mathbb { Z } } f _ { n } \left ( \frac { r } { \rho } \right ) ^ { | n | } \mathrm e ^ { \mathrm i n \theta },$$

↦

where f n denotes the n -th Fourier coefficient of the function θ → f ( ρ e i θ ). We next represent u ( N ) as an infinite series using { f n } n ∈ Z . The concrete expression (4.10) of coefficients { Q k } N k =1 yields that

$$\sl i t e s e r s i u n \{ f _ { n } \} _ { n \in \mathbb { Z } }. \text{ The concrete expression } ( 4. 1 0 ) \text{ of } \text{ coefficients } \{ Q _ { k } \} _ { k = 1 } ^ { n } \text{ yields } \\ u ^ { ( N ) } ( r e ^ { i \theta } ) & = \sum _ { k = 1 } ^ { N } \left [ \frac { 1 } { [ G ^ { \top } G ] _ { k } } \sum _ { j = 1 } ^ { N } \psi _ { k } ( 1, \theta _ { j } ) f ( \rho e ^ { i \theta _ { j } } ) \right ] \psi _ { k } ( r / \rho, \theta ) \\ & = \sum _ { k = 1 } ^ { N } \frac { 1 } { [ G ^ { \top } G ] _ { k } } \sum _ { j = 1 } ^ { N } \psi _ { k } ( 1, \theta _ { j } ) \left ( \sum _ { n \in \mathbb { Z } } f _ { n } e ^ { i n \theta _ { j } } \right ) \psi _ { k } ( r / \rho, \theta ) \\ & = \sum _ { n \in \mathbb { Z } } f _ { n } \sum _ { k = 1 } ^ { N } \frac { 1 } { [ G ^ { \top } G ] _ { k } } \left ( \sum _ { j = 1 } ^ { N } \psi _ { k } ( 1, \theta _ { j } ) e ^ { i n \theta _ { j } } \right ) \psi _ { k } ( r / \rho, \theta ). \\ \text{a 4.3.} \text{ It holds that }$$

Lemma 4.3. It holds that

$$\alpha _ { n } ^ { ( N ) } ( r / \rho, \theta ) & \coloneqq \sum _ { k = 1 } ^ { N } \left ( \frac { 1 } { [ G ^ { \top } G ] _ { k } } \sum _ { j = 1 } ^ { N } \psi _ { k } ( 1, \theta _ { j } ) e ^ { i n \theta _ { j } } \right ) \psi _ { k } ( r / \rho, \theta ) \\ & \quad = \begin{cases} 1 & \text{if $n\equiv 0$,} \\ \frac { 1 } { 1 + \kappa ^ { N - 2 m } } \left ( w ^ { m } + \kappa ^ { N - 2 m } \overline { w } ^ { N - m } \right ) & \text{if $n\equiv m$,} \\ \frac { 1 } { 1 + \kappa ^ { N - 2 m } } \left ( \overline { w } ^ { m } + \kappa ^ { N - 2 m } w ^ { N - m } \right ) & \text{if $n\equiv -m$} \end{cases} \\ = 1, \dots, ( N - 1 ) / 2, \, \text{where $we$ set $w\coloneqq (r/\rho)e^{i\theta}$.}$$

with m = 1 , . . . , ( N -1) / 2 , where we set w := ( r/ρ )e i θ .

Proof. We first compute ∑ N j =1 ψ k (1 , θ j )e i nθ j . The elementary relation (4.1) and expression (4.3) yield that

$$\sum _ { j = 1 } ^ { N } \psi _ { 1 } ( 1, \theta _ { j } ) = \sqrt { N } \sum _ { j = 1 } ^ { N } \mathrm e ^ { i n \theta _ { j } } = \begin{cases} N ^ { 3 / 2 } & \text{if $n\equiv 0$,} \\ 0 & \text{otherwise} \end{cases}$$

for k = 1,

$$\text{for } k = 1, \\ & \sum _ { j = 1 } ^ { N } \psi _ { 2 m } ( 1, \theta _ { j } ) e ^ { i n \theta _ { j } } \\ = \sqrt { \frac { N } { 2 } } \sum _ { j = 1 } ^ { N } \left ( \frac { e ^ { i ( m + n ) \theta _ { j } } - e ^ { - i ( m - n ) \theta _ { j } } } { 2 i } - \kappa ^ { N - 2 m } \frac { e ^ { i ( N - ( m - n ) ) \theta _ { j } } - e ^ { - i ( N - ( m + n ) ) \theta _ { j } } } { 2 i } \right ) \\ = \text{i} \left ( \frac { N } { 2 } \right ) ^ { 3 / 2 } \times \begin{cases} 1 + \kappa ^ { N - 2 m } & \text{if } n \equiv m, \\ - ( 1 + \kappa ^ { N - 2 m } ) & \text{if } n \equiv - m, \\ 0 & \text{otherwise} \end{cases}$$

for k = 2 m with m = 1 , . . . , ( N -1) / 2, and

$$& \text{for $k=2m$ with $m=1,\dots,(N-1)/2$, and} \\ & \quad \sum _ { j = 1 } ^ { N } \psi _ { 2 m + 1 } ( 1, \theta _ { j } ) e ^ { i n \theta _ { j } } \\ & \quad = \sqrt { \frac { N } { 2 } } \sum _ { j = 1 } ^ { N } \left ( \frac { e ^ { i ( m + n ) \theta _ { j } } + e ^ { - i ( m - n ) \theta _ { j } } } { 2 } + \kappa ^ { N - 2 m } \frac { e ^ { i ( N - ( m - n ) ) \theta _ { j } } + e ^ { - i ( N - ( m + n ) ) \theta _ { j } } } { 2 } \right ) \\ & \quad = \left ( \frac { N } { 2 } \right ) ^ { 3 / 2 } \times \begin{cases} 1 + \kappa ^ { N - 2 m } & \text{if $(n\equiv m) \vee (n\equiv -m)$,} \\ 0 & \text{otherwise} \end{cases} \\ & \text{for $k=2m+1$ with $m=1,\dots,(N-1)/2$. Therefore, it follows from the above result and Lemma 4.2$}$$

for k = 2 m +1 with m = 1 , . . . , ( N -1) / 2. Therefore, it follows from the above result and Lemma 4.2 that

$$\overset { N } { \underset { k = 1 } { \sum } } \frac { 1 } { [ { \mathbf G } ^ { \top } { \mathbf G } ] _ { k } } \left ( \overset { N } { \underset { j = 1 } { \sum } } \psi _ { k } ( 1, \theta _ { j } ) e ^ { i n \theta _ { j } } \right ) \psi _ { k } ( r / \rho, \theta ) = \frac { 1 } { [ { \mathbf G } ^ { \top } { \mathbf G } ] _ { 1 } } \cdot N ^ { 3 / 2 } \psi _ { 1 } ( r / \rho, \theta ) = 1$$

for n ≡ 0,

$$\text{for } n \equiv 0, \\ & \sum _ { k = 1 } ^ { N } \frac { 1 } { [ G ^ { \top } G ] _ { k } } \left ( \sum _ { j = 1 } ^ { N } \psi _ { k } ( 1, \theta _ { j } ) e ^ { i n \theta _ { j } } \right ) \psi _ { k } ( r / \rho, \theta ) \\ & = \left ( \frac { N } { 2 } \right ) ^ { 3 / 2 } ( 1 + \kappa ^ { N - 2 m } ) \left ( \frac { \mathbf i } { [ G ^ { \top } G ] _ { 2 m } } \psi _ { 2 m } ( r / \rho, \theta ) + \frac { 1 } { [ G ^ { \top } G ] _ { 2 m + 1 } } \psi _ { 2 m + 1 } ( r / \rho, \theta ) \right ) \\ & = \frac { 1 } { 1 + \kappa ^ { N - 2 m } } \left ( w ^ { m } + \kappa ^ { N - 2 m } \overline { w } ^ { N - m } \right ) \\ \text{for } n \equiv m, \text{ and }$$

for n ≡ m , and

$$\text{for } n \equiv m, \text{ and } \\ \sum _ { k = 1 } ^ { N } \frac { 1 } { [ G ^ { \top } G ] _ { k } } \left ( \sum _ { j = 1 } ^ { N } \psi _ { k } ( 1, \theta _ { j } ) e ^ { i n \theta _ { j } } \right ) \psi _ { k } ( r / \rho, \theta ) \\ = \left ( \frac { N } { 2 } \right ) ^ { 3 / 2 } ( 1 + \kappa ^ { N - 2 m } ) \left ( \frac { - 1 } { [ G ^ { \top } G ] _ { 2 m } } \psi _ { 2 m } ( r / \rho, \theta ) + \frac { 1 } { [ G ^ { \top } G ] _ { 2 m + 1 } } \psi _ { 2 m + 1 } ( r / \rho, \theta ) \right ) \\ = \frac { 1 } { 1 + \kappa ^ { N - 2 m } } \left ( \overline { w } ^ { m } + \kappa ^ { N - 2 m } w ^ { N - m } \right ) \\ \text{Using the functions } \varpi _ { \cdot } ^ { ( N ) } \text{ defined in Lemma } 4. 3. \text{ we can represent } u ^ { ( N ) } ( r e ^ { i \theta } ) \text{ as}$$

for n ≡ -m .

Using the functions φ ( N ) n defined in Lemma 4.3, we can represent u ( N ) ( r e i θ ) as

$$u ^ { ( N ) } ( r \mathrm e ^ { \mathrm i \theta } ) = \sum _ { n \in \mathbb { Z } } f _ { n } \varphi _ { n } ^ { ( N ) } ( r / \rho, \theta ).$$

Since u and u ( N ) are harmonic in B ρ and continuous on B ρ , the maximum principle gives the following bound:

$$\| u - u ^ { ( N ) } \| _ { L ^ { \infty } ( \mathbb { B } _ { \rho } ) } \leq \sum _ { n \in \mathbb { Z } } | f _ { n } | g _ { n } ^ { ( N ) },$$

$$\| u - u ^ { ( N ) } \| _ { L ^ { \infty } ( \mathbb { B } _ { \rho } ) } \leq \sum _ { n \in \mathbb { Z } } | f _ { n } | g _ { n } ^ { ( N ) },$$

where

$$g _ { n } ^ { ( N ) } \coloneqq \sup _ { \theta \in \mathbb { R } } \left | e ^ { i n \theta } - \varphi _ { n } ^ { ( N ) } ( 1, \theta ) \right |.$$

To derive the desired error estimate, the local and global estimates on g ( N ) n are required.

$$\text{Lemma } 4. 4. \quad \text{ (i) } \ F o r \ a n y \ n \in \mathbb { N }, \ i t \ h o l d s \ t h a t$$

$$g _ { n } ^ { ( N ) } \leq 2.$$

## (ii) It holds that

$$g _ { 0 } ^ { ( 0 ) } = 0.$$

Moreover, for n = 1 , . . . , ( N -1) / 2 , it holds that

$$g _ { n } ^ { ( N ) } \leq 2 \kappa ^ { N - 2 n }.$$

Proof. (i) Lemma 4.3 yields that

$$\cdot \vee \nu \text{mm} \colon \infty \nu \text{raw} \\ \varphi _ { n } ^ { ( N ) } ( 1, \theta ) = \begin{cases} 1 & \text{if $n\equiv 0$,} \\ \frac { 1 } { 1 + \kappa ^ { N - 2 m } } \left ( e ^ { i m \theta } + \kappa ^ { N - 2 m } e ^ { - i ( N - m ) \theta } \right ) & \text{if $n\equiv m$,} \\ \frac { 1 } { 1 + \kappa ^ { N - 2 m } } \left ( e ^ { - i m \theta } + \kappa ^ { N - 2 m } e ^ { i ( N - m ) \theta } \right ) & \text{if $n\equiv -m$} \\ \cdot \quad \nu \text{raw} \cdot \infty \nu \colon \infty \nu \colon \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \, ( N ) \, \ast \, \omega \, \ast \, e \end{cases}$$

with m = 1 , . . . , ( N -1) / 2. Therefore, the triangle inequality implies that φ ( N ) n (1 , θ ) ≤ 1 for any n ∈ Z and θ ∈ R . We thus obtain

$$g _ { n } ^ { ( N ) } \leq \sup _ { \theta \in \mathbb { R } } \left ( | \mathrm e ^ { \mathrm i n \theta } | + | \varphi _ { n } ^ { ( N ) } ( 1, \theta ) | \right ) \leq 2.$$

(ii) Since φ ( N ) 0 (1 , θ ) = 1 as shown in the above, we readily see that g ( N ) 0 = 0 holds. For n = 1 , . . . , ( N -1) / 2, we have

$$\Gamma ^ { \prime \prime } - \varphi _ { n } ^ { ( N ) } ( 1, \theta ) & = e ^ { i n \theta } - \frac { 1 } { 1 + \kappa ^ { N - 2 n } } \left ( e ^ { i n \theta } + \kappa ^ { N - 2 n } e ^ { - i ( N - n ) \theta } \right ) \\ & = \frac { \kappa ^ { N - 2 n } } { 1 + \kappa ^ { N - 2 n } } \left ( e ^ { i n \theta } - e ^ { - i ( N - n ) \theta } \right ). \\ \intertext { a t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d j u t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d I t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d l i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t j u t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i ni d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i in d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i, i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d it i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i c t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d e x \, \quad \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \gamma } \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \amma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \\ \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \implement \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \Gamma \paramma \Gamma \Gamma \Gamma \$$

Hence, the triangle inequality yields that

$$\left | \text{e} ^ { i n \theta } - \varphi _ { n } ^ { ( N ) } ( 1, \theta ) \right | \leq \frac { 2 \kappa ^ { N - 2 n } } { 1 + \kappa ^ { N - 2 n } } \leq 2 \kappa ^ { N - 2 n },$$

which is the desired estimate.

□

Proof of Theorem 1.3. (i), (ii) Note that φ ( N ) -n = φ ( N ) n holds, which implies that g ( N ) -n = g ( N ) n . Based on this relation, we decompose the right-hand side of (4.11) as follows:

$$\sum _ { n \in \mathbb { Z } } | f _ { n } | g _ { n } ^ { ( N ) } = | f _ { 0 } | g _ { 0 } ^ { ( N ) } + \sum _ { n = 1 } ^ { \lfloor N / 3 \rfloor - 1 } ( | f _ { n } | + | f _ { - n } | ) g _ { n } ^ { ( N ) } + \sum _ { n = \lfloor N / 3 \rfloor } ^ { \infty } ( | f _ { n } | + | f _ { - n } | ) g _ { n } ^ { ( N ) } \equiv I + \amalg \amalg.$$

Since g ( N ) 0 = 0 by Lemma 4.4(ii), we have I = 0. Concerning II, since | f n | ≤ ∥ f ∥ L ∞ ( ∂ B ρ ) holds, Lemma 4.4(ii) implies that

$$\text{II} \leq 4 \| f \| _ { L ^ { \infty } ( \partial \mathbb { B } _ { \rho } ) } \sum _ { n = 1 } ^ { \lfloor N / 3 \rfloor - 1 } \kappa ^ { N - 2 n } \leq 4 \| f \| _ { L ^ { \infty } ( \partial \mathbb { B } _ { \rho } ) } \kappa ^ { N } \cdot \frac { \kappa ^ { - 2 \lfloor N / 3 \rfloor } } { \kappa ^ { - 2 } - 1 } \leq \frac { 4 \| f \| _ { L ^ { \infty } ( \partial \mathbb { B } _ { \rho } ) } } { \kappa ^ { - 2 } - 1 } \kappa ^ { \frac { N } { 3 } },$$

Lemma 4.4(i) yields that

$$\text{III} \leq 2 \sum _ { n = \lfloor N / 3 \rfloor } ^ { \infty } ( | f _ { n } | + | f _ { - n } | ) = \begin{cases} o ( 1 ) & \text{for case (i)}, \\ O ( N ^ { - \alpha + 1 } ) & \text{for case (ii)}. \end{cases}$$

Summarizing the above, we complete the proof of (i) and (ii).

(iii) We decompose the right-hand side of (4.11) as

$$\sum _ { n \in \mathbb { Z } } | f _ { n } | g _ { n } ^ { ( N ) } = | f _ { 0 } | g _ { 0 } ^ { ( N ) } + \sum _ { n = 1 } ^ { ( N - 1 ) / 2 } ( | f _ { n } | + | f _ { - n } | ) g _ { n } ^ { ( N ) } + \sum _ { n = ( N + 1 ) / 2 } ^ { \infty } ( | f _ { n } | + | f _ { - n } | ) g _ { n } ^ { ( N ) } \equiv I + \text{II} + \text{III}.$$

Figure 1. Plots of the L ∞ -norms of errors and the condition numbers. (Left) L ∞ -error; (Right) condition number when Ω is a disk region.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00212v1-figure-32ebaae2-d34c-4994-9411-6b0b0f9fbb20.jpg)

Since g (0) 0 = 0 by Lemma 4.4(ii), we have I = 0. The assumption | f n | = O( a | n | ) ( | n | →∞ ) assures the existence of positive constant C satisfying | f n | ≤ Ca | n | for all n ∈ Z . Combining this relation with Lemma 4.4(ii), we have

$$\text{$\text{$columns $k$} } \text{$\text{$\text{$\text{$k$} } } / \text{$\text{$\text{$k$} } } \text{$\text{$\text{$k$} } } \text{$\text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\quad \text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$.} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \quad \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$k$} } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$j$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$,} } \text{$k$} } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$j$} } } \text{\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \text{$\text{$k$} } } \end{$$

where we have used

$$\sum _ { n = 1 } ^ { m } \tau ^ { n } \leq \begin{cases} \frac { \tau ^ { m + 1 } } { \tau - 1 } & \text{if $\tau>1$}, \\ m & \text{if $\tau=1$}, \\ \frac { 1 } { 1 - \tau } & \text{if $0<\tau<1$}. \end{cases} \\ \omega + \cdots \sim \int \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimitimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { N = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolimits _ { n = 1 } ^ { n } \nolipseq \nolimseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nulseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \ncolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolorseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nocolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nvolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nalseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nolseq \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodsc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodcss \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \noddesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodisc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodeless \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodead \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodDesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \nodesc \n$$

Concerning III, the assumption | f n | = O( a | n | ) ( n →∞ ) and Lemma 4.4 imply that

$$\text{III} \leq 4 C \sum _ { n = ( N + 1 ) / 2 } ^ { \infty } a ^ { n } = \frac { 4 C } { 1 - a } a ^ { ( N + 1 ) / 2 }.$$

Combining the above estimates, we complete the proof of (iii).

□

## 5. Numerical experiments

In this section, we show several results of numerical experiments that compare the original MFS and MFS-QR.

5.1. Disk. We first consider the problem (1.3) with boundary value f ( x, y ) = x y 2 3 . We set parameters as ρ := 1 and R := 1 1 or . R := 1 5. . We change N as N = 2 l + 1 with l = 1 2 , , . . . , l max := 150, and compute condition numbers and approximation errors. Figure 1 shows the result of numerical experiments. It can be observed that the L ∞ -norms of approximation errors decay exponentially with respect to N for both the DSM and the DSM-QR. The condition number increases exponentially for the DSM, while it is of O(1) for the DSM-QR.

5.2. Jordan region. We next deal with the case where Ω is a Jordan region. In the previous study [2], it was pointed out that the MFS-QR reduces the ill-conditioning but cannot remove it completely when the problem region Ω is not a disk. In this subsection, we develop an analogue of the MFS-QR for Jordan region with analytic boundary using conformal mappings.

Since Ω is a Jordan region, the Riemann mapping theorem assures the existence of a conformal map Ψ: B 1 → Ω. Then, the Osgood-Carath´ eodory theorem extends Ψ to a homeomorphism of Ω onto B 1 , denoted by the same symbol Ψ. Since ∂ Ω is analytic, we can extend Ψ conformally to B τ with some τ &gt; 1 by the Schwarz reflection principle. As studied in [11], for the MFS, if we arrange the collocation points { z j } N j =1 ∈ ∂ Ω and { ζ k } N k =1 ⊂ C \ Ω by

$$j$$

$$z _ { j } & \coloneqq \Psi ( \omega ^ { j } ), \quad j = 1, 2, \dots, N, \\ \zeta _ { k } & \coloneqq \Psi ( R \omega ^ { k } ), \quad k = 1, 2, \dots, N$$

with R ∈ (1 , τ ) and ω := exp(2 π /N i ), an approximate solution uniquely exists for sufficiently large N , and an approximation error decays exponentially for analytic boundary data. The essense of the proof is to decompose the basis function into the disk part and the perturbation part as follows:

$$\log | \Psi ( r e ^ { i \theta } ) - \Psi ( R e ^ { i \phi } ) | = \log \left | \frac { \Psi ( r e ^ { i \theta } ) - \Psi ( R e ^ { i \phi } ) } { r e ^ { i \theta } - R e ^ { i \phi } } \right | + \log | r e ^ { i \theta } - R e ^ { i \phi } |.$$

Concerning the DSM, the author developed mathematical theory in [20]. The similar idea plays an essential role, and the basis function is decomposed as

$$- \text{Re} \left ( \frac { \nu ( \Psi ( R e ^ { i \phi } ) ) } { \Psi ( r e ^ { i \theta } ) - \Psi ( R e ^ { i \phi } ) } \right ) = - \text{Re} \left ( \frac { \nu ( \Psi ( R e ^ { i \phi } ) ) } { \Psi ( r e ^ { i \theta } ) - \Psi ( R e ^ { i \phi } ) } - \frac { \nu ( R e ^ { i \phi } ) } { r e ^ { i \theta } - R e ^ { i \phi } } \right ) - \text{Re} \left ( \frac { \nu ( R e ^ { i \phi } ) } { r e ^ { i \theta } - R e ^ { i \phi } } \right ),$$

where ν R ( e i ϕ ) denotes the unit outward normal vector to the circle ∂ B R at the point R e i ϕ , whereas ν (Ψ( R e i ϕ )) does the one to the Jordan curve Ψ( ∂ B R ) at Ψ( R e i ϕ ). The first-term in the right-hand side can be regarded as a good perturbation in some sense. The second-term in the right-hand side is the basis function for the case of disk, so it can be modified via the DSM-QR. Hence, it is natural to imagine that the set of functions

$$\dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots$$

would be a new basis function, which removes the ill-conditioning completely. Although there is no mathematical theory concerning the above idea, we below show several numerical experiments that exemplify the effectiveness of the above strategy.

Example 5.1 . We first show the result where the conformal mapping is given by a polynomial:

$$\Psi ( z ) \coloneqq \frac { 4 } { 5 } z + \frac { 1 } { 1 0 } z ^ { 5 }.$$

We use the same boundary value f ( x, y ) = x y 2 3 , and set parameters as R := 1 05, . N = 2 l +1, and l max := 500. Figure 3 demonstrates results of numerical experiments. As in the case of disk, we can observe that the approximation error decays exponentially for both the DSM and the DSM-QR. Concerning the condition number, it increases exponentially for the DSM, while it is of O(1) for the DSM-QR. Namely, the aforementioned strategy works well.

Example 5.2 . To strengthen the validity of the above idea, another numerical experiment will be performed. We consider the Joukowski-type transformation, defined as

$$\Psi ( z ) \coloneqq w + \frac { 1 } { 2 w }, \ \ w \coloneqq z - \frac { 1 } { 5 } + \frac { i } { 5 }.$$

We use the same boundary value f ( x, y ) = x y 2 3 , and set parameters as R := 1 1, . N = 2 l + 1, and l max := 500. In this case, the curvature of boundary ∂ Ω is very large at the right end, and it is generally difficult to obtain a highly accurate approximate solution. However, as shown in Figure 5, DSM and DSM-QR show exponential decay of approximation error. This indicates that DSM and

Figure 2. Arrangements of the collocation points, the singular points, and the dipole moments in Example 5.1.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00212v1-figure-db845142-4393-4066-8df2-d214f3e6e715.jpg)

Figure 3. Plots of the L ∞ -norms of errors and the condition numbers in Example 5.1. (Left) L ∞ -error; (Right) condition number.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00212v1-figure-e77fd5ae-a3dd-4f51-a44e-4f6b64484214.jpg)

DSM-QR are even more powerful in regions with complex geometry. Moreover, even in this case, the ill-conditionality is completely eliminated in the DSM-QR, which truly demonstrates the high usefulness of the DSM-QR.

## 6. Concluding remarks

This paper puts forth a new approach, designated as the DSM-QR, which draws inspiration from the MFS-QR. It then proceeds to undertake a comprehensive mathematical analysis of this proposed approach within the context of the disk region. Theorems 1.1 and 1.3 demonstrate that the approximate solution exists uniquely and that the approximation error decays in accordance with the expected order in the context of the MFS. A noteworthy outcome is the mathematical proof presented in Theorem 1.2, which demonstrates that the condition number is of O(1). Moreover, despite the absence

Figure 4. Arrangements of the collocation points, the singular points, and the dipole moments in Example 5.2.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00212v1-figure-e0cc5677-4c0e-4d6b-b89b-6656cf02b286.jpg)

Figure 5. Plots of the L ∞ -norms of errors and the condition numbers in Example 5.2. (Left) L ∞ -error; (Right) condition number.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00212v1-figure-5223c769-4b1b-46b0-98fc-7cfca9705d6d.jpg)

of a mathematical analysis, we were able to extend the DSM-QR to the general Jordan region through the use of conformal mapping. In a previous study [3], a method based on singular value decomposition was shown to be effective in removing ill-conditionality in star-shaped domains. However, the main feature of the method presented in this paper is that it is based on conformal mapping, which allows for domain shapes to be unrestricted.

Further research should be conducted in the following areas: Firstly, the mathematical justification of the extension to Jordan regions should be considered. This can be achieved by applying the mathematical analysis methods of the DSM in the Jordan domain [20]. An additional avenue for exploration is the expansion of the problem region in terms of its dimensionality. This is a direction that has been pursued by the MFS-QR and the MFS-SVD as well, but thus far, these methods have only demonstrated the ability to remove ill-conditionality in planar regions. They have not yet succeeded in extending this approach to domains of greater dimensionality, such as three-dimensional regions. Moreover, it would be advantageous to consider extensions to problems in multiply-connected domains, even if they are in two-dimensional regions, for practical applications. Finally, since the MFS-QR has

been developed for the Helmholtz equation [1], it is crucial to investigate whether the DSM-QR can be extended to other equations as well.

Acknowledgements. This work was supported by JSPS KAKENHI Grant Numbers JP18K13455, JP22K03425, JP22K18677, JP23H00086.

## References

- [1] P. R. S. Antunes , A numerical algorithm to reduce ill-conditioning in meshless methods for the Helmholtz equation , Numer. Algorithms, 79 (2018), pp. 879-897.
- [2] , Reducing the ill conditioning in the method of fundamental solutions , Adv. Comput. Math., 44 (2018), pp. 351-365.
- [3] , A well-conditioned method of fundamental solutions for Laplace equation , Numer. Algorithms, 91 (2022), pp. 1381-1405.
- [4] A. H. Barnett and T. Betcke , Stability and convergence of the method of fundamental solutions for Helmholtz problems on analytic domains , J. Comput. Phys., 227 (2008), pp. 7003-7026.
- [5] H. Brand´ en, S. Holmgren, and P. Sundqvist , Discrete fundamental solution preconditioning for hyperbolic systems of PDE , J. Sci. Comput., 30 (2007), pp. 35-60.
- [6] H. Brand´ en and P. Sundqvist , Preconditioners based on fundamental solutions , BIT, 45 (2005), pp. 481-494.
- [7] S.-I. Ei, H. Ochiai, and Y. Tanaka , Method of fundamental solutions for Neumann problems of the modified Helmholtz equation in disk domains , J. Comput. Appl. Math., 402 (2022), pp. Paper No. 113795, 27.
- [8] G. B. Folland , Introduction to partial differential equations , Princeton University Press, Princeton, NJ, second ed., 1995.
- [9] B. Fornberg and C. Piret , A stable algorithm for flat radial basis functions on a sphere , SIAM J. Sci. Comput., 30 (2007/08), pp. 60-80.
- [10] M. Katsurada , A mathematical study of the charge simulation method. II , J. Fac. Sci. Univ. Tokyo Sect. IA Math., 36 (1989), pp. 135-162.
- [11] , Asymptotic error analysis of the charge simulation method in a Jordan region with an analytic boundary , J. Fac. Sci. Univ. Tokyo Sect. IA Math., 37 (1990), pp. 635-657.
- [12] , Charge simulation method using exterior mapping functions , Japan J. Indust. Appl. Math., 11 (1994), pp. 47-61.
- [13] M. Katsurada and H. Okamoto , A mathematical study of the charge simulation method. I , J. Fac. Sci. Univ. Tokyo Sect. IA Math., 35 (1988), pp. 507-518.
- [14] M. Katsurada and H. Okamoto , The collocation points of the fundamental solution method for the potential problem , Comput. Math. Appl., 31 (1996), pp. 123-137.
- [15] Y. Kazashi and M. Sugihara , Error estimation for the invariant scheme of charge simulation method on a disc with scattered points , Jpn. J. Ind. Appl. Math., 31 (2014), pp. 305-321.
- [16] X. Li , On convergence of the method of fundamental solutions for solving the Dirichlet problem of Poisson's equation , Adv. Comput. Math., 23 (2005), pp. 265-277.
- [17] , Convergence of the method of fundamental solutions for Poisson's equation on the unit sphere , Adv. Comput. Math., 28 (2008), pp. 269-282.
- [18] Z.-C. Li , The method of fundamental solutions for annular shaped domains , J. Comput. Appl. Math., 228 (2009), pp. 355-372.

[19]

K. Murota

,

Comparison of conventional and 'invariant' schemes of fundamental solutions method for annular domains

,

Japan J. Indust. Appl. Math., 12 (1995), pp. 61-85.

- [20] K. Sakakibara , Analysis of the dipole simulation method for two-dimensional Dirichlet problems in Jordan regions with analytic boundaries , BIT, 56 (2016), pp. 1369-1400.
- [21] , Asymptotic analysis of the conventional and invariant schemes for the method of fundamental solutions applied to potential problems in doubly-connected regions , Jpn. J. Ind. Appl. Math., 34 (2017), pp. 177-228.
- [22] , Method of fundamental solutions for biharmonic equation based on Almansi-type decomposition , Appl. Math., 62 (2017), pp. 297-317.
- [23] Y.-S. Smyrlis , The method of fundamental solutions: a weighted least-squares approach , BIT, 46 (2006), pp. 163194.
- (K. Sakakibara) Faculty of Mathematics and Physics, Institute of Science and Engineering, Kanazawa University, Kakuma-machi, Kanazawa-shi, Ishikawa 920-1192, Japan; RIKEN iTHEMS, 2-1 Hirosawa, Wakoshi, Saitama 351-0198, Japan

Email address : ksakaki@se.kanazawa-u.ac.jp
Koya Sakakibara
Well-conditioned dipole-type method of fundamental solutions: derivation and its mathematical analysis
In this paper, we examine the dipole-type method of fundamental solutions,
which can be conceptualized as a discretization of the "singularity-removed"
double-layer potential. We present a method for removing the
ill-conditionality, which was previously considered a significant challenge,
and provide a mathematical analysis in the context of disk regions. Moreover,
we extend the proposed method to the general Jordan region using conformal
mapping, demonstrating the efficacy of the proposed method through numerical
experiments.