Parallel implementation of Thomas algorithm for the 2D heat equation

Journal of Mathematics Mechanics and Computer Science, 2019

Full text Open access: Gold Metadata (CSL-JSON)
DOI 10.26577/jmmcs-2019-3-24
Data confidence 70 medium
How this score is calculated

It is built from checkable properties of the record itself: whether the article has a DOI, how many independent sources confirmed its data, whether the required fields are filled in, and how recent the latest confirmation is. Disagreement between sources and a missing identifier lower it. A value of 80 or above means the record was published without further checking; 60 to 79 means it is published at reduced confidence, with a manual check advisable. A record verified by an editor by hand is never pushed below 90 by a recalculation. This scores our data about the article, not the article itself.

Abstract

In this paper was considered a parallel implementation of the Thomas algorithm for the 2D heat equation.MPI was chosen as the technology for parallelization.The numerical solution of the twodimensional heat conduction problem was solved using the two step iteration process of alternating direction implicit method (ADI).The Thomas algorithm is simple to implement a sequential program, but difficult to parallelize due to dependent data transfers.When applying this method to solve the 2D heat equation, there is a need to implement Thomas algorithm along each direction x-axis and y-axis.It was implemented the parallelization of this problem using the Yanenko method using 1D and 2D data decomposition.In particular, in 2D data decomposition, the Yanenko method was used along each x-axis and y-axis direction.In the article the speedup and efficiency of parallel programs using 1D and 2D data decomposition were shown in the form of tables and graphs.The presented algorithm was tested on a cluster of the computing center of Novosibirsk State University for a different number of points in the computational domain (from 512x512 to 4096x4096).The obtained test results are presented and analyzed, on the basis of which the features of the used decompositions are described.