We describe a GPUand multicore-oriented implementation technique for a key component of finite element based simulation toolkits for partial differential equations on unstructured grids: Geometric Multigrid solvers. We use efficient sparse matrix-vector multiplications throughout the solver pipeline: within the coarse-grid solver, smoothers and even grid transfers. Our implementation can handle...

Geologic survey procedures accumulate large volumes of structured and unstructured data. Fully exploiting the knowledge and information that are included in geological big data and improving the accessibility of large volumes of data are important endeavors. In this paper, which is based on the architecture of the geological survey information cloud-computing platform (GSICCP) and big-data-rela...

We propose a multigrid method for solving large-scale sparse linear systems arising from discretizations of partial differential equations, such as those from finite element and generalized finite difference (GFD) methods. Our proposed method has the following two characteristics. First, we introduce a hybrid geometric+algebraic multigrid method, or HyGA, to leverage the rigor, accuracy and eff...

In this paper we describe general software utilities for performing unstructured sparse matrix-vector multiplications on distributed-memory message-passing computers. The matrix-vector multiply comprises an important kernel in the solution of large sparse linear systems by iterative methods. Our focus is to present the data structures and communication parameters necessary for these utilities f...

A distance field is a representation of the closest distance from a point to a given surface. Distance fields are widely used in applications ranging from computer vision, physics and computer graphics and have been the subject of research of many authors in the last decade. Most of the methods for computing distance fields are devoted to Cartesian grids while little attention has been paid to ...

A new Navier-Stokes algorithm for use on unstructured triangular meshes is presented. Spatial discretization of the governing equations is achieved using a finite-element Galerkin approximation, which can be shown to be equivalent to a finite-volume approximation for regular equilateral triangular meshes. Integration to steady state is performed using a multistage time-stepping scheme, and conv...

This paper deals with a high-order accurate discontinuous finite element method for the numerical solution of the Euler equations. The method combines two key ideas which are at the basis of the finite volume and of the finite element method, the physics of wave propagation being accounted for by means of Riemann problems and accuracy being obtained by means of high-order polynomial approximati...

The simulation of fusion plasmas in realistic magnetic configurations and tokamak geometries still requires the development advanced numerical algorithms owing to complexity problem. In this context, we propose a Hybrid Discontinuous Galerkin (HDG) method solve 2D transport fluid equations wall geometries. This high-order solver can handle equilibrium free structured unstructured meshes allowin...