arxiv: v3 [math.na] 1 Jan 2015

Transcription

1 On Solving Pentadiagonal Linear Systems via Transformations arxiv: v3 [mathna] 1 Jan 2015 A A KARAWIA Computer Science Unit Deanship of Educational Services Qassim University POBox 6595 Buraidah Saudi Arabia kraoieh@quedusa Abstract Many authors studied numerical algorithms for solving the linear systems of the pentadiagonal type The well-known Fast Pentadiagonal System Solver algorithm is an example of such algorithms The current article is describes new numerical and symbolic algorithms for solving pentadiagonal linear systems via transformations New algorithms are natural generalization of the work presented in [Moawwad El-Mikkawy and Faiz Atlan Algorithms for Solving Linear Systems of Equations of Tridiagonal Type via Transformations Applied Mathematics ] The symbolic algorithms remove the cases where the numerical algorithms fail The computational cost of our algorithms is given Some examples are given in order to illustrate the effectiveness of the proposed algorithms All of the experiments are performed on a computer with the aid of programs written in MATLAB Keywords:Pentadiagonal matrix; Backward pentadiagonal; Permutation matrix; Linear systems; Algorithm; MATLAB AMS Subject Classification:15A15; 15A23; 68W30; 11Y05; 33F10; F21; G10 1 Introduction The pentadiagonal linear systems denoted by (PLS) take the forms: PX = Y (11) where P is n by n pentadiagonal matrix given by d 1 a 1 b c 2 d 2 a 2 b e 3 c 3 d 3 a 3 b e 4 c 4 d 4 a 4 b P = n 4 (12) 0 0 e n 2 c n 2 d n 2 a n 2 b n e n 1 c n 1 d n 1 a n e n c n d n and X = (x 1 x 2 x n ) t Y = (y 1 y 2 y n ) t are vectors of length n This kind of linear systems is well known in the literature [1-11] and often these types of linear systems are widely used in areas of science and engineering for example in numerical solution of ordinary and partial differential equations (ODE and PDE) interpolation problems boundary value problems (BVP) parallel Home Address: Mathematics Department Faculty of Science Mansoura University Mansoura Egypt E- mail:abibka@mansedueg

2 computing Physics matrix algebra[4-13] In [7] the author presented an efficient computational algorithm for solving periodic pentadiagonal linear systems The algorithm is based on the LU factorization of the periodic pentadiagonal matrix New algorithms are described for solving periodic pentadiagonal linear systems based on the use of any pentadiagonal linear solver and the author described a symbolic algorithm for solving pentadiagonal linear systems[8] In [9] the authors discussed the general nonsymmetric problem and proposed an algorithm for solving nonsymmetric penta-diagonal Toeplitz linear systems A fast algorithm for solving a large system with a symmetric Toeplitz penta-diagonal coefficient matrix is presented[10] This efficient method is based on the idea of a system perturbation followed by corrections and is competitive with standard methods In [11] the authors described an efficient computational algorithm and symbolic algorithm for solving nearly pentadiagonal linear systems based on the LU factorization of the nearly pentadiagonal matrix In this paper we show that more efficient algorithms are derived via transformations that can be regarded as a natural generalization of the efficient algorithms in [14] The current paper is organized as follows: In Section 2 new numerical algorithms for solving a pentadiagonal linear system are presented New symbolic algorithms for solving a pentadiagonal linear system are constructed in Section 3 In Section 4 three illustrative examples are presented Conclusions of the work are given in Section 5 2 Numerical Algorithms for Solving PLS In this section we shall focus on the construction of new numerical algorithms for computing the solution of pentadiagonal linear system For this purpose it is convenient to give five vectors α = (α 1 α 2 α n 1 ) β = (β 1 β 2 β n 2 ) Z = (z 1 z 2 z n ) γ = (γ 2 γ 3 γ n ) and µ = ( µ n ) where α i = { a1 i = 1 a i β i 1γ i i = 23n 1 (21) and β i = b i i = 23n 1 (22) y 1 i = 1 y z i = 2 z 1γ 2 i = 2 (23) y i z i 2e i z i 1γ i i = 34n { γ i = c 2 i = 2 c i α i 2 e i i = 34n d 1 i = 1 = d 2 α 1 γ 2 i = 2 d i β i 2 e i α i 1 γ i i = 34n By using the vectors α β Z γ and µ together with the suitable elementary row operations we see that the system (11) may be transformed to the equivalent linear system: 1 α 1 β α 2 β α 3 β α 4 β α n 2 β n α n x 1 x 2 x 3 x 4 x n 2 x n 1 x n = z 1 z 2 z 3 z 4 z n 2 z n 1 z n (24) (25) (26)

3 The transformed system (26) is easy to solve by a backward substitution Consequently the PLS (11) can be solved using the following algorithm: Algorithm 21 First numerical algorithm for solving pentadiagonal linear system To find the solution of PLS (11) using the transformed system (26) we may proceed as follows: INPUT order of the matrix n and the components d i a i b i c i e i f i i = 12n(a n = b n = b n 1 = c 1 = e 1 = e 2 = 0) OUTPUT The solution vector x = (x 1 x 2 x n ) t Step 1: Use DETGPENTA algorithm [13] to check the non-singularity of the coefficient matrix of the system (13) Step 2: If det(p) = 0 then Exit and Print Message ( No solutions ) end if Step 3: Set = d 1 α 1 = a1 β 1 = b1 and z 1 = y1 Step 4: Set γ 2 = c 2 = d 2 α 1 γ 2 α 2 = a2 β1γ2 β 2 = b2 and z 2 = y2 z1γ2 Step 5: For i=34n-2 do γ i = c i α i 2 e i = d i β i 2 e i α i 1 γ i α i = ai βi 1γi β i = bi z i = yi zi 2ei zi 1γi End do γ n 1 = c n 1 α n 3 e n 1 µ n 1 = d n 1 β n 3 e n 1 α n 2 γ n 1 α n 1 = an 1 βn 2γn 1 µ n 1 γ n = c n α n 2 e n µ n = d n β n 2 e n α n 1 γ n z n 1 = yn 1 zn 2en 1 zn 2γn 1 µ n 1 z n = yn zn 1en zn 1γn µ n Step 6: Compute the solution vector X = (x 1 x 2 x n ) t using x n = z n x n 1 = z n 1 α n 1 x n For i=n-2 n-3 1 do x i = z i α i x i+1 β i x i+2 End do The numerical Algorithm 21 will be referred to as PTRANS-I algorithm The computational cost of PTRANS-I algorithm is 19n 29 operations The conditions 0i = 12n are sufficient for its validity In a similar manner we may consider five vectors σ = (σ 2 σ 3 σ n ) φ = (φ 3 φ 4 φ n ) W = (w 1 w 2 w n ) ρ = (ρ 1 ρ 2 ρ n 1 ) and ψ = (ψ 1 ψ 2 ) where σ i = { cn c i φ i+1ρ i i = n i = n 1n 22 (27) φ i = e i i = nn 13 (28) y n i = n y n 1 w nρ n 1 w i = 1 i = n 1 y i w i+2b i w i+1ρ i i = n 2n 31 (29) { an 1 i = n 1 ρ i = a i σ i+2 b i i = n 2n 31 (210)

4 and d n = d n 1 σ n ρ n 1 d i φ i+2 b i σ i+1 ρ i i = n i = n 1 i = n 2n 31 (211) Now we will present another algorithm for solving PLS As in PTRANS-I algorithm by using the vectors σ φ W ρ and ψ together with the suitable elementary row operations we see that the system (11) may be transformed to the equivalent linear system: σ φ 3 σ φ 4 σ φ n 2 σ n φ n 1 σ n φ n σ n 1 x 1 x 2 x 3 x 4 x n 2 x n 1 x n = w 1 w 2 w 3 w 4 w n 2 w n 1 w n (212) The transformed system (212) is easy to solve by a forward substitution Consequently the PLS (11) can be solved using the following algorithm: Algorithm 22 Second numerical algorithm for solving pentadiagonal linear system To find the solution of PLS (11) using the transformed system (212) we may proceed as follows: INPUT order of the matrix n and the components d i a i b i c i e i f i i = 12n(a n = b n = b n 1 = c 1 = e 1 = e 2 = 0) OUTPUT The solution vector x = (x 1 x 2 x n ) t Step 1: Use DETGPENTA algorithm [13] to check the non-singularity of the coefficient matrix of the system (13) Step 2: If det(p) = 0 then Exit and Print Message ( No solutions ) end if Step 3: Set = d n σ n = cn φ n = en and w n = yn Step 4: Set ρ n 1 = a n 1 1 = d n 1 σ n ρ n 1 σ n 1 = cn 1 φnρn 1 y n 1 w nρ n 1 1 Step 5: For i=n-2 n-3 3 do ρ i = a i σ i+2 b i = d i φ i+2 b i σ i+1 ρ i σ i = ci φi+1ρi φ i = ei w i = yi wi+2bi wi+1ρi End do ρ 2 = a 2 σ 4 b 2 ψ 2 = d 2 φ 4 b 2 σ 3 ρ 2 σ 2 = c2 φ4ρ2 ψ 2 ρ 1 = a 1 σ 3 b 1 ψ 1 = d 1 φ 3 b 1 σ 2 ρ 1 w 2 = y2 w4b2 w3ρ2 ψ 2 w 1 = y1 w3b1 w2ρ1 ψ 1 Step 6: Compute the solution vector X = (x 1 x 2 x n ) t using x = w 1 x 2 = w 2 σ 2 x 1 For i=3 4 n do 1 φ n 1 = en 1 1 and w n 1 =

5 x i = w i σ i x i 1 φ i x i 2 End do The numerical Algorithm 22 will be referred to as PTRANS-II algorithm The computational cost of PTRANS-II algorithm is 19n 29 operations Also the conditions 0i = 12n are sufficient for its validity If = 0 or = 0 for any i {12n} then PTRANS-I and PTRANS-II algorithm fail to solve pentadiagonal linear systems respectively So in the next section we developed two symbolic algorithms in order to remove the cases where the numerical algorithms fail The parameter p in the following symbolic algorithms is just a symbolic name It is a dummy argument and its actual value is zero 3 Symbolic Algorithms for Solving PLS In this section we shall focus on the construction of new symbolic algorithms for computing the solution of pentadiagonal linear systems The following algorithm is a symbolic version of PTRANS-I algorithm: Algorithm 31 First symbolic algorithm for solving pentadiagonal linear system To find the solution of PLS (11) using the transformed system (26) we may proceed as follows: INPUT order of the matrix n and the components d i a i b i c i e i f i i = 12n(a n = b n = b n 1 = c 1 = e 1 = e 2 = 0) OUTPUT The solution vector x = (x 1 x 2 x n ) t Step 1: Use DETGPENTA algorithm [13] to check the non-singularity of the coefficient matrix of the system (13) Step 2: If det(p) = 0 then Exit and Print Message ( No solutions ) end if Step 3: Set = d 1 If = 0 then = p end if Step 4: Set α 1 = a1 β 1 = b1 z 1 = y1 and γ 2 = c 2 Step 5: Set = d 2 α 1 γ 2 If = 0 then = p end if Step 6: Set α 2 = a2 β1γ2 β 2 = b2 and z 2 = y2 z1γ2 Step 7: For i=34n-2 do γ i = c i α i 2 e i = d i β i 2 e i α i 1 γ i If = 0 then = p end if α i = ai βi 1γi β i = bi z i = yi zi 2ei zi 1γi End do γ n 1 = c n 1 α n 3 e n 1 µ n 1 = d n 1 β n 3 e n 1 α n 2 γ n 1 If µ n 1 = 0 then µ n 1 = p end if α n 1 = an 1 βn 2γn 1 µ n 1 γ n = c n α n 2 e n µ n = d n β n 2 e n α n 1 γ n If µ n = 0 then µ n = p end if z n 1 = yn 1 zn 2en 1 zn 2γn 1 µ n 1 z n = yn zn 1en zn 1γn µ n Step 8: Compute the solution vector X = (x 1 x 2 x n ) t using x n = z n x n 1 = z n 1 α n 1 x n For i=n-2 n-3 1 do x i = z i α i x i+1 β i x i+2 End do Step 9: Substitute p = 0 in all expressions of the solution vector x i i = 12n The symbolic Algorithm 31 will be referred to as SPTRANS-I algorithm Now we are going to give the symbolic version of PTRANS-II algorithm:

6 Algorithm 32 Second symbolic algorithm for solving pentadiagonal linear system To find the solution of PLS (11) using the transformed system (212) we may proceed as follows: INPUT order of the matrix n and the components d i a i b i c i e i f i i = 12n(a n = b n = b n 1 = c 1 = e 1 = e 2 = 0) OUTPUT The solution vector x = (x 1 x 2 x n ) t Step 1: Use DETGPENTA algorithm [13] to check the non-singularity of the coefficient matrix of the system (13) Step 2: If det(p) = 0 then Exit and Print Message ( No solutions ) end if Step 3: Set = d n If = 0 then = p end if Step 4: σ n = cn φ n = en w n = yn and ρ n 1 = a n 1 Step 5: Set 1 = d n 1 σ n ρ n 1 If 1 = 0 then 1 = p end if Step 6: σ n 1 = cn 1 φnρn 1 1 Step 7: For i=n-2 n-3 3 do ρ i = a i σ i+2 b i = d i φ i+2 b i σ i+1 ρ i If = 0 then = p end if σ i = ci φi+1ρi φ i = ei φ n 1 = en 1 1 and w n 1 = yn 1 wnρn 1 1 w i = yi wi+2bi wi+1ρi End do ρ 2 = a 2 σ 4 b 2 ψ 2 = d 2 φ 4 b 2 σ 3 ρ 2 If ψ 2 = 0 then ψ 2 = p end if σ 2 = c2 φ4ρ2 ψ 2 ρ 1 = a 1 σ 3 b 1 ψ 1 = d 1 φ 3 b 1 σ 2 ρ 1 If ψ 1 = 0 then = p end if w 2 = y2 w4b2 w3ρ2 ψ 2 w 1 = y1 w3b1 w2ρ1 ψ 1 Step 8: Compute the solution vector X = (x 1 x 2 x n ) t using x = w 1 x 2 = w 2 σ 2 x 1 For i=3 4 n do x i = w i σ i x i 1 φ i x i 2 End do Step 9: Substitute p = 0 in all expressions of the solution vector x i i = 12n The symbolic Algorithm 32 will be referred to as SPTRANS-II algorithm Corollary 31 (generalization version of Corollary 21 in [14]) Let ˆP be the backward matrix of the pentadiagonal matrix P in (12) and given by: ˆP = 0 0 b 1 a 1 d b 2 a 2 d 2 c b 3 a 3 d 3 c 3 e b 4 a 4 d 4 c 4 e 4 0 b n 2 a n 2 d n 2 c n 2 e n a n 1 d n 1 c n 1 e n d n c n e n 0 0 n 4 (31)

7 Then the backward pentadiagonal linear system ˆPV = Y V = (v 1 v 2 v n ) t (32) has the solution: v i = x n i+1 i = 12 n where j is the floor function of j and X = (x 1 x 2 x n ) t is the solution vector of the linear system (11) Proof: Consider the n n permutation matrix M defined by: M = (33) For this matrix we have: Since M 1 = M (34) ˆP = PM (35) Then using (34) and (35) the result follows Corollary 32 (generalization version of Corollary 22 in [14]) The determinants of the coefficient matrices P and ˆP in (12) and (31) are given respectively by: and det(p) = det(ˆp) = ( 1) n(n 1) 2 n = i=1 n i=1 n (36) i=1 = ( 1) n(n 1) 2 where µ n and ψ 1 ψ 2 satisfy (25) and (211) respectively Proof: Using (26) (212) and (35) the result follows 4 ILLUSTRATIVE EXAMPLES n (37) In this section we are going to give three examples for the sake of illustration All experiments performed in MATLAB R2014a with an Intel(R) Core(TM) i7-4700mq CPU@240GHz 240 GHz Example 41(Case I: 0 and 0 for all i) Find the solution of the following pentadiagonal linear system of size x x x x x x x x x x 10 i=1 = Solution: We have n = 10 d = ( ) t a = ( ) t b = ( ) t c = ( ) t e = ( ) t and y = ( ) t i)- Applying the PTRANS-I algorithm it yields (41)

8 µ = ( )t det(p) = 10 i=1 = PTRANS-I(ndabcey)=( ) t ii)- Applying the PTRANS-II algorithm it yields ψ = ( )t det(p) = 10 i=1 = PTRANS-II(ndabcey)=( ) t Example 42(Case II: = 0 and = 0 for some i) Find the solution of the following pentadiagonal linear system of size x x x 3 = 3 9 (42) x 4 6 Solution: We have n = 4 d = (3 2 13) t a = (275) t b = (11) t c = (0 322) t e = (0031) t and y = (6396) t The numerical algorithms PTRANS-I and PTRANS-II fail to solve the pentadiagonal linear system (42) since = 0 i)- Applying the SPTRANS-I algorithm it yields µ = (3p 2 8p 21 p ) t det(p) = ( 4 i=1 ) p=0 = 126 SPTRANS-I(ndabcey)=(( (25p 42) (16p 42) 21 (8p 21) 21(p 2) 2(8p 21) (9p 21) (8p 21) )t ) p=0 = (1111) t ii)- Applying the SPTRANS-II algorithm it yields ψ = ( )t det(p) = 4 i=1 = 126 SPTRANS-II(ndabcey)=(1111) t Example 43 We consider the following n n pentadiagonal linear system in order to demonstrate the efficiency of algorithms 31 and x 1 x 2 x 3 x 4 x n 2 x n 1 x n = It can be verified that the exact solution is x = (111) t In Table 1 we give a comparison of the absolute error and running time between PTRANS-I PTRANS-II algorithms Algorithm 3[8] and A\b function in Matlab to compute x for different sizes

9 Table1 n x x and CPU time(s) PTRANS-I PTRANS-II Algorithm 3[8] A\b(MATLAB) Table1indicatesthat the valueofabsoluteerrorforlargevaluesofnforptrans-iiislessthan the other algorithms Also note that from Table 1 it is clear that the value of the runing time for Algorithm 3[8] is less than the other algorithms for large sizes Acknowledgements The author wishes to thank anonymous referees for useful comments that enhanced the quality of this paper The author is grateful to Prof Dr MEA El-Mikkawy for providing him with the references [14] and [15] 5 CONCLUSION There are many numerical algorithms in current use for solving linear systems of pentadiagonal type All numerical algorithms including the PTRANS-I and PTRANS-II algorithms of the current paper fail to solve the pentadiagonallinear system if = 0 and = 0 for any i {12n} The symbolic algorithms SPTRANS-I and SPTRANS-II of the current paper are constructed in order to remove the cases where the numerical algorithms fail From some numerical examples we have learned that SPTRANS-II algorithm works as well as Algorithm 3[8] and (A\y)MATLAB algorithms Hence it may become a useful tool for solving linear systems of pentadiagonal type References [1] TA Davis Direct Methods for Sparse Linear Systems SIAM Philadelphia 2006 [2] JW Demmel Applied Numerical Linear Algebra SIAM 1997 [3] GH Golub CF Van Loan Matrix Computations third ed The Johns Hopkins University Press Baltimore and London 1996 [4] J P Killingbeck G Jolicard The Folding Algorithm for Pentadiagonal Matrices Physics Letters A 166 (1992) [5] J Dongarra A Sameh On Some Parallel Banded Solvers Parallel Computing 1 (1984) [6] U Meier A Parallel Partition Method for Solving Banded Systems of Linear Equations Parallel Computing 2 (1985) [7] A A Karawia A computational algorithm for solving periodic pentadiagonal linear systems Applied Mathematics and Computation 174(2006) [8] T Sogabe New algorithms for solving periodic tridiagonal and periodic pentadiagonal linear systems Applied Mathematics and Computation 202(2008) [9] L Du T Sogabe S Zhang An algorithm for solving nonsymmetric penta-diagonal Toeplitz linear systemsoriginal Applied Mathematics and Computation 244(2014)10-15 [10] SS Nemani A fast algorithm for solving Toeplitz penta-diagonal systemsoriginal Applied Mathematics and Computation 215(2010) [11] A A Karawia and S A El-Shehawy Nearly Pentadiagonal Linear Systems Journal of Natural Sciences and Mathematics Qassim University 5(2012)

10 [12] J M McNally Fast Algorithm for Solving Diagonally Dominant Symmetric Pentadiagonal Toeplitz Systems Journal of Computational and Applied Mathematics 234 (2010) [13] T Sogabe A Fast Numerical Algorithm for the Determinant of a Pentadiagonal Matrix Applied Mathematics and Computation 196 (2008) [14] M El-Mikkawy F Atlan Algorithms for Solving Linear Systems of Equations of Tridiagonal Type via Transformations Applied Mathematics 5(2014) [15] M El-Mikkawy A Fast and Reliable Algorithm for Evaluating nth Order Pentadiagonal Determinants Applied Mathematics and Computation 202 (1) (2008)