Abstract:In recent years, matrix‐free conjugate gradient (CG) solvers with graphics processing unit (GPU) acceleration have been effectively used to reduce the execution timings of finite element method (FEM)‐based engineering simulations. However, there is not much in the literature that discusses the application of matrix‐free CG solvers for elastoplasticity. The primary challenge is the presence of both elastic and plastic states, which introduce branching issues in the parallel computation of sparse matrix‐vector multiplication (SpMV). The current work proposes an efficient split kernel strategy for elastoplasticity that segregates elements in elastic and plastic zones and allows the application of standard GPU‐based matrix‐free SpMV strategies in the literature to elastoplastic simulation. The proposed strategy avoids branching issues, facilitates coalesced memory access, allows efficient usage of on‐chip memory, and uses minimum storage to realize large‐scale simulation on a GPU with limited memory. We also propose matrix‐free SpMV strategies for GPU implementation based on an element‐by‐element technique that makes efficient use of memory hierarchy available in a GPU. The computational efficiency of the proposed strategies is demonstrated over three large‐scale benchmark examples from elastoplasticity, and the performance results are compared with GPU optimized node‐based and degrees of freedom (DOF)‐based matrix‐free strategies. The results show that the splitting of computation is most suitable for low plasticity problems, where most of the matrix‐free strategies show significant speedups with respect to single kernel strategy for elastoplasticity. The proposed element‐by‐element strategy using node‐wise thread allocation is found to have the best performance, achieving up to 7.16× and 3.35× speedup over node‐based and DOF‐based strategy, respectively. For problems with higher amounts of plasticity, the proposed strategies perform slightly better and achieve modest speedups due to advantages in the initial stages of Newton iterations. In terms of wall‐clock timings, the proposed matrix‐free solver is found to be up to 2× faster than GPU optimized assembly‐based elastoplasticity solver.

Enhancing data locality of the conjugate gradient method for high-order matrix-free finite-element implementations

Parallel Conjugate Gradient Method Based on Spline Difference Method for the One-Dimensional Heat Equation

Performant low-order matrix-free finite element kernels on GPU architectures

High-efficiency Improved Symmetric Successive Over-Relaxation Preconditioned Conjugate Gradient Method for Solving Large-Scale Finite Element Linear Equations

Code Generation and Performance Engineering for Matrix-Free Finite Element Methods on Hybrid Tetrahedral Grids

Optimizing Multi-Grid Preconditioned Conjugate Gradient Method on Multi-Cores

Matrix-free multigrid block-preconditioners for higher order Discontinuous Galerkin discretisations

Development of GPU‐based matrix‐free strategies for large‐scale elastoplasticity analysis using conjugate gradient solver

A Hierarchical Grid Algorithm for Accelerating High-Performance Conjugate Gradient Benchmark on Sunway Many-Core Processor

Fast hardware-aware matrix-free algorithm for higher-order finite-element discretized matrix multivector products on distributed systems

Multilevel Interior Penalty Methods on GPUs

Fast hardware-aware matrix-free algorithms for higher-order finite-element discretized matrix multivector products on distributed systems

High-Performance Implementation of Matrix-Free Runge-Kutta Discontinuous Galerkin Method for Euler Equations

Application of Diagonally Perturbed Incomplete Factorization Preconditioned Conjugate Gradient Algorithms for Edge Finite-Element Analysis of Helmholtz Equations

A Cooperative Conjugate Gradient Method for Linear Systems Permitting Multithread Implementation of Low Complexity

Manycore Parallel Computing for a Hybridizable Discontinuous Galerkin Nested Multigrid Method

Fundamental Data Structures for Matrix-Free Finite Elements on Hybrid Tetrahedral Grids

Deep Learning-Enhanced Preconditioning for Efficient Conjugate Gradient Solvers in Large-Scale PDE Systems

Efficient parallelization strategy for real-time FE simulations

Neural operators meet conjugate gradients: The FCG-NO method for efficient PDE solving

Exploiting mesh structure to improve multigrid performance for saddle-point problems