Source author record

Mohammad Shafaet Islam

Mohammad Shafaet Islam appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Numerical Analysis Mathematical Software

Catalog footprint

What is connected

2works

3topics

1close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

A data driven heuristic for rapid convergence of Scheduled Relaxation Jacobi schemes

The Scheduled Relaxation Jacobi (SRJ) method is a viable candidate as a high performance linear solver for elliptic partial differential equations (PDEs). The method greatly improves the convergence of the standard Jacobi iteration by applying a sequence of $M$ well-chosen overrelaxation and underrelaxation factors in each cycle of the algorithm to effectively attenuate the solution error. In previous work, optimal SRJ schemes (sets of relaxation factors) have been derived to accelerate convergence for specific discretizations of elliptic PDEs. In this work, we develop a family of SRJ schemes which can be applied to solve elliptic PDEs regardless of the specific discretization employed. To achieve favorable convergence, we train an algorithm to select which scheme in this family to apply at each cycle of the linear solve process, based on convergence data collected from applying these schemes to the one-dimensional Poisson equation. The automatic selection heuristic that is developed based on this limited data is found to provide good convergence for a wide range of problems.

preprint2020arXiv

Hierarchical Jacobi Iteration for Structured Matrices on GPUs using Shared Memory

High fidelity scientific simulations modeling physical phenomena typically require solving large linear systems of equations which result from discretization of a partial differential equation (PDE) by some numerical method. This step often takes a vast amount of computational time to complete, and therefore presents a bottleneck in simulation work. Solving these linear systems efficiently requires the use of massively parallel hardware with high computational throughput, as well as the development of algorithms which respect the memory hierarchy of these hardware architectures to achieve high memory bandwidth. In this paper, we present an algorithm to accelerate Jacobi iteration for solving structured problems on graphics processing units (GPUs) using a hierarchical approach in which multiple iterations are performed within on-chip shared memory every cycle. A domain decomposition style procedure is adopted in which the problem domain is partitioned into subdomains whose data is copied to the shared memory of each GPU block. Jacobi iterations are performed internally within each block's shared memory, avoiding the need to perform expensive global memory accesses every step. We test our algorithm on the linear systems arising from discretization of Poisson's equation in 1D and 2D, and observe speedup in convergence using our shared memory approach compared to a traditional Jacobi implementation which only uses global memory on the GPU. We observe a x8 speedup in convergence in the 1D problem and a nearly x6 speedup in the 2D case from the use of shared memory compared to a conventional GPU approach.