Parallel tridiagonal matrix inversion with a hybrid multigrid–Thomas algorithm method

Author(s):  
J.T. Parker ◽  
P.A. Hill ◽  
D. Dickinson ◽  
B.D. Dudson
2008 ◽  
Vol 227 (6) ◽  
pp. 3174-3190 ◽  
Author(s):  
Dan Erik Petersen ◽  
Hans Henrik B. Sørensen ◽  
Per Christian Hansen ◽  
Stig Skelboe ◽  
Kurt Stokbro

2019 ◽  
Vol 68 (7) ◽  
pp. 6272-6285 ◽  
Author(s):  
Chuan Zhang ◽  
Xiao Liang ◽  
Zhizhen Wu ◽  
Feng Wang ◽  
Shunqing Zhang ◽  
...  

Author(s):  
Natalia M. Gavrilova ◽  
Yuri A. Plotonenko ◽  
Andrey A. STUPNIKOV

One of the most important ways of improving the speed of complex task solving is employing a multiprocessor computational system. This paper describes the experience of software development for research management and solving educational problems using parallel computing technologies. The authors describe approaches to computation parallelization using a multiprocessor system with shared memory within a task of finding a numerical root of a system of linear equations with a tridiagonal coefficient matrix that appears when solving a boundary problem for a partial differential equation of parabolic type, the heat equation. Additionally, the approaches to parallelization implementation of the tridiagonal matrix method for the heat equation in the two-dimensional case within a numerical root-finding algorithm using the alternating-direction implicit method for a multiprocessor system with shared memory are described. A finite-difference method of variable directions is used to find a numerical root of a heat equation in the two-dimensional case. Sequential and parallel algorithms (two-sided Thomas algorithm and multithread horizontal block Thomas algorithm) that fit an execution on computational systems with shared memory have been used to implement a tridiagonal matrix method. Two parallel computation organization technologies for computational systems with shared memory have been used for computation parallelization: one based on OpenMP technology and one using .NET framework facilities. The parallelization process and load balance serving have been performed by means of the environment in the first case as manual operation of threads parallelization process is allowed in the latter one. As an assessment of the described approach performance, the calculation time for sequential and parallel algorithms is given depending on the task’s size and the number of threads used. Comparison of the considered parallelization algorithms and implementation technologies is performed based on the analysis of the resulting acceleration. This paper shows that total computation time is several times smaller and calculation acceleration is several times bigger when using Thread instead of OpenMP. An application has been developed that allows obtaining a visual result of modelling of process of temperature propagation in the study area using parallel calculation technologies in real time.


1962 ◽  
Vol 5 (8) ◽  
pp. 438-439
Author(s):  
Richard George
Keyword(s):  

2020 ◽  
Vol 12 (11) ◽  
pp. 1747 ◽  
Author(s):  
Yin Zhang ◽  
Qiping Zhang ◽  
Yongchao Zhang ◽  
Jifang Pei ◽  
Yulin Huang ◽  
...  

Deconvolution methods can be used to improve the azimuth resolution in airborne radar imaging. Due to the sparsity of targets in airborne radar imaging, an L 1 regularization problem usually needs to be solved. Recently, the Split Bregman algorithm (SBA) has been widely used to solve L 1 regularization problems. However, due to the high computational complexity of matrix inversion, the efficiency of the traditional SBA is low, which seriously restricts its real-time performance in airborne radar imaging. To overcome this disadvantage, a fast split Bregman algorithm (FSBA) is proposed in this paper to achieve real-time imaging with an airborne radar. Firstly, under the regularization framework, the problem of azimuth resolution improvement can be converted into an L 1 regularization problem. Then, the L 1 regularization problem can be solved with the proposed FSBA. By utilizing the low displacement rank features of Toeplitz matrix, the proposed FSBA is able to realize fast matrix inversion by using a Gohberg–Semencul (GS) representation. Through simulated and real data processing experiments, we prove that the proposed FSBA significantly improves the resolution, compared with the Wiener filtering (WF), truncated singular value decomposition (TSVD), Tikhonov regularization (REGU), Richardson–Lucy (RL), iterative adaptive approach (IAA) algorithms. The computational advantage of FSBA increases with the increase of echo dimension. Its computational efficiency is 51 times and 77 times of the traditional SBA, respectively, for echoes with dimensions of 218 × 400 and 400 × 400 , optimizing both the image quality and computing time. In addition, for a specific hardware platform, the proposed FSBA can process echo of greater dimensions than traditional SBA. Furthermore, the proposed FSBA causes little performance degradation, when compared with the traditional SBA.


Sign in / Sign up

Export Citation Format

Share Document