Next Article in Journal
About the Resolvent Kernel of Neutral Linear Fractional System with Distributed Delays
Next Article in Special Issue
Dual Variational Formulations for a Large Class of Non-Convex Models in the Calculus of Variations
Previous Article in Journal
Anti-Disturbance Fault-Tolerant Constrained Consensus for Time-Delay Faulty Multi-Agent Systems with Semi-Markov Switching Topology
Previous Article in Special Issue
One-Step Suicide Substrate Inactivation Kinetics of a Ping-Pong Reaction with One Substrate Undergoing Disproportionation: A Theoretical Approach with Approximate Solutions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solving Poisson Equations by the MN-Curve Approach

Department of Data Science, Providence University, Shalu, Taichung 43301, Taiwan
Mathematics 2022, 10(23), 4582; https://doi.org/10.3390/math10234582
Submission received: 25 October 2022 / Revised: 21 November 2022 / Accepted: 28 November 2022 / Published: 2 December 2022
(This article belongs to the Special Issue Numerical Analysis and Scientific Computing II)

Abstract

:
In this paper, we adopt the choice theory of the shape parameters contained in the smooth radial basis functions to solve Poisson equations. Luh’s choice theory, based on harmonic analysis, is mathematically complicated and applies only to function interpolation. Here, we aim at presenting an easily accessible approach to solving differential equations with the choice theory which proves to be very successful, not only by its easy accessibility but also by its striking accuracy and efficiency. Our emphases are on the highly reliable prediction of the optimal value of the shape parameter and the extremely small approximation errors of the numerical solutions to the differential equations. We hope that our approach can be accepted by both mathematicians and non-mathematicians.
MSC:
31A30; 35J05; 35J25; 35J67; 35Q40; 35Q70; 65D05; 65L10; 65N35

1. Introduction

Here, we focus on the generalized multiquadrics
ϕ ( x ) : = ( 1 ) β / 2 ( c 2 + x 2 ) β / 2 , β R \ 2 N 0 , c > 0 , x R d ,
where β / 2 denotes the smallest integer greater than or equal to β / 2 . These are the most popular radial basis functions (RBFs) and are frequently used in the collocation method of solving partial differential equations. The choice of the shape parameter c contained in ϕ ( x ) has been considered to be a hard and important question in this field for decades and often leads to giving up this approach. Hitherto, there is no strict theory about its optimal choice when dealing with partial differential equations (PDEs). Although Luh’s theory, called the c theory by E. Kansa, can predict its optimal value almost exactly, it applies mainly to function interpolation only and involves the complicated theory of harmonic analysis. Scientists, especially non-mathematicians, still do not know how to choose it when solving PDEs with RBFs.
As the inventor of the choice theory, the author knows that the theory applies to PDEs as well, maybe with a moderate search when necessary. The main reason is that collocation is, in spirit, a kind of interpolation. Moreover, Dirichlet conditions do offer interpolation points on the boundary. As can be seen in Luh [1,2], when c is chosen according to the MN curves, the accuracy of the function approximation is incredibly good, both in theory and practice. It is not hard to imagine that the combination of the c theory and collocation may lead to very good results in the field of numerical PDEs.
In this paper, we follow Kansa’s route [3,4,5] to make collocation for solving Poisson equations but in a totally different way of choosing c. Basically, we discard the traditional time-consuming trial-and-error search and adopt the theoretically predicted optimal value of c. It greatly increases the time efficiency of finding a suitable c because it avoids the traditional way of solving a linear system for each trial of the c value. Experiments show that such a c does produce a very good numerical solution to the PDE, even if the value of c is not exactly the experimentally optimal one. If one insists on finding the experimentally optimal value, it can be achieved by a moderate search. The stopping criterion is totally different from Kansa’s approach and perhaps has never appeared in the literature. Our stopping criterion proves to be very reliable and does lead to the experimentally optimal value of c. The highly precise prediction greatly relieves the trouble of finding a suitable c value when using RBFs to solve PDEs.
For readers who are interested in the historical development and traditional methods of finding a suitable c value, we refer them to Kansa [3,4,5], Powell [6], Madych [7], Rippa [8], Fasshauer [9], and Huang et al. [10,11], although they are not needed in this paper.

2. Materials and Methods

2.1. Fundamental Theory

We revisit the c theory introduced in Luh [12]. Some basic definitions are cited as follows.
Definition 1.
For any positive number σ,
E σ : = f L 2 ( R d ) : | f ^ ( ξ ) | 2 e | ξ | 2 / σ d ξ <
where f ^ denotes the Fourier transform of f. For each f E σ , its norm is
f E σ : = | f ^ ( ξ ) | 2 e | ξ | 2 / σ d ξ 1 / 2 .
Here, f L 2 ( R d ) means that f 2 is integrable. In our theory, all functions f E σ can be interpolated by a function s ( x ) of the form
s ( x ) : = i = 1 N c i ϕ ( x x i ) + p ( x ) ,
where p ( x ) P m 1 , the space of polynomials of degree less than or equal to m 1 in R d , and X = { x 1 , , x N } is the set of centers (interpolation points). For m = 0 , P m 1 : = { 0 } . The integer m : = β / 2 . A well-established theory of RBF (radial basis functions) asserts that any set of data points ( x 1 , y 1 ) , , ( x N , y N ) , x i R d , y i R 1 for all i, can be interpolated by s ( x ) if i = 1 N c i p j ( x i ) = 0 , for j = 1 , , Q , where { p 1 , , p Q } is a basis of P m 1 . The coefficients c 1 , , c N and b 1 , , b Q in p ( x ) : = i = 1 Q b i p i ( x ) exist and are unique. For β = 1 , the polynomial p ( x ) can be dropped. For β = 1 , p ( x ) is a zero polynomial.
For any set Ω R d and any set X = { x 1 , , x N } of centers contained in Ω , the fill distance is defined by
δ ( Ω , X ) : = sup x Ω inf i = 1 , , N x x i ,
which measures the spacing of the centers in Ω . The smaller δ is, the more centers are needed. If one investigates Luh [1,2] carefully, it will be found that the parameter δ was defined there in a very different and complicated way. Here, we relax the strict requirement and let it be fill distance. It will greatly ease the pain of using the approach in those two papers, both for mathematicians and non-mathematicians. The reason we can do so is that the two kinds of δ ’s are in spirit the same thing. We also relax the requirement that the domain should be a simplex. Instead, we allow it be a cube in any dimensions. This will also greatly make things easier. In fact, the domain can be of any shape in our unorthodox approach. In this paper, Ω denotes the function domain. We require that the diameter r of Ω satisfies b 0 / 2 r b 0 for some parameter b 0 which can be arbitrarily determined by us. The number b 0 plays a crucial role in the error bound of the function interpolation.
In order to bound | f ( x ) s ( x ) | , we need the following definition.
Definition 2.
Let d and β be as in (1). The numbers ρ and Δ 0 are defined as follows.
(a)
Suppose β < d 3 . Let t = ( d β 3 ) / 2 . Then,
(i)
If β < 0 , ρ = ( 3 + t ) / 3 a n d Δ 0 = ( 2 + t ) ( 1 + t ) 3 ρ 2 ;
(ii)
If β > 0 , ρ = 1 + t 2 β / 2 + 3 a n d Δ 0 = ( 2 m + 2 + t ) ( 2 m + 1 + t ) ( 2 m + 3 ) ρ 2 m + 2 where m = β / 2 .
(b)
Suppose d 3 β < d 1 . Then, ρ = 1 and Δ 0 = 1 .
(c)
Suppose β d 1 . Let t = ( d β 3 ) / 2 . Then,
ρ = 1 a n d Δ 0 = 1 ( 2 m + 2 ) ( 2 m + 1 ) ( 2 m t + 3 ) w h e r e m = β / 2 .
The bound of | f ( x ) s ( x ) | is a very complicated expression and involves ρ and Δ 0 . We only extract its essential part and denote it by M N ( c ) , called MN function, whose graph is called MN curve. We have
| f ( x ) s ( x ) | M N ( c ) f E σ .
Here, c is just the shape parameter introduced in (1). The optimal choice of c should minimize the value of M N ( c ) .
The definition of the MN function is also complicated. There are three cases.
Case 1: β < 0 , | d + β | 1 and d + β + 1 0 . Let f E σ and ϕ ( x ) be as in (1). For any fixed fill distance δ satisfying 0 < δ b 0 / 2 , the optimal value of c in [ 24 ρ δ , ) is the number minimizing
M N ( c ) : = 8 ρ c ( β d 1 ) / 4 ( ξ * ) ( d + β + 1 ) / 2 e c ξ * ( ξ * ) 2 / σ 1 / 2 ( 2 3 ) c / ( 24 ρ δ ) if c [ 24 ρ δ , 12 b 0 ρ ) , 2 3 b 0 c ( β d + 1 ) / 4 ( ξ * ) ( d + β + 1 ) / 2 e c ξ * ( ξ * ) 2 / σ 1 / 2 ( 2 3 ) b 0 / ( 2 δ ) if c [ 12 b 0 ρ , ) ,
where
ξ * = c σ + c 2 σ 2 + 4 σ ( d + β + 1 ) 4 .
Case 2: β = 1 and d = 1 . Let f E σ and ϕ ( x ) be as in (1). For any fill distance δ satisfying 0 < δ b 0 / 2 , the optimal value of c in [ 24 ρ δ , ) is the number minimizing
M N ( c ) : = 8 ρ c ( β 1 ) / 2 1 ln 2 + 2 3 M ( c ) 1 / 2 ( 2 3 ) c / ( 24 ρ δ ) if c [ 24 ρ δ , 12 b 0 ρ ) , 2 3 b 0 c β / 2 1 ln 2 + 2 3 M ( c ) 1 / 2 ( 2 3 ) b 0 / ( 2 δ ) if c [ 12 b 0 ρ , ) ,
where
M ( c ) : = e 1 1 / ( c 2 σ ) if 0 < c 2 3 σ , g ( c σ + c 2 σ 2 + 4 σ 4 ) if 2 3 σ < c ,
g being defined by g ( ξ ) : = c ξ e c ξ ξ 2 / σ .
Case 3: β > 0 and d 1 . Let f E σ and ϕ ( x ) be as in (1). For any fixed fill distance δ satisfying 0 < δ b 0 / 2 , the optimal value of c in [ 24 ρ δ , ) is the number minimizing
M N ( c ) : = 8 ρ c ( β d 1 ) / 4 ( ξ * ) ( 1 + β + d ) / 2 e c ξ * e ( ξ * ) 2 / σ 1 / 2 ( 2 3 ) c / ( 24 ρ δ ) if c [ 24 ρ δ , 12 b 0 ρ ) , 2 3 b 0 c ( 1 + β d ) / 4 ( ξ * ) ( 1 + β + d ) / 2 e c ξ * e ( ξ * ) 2 / σ 1 / 2 ( 2 3 ) b 0 / ( 2 δ ) if c [ 12 b 0 ρ , ) ,
where
ξ * = c σ + c 2 σ 2 + 4 σ ( 1 + β + d ) 4 .
The requirement c 24 ρ δ can be ignored due to trivial reasons. Moreover, the extreme cases when c tends to 0 or infinity are not important because the optimal value of c is always a finite positive number. For detailed explanations, we refer the reader to Luh [12].
The reader should not get scared by the complicated definition of the MN function because what we need is mainly its graph which can be easily sketched by Mathematica. This function was so named by the author just to honor Professors W. R. Madych and S. A. Nelson for their outstanding contributions to the development of radial basis functions.
Seemingly, the MN-curve theory works only for function interpolation. In fact, it works for the collocation method of solving PDEs as well. The main reason is that collocation is in spirit a kind of interpolation. Combining collocation and the MN-curve theory leads to very good results, as we shall see in the next section.

2.2. Poisson Equations

The partial differential equations we deal with are of the form
u x x ( x , y ) + u y y ( x , y ) = f ( x , y ) for ( x , y ) Ω \ Ω , u ( x , y ) = g ( x , y ) for ( x , y ) Ω ,
where Ω is the domain with boundary Ω , and f , g are given functions. The reason we choose the Dirichlet condition as the boundary condition is that this setting is closer to function interpolation. For simplicity, we let Ω be a square.
The numerical solution to Equation (3) will be denoted by u ^ which is in the form of s ( x ) introduced in Equation (2). Although we cannot interpolate u with u ^ in the interior of the domain Ω , we require that u ^ satisfies the differentiation property defined by Equation (3) at some interior points, called collocation points. Then, u ^ serves as an approximation to u. As for the boundary condition, because g is explicitly given, we require that u ^ interpolates u at some boundary points. In spirit, the function u ^ is a kind of interpolator of u, like s ( x ) in Equation (2). The coefficients c i in Equation (2) can then be obtained by solving a system of linear equations. This technique of establishing a numerical solution u ^ is called the collocation method. A typical example of applying this approach to solving Poisson equations can be found in Kansa [3,4,5]. The choice of the shape parameter c can then be determined by our MN-curve theory. The technique of solving (3) by the MN-curve theory will be introduced in detail in our experiments which can be easily understood, even if the reader has not learned the collocation method.

3. Results

3.1. One-Dimensional Experiment

Although we are interested mainly in two-dimensional problems, as a prelude, a one-dimensional problem is illustrated and tested so that the reader can grasp the central idea and obtain a simple understanding of our approach.
In this experiment, the solution function is u ( x ) = e ( σ / 2.1 ) x 2 where σ = 1 . It satisfies the equations
u x x ( x ) = e ( 1 / 2.1 ) x 2 [ ( 2 / 2.1 ) 2 x 2 2 / 2.1 ] , u ( 0 ) = 1 , u ( 10 ) = e 100 / 2.1
in the domain [ 0 , 10 ] . Then, u E σ , σ = 1 , and the MN curves of Case 2 for interpolation apply if we choose β = 1 . The reason we adopt the inverse multiquadrics is that their programming is easier because the polynomial p ( x ) in (2) disappears. For β > 0 , where p ( x ) exists, the way of dealing with interpolation for E σ functions can be seen in Luh [2]. We offer six MN curves in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6, which serve as the essential error bounds for the function interpolations. The number b 0 , which greatly affects the MN curves, denotes the diameter of the interpolation domain. In these figures, it is easily seen that as the fill distance δ decreases, i.e., the number of data points increases, the optimal values of c move to 120 and are fixed there at last. The empirical results, as shown in Luh [12], show that one should choose c = 120 to make the approximation.
Now, we let
u ^ ( x ) = j = 1 201 λ j ϕ ( x x j ) ,
where ϕ ( x ) = 1 / 120 2 + | x | 2 , and require that u ^ satisfies
u ^ x x ( x j ) = e ( 1 / 2.1 ) x j 2 [ ( 2 / 2.1 ) 2 x j 2 2 / 2.1 ] for j = 2 , , 200 ,
and
u ^ ( 0 ) = 1 , u ^ ( 10 ) = e 100 / 2.1 ,
where x 1 = 0 , x 201 = 10 , and x j = 0.05 ( j 1 ) for j = 1 , , 201 .
This is a standard collocation setting. The function u ^ ( x ) is used to approximate the exact solution u ( x ) of the differential equation. After solving the linear equations for λ j s , we test | u ( x ) u ^ ( x ) | at 400 test points z 1 , , z 400 evenly spaced in [ 0 , 10 ] and find its root-mean-square error (RMS)
R M S = 1 400 i = 1 400 | u ( z i ) u ^ ( z i ) | 2 1 / 2 = 1.25 × 10 83 .
The condition number of the linear system is 4.4 × 10 643 . With the arbitrarily precise computer software Mathematica, we kept 800 effective digits to the right of the decimal point for each step of the calculation, successfully overcoming the problem of ill-conditioning. The computer time for solving the linear system is less than one second. We did not test smaller fill distances δ s and different c s because the RMS is already satisfactory.
Note that, here, u ^ ( x ) is in the same form as the interpolating function defined in (2) where m = 0 because β = 1 . Although the collocation points are not exactly the interpolation points, they are in spirit like the interpolation points. It is not surprising that | u ( x ) u ^ ( x ) | is very small.

3.2. Two-Dimensional Experiment

Here, the solution function is u ( x , y ) = e ( σ / 2.1 ) ( x 2 + y 2 ) where σ = 10 36 . The domain is a large square with vertices ( 0 , 0 ) , ( 10 16 , 0 ) , ( 10 16 , 10 16 ) , and ( 0 , 10 16 ) . The function u ( x , y ) satisfies
u x x ( x , y ) + u y y ( x , y ) = ( 2 σ / 2.1 ) e ( σ / 2.1 ) ( x 2 + y 2 ) [ 2 ( 2 σ / 2.1 ) ( x 2 + y 2 ) ]
for ( x , y ) in the interior of the domain Ω = { ( x , y ) | 0 x 10 16 , 0 y 10 16 } and
u ( x , y ) = e ( σ / 2.1 ) ( x 2 + y 2 )
for ( x , y ) on the boundary Ω .
It can be easily checked that u E σ where σ = 10 36 and Case 1 of the MN curves applies if we choose β = 1 . Five MN curves are shown in Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11.
All these curves show that one should choose c = 7000 × 2 × 10 14 0.99 × 10 18 as the shape parameter in ϕ ( x , y ) = 1 / c 2 + x 2 + y 2 . We let u ^ ( x , y ) : = j = 1 N d λ j ϕ ( x x j , y y j ) be the approximate solution of the Poisson equation satisfying (5) and (6), and require that it satisfy
u ^ x x ( x , y ) + u ^ y y ( x , y ) = ( 2 σ / 2.1 ) e ( σ / 2.1 ) ( x 2 + y 2 ) [ 2 ( 2 σ / 2.1 ) ( x 2 + y 2 ) ]
for ( x , y ) = ( x j , y j ) where ( x j , y j ) belong to the interior of Ω , i.e., Ω \ Ω , and N d denotes the number of data points used. Moreover,
u ^ ( x , y ) = e ( σ / 2.1 ) ( x 2 + y 2 )
for ( x , y ) = ( x j , y j ) where ( x j , y j ) belong to the boundary Ω .
A grid of 41 × 41 is adopted. Hence, there are 1681 data points ( x j , y j ) altogether. Among them, 160 are boundary points where the Dirichlet condition occurs. Thus, the fill distance is δ = 1.25 2 E14. The reason we choose such a delta is that it is close to the one in Figure 7. When applying the MN curves, we consider all the 1681 data points to be the interpolation points, even though it is not theoretically rigorous. As explained in Luh [12], it is supposed to work well. However, something important must be pointed out. Although MN curves can be used to predict almost exactly the optimal value of c for function interpolation, a moderate search may be needed if this approach is used in a non-rigorous way. We begin with the theoretically predicted optimal value c = 7000 2 × 10 14 and test two nearby values, one larger and the other smaller. Then, we check the RMS on the boundary for each c and choose the direction which makes the RMS smaller. Continuing to choose c in this direction, we stop when the RMS values begin to grow. Our experiment shows that not many steps are needed, and the finally obtained c does produce the best result.
The experimental results are presented in Table 1. Here, R M S , N d , N t , and C O N D denote the root-mean-square error, number of data points, number of test points, and condition number of the linear system, respectively. The optimal value of c is marked with the symbol *. We use R M S b d y to denote the root-mean-square error of the approximation on the boundary, generated by 800 test points located on the boundary. In the entire domain Ω , 6400 test points were used to generate the R M S s. The most time-consuming command of solving the system of linear equations took about 30 min for each c. Although we adopted 1200 effective digits for each step of the calculation, it still worked with an acceptable time efficiency.
Note that the optimal value of c is 800 2 × 10 14 which coincides with the value chosen by our stopping criterion based on R M S b d y . This optimal value was obtained without much effort. The most time-consuming step is solving a linear system for each trial of the c value whenever the condition number is very large. Fortunately, it took only 30 min in our experiment. Obviously, we could have obtained better RMS values by increasing the number of data points, whereas we did not do so because the approximation was already quite good. Moreover, it can be seen that the theoretical value 7000 2 × 10 14 is not the same as the experimental value. The main reason is that we relaxed two strict requirements in Luh [1,2], as explained in the paragraph preceding Definition 2.2. Another important reason is that, according to the strict theory in [1,2], the MN functions are defined only for function interpolation, not for the collocation method of solving PDEs.
In this paper, all the calculations were not performed by double precision. Instead, they were made by the arbitrarily precise computer software Mathematica. In the MN-curve approach, the condition number of the linear system may be very large. If double precision is adopted, the final results may be meaningless. For example, in Table 1, when c = 8000 2 × 10 14 , the condition number is 2.2 × 10 730 . One has to adopt at least 750 effective digits to the right of the decimal point for each step of the calculation. In order to increase the reader’s confidence, we adopted 1200 digits at the cost of spending more time. Hence, our calculations should be reliable. This is about the stability. As for the convergence, in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11, it can be clearly seen that as the fill distance δ decreases, the value of MN(c), which denotes the essential error bound of the approximation, decreases rapidly. It fully reflects a salient characteristic of the RBF approach. In our experiments, we did not test different fill distances because our primary concern is finding a good c value, not the convergence rate.
Seemingly, it is a limitation of our approach that one has to adopt Mathematica rather than the widely used Matlab. In fact, it probably can be considered to be a limitation of the entire RBF approach. In Madych [7], an incomplete experiment was presented. Meaningless results appear in that experiment whenever the condition number is greater than 10 16 . What Madych used was Matlab and a double-precision scheme. In the paper, Madych said that he did not know how to overcome this trouble. By virtue of Mathematica, we can now handle it successfully. The severe ill-conditioning of the RBF approach seems to be an inherent problem which we have to face.

4. Discussion

The author believes that only a very small part of the PDE problems has been dealt with in this paper. For future research, we should try to explore different boundary conditions, different domains, and various PDEs. These will be an endless work.
Another paper by Oaxaca-Adams et al. [13] may be worth mentioning. In [13], a curve is also used to obtain the optimal point of interest, and in a similar way, a numerical approach is used that allows numerically solving a type of systems of interest in applied sciences.
In physics, as said by Kansa, many numerical solutions to PDEs are not bad, but truly good solutions are rarely seen. E. Kansa invented the collocation method and opened a new route to solving them. The combination of the c theory and collocation does produce very good results, as shown in our experiments. Maybe this is just a starting point. We are still facing a huge challenge and have a lot of work to conduct in the future.

Funding

This research was funded by National Science and Technology Council through the MOST project 107-2115-M-126-005. The APC was included.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Luh, L.-T. The mystery of the shape parameter IV. Eng. Anal. Bound. Elem. 2014, 48, 24–31. [Google Scholar] [CrossRef] [Green Version]
  2. Luh, L.-T. The mystery of the shape parameter III. Appl. Comput. Harmon. Anal. 2016, 40, 186–199. [Google Scholar] [CrossRef]
  3. Kansa, E.J. Multiquadrics—A scattered data approximation scheme with applications to computational fluid dynamics I: Surface approximations and partial derivative estimates. Comput. Math. Appl. 1990, 19, 127–145. [Google Scholar] [CrossRef] [Green Version]
  4. Kansa, E.J. Multiquadrics-a scattered data approximation scheme with applications to computational fluid dynamics II: Solutions to parabolic, hyperbolic, and elliptic partial differential equations. Comput. Math. Appl. 1990, 19, 147–161. [Google Scholar] [CrossRef] [Green Version]
  5. Kansa, E.J.; Holoborodko, P. On the ill-conditioned nature of C RBF strong collocation. Eng. Anal. Bound. Elem. 2017, 78, 26–30. [Google Scholar] [CrossRef]
  6. Powell, M.J.D. Univariate Multiquadric Interpolation: Some Recent Results. In Curves and Surfaces; Academic Press: San Diego, CA, USA, 1991; pp. 371–382. [Google Scholar]
  7. Madych, W.R. Miscellaneous error bounds for multiquadric and related interpolators. Comput. Math. Appl. 1992, 24, 121–138. [Google Scholar] [CrossRef] [Green Version]
  8. Rippa, S. An algorithm for selecting a good value for the parameter c in a radial basis function interpolation. Adv. Comput. Math. 1999, 11, 193–210. [Google Scholar] [CrossRef]
  9. Fasshauer, G. Meshfree Approximation Methods with MATLAB; World Scientific Publishers: Singapore, 2007. [Google Scholar]
  10. Huang, C.-S.; Lee, C.-F.; Cheng, A.H.D. Error estimate, optimal shape parameter, and high precision computation of multiquadric collocation method. Eng. Anal. Bound Elem. 2007, 31, 615–623. [Google Scholar] [CrossRef]
  11. Huang, C.-S.; Yen, H.-D.; Cheng, A.H.D. On the increasingly flat radial basis function and optimal shape parameter for the solution of elliptic PDEs. Eng. Anal. Bound Elem. 2010, 34, 802–809. [Google Scholar] [CrossRef]
  12. Luh, L.-T. The choice of the shape parameter-a friendly approach. Eng. Anal. Bound. Elem. 2019, 98, 103–109. [Google Scholar] [CrossRef]
  13. Oaxaca-Adams, G.; Villafuerte-Segura, R.; Aguirre-Hernández, B. On nonfragility of controllers for time delay systems: A numerical approach. J. Frankl. Inst. 2021, 358, 4671–4686. [Google Scholar] [CrossRef]
Figure 1. Here, d = 1 , β = 1 , b 0 = 10 , and σ = 1 .
Figure 1. Here, d = 1 , β = 1 , b 0 = 10 , and σ = 1 .
Mathematics 10 04582 g001
Figure 2. Here, d = 1 , β = 1 , b 0 = 10 , and σ = 1 .
Figure 2. Here, d = 1 , β = 1 , b 0 = 10 , and σ = 1 .
Mathematics 10 04582 g002
Figure 3. Here, d = 1 , β = 1 , b 0 = 10 , and σ = 1 .
Figure 3. Here, d = 1 , β = 1 , b 0 = 10 , and σ = 1 .
Mathematics 10 04582 g003
Figure 4. Here, d = 1 , β = 1 , b 0 = 10 , and σ = 1 .
Figure 4. Here, d = 1 , β = 1 , b 0 = 10 , and σ = 1 .
Mathematics 10 04582 g004
Figure 5. Here, d = 1 , β = 1 , b 0 = 10 , and σ = 1 .
Figure 5. Here, d = 1 , β = 1 , b 0 = 10 , and σ = 1 .
Mathematics 10 04582 g005
Figure 6. Here, d = 1 , β = 1 , b 0 = 10 , and σ = 1 .
Figure 6. Here, d = 1 , β = 1 , b 0 = 10 , and σ = 1 .
Mathematics 10 04582 g006
Figure 7. Here, d = 2 , β = 1 , b 0 = 2 E16, and σ = 1 E 36 .
Figure 7. Here, d = 2 , β = 1 , b 0 = 2 E16, and σ = 1 E 36 .
Mathematics 10 04582 g007
Figure 8. Here, d = 2 , β = 1 , b 0 = 2 E16, and σ = 1 E 36 .
Figure 8. Here, d = 2 , β = 1 , b 0 = 2 E16, and σ = 1 E 36 .
Mathematics 10 04582 g008
Figure 9. Here, d = 2 , β = 1 , b 0 = 2 E16, and σ = 1 E 36 .
Figure 9. Here, d = 2 , β = 1 , b 0 = 2 E16, and σ = 1 E 36 .
Mathematics 10 04582 g009
Figure 10. Here, d = 2 , β = 1 , b 0 = 2 E16, and σ = 1 E 36 .
Figure 10. Here, d = 2 , β = 1 , b 0 = 2 E16, and σ = 1 E 36 .
Mathematics 10 04582 g010
Figure 11. Here, d = 2 , β = 1 , b 0 = 2 E16, and σ = 1 E 36 .
Figure 11. Here, d = 2 , β = 1 , b 0 = 2 E16, and σ = 1 E 36 .
Mathematics 10 04582 g011
Table 1. δ = 2 × 1.25 × 10 14 , b 0 = 2 × 10 16 , N d = 1681 , N t = 6400 .
Table 1. δ = 2 × 1.25 × 10 14 , b 0 = 2 × 10 16 , N d = 1681 , N t = 6400 .
c 300 2 × 10 14 400 2 × 10 14 500 2 × 10 14 600 2 × 10 14 700 2 × 10 14
R M S 5.2 × 10 134 3.7 × 10 139 1.5 × 10 143 2.1 × 10 145 5.8 × 10 147
C O N D 1.1 × 10 505 5.8 × 10 524 1.2 × 10 540 3.9 × 10 552 1.5 × 10 563
R M S b d y 1.07 × 10 143 7.2 × 10 149 3.2 × 10 153 4.1 × 10 155 1.1 × 10 156
c 800 2 × 10 14 * 900 2 × 10 14 1000 2 × 10 14 1100 2 × 10 14 1200 2 × 10 14
R M S 3.6 × 10 148 7.9 × 10 148 8.2 × 10 146 5.2 × 10 144 6.0 × 10 143
C O N D 2.1 × 10 572 2.6 × 10 580 4.4 × 10 587 1.5 × 10 594 1.4 × 10 600
R M S b d y 7.2 × 10 158 1.7 × 10 157 1.6 × 10 155 1.0 × 10 153 1.2 × 10 152
c 3000 2 × 10 14 5000 2 × 10 14 7000 2 × 10 14 8000 2 × 10 14
R M S 2.3 × 10 133 7.1 × 10 127 1.1 × 10 122 5.4 × 10 121
C O N D 1.1 × 10 663 1.2 × 10 698 1.5 × 10 721 2.2 × 10 730
R M S b d y 4.4 × 10 143 1.8 × 10 136 5.1 × 10 132 3.2 × 10 130
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Luh, L.-T. Solving Poisson Equations by the MN-Curve Approach. Mathematics 2022, 10, 4582. https://doi.org/10.3390/math10234582

AMA Style

Luh L-T. Solving Poisson Equations by the MN-Curve Approach. Mathematics. 2022; 10(23):4582. https://doi.org/10.3390/math10234582

Chicago/Turabian Style

Luh, Lin-Tian. 2022. "Solving Poisson Equations by the MN-Curve Approach" Mathematics 10, no. 23: 4582. https://doi.org/10.3390/math10234582

APA Style

Luh, L. -T. (2022). Solving Poisson Equations by the MN-Curve Approach. Mathematics, 10(23), 4582. https://doi.org/10.3390/math10234582

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop