Solution of Strictly Diagonal Dominant Tridiagonal Systems on Vector Computers
ثبت نشده
چکیده
24 restriction. As and , the minimum of δ (l) is achieved on the boundary of R. It is not difficult to see that this minimum is in x = y = , independently of the values of the greek letters. With this, we can see that: Now, if , we have: which allows to finish the proof: x y , () R ∈ 0 δ l 1 − () 1 − , [ ] 0 δ l 1 − () 1 − , [ ] × = x ∂ ∂δ l () 0 ≠ y ∂ ∂δ l () 0 ≠ δ l 1 − () 1 − δ l () δ l 1 − () a i 2 − a i 1 c i 2 − δ l 1 − () 1 − − c i 2 + c i 1 a i 2 + δ l 1 − () 1 − − + ≥ ∇ () = c i 2 ± a i 2 ± + δ l 1 − () 1 − = a i 2 − 1 c i 2 − δ l 1 − () 1 − − δ l 1 − () 1 − ≤ c i 2 + 1 a i 2 + δ l 1 − () 1 − − δ l 1 − () 1 − ≤ , ∇ () δ l 1 − δ l 1 − 1 − a i δ l 1 − 1 − c i + δ l 1 − 2 a i c i + δ l 1 − 3 ≥ ≥ ≥ 23 Also, the values of d' i-1 , d' i and d' i+1 can be found easily by applying transformations to the corresponding equations. With this, the diagonal dominance of equation i after the transformations of step l is exactly: Now, we want to bound δ (l) in terms of δ (l-1). We do that by finding the minimum for δ (l). This minimum bounds the diagonal dominance of the system at step l with respect to the equations at a given step l-1. Note that the worst case for δ (l) (δ (l) minimum) is when. All other cases do not cause so many cancellations as this case. Taken this into account, we know that in the worst case and. So, defining and we have where greek letters are defined …
منابع مشابه
OPM: a parallel method to solve banded systems of equations
studied in this report. We can see, in general, that low accuracies will require less overlapping and, in those cases, the algorithm will run faster. As the diagonal dominance increases, the amount of overlapping will be lower and the timings for different accuracies within the same diagonal dominance will get closer and closer. In this report we presented the Overlapped Partitions Method, a me...
متن کاملIncomplete Cyclic Reduction of Banded and Strictly Diagonally Dominant Linear Systems
The ScaLAPACK library contains a pair of routines for solving banded linear systems which are strictly diagonally dominant by rows. Mathematically, the algorithm is complete block cyclic reduction corresponding to a particular block partitioning of the system. In this paper we extend Heller’s analysis of incomplete cyclic reduction for block tridiagonal systems to the ScaLAPACK case. We obtain ...
متن کاملA Parallel Two-Level Hybrid Method for Diagonal Dominant Tridiagonal Systems
A new method, namely the Parallel Two-Level Hybrid (PTH) method, is developed to solve tridiagonal systems on parallel computers. PTH is designed based on Parallel Diagonal Dominant (PDD) algorithm. Like PDD, PTH is highly scalable. It provides accurate solutions when PDD may not be applicable and maintains a near PDD performance when the underlying machine ensemble size is large. By controllin...
متن کاملSome aspects of the cyclic reduction algorithm for block tridiagonal linear systems
The solution of a general block tridiagonal linear system by a cyclic odd-even reduction algorithm is considered. Under conditions of diagonal dominance, norms describing the off-diagonal blocks relative to the diagonal blocks decrease quadratically with each reduction. This allows early termination of the reduction when an approximate solution is desired. The algorithm is well-suited for paral...
متن کاملDirect and Iterative Methods for Block Tridiagonal Linear Systems
Block tridiagonal systems of linear equations occur frequently in scientific computations, often forming the core of more complicated problems. Numerical methods for solution of such systems are studied with emphasis on efficient methods for a vector computer. A convergence theory for direct methods under conditions of block diagonal dominance is developed, demonstrating stability, convergence ...
متن کاملOptimizing OPM for tridiagonal systems on vector uniprocessors
In the paper, we compare the optimum OPM to the Cyclic Reduction and the Tricyclic Reduction algorithms. In order to compare OPM with CR and TR we build a model for each of these two methods, too. The models built for the three algorithms show to be accurate (mean quadratic error lower than a 5% for the cases tested). For this reason the comparison of the three methods is made with their models...
متن کامل