Sparse Matrix-Vector Multiplication for Finite Element Method Matrices on FPGAs

Author(s):  
Yousef El-Kurdi ◽  
Warren Gross ◽  
Dennis Giannacopoulos
2008 ◽  
Vol 178 (8) ◽  
pp. 558-570 ◽  
Author(s):  
Yousef Elkurdi ◽  
David Fernández ◽  
Evgueni Souleimanov ◽  
Dennis Giannacopoulos ◽  
Warren J. Gross

Author(s):  
Vikalp Mishra ◽  
Krishnan Suresh

A serious computational bottle-neck in finite element analysis today is the solution of the underlying system of equations. To alleviate this problem, researchers have proposed the use of graphics programmable units (GPU) for fast iterative solution of such equations. Indeed, researchers have shown that a GPU-implementation of a double-precision sparse-matrix-vector multiplication (that underlies all iterative methods) is approximately an order of magnitude faster than that of an optimized CPU implementation. Unfortunately, fast matrix-vector multiplication alone is insufficient… a good preconditioner is necessary for rapid convergence. Furthermore, most modern preconditioners, such as incomplete Cholesky, are expensive to compute, and cannot be easily ported to the GPU. In this paper, we propose a special class of preconditioners for the analysis of thin structures, such as beams and plates. The proposed preconditioners are developed by combining the multi-grid method, with recently developed dual-representation method for thin structures. It is shown, that these preconditioners are computationally inexpensive, perform better than standard pre-conditioners, and can be easily ported to the GPU.


2016 ◽  
Vol 33 (8) ◽  
pp. 2339-2355 ◽  
Author(s):  
Yunfei Liu ◽  
Jun Lv ◽  
Xiaowei Gao

Purpose The purpose of this paper is to introduce a new method called simultaneous elimination and back-substitution method (SEBSM) to solve a system of linear equations as a new finite element method (FEM) solver. Design/methodology/approach In this paper, a new technique assembling the global stiffness matrix will be proposed and meanwhile the direct method SEBSM will be applied to solve the equations formed in FEM. Findings The SEBSM solver for FEM with the present assembling technique has distinct advantages in both computational time and memory space occupation over the conventional methods, such as the Gauss elimination and LU decomposition methods. Originality/value The developed solver requires less memory space no matter the coefficient matrix is a typical sparse matrix or not, and it is applicable to both symmetric and unsymmetrical linear systems of equations. The processes of assembling matrix and dealing with constraints are straightforward, so it is convenient for coding. Compared to the previous solvers, the proposed solver has favorable universality and good performances.


Sign in / Sign up

Export Citation Format

Share Document