Next Article in Journal
Linear Optimization of Polynomial Rational Functions: Applications for Positivity Analysis
Next Article in Special Issue
Oscillation Criteria of Higher-order Neutral Differential Equations with Several Deviating Arguments
Previous Article in Journal
An Extended EDAS Method with a Single-Valued Complex Neutrosophic Set and Its Application in Green Supplier Selection
Previous Article in Special Issue
Weighted Fractional Iyengar Type Inequalities in the Caputo Direction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Alternating Asymmetric Iterative Algorithm Based on Domain Decomposition for 3D Poisson Problem

School of Mathematics and Statistics, Ningxia University, Yinchuan 750021, Ningxia, China
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(2), 281; https://doi.org/10.3390/math8020281
Submission received: 2 January 2020 / Revised: 31 January 2020 / Accepted: 15 February 2020 / Published: 19 February 2020
(This article belongs to the Special Issue Multivariate Approximation for solving ODE and PDE)

Abstract

:
Poisson equation is a widely used partial differential equation. It is very important to study its numerical solution. Based on the strategy of domain decomposition, the alternating asymmetric iterative algorithm for 3D Poisson equation is provided. The solution domain is divided into several sub-domains, and eight asymmetric iterative schemes with the relaxation factor for 3D Poisson equation are constructed. When the numbers of iteration are odd or even, the computational process of the presented iterative algorithm are proposed respectively. In the calculation of the inner interfaces, the group explicit method is used, which makes the algorithm to be performed fast and in parallel, and avoids the difficulty of solving large-scale linear equations. Furthermore, the convergence of the algorithm is analyzed theoretically. Finally, by comparing with the numerical experimental results of Jacobi and Gauss Seidel iterative algorithms, it is shown that the alternating asymmetric iterative algorithm based on domain decomposition has shorter computation time, fewer iteration numbers and good parallelism.

1. Introduction

Poisson equation is an elliptic partial differential equation, which frequently appears in many fields such as fluid dynamics, heat transfer, electromagnetics, acoustics, electrostatics mechanical engineering and so on. Many researches on studding the numerical techniques to approximate the solution of Poisson equation have been made in the past few decades. The application of finite difference methods for solving Poisson equation will normally lead to a large, block, and sparse system of equations. Direct methods and iterative methods [1] are normally considered as common approaches for solving such system of equations. Several high precision multigrid and compact difference methods are given in [2,3,4,5,6]. Romao et al. [7,8] provides the Galerkin and least-squares finite element methods in the solution of 3D Poisson equation. In [9,10], the Haar wavelet methods are given. Speyer et al. [11] provide a preconditioned bi-conjugate gradient stabilized method which is efficient, albeit nonmonotonic and convergent.
With the continuous improvement of computer hardware, people are more and more focused on solving large-scale scientific and engineering problems quickly and efficiently on parallel computers. Therefore, people wish to find some direct methods and iterative methods, which have the characteristic of much better solving elliptic equations and easier parallel implementation. In recent years, parallel algorithms are also constantly emerging. Several new parallel methods of direct solution are proposed. P. Valero-Lara and A. Pinelli et al. [12] provide the implementation of a fast solver based on a block cyclic reduction algorithm for the linear systems of a three dimensional separable elliptic problem. And they also study on the parallel characteristics of an algorithm for the direct solution of linear systems with a block-tridiagonal coefficient matrix (BLKTRI problem) [13]. C. P. Stone et al. [14] analyze the performance of a block tridiagonal benchmark. Many authors have given the implementation of scalar tridiagonal solver on GPUs [15,16,17]. Y. Zhang [16] also illustrates several methods to solve tridiagonal systems on GPUs. Because of the direct method to solve large-scale sparse and block diagonal equation systems, when the coefficient matrix is close to singularity, the calculation will often stop or make mistakes. So people also find iterative methods which can be solved by constructing some efficient iterative schemes to approximate the problem itself, so that the iteration can reach a certain accuracy. In [18,19,20,21], a class of efficient parallel finite difference iterative algorithms for Poisson equation were also proposed.
In addition, the domain decomposition method [22] is also a powerful tool for parallel implementation, which studies parallelization from the model level of physical problems. This kind of method can decompose scale problem into small-scale problem and solve serial problem into parallel problem. The explicit-implicit domain decomposition method is proposed by Kuznetsov [23]. Because the numerical boundary conditions on the internal boundary are often not the same as those of the original mathematical model or the corresponding physical problems, different methods to obtain the internal boundary information form different explicit-implicit domain decomposition (EIDD) methods. This leads to the idea of parallel implementation for iterative method based on domain decomposition. In [24], the authors have proposed a kind of finite difference parallel iterative algorithm for two-dimensional Poisson problem, and verified its efficiency and accuracy.
This paper extends to the study of the domain decomposition method for three-dimensional Poisson problem. Several finite difference asymmetric iterative schemes are constructed, and each asymmetric iterative schemes are used to solve the sub-domains alternatively and in parallel; in the processing of inner interfaces, group explicit (GE) method [25,26] is used. The calculation on the whole solution domain is explicit but using the implicit iterative schemes, which greatly avoids the difficulty of solving linear equations and improves the calculation speed and accuracy. When the number of iteration is odd or even, the iterative process of the presented algorithm is given respectively, and a kind of efficient iterative algorithm is established based on domain decomposition for solving three-dimensional Poisson equation.
This paper is outlined as follows. In Section 2, we present several asymmetric iterative schemes. Section 3 gives the alternating asymmetric iterative algorithm. And the convergence and the optimal relaxation factor are obtained in Section 4. In Section 5, we perform the numerical experiments to examine the presented algorithm. Finally we give the conclusion of this paper in Section 6.

2. Asymmetric Iterative Schemes

Consider the three-dimensional Poisson problem,
2 u x 2 + 2 u y 2 + 2 u z 2 = f ( x , y , z ) , ( x , y , z ) Ω ,
with the boundary condition,
u ( x , y , z ) = g ( x , y , z ) , ( x , y , z ) Ω .
where Ω = [ 0 , L ] × [ 0 , M ] × [ 0 , K ] , and Ω is the boundary of the domain Ω . We divide the solution domain Ω into uniform grid, the space step h x = L / l in x direction, h y = M / m in y direction and h z = K / s in z direction. For implicity, the space steps are assumed equal that h x = h y = h z = h . Denote x i = i h , i = 0 , 1 , , l ; y j = j h , j = 0 , 1 , , m ; z k = k h , k = 0 , 1 , , s ; u i , j , k ( n ) as numerical solution on the nth iteration level at the grid node ( x i , y j , z k ) . We can give the classical difference discretization in Equation (3) for the 3D Poisson Equation (1),
u i + 1 , j , k 2 u i , j , k + u i 1 , j , k h 2 + u i , j + 1 , k 2 u i , j , k + u i , j 1 , k h 2 + u i , j , k + 1 2 u i , j , k + u i , j , k 1 h 2 = f i , j , k ,
namely,
u i , j , k 1 6 ( u i + 1 , j , k + u i 1 , j , k + u i , j + 1 , k + u i , j 1 , k + u i , j , k + 1 + u i , j , k 1 h 2 f i , j , k ) = 0 .
Then we construct eight asymmetric iterative schemes by the difference operator L with the relaxation factor ω as follows,
L 1 u i , j , k ( n + 1 ) = u i , j , k ( n + 1 ) 1 6 ω ( u i + 1 , j , k ( n + 1 ) + u i , j + 1 , k ( n + 1 ) + u i , j , k + 1 ( n + 1 ) ) + u i 1 , j , k ( n ) + u i , j 1 , k ( n ) + u i , j , k 1 ( n ) + ( 1 ω ) ( u i + 1 , j , k ( n ) + u i , j + 1 , k ( n ) + u i , j , k + 1 ( n ) ) h 2 f i , j , k ,
L 2 u i , j , k ( n + 1 ) = u i , j , k ( n + 1 ) 1 6 ω ( u i 1 , j , k ( n + 1 ) + u i , j + 1 , k ( n + 1 ) + u i , j , k + 1 ( n + 1 ) ) + u i + 1 , j , k ( n ) + u i , j 1 , k ( n ) + u i , j , k 1 ( n ) + ( 1 ω ) ( u i 1 , j , k ( n ) + u i , j + 1 , k ( n ) + u i , j , k + 1 ( n ) ) h 2 f i , j , k ,
L 3 u i , j , k ( n + 1 ) = u i , j , k ( n + 1 ) 1 6 ω ( u i 1 , j , k ( n + 1 ) + u i , j 1 , k ( n + 1 ) + u i , j , k + 1 ( n + 1 ) ) + u i + 1 , j , k ( n ) + u i , j + 1 , k ( n ) + u i , j , k 1 ( n ) + ( 1 ω ) ( u i 1 , j , k ( n ) + u i , j 1 , k ( n ) + u i , j , k + 1 ( n ) ) h 2 f i , j , k ,
L 4 u i , j , k ( n + 1 ) = u i , j , k ( n + 1 ) 1 6 ω ( u i + 1 , j , k ( n + 1 ) + u i , j 1 , k ( n + 1 ) + u i , j , k + 1 ( n + 1 ) ) + u i 1 , j , k ( n ) + u i , j + 1 , k ( n ) + u i , j , k 1 ( n ) + ( 1 ω ) ( u i + 1 , j , k ( n ) + u i , j 1 , k ( n ) + u i , j , k + 1 ( n ) ) h 2 f i , j , k ,
L 5 u i , j , k ( n + 1 ) = u i , j , k ( n + 1 ) 1 6 ω ( u i + 1 , j , k ( n + 1 ) + u i , j + 1 , k ( n + 1 ) + u i , j , k 1 ( n + 1 ) ) + u i 1 , j , k ( n ) + u i , j 1 , k ( n ) + u i , j , k + 1 ( n ) + ( 1 ω ) ( u i + 1 , j , k ( n ) + u i , j + 1 , k ( n ) + u i , j , k 1 ( n ) ) h 2 f i , j , k ,
L 6 u i , j , k ( n + 1 ) = u i , j , k ( n + 1 ) 1 6 ω ( u i 1 , j , k ( n + 1 ) + u i , j + 1 , k ( n + 1 ) + u i , j , k 1 ( n + 1 ) ) + u i + 1 , j , k ( n ) + u i , j 1 , k ( n ) + u i , j , k + 1 ( n ) + ( 1 ω ) ( u i 1 , j , k ( n ) + u i , j + 1 , k ( n ) + u i , j , k 1 ( n ) ) h 2 f i , j , k ,
L 7 u i , j , k ( n + 1 ) = u i , j , k ( n + 1 ) 1 6 ω ( u i 1 , j , k ( n + 1 ) + u i , j 1 , k ( n + 1 ) + u i , j , k 1 ( n + 1 ) ) + u i + 1 , j , k ( n ) + u i , j + 1 , k ( n ) + u i , j , k + 1 ( n ) + ( 1 ω ) ( u i 1 , j , k ( n ) + u i , j 1 , k ( n ) + u i , j , k 1 ( n ) ) h 2 f i , j , k ,
L 8 u i , j , k ( n + 1 ) = u i , j , k ( n + 1 ) 1 6 ω ( u i + 1 , j , k ( n + 1 ) + u i , j 1 , k ( n + 1 ) + u i , j , k 1 ( n + 1 ) ) + u i 1 , j , k ( n ) + u i , j + 1 , k ( n ) + u i , j , k + 1 ( n ) + ( 1 ω ) ( u i + 1 , j , k ( n ) + u i , j 1 , k ( n ) + u i , j , k 1 ( n ) ) h 2 f i , j , k ,
Figure 1 represents the distribution of unknown solution at the ( n + 1 ) th iteration level for the eight asymmetric iterative schemes (5)–(12).

3. Alternating Asymmetric Iterative Algorithm Based on Domain Decomposition

3.1. The Domain Decomposition

We can divide the 3D solution domain Ω into multi-subdomains. For simplicity, we use six grid planes x = p , x = p + 1 , y = q , y = q + 1 , z = r , z = r + 1 to discrete the solution domain Ω into eight sub-domains, and note Ω i , i = 1 , 2 , , 8 as subsets of grid points, while p , q , r are positive integers with p [ 1 , l ] , q [ 1 , m ] , r [ 1 , s ] . Denote π i is the interfaces of the sub-domain Ω i . The sorting order of sub-domains is as follows: the subspaces above z = r + 1 are sorted anticlockwise starting from the upper right sub-domain of the inner layer, and the sub-domains under z = r are sorted anticlockwise starting from the lower right subspace (as shown in Figure 2). The specific description is as follows:
Ω 1 : ( x , y , z ) | x = l 1 , l 2 , , p + 1 ; y = m 1 , m 2 , , q + 1 ; z = s 1 , s 2 , , r + 1 , Ω 2 : ( x , y , z ) | x = 1 , 2 , , p ; y = m 1 , m 2 , , q + 1 ; z = s 1 , s 2 , , r + 1 , Ω 3 : ( x , y , z ) | x = 1 , 2 , , p ; y = 1 , 2 , , q ; z = s 1 , s 2 , , r + 1 , Ω 4 : ( x , y , z ) | x = l 1 , l 2 , , p + 1 ; y = 1 , 2 , , q ; z = s 1 , s 2 , , r + 1 , Ω 5 : ( x , y , z ) | x = l 1 , l 2 , , p + 1 ; y = m 1 , m 2 , , q + 1 ; z = 1 , 2 , , r , Ω 6 : ( x , y , z ) | x = 1 , 2 , , p ; y = m 1 , m 2 , , q + 1 ; z = 1 , 2 , , r , Ω 7 : ( x , y , z ) | x = 1 , 2 , , p ; y = 1 , 2 , , q ; z = 1 , 2 , , r , Ω 8 : ( x , y , z ) | x = l 1 , l 2 , , p + 1 ; y = 1 , 2 , , q ; z = 1 , 2 , , r .

3.2. Algorithm Implementation

In this subsection, we provide a new alternating asymmetric iterative (AAI) algorithm based on domain decomposition for 3D Poisson problem (1) and (2). We give different computational processes in each sub-domains at the odd iteration layers and even iteration layers respectively, and use the asymmetric iterative schemes alternatively. The Group Explicit (GE) method is used to solve the inner interfaces, which makes the algorithm to be computed fast and in parallel, and avoids the difficulty of solving large-scale linear equations.

3.2.1. Implementation of Odd Level Iteration

When the iteration number are odd, namely, n = 2 a + 1 , a = 0 , 1 , , we solve the grid nodes from the boundaries to the inner interfaces step by step, that is, using the asymmetric iterative schemes (5)–(12) to solve the grid points in Ω i , i = 1 , 2 , , 8 respectively:
L 1 u i , j , k ( 2 a + 1 ) = 0 , ( i , j , k ) Ω 1 , L 2 u i , j , k ( 2 a + 1 ) = 0 , ( i , j , k ) Ω 2 , L 3 u i , j , k ( 2 a + 1 ) = 0 , ( i , j , k ) Ω 3 , L 4 u i , j , k ( 2 a + 1 ) = 0 , ( i , j , k ) Ω 4 , L 5 u i , j , k ( 2 a + 1 ) = 0 , ( i , j , k ) Ω 5 , L 6 u i , j , k ( 2 a + 1 ) = 0 , ( i , j , k ) Ω 6 , L 7 u i , j , k ( 2 a + 1 ) = 0 , ( i , j , k ) Ω 7 , L 8 u i , j , k ( 2 a + 1 ) = 0 , ( i , j , k ) Ω 8 .
Obviously, the numerical solution can be obtained independently in parallel when the iteration numbers are odd, which saves a lot of computational time compared with the full-implicit iteration case. In addition, although the asymmetric iterative schemes are implicit, the computational process can be transformed into explicit, which can obviously improve the calculation speed and avoid solving large and complex linear equations.

3.2.2. Implementation of Even Level Iteration

When the iteration number is even, namely, n = 2 a + 2 , a = 0 , 1 , , we calculate the numerical solution from the inner interfaces to the boundaries step by step. Where the computational process of the interfaces π = π i , i = 1 , 2 , , 8 includes three parts:
  • Interfaces I (namely, the grid nodes at the center of the domain Ω ) (shown in Figure 3): ( p , q , r + 1 ) , ( p + 1 , q , r + 1 ) , ( p + 1 , q + 1 , r + 1 ) , ( p , q + 1 , r + 1 ) , ( p , q , r ) , ( p + 1 , q , r ) , ( p + 1 , q + 1 , r ) , ( p , q + 1 , r ) .
  • Interfaces II: the interface lines except Interfaces I (shown in Figure 3).
  • Interfaces III: the interfaces except Interfaces I and II, namely,
    π \ I n t e r f a c e s I I n t e r f a c e s I I .
Therefore, it can be seen that the interfaces π = Interfaces I ∪ Interfaces II ∪ Interfaces III. When the interfaces π are solved, the inner grid nodes in the domain Ω can be solved in order like the odd case of iteration numbers in Section 3.2.1. We give the computational procedures in detail as follows.
(1) The solution to Interfaces I
We use the asymmetric iterative schemes (5)–(12) to solve the Interfaces I, then the following linear equations can be obtained:
M 1 u p + 1 , q + 1 , r + 1 ( 2 a + 2 ) u p , q + 1 , r + 1 ( 2 a + 2 ) u p , q , r + 1 ( 2 a + 2 ) u p + 1 , q , r + 1 ( 2 a + 2 ) u p + 1 , q + 1 , r ( 2 a + 2 ) u p , q + 1 , r ( 2 a + 2 ) u p , q , r ( 2 a + 2 ) u p + 1 , q , r ( 2 a + 2 ) = N 1 ,
where matrices M 1 and N 1 are represented as bellow,
M 1 = 6 ω ω ω ω 6 ω ω ω 6 ω ω ω ω 6 ω ω 6 ω ω ω ω 6 ω ω ω 6 ω ω ω ω 6 , N 1 = e 1 e 2 e 3 e 4 e 5 e 6 e 7 e 8 ,
e 1 = u p + 2 , q + 1 , r + 1 ( 2 a + 1 ) + u p + 1 , q + 2 , r + 1 ( 2 a + 1 ) + u p + 1 , q + 1 , r + 2 ( 2 a + 1 ) + ( 1 ω ) ( u p , q + 1 , r + 1 ( 2 a + 1 ) + u p + 1 , q , r + 1 ( 2 a + 1 ) + u p + 1 , q + 1 , r ( 2 a + 1 ) ) h 2 f p + 1 , q + 1 , r + 1 , e 2 = u p 1 , q + 1 , r + 1 ( 2 a + 1 ) + u p , q + 2 , r + 1 ( 2 a + 1 ) + u p , q + 1 , r + 2 ( 2 a + 1 ) + ( 1 ω ) ( u p + 1 , q + 1 , r + 1 ( 2 a + 1 ) + u p , q , r + 1 ( 2 a + 1 ) + u p , q + 1 , r ( 2 a + 1 ) ) h 2 f p , q + 1 , r + 1 , e 3 = u p 1 , q , r + 1 ( 2 a + 1 ) + u p , q 1 , r + 1 ( 2 a + 1 ) + u p , q , r + 2 ( 2 a + 1 ) + ( 1 ω ) ( u p + 1 , q , r + 1 ( 2 a + 1 ) + u p , q + 1 , r + 1 ( 2 a + 1 ) + u p , q , r ( 2 a + 1 ) ) h 2 f p , q , r + 1 , e 4 = u p + 2 , q , r + 1 ( 2 a + 1 ) + u p + 1 , q 1 , r + 1 ( 2 a + 1 ) + u p + 1 , q , r + 2 ( 2 a + 1 ) + ( 1 ω ) ( u p , q , r + 1 ( 2 a + 1 ) + u p + 1 , q + 1 , r + 1 ( 2 a + 1 ) + u p + 1 , q , r ( 2 a + 1 ) ) h 2 f p + 1 , q , r + 1 , e 5 = u p + 2 , q + 1 , r ( 2 a + 1 ) + u p + 1 , q + 2 , r ( 2 a + 1 ) + u p + 1 , q + 1 , r 1 ( 2 a + 1 ) + ( 1 ω ) ( u p , q + 1 , r ( 2 a + 1 ) + u p + 1 , q , r ( 2 a + 1 ) + u p + 1 , q + 1 , r + 1 ( 2 a + 1 ) ) h 2 f p + 1 , q + 1 , r , e 6 = u p 1 , q + 1 , r ( 2 a + 1 ) + u p , q + 2 , r ( 2 a + 1 ) + u p , q + 1 , r 1 ( 2 a + 1 ) + ( 1 ω ) ( u p + 1 , q + 1 , r ( 2 a + 1 ) + u p , q , r ( 2 a + 1 ) + u p , q + 1 , r + 1 ( 2 a + 1 ) ) h 2 f p , q + 1 , r , e 7 = u p 1 , q , r ( 2 a + 1 ) + u p , q 1 , r ( 2 a + 1 ) + u p , q , r 1 ( 2 a + 1 ) + ( 1 ω ) ( u p + 1 , q , r ( 2 a + 1 ) + u p , q + 1 , r ( 2 a + 1 ) + u p , q , r + 1 ( 2 a + 1 ) ) h 2 f p , q , r , e 8 = u p + 2 , q , r ( 2 a + 1 ) + u p + 1 , q 1 , r ( 2 a + 1 ) + u p + 1 , q , r 1 ( 2 a + 1 ) + ( 1 ω ) ( u p , q , r ( 2 a + 1 ) + u p + 1 , q + 1 , r ( 2 a + 1 ) + u p + 1 , q , r + 1 ( 2 a + 1 ) ) h 2 f p + 1 , q , r .
Then we just solve the above eight-order sparse linear Equation (14) to obtain the numerical solution of the i n t e r f a c e I .
(2) The solution to Interfaces II
The computational procedure of the Interfaces I is depending on the use of GE method based on eight points per group. Similarly, we use the GE method based on four points per group to solve the Interfaces II between the inner boundaries of eight subspaces Ω i , i = 1 , 2 , , 8 . Take one group of the Interfaces II for example to illustrate the order of the solution process. Figure 3 gives the direction of the iteration computation.
Using the asymmetric iterative schemes (6), (7), (10), (11) to solve the grid nodes ( i , q , r ) , ( i , q , r + 1 ) , ( i , q + 1 , r + 1 ) , ( i , q + 1 , r ) i = p + 2 , p + 3 , , l 1 (shown in Figure 3), then we can provide the following fourth-order linear equations:
6 ω ω ω 6 ω ω 6 ω ω ω 6 u i , q , r ( 2 a + 2 ) u i , q , r + 1 ( 2 a + 2 ) u i , q + 1 , r + 1 ( 2 a + 2 ) u i , q + 1 , r ( 2 a + 2 ) = d 1 d 2 d 3 d 4 ,
where
d 1 = ω u i 1 , q , r ( 2 a + 2 ) + u i + 1 , q , r ( 2 a + 1 ) + u i , q 1 , r ( 2 a + 1 ) + u i , q , r 1 ( 2 a + 1 ) + ( 1 ω ) ( u i 1 , q , r ( 2 a + 1 ) + u i , q + 1 , r ( 2 a + 1 ) + u i , q , r + 1 ( 2 a + 1 ) ) , d 2 = ω u i 1 , q , r + 1 ( 2 a + 2 ) + u i + 1 , q , r + 1 ( 2 a + 1 ) + u i , q 1 , r + 1 ( 2 a + 1 ) + u i , q , r + 2 ( 2 a + 1 ) + ( 1 ω ) ( u i 1 , q , r + 1 ( 2 a + 1 ) + u i , q + 1 , r + 1 ( 2 a + 1 ) + u i , q , r ( 2 a + 1 ) ) , d 3 = ω u i 1 , q + 1 , r + 1 ( 2 a + 2 ) + u i + 1 , q + 1 , r + 1 ( 2 a + 1 ) + u i , q + 2 , r + 1 ( 2 a + 1 ) + u i , q + 1 , r + 2 ( 2 a + 1 ) + ( 1 ω ) ( u i 1 , q + 1 , r + 1 ( 2 a + 1 ) + u i , q , r + 1 ( 2 a + 1 ) + u i , q + 1 , r ( 2 a + 1 ) ) , d 4 = ω u i 1 , q + 1 , r ( 2 a + 2 ) + u i + 1 , q + 1 , r ( 2 a + 1 ) + u i , q + 2 , r ( 2 a + 1 ) + u i , q + 1 , r 1 ( 2 a + 1 ) + ( 1 ω ) ( u i 1 , q + 1 , r ( 2 a + 1 ) + u i , q , r ( 2 a + 1 ) + u i , q + 1 , r + 1 ( 2 a + 1 ) ) .
Then the numerical solution of such a set of inner boundary points can be calculated quickly only by solving the fourth order sparse linear Equation (15). In the same way, the points on the other five groups of inner boundary lines are also calculated by Group Explicit method, and we will not represent them one by one here.
(3) The solution to Interfaces III
It can be seen from the calculation process of Interfaces I and II we solve the Interfaces III just depending on the group explicit method based on two points a group. The specific calculation process is the same as above (1) and (2), and we do not repeat it.
Finally, taking the above results as the interface conditions, we use the asymmetric iterative schemes different from the schemes at the odd levels to solve the inner points on Ω i , i = 1 , 2 , , 8 .
L 7 u i , j , k ( 2 a + 2 ) = 0 , ( i , j , k ) Ω 1 \ π 1 , L 8 u i , j , k ( 2 a + 2 ) = 0 , ( i , j , k ) Ω 2 \ π 2 , L 5 u i , j , k ( 2 a + 2 ) = 0 , ( i , j , k ) Ω 3 \ π 3 , L 6 u i , j , k ( 2 a + 2 ) = 0 , ( i , j , k ) Ω 4 \ π 4 , L 3 u i , j , k ( 2 a + 2 ) = 0 , ( i , j , k ) Ω 5 \ π 5 , L 4 u i , j , k ( 2 a + 2 ) = 0 , ( i , j , k ) Ω 6 \ π 6 , L 1 u i , j , k ( 2 a + 2 ) = 0 , ( i , j , k ) Ω 7 \ π 7 , L 2 u i , j , k ( 2 a + 2 ) = 0 , ( i , j , k ) Ω 8 \ π 8 .
Through the above specific implementation process of the domain decomposition iteration algorithm, we can see that the calculation of numerical solution can transform implicit iteration to explicit calculation no matter on the odd or even iteration layer. Combining with the domain decomposition method, the alternating asymmetric iterative (AAI) algorithm is well performed in parallel.

4. The Algorithm Convergence

In the last section, we propose a new AAI algorithm based domain decomposition for solving the three-dimensional Poisson problem (1)–(2), which can be written in the following matrix form,
u ( n + 1 ) = T ω u ( n ) + b ,
while T ω is the iterative matrix of AAI algorithm and b is the right-term. Then we can give the following theorem:
Theorem 1.
The sufficient and necessary conditions for the convergence of AAI algorithm are as follows:
ρ ( T ω ) < 1 ,
where ρ ( T ω ) is the corresponding spectral radius of the iterative matrix T ω .
Consider the eigenvalue problem of Equation (17):
T ω x = λ x ,
Due to the asymmetry of the schemes (5) and (12) in the calculation direction, we take one of the iteration scheme (5) as an example,
λ u i , j , k 1 6 [ ( λ ω + 1 ω ) u i , j , k + ( λ ω + 1 ω ) u i , j 1 , k + ( λ ω + 1 ω ) u i , j , k 1 + u i + 1 , j , k + u i , j + 1 , k + u i , j , k + 1 ] = 0 .
Firstly, we give the relationship between the eigenvalues μ of Jacobi iterative matrix B and the eigenvalues λ of the AAI iterative matrix T ω . Let V i , j , k be the eigenvectors of the Jacobi iterative matrix, then
u i , j , k = [ ± ( λ ω + 1 ω ) 1 2 ] i + j + k V i , j , k .
Taking Equation (21) into Equation (20), we can obtain,
μ V i , j , k 1 6 [ V i 1 , j , k + V i , j 1 , k + V i , j , k 1 + V i + 1 , j , k + V i , j + 1 , k + V i , j , k + 1 ] = 0 .
where
μ = ± λ ( λ ω + 1 ω ) 1 2 .
If λ is the eigenvalue of the matrix T ω , then
μ 2 ( λ ω + 1 ω ) = λ 2 .
Equation (24) determines that μ is eigenvalue of the matrix B, which is Jacobi iteration matrix of Poisson equation. On the contrary, if μ is eigenvalue of the matrix B, it can be determined only if there is a relationship between the eigenvalues λ of Jacobi iteration matrix and the eigenvalues μ of the given iteration matrix in the Equation (24).
In particular, it is shown that the iterative schemes (5)–(12) are Gauss-Seidel iterative schemes in fact when ω = 1 . Then the presented AAI algorithm has obvious convergence since λ = μ 2 < 1 .
Second, we discuss the changes of ρ ( T ω ) about ω .
From Equation (24), we can see that the eigenvalue λ depends on the relaxation factor ω and the eigenvalue μ of Jacobi iteration matrix. Suppose 0 μ 1 , 0 < ω < 2 , the two eigenvalues are obtained by Equation (24):
λ 1 ( ω , μ ) = μ 2 ω 2 + μ ( μ ω 2 ) 2 ( ω 1 ) ,
λ 2 ( ω , μ ) = μ 2 ω 2 μ ( μ ω 2 ) 2 ( ω 1 ) ,
Define
M ( ω , μ ) = max { | λ 1 ( ω , μ ) | , | λ 2 ( ω , μ ) | } ,
by the discriminant equaling to zero, namely,
Δ = μ 4 ω 2 4 μ 2 ( ω 1 ) = 0 .
Then the root of the Equation (24) is
ω μ = 2 ( 1 1 μ 2 ) μ 2 , 0 < ω μ < 2 .
When 0 < ω < ω μ , Δ > 0 , we can get
λ 1 ( ω , μ ) > λ 2 ( ω , μ ) > 0 .
When ω μ < ω < 2 , the eigenvalue λ 1 ( ω , μ ) and λ 2 ( ω , μ ) are conjugate complex, therefore
| λ 1 ( ω , μ ) | = | λ 2 ( ω , μ ) | = μ 2 ( ω 1 ) .
Due to Equations (30) and (31), we can give
M ( ω , μ ) = λ 1 ( ω , μ ) , 0 < ω < ω μ , μ 2 ( ω 1 ) , ω μ < ω < 2 .
It is obviously seen that
M ( ω , μ ) < 1 .
In fact, if ω μ < ω < 2 , Equation (33) is ture clearly; Otherwise 0 < ω < ω μ , and
M ( ω , μ ) = λ 1 ( ω , μ ) < μ 2 ω 2 + μ ( μ ω 2 ) 2 μ ω + 1 , ( 0 < ω < ω μ ) = μ < 1 .
Therefore, ρ ( T ω ) < 1 , Equation (18) is proved and the presented AAI algorithm is convergent.
Obviously, the spectrum radius ρ ( T ω ) of the presented iterative matrix depends on the relaxation factor ω , so choosing approximate ω is important to the number of iterations and the convergence rate.
Since the optimal relaxation factor ω o p t is obtained for 2D Poisson problem in [25], we can also provide the same computation for 3D case. When
ω = ω o p t = 2 1 + 1 ρ ( B ) 2 + ε , ( ε > 0 ) ,
ρ ( T ω o p t ) obtains the minimum
ρ ( T ω o p t ) = ( 1 1 ρ ( B ) 2 ) 2 + ε ρ ( B ) 2 , ( ε > 0 ) ,
where ρ ( B ) is the spectrum radius of Jacobi iterative matrix, and ε is a positive, sufficiently small number. The optimal relaxation factor ω o p t can be theoretically evaluated by Equation (36).

5. Numerical Experiments

In order to confirm the effectiveness of the AAI algorithm, the following experiments are carried out.
The initial iterative values u i , j , k ( 0 ) = 0 ( i = 1 , 2 , , l 1 ; j = 1 , 2 , , m 1 ; k = 1 , 2 , , s 1 ) is given.
(1) Consider the 3D Laplace equation
2 u x 2 + 2 u y 2 + 2 u z 2 = 0 , ( x , y , z ) [ 0 , 1 ] 3 ,
with the boundary condition,
u ( 0 , y , z ) = s i n ( y + z ) ; u ( 1 , y , z ) = e x p ( 2 ) s i n ( y + z ) ; u ( x , 0 , z ) = e x p ( 2 x ) s i n ( z ) ; u ( x , 1 , z ) = e x p ( 2 x ) s i n ( 1 + z ) ; u ( x , y , 0 ) = e x p ( 2 x ) s i n ( y ) ; u ( x , y , 1 ) = e x p ( 2 x ) s i n ( y + 1 ) .
The exact solution of the 3D Poisson problem (37)–(38) is u ( x , y , z ) = e x p ( 2 x ) s i n ( y + z ) . Let u ( x i , y j , z k ) be the exact solution and u i , j , k ( n ) the nth iterative solution, the errors are calculated in L -norm as:
E ( n ) , h = max i , j , k ( e h ( n ) ( i , j , k ) ) = max i , j , k | u ( x i , y j , z k ) u i , j , k ( n ) | .
Moreover, the rate of convergence in space is calculated by
R a t e o f c o n v e r g e n c e l o g ( E , h 1 / E , h 2 ) l o g ( h 1 / h 2 ) .
where h 1 , h 2 are the space steps.
Table 1 gives the errors E of the presented alternating asymmetric iteration algorithm based on domain decomposition for the 3D Laplace problem (37)–(38) with different values of ω when l = m = s = 31 , h = 1 / 30 , n = 150 , we can obviously see that the errors is relatively smaller when the relaxation factor ω is about 1.9 . we further see that the errors get the minimum when ω is about 1.82 shown in Figure 4a, which is match with the result of Equation (35). Figure 4b performs the errors with z = 0.5 when ω = 1.82 , which illustrate the effectiveness of the AAI algorithm.
From Table 2 and Table 3, we can see the iteration numbers of the AAI algorithm is the least during the Jacobi, Gauss-Seidel iterative methods under some error controls when h = 1 / 30 , 1 / 50 . In addition, the AAI algorithm obtains shorter times than the Jacobi and the Gauss-Seidel methods when the number of the grid nodes is in increasing. Table 4 gives the convergence rates and errors of the AAI algorithm. In computation of the rates of convergence in space, the spatial steps are taken as h = 1 / ( 16 + 8 d ) , d = 0 , 1 , , 4 . We can see the rates is of order 2 in space and the errors can up to 10 5 .
Figure 5 provides the errors, relative errors and numerical solutions at z = 1 / 3 , 2 / 3 when l = m = s = 31 , h = 1 / 30 , n = 150 , ω = 1.82 , which shows the AAI algorithm is effect and accurate.
(2) Consider the 3D Poisson equation
2 u x 2 + 2 u y 2 + 2 u z 2 = 3 s i n ( x ) s i n ( y ) s i n ( z ) , ( x , y , z ) [ 0 , 1 ] 3 ,
with the boundary condition,
u ( 0 , y , z ) = 0 ; u ( 1 , y , z ) = s i n ( 1 ) s i n ( y ) s i n ( z ) ; u ( x , 0 , z ) = 0 ; u ( x , 1 , z ) = s i n ( x ) s i n ( 1 ) s i n ( z ) ; u ( x , y , 0 ) = 0 ; u ( x , y , 1 ) = s i n ( x ) s i n ( y ) s i n ( 1 ) .
The exact solution of the 3D Poisson problem (40)–(41) is u ( x , y , z ) = s i n ( x ) s i n ( y ) s i n ( z ) .
Table 5 gives the errors E of AAI algorithm for the problem 2 with the different values of ω when l = m = s = 31 , h = 1 / 30 , n = 150 , and Figure 6 shows the errors get nearly the minimum 2.2319 × 10 4 while ω = 1.76 . which show the effect of ω to the AAI algorithm. Table 6, Table 7 and Table 8 give the iteration numbers and times under some error controls, which also show obviously the presented algorithm has smaller iteration numbers than the Jacobi and Gauss-Seidel methods. The computational times are shorter with the grid points increasing.
Figure 7 shows the errors, relative errors and numerical solutions of the AAI iterative algorithm based on domain decomposition for the problem 2 at z = 1 / 3 , 2 / 3 when l = m = s = 51 , h = 1 / 50 , n = 150 , ω = 1.96 . All of the numerical experiments examine the effectiveness and accuracy of the presented AAI algorithm.
Since the presented AAI algorithm based on domain decomposition is constructed by the asymmetrical iterative schemes and GE method, which has interior parallelism and is easy to be implemented. During the process of the implementation of the AAI algorithm, there’s no need to solve the large-scale sparse block tridiagonal matrices. When the iteration numbers are odd or even, we just to solve the 3D problems by the constructed iterative schemes, which are computed independently. Once the interfaces are solved by GE method, the other grid points can also be solved by the constructed iterative schemes directly and in parallel. Therefore, the key of parallel implementation lies in information transfer and time cost of the inner boundary. In fact, the whole computation is explicit but convergent, which save the most of the consuming time.
So In this paper, we use the Matlab software to implement the presented AAI algorithm. We extend this idea to solve the other time-dependent high-dimensional problems, and compare the times, speedup, caches and so on. The detailed parallel implementation and performance analysis are provided in [27].

6. Conclusions

In this paper, we provide a new alternating asymmetric iterative (AAI) algorithm for 3D Poisson problem based on domain decomposition. We use several asymmetrical iterative schemes to solve the sub-domains respectively. Meanwhile the asymmetrical iterative schemes are alternatively used on odd and even iteration levels to improve the accuracy. Moreover, we give the convergence of the algorithm and the optimal relaxation factor. Finally, several numerical experiments are taken to examine the effectiveness and accuracy of the presented algorithm.
The study will be extended to other high-dimensional diffusion problems and wave problems and so on, and also can be used to solve on more multi-subdomains, and the corresponding new algorithms will be designed. we will report these soon.

Author Contributions

Conceptualization, and Methodology and Writing—original draft preparation, Q.X.; Formal analysis and Writing—review, Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

The research of the first author was partially supported by the Natural Science Foundations of China (No. 61662059), the Major Innovation Projects for Building First-calss Universities in China’s Western Region (No. ZKZD2017009), the Natural Science Foundations of Ningxia Province (No. NZ2018AAC03026), and the Fourth Batch of the Ningxia Youth Talents Supporting Program (No.TJGC2019012). The research of the second author was partially supported by the Natural Science Foundations of China (No. 11501313), the Natural Science Foundations of Ningxia Province (No. 2019AAC02001), the Project funded by the China Postdoctoral Science Foundation (No. 2017M621343), and the Third Batch of the Ningxia Youth Talents Supporting Program (No. TJGC2018037).

Acknowledgments

The authors thank the Reviewers and the Editors for their valuable comments and suggestions on an earlier version of this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Saad, Y. Iterative Methods for Sparse Linear Systems, 2nd ed.; SIAM: Philadelphia, PA, USA, 2003. [Google Scholar]
  2. Ge, Y.; Cao, F.; Zhang, J. A transformation-free HOC scheme and multigrid method for solving the 3D Poisson equation on nonuniform grids. J. Comput. Phys. 2013, 234, 199–216. [Google Scholar] [CrossRef]
  3. Abide, S.; Zeghmati, B. Multigrid defect correction and fourth-order compact scheme for Poisson equation. Comput. Math. Appl. 2017, 73, 1433–1444. [Google Scholar] [CrossRef] [Green Version]
  4. Fen, H.; Zhang, B.L.; Liu, Y. Mathematics Stencil of the finite difference method for Poisson equation and its application. China Sci. 2005, 35, 901–909. (In Chinese) [Google Scholar]
  5. Spotz, W.F.; Carey, G.F. A high-order compact formulation for the 3D Poisson equation. Numer. Methods Part. Diff. Equ. 1996, 12, 235–243. [Google Scholar] [CrossRef]
  6. Gupta, M.M.; Kouatchou, J. Symbolic derivation of finite difference approximations for the three dimensional Poisson equation. Numer. Methods Part. Diff. Equ. 1998, 14, 593–606. [Google Scholar] [CrossRef]
  7. Romao, E.C.; Campos, M.D.; Moura, L.F.M. Application of the galerkin and least-squares finite element methods in the solution of 3D Poisson and Helmholtz equations. Comput. Math. Appl. 2011, 62, 4288–4299. [Google Scholar] [CrossRef] [Green Version]
  8. Nintcheu Fata, S. Semi-analytic treatment of the three-dimensional Poisson equation via a Galerkin BIE method. J. Comput. Appl. Math. 2011, 236, 1216–1225. [Google Scholar] [CrossRef] [Green Version]
  9. Shi, Z.; Cao, Y.Y.; Chen, Q.J. Solving 2D and 3D Poisson equations and biharmonic equations by the Haar wavelet method. Appl. Math. Model. 2012, 36, 5143–5161. [Google Scholar] [CrossRef]
  10. Shi, Z.; Cao, Y.Y. A spectral collocation method based on Haar wavelets for Poisson equations and biharmonic equations. Appl. Math. Model. 2011, 54, 2858–2868. [Google Scholar] [CrossRef]
  11. Speyer, G.; Vasileska, D.; Goodnick, S.M. Efficient Poisson equation solvers for large scale 3D simulations. In Proceedings of the 2001 International Conference on Modeling and Simulation of Microsystems-MSM, Hilton Head Island, SC, USA, 19–21 March 2001; pp. 23–26. [Google Scholar]
  12. Pedro, V.-L.; Alfredo, P.; Manuel, P.-M. Fast finite difference Poisson solvers on heterogeneous architectures. Comput. Phys. Commun. 2014, 185, 1265–1272. [Google Scholar]
  13. Pedro, V.-L.; Alfredo, P.; Julien, F.; Manuel, P.-M. Block tridiagonal solvers on heterogeneous architectures. In Proceedings of the 2012 IEEE 10th International Symposium on Parallel and Distributed Processing with Applications, Leganes, Spain, 10–13 July 2012; pp. 609–617. [Google Scholar]
  14. Stone, C.P.; Duque, E.P.N.; Zhang, Y.; Car, D.; Owens, J.D.; Davis, R.L. 20th AIAA Computational Fluid Dynamics Conference. AIAA-306 2011, 3221, 307. [Google Scholar]
  15. Yang, W.D.; Li, K.L.; Li, K.Q. A parallel solving method for block-tridiagonal equations on CPUCGPU heterogeneous computing systems. J. Supercomput. 2017, 73, 1760–1781. [Google Scholar] [CrossRef]
  16. Zhang, Y.; Cohen, J.; Owens, J.D. Proceedings of the 15th ACM SIGPLAN Symposium on Principles and Practice of Parallel Programming; ACM: New York, NY, USA, 2010; p. 127. [Google Scholar]
  17. Davidson, A.; Zhang, Y.; Owens, J.D. IEEE International Parallel and Distributed Processing Symposium; IEEE: Piscataway, NJ, USA, 2011; p. 956. [Google Scholar]
  18. Zhang, B.L.; Yuan, G.X.; Liu, X.P. Parallel Finite Difference Methods for Partial Differential Equations; Science Press: Beijing, China, 1994; pp. 7–67. [Google Scholar]
  19. Mohanty, R.K. The samrt-BLAGE algorithm for singularly perturbed 2D elliptic partial differential equations. Appl. Math. Comput. 2007, 190, 321–331. [Google Scholar]
  20. Tavakoli, R.; Davami, P. A new parallel Gauss-Seidel method based on alternating group explicit method and domain decomposition method. Appl. Math. Comput. 2007, 188, 713–719. [Google Scholar] [CrossRef]
  21. Abdullah, A.R.; Ali, N.M. The comparative study of parallel alternating-type iterative methods. Appl. Math. Comput. 1996, 74, 331–344. [Google Scholar]
  22. Keyes, D.E.; Gropp, W.D. A comparison of domain decomposition on techniques for elliptic partial differential equations and their parallel implementation. SIAM J. Sci. Stat. Comput. 1987, 8, s166–s202. [Google Scholar] [CrossRef]
  23. Kuznetsov, Y. New algorithm for approximate realization of implicit difference scheme. Soviet J. Numer. Anal. Math. Model. 1988, 3, 99–114. [Google Scholar] [CrossRef]
  24. Xu, Q.Y.; Wang, W.Q. Wenqiang, W. A new parallel iterative algorithm for solving 2D Poisson equation. Numer. Methods Part. Differ. Eqs. 2011, 27, 829–853. [Google Scholar] [CrossRef]
  25. Ng, K.F.; Mohd Ali, N.H. Performance analysis of explicit group parallel algorithms for distributed memory multicomputer. Parallel Comput. 2008, 34, 427–440. [Google Scholar] [CrossRef]
  26. Evans, D. Alternating group explicit method for the diffusion equations. Appl. Math. Model. 1985, 9, 201–206. [Google Scholar] [CrossRef]
  27. Xu, Q.Y.; Liu, Z.Y. On the parallel implementation of several iterative algorithms for the 3D Poisson and Diffusion problems. 2020; in preparation. [Google Scholar]
Figure 1. The asymmetric iterative schemes (5)–(12) for the 3D Poisson equation with the relaxation factor ω .
Figure 1. The asymmetric iterative schemes (5)–(12) for the 3D Poisson equation with the relaxation factor ω .
Mathematics 08 00281 g001
Figure 2. The solution domain is divided into eight sub-domains.
Figure 2. The solution domain is divided into eight sub-domains.
Mathematics 08 00281 g002
Figure 3. The computation of the interfaces I , I I , I I I .
Figure 3. The computation of the interfaces I , I I , I I I .
Mathematics 08 00281 g003
Figure 4. The E with the different ω and the errors at z = 0.5 when l = m = s = 31 , h = 1 / 30 , n = 150 . (a) E ; (b) the errors at z = 0.5 when ω = 1.82 .
Figure 4. The E with the different ω and the errors at z = 0.5 when l = m = s = 31 , h = 1 / 30 , n = 150 . (a) E ; (b) the errors at z = 0.5 when ω = 1.82 .
Mathematics 08 00281 g004
Figure 5. The errors, relative errors and numerical solutions of the AAI algorithm at z = 1 / 3 , 2 / 3 when l = m = s = 51 , h = 1 / 50 , n = 150 , ω = 1.94 . (a) The errors at z = 1 / 3 ; (b) The errors at z = 2 / 3 ; (c) The relative errors at z = 1 / 3 ; (d) The relative errors at z = 2 / 3 ; (e) The numerical solution at z = 1 / 3 ; (f) The numerical solution at z = 2 / 3 .
Figure 5. The errors, relative errors and numerical solutions of the AAI algorithm at z = 1 / 3 , 2 / 3 when l = m = s = 51 , h = 1 / 50 , n = 150 , ω = 1.94 . (a) The errors at z = 1 / 3 ; (b) The errors at z = 2 / 3 ; (c) The relative errors at z = 1 / 3 ; (d) The relative errors at z = 2 / 3 ; (e) The numerical solution at z = 1 / 3 ; (f) The numerical solution at z = 2 / 3 .
Mathematics 08 00281 g005
Figure 6. The E with the different ω and the errors at z = 0.5 when l = m = s = 31 , h = 1 / 30 , n = 150 . (a) The errors E with the different ω ; (b) The errors at z = 0.5 when ω = 1.76 .
Figure 6. The E with the different ω and the errors at z = 0.5 when l = m = s = 31 , h = 1 / 30 , n = 150 . (a) The errors E with the different ω ; (b) The errors at z = 0.5 when ω = 1.76 .
Mathematics 08 00281 g006
Figure 7. The errors and numerical solutions of the AAI algorithm when l = m = s = 51 , h = 1 / 50 , n = 150 , ω = 1.96 . (a) The errors at z = 1 / 3 ; (b) The errors at z = 2 / 3 ; (c) The relative errors at z = 1 / 3 ; (d) The relative errors at z = 2 / 3 ; (e) The numerical solution at z = 1 / 3 ; (f) The numerical solution at z = 2 / 3 .
Figure 7. The errors and numerical solutions of the AAI algorithm when l = m = s = 51 , h = 1 / 50 , n = 150 , ω = 1.96 . (a) The errors at z = 1 / 3 ; (b) The errors at z = 2 / 3 ; (c) The relative errors at z = 1 / 3 ; (d) The relative errors at z = 2 / 3 ; (e) The numerical solution at z = 1 / 3 ; (f) The numerical solution at z = 2 / 3 .
Mathematics 08 00281 g007
Table 1. The errors E of alternating asymmetric iterative (AAI) algorithm based on domain decomposition with the different values of ω when l = m = s = 31 , h = 1 / 30 , n = 150 .
Table 1. The errors E of alternating asymmetric iterative (AAI) algorithm based on domain decomposition with the different values of ω when l = m = s = 31 , h = 1 / 30 , n = 150 .
ω E ω E
1.04.4434 × 10 1 1.63.9316 × 10 2
1.13.6971 × 10 1 1.71.0478 × 10 2
1.22.9383 × 10 1 1.84.8059 × 10 4
1.32.1896 × 10 1 1.92.0334 × 10 4
1.41.4838 × 10 1 2.0errors
1.58.6993 × 10 2
Table 2. The iteration numbers and time(s) of Jacobi, Gauss-Seidel and the AAI algorithm under the different error controls when l = m = s = 31 , h = 1 / 30 .
Table 2. The iteration numbers and time(s) of Jacobi, Gauss-Seidel and the AAI algorithm under the different error controls when l = m = s = 31 , h = 1 / 30 .
E 10 3 10 4 10 5
NumbersTime(s)NumbersTime(s)NumbersTime(s)
Jacobi9873.390614014.421917745.0000
Gauss-Seidel4912.42196983.14068853.5469
AAI ( ω = 1.82 ) 936.89061239.640614510.1093
Table 3. The iteration numbers and time(s) of Jacobi, Gauss-Seidel and the AAI algorithm under the different error controls when l = m = s = 51 , h = 1 / 50 .
Table 3. The iteration numbers and time(s) of Jacobi, Gauss-Seidel and the AAI algorithm under the different error controls when l = m = s = 51 , h = 1 / 50 .
E 10 3 10 4 10 5
NumbersTime(s)NumbersTime(s)NumbersTime(s)
Jacobi274521.5156390534.5000501837.9219
Gauss-Seidel136817.2344194820.5938250535.4688
AAI ( ω = 1.94 ) 1159.437520918.343728322.1250
Table 4. The convergence rates and errors E of the presented iterative algorithm when ω = 1.95 , n = 3000 .
Table 4. The convergence rates and errors E of the presented iterative algorithm when ω = 1.95 , n = 3000 .
Numbers 16 × 16 × 16 24 × 24 × 24 32 × 32 × 32 40 × 40 × 40 48 × 48 × 48
Rates1.99091.99061.99561.9979
Errors1.4277 × 10 4 6.3687 × 10 5 3.5920 × 10 5 2.3011 × 10 5 1.5986 × 10 5
Table 5. The errors E of AAI algorithm based on domain decomposition with the different values of ω when l = m = s = 31 , h = 1 / 30 , n = 150 .
Table 5. The errors E of AAI algorithm based on domain decomposition with the different values of ω when l = m = s = 31 , h = 1 / 30 , n = 150 .
ω E ω E
1.04.2402 × 10 2 1.63.6712 × 10 3
1.13.4981 × 10 2 1.79.4888 × 10 4
1.22.7756 × 10 2 1.83.9125 × 10 4
1.32.0657 × 10 2 1.95.8138 × 10 4
1.41.3986 × 10 2 2.0errors
1.58.1573 × 10 3
Table 6. The iteration numbers and time(s) of Jacobi, Gauss-Seidel and the AAI algorithm under the different error controls when l = m = s = 31 , h = 1 / 30 .
Table 6. The iteration numbers and time(s) of Jacobi, Gauss-Seidel and the AAI algorithm under the different error controls when l = m = s = 31 , h = 1 / 30 .
E 10 3 10 4
NumbersTime(s)NumbersTime(s)
Jacobi5572.67199763.0781
Gauss-Seidel2880.89064981.4688
AAI ( ω = 1.96 ) 295.671913513.9531
Table 7. The iteration numbers and time(s) of Jacobi, Gauss-Seidel and the AAI algorithm under the different error controls when l = m = s = 51 , h = 1 / 50 .
Table 7. The iteration numbers and time(s) of Jacobi, Gauss-Seidel and the AAI algorithm under the different error controls when l = m = s = 51 , h = 1 / 50 .
E 10 3 10 4
NumbersTime(s)NumbersTime(s)
Jacobi154815.7031271224.1875
Gauss-Seidel7908.3594137215.2500
AAI ( ω = 1.96 ) 5114.546918130.2344
Table 8. The iteration numbers and time(s) of Jacobi, Gauss-Seidel and the AAI algorithm under the different error controls when l = m = s = 71 , h = 1 / 70 .
Table 8. The iteration numbers and time(s) of Jacobi, Gauss-Seidel and the AAI algorithm under the different error controls when l = m = s = 71 , h = 1 / 70 .
E 10 3 10 4
NumbersTime(s)NumbersTime(s)
Jacobi303498.54685317176.6562
Gauss-Seidel153853.5156268192.9531
AAI ( ω = 1.96 ) 9131.718720154.2968

Share and Cite

MDPI and ACS Style

Xu, Q.; Liu, Z. Alternating Asymmetric Iterative Algorithm Based on Domain Decomposition for 3D Poisson Problem. Mathematics 2020, 8, 281. https://doi.org/10.3390/math8020281

AMA Style

Xu Q, Liu Z. Alternating Asymmetric Iterative Algorithm Based on Domain Decomposition for 3D Poisson Problem. Mathematics. 2020; 8(2):281. https://doi.org/10.3390/math8020281

Chicago/Turabian Style

Xu, Qiuyan, and Zhiyong Liu. 2020. "Alternating Asymmetric Iterative Algorithm Based on Domain Decomposition for 3D Poisson Problem" Mathematics 8, no. 2: 281. https://doi.org/10.3390/math8020281

APA Style

Xu, Q., & Liu, Z. (2020). Alternating Asymmetric Iterative Algorithm Based on Domain Decomposition for 3D Poisson Problem. Mathematics, 8(2), 281. https://doi.org/10.3390/math8020281

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