This paper describes a new parallel, scalable and robust finite element based solver for the
first-order Stokes momentum balance equations for ice flow. The solver, known as …

We present the full approximation scheme constraint decomposition (FASCD) multilevel
method for solving variational inequalities (VIs). FASCD is a joint extension of both the full …

The aim of this paper is to propose a multigrid method to obtain the numerical solution of the
one‐dimensional nonlinear sine‐Gordon equation. The finite difference equations at all …

Nonlinear multigrid methods such as the Full Approximation Scheme (FAS) and Newton-
multigrid (Newton-MG) are well established as fast solvers for nonlinear PDEs of elliptic and …

Abstract The Truncated Nonsmooth Newton Multigrid method is a robust and efficient
solution method for a wide range of block-separable convex minimization problems, typically …

Numerical glacier and ice-sheet models compute evolving ice geometry and velocity fields
using various stress-balance approximations and boundary conditions. At high spatial …

The isothermal, non-sliding shallow-ice approximation, combined with mass conservation, is
a fundamental model for ice-sheet and glacier flow. It determines the ice extent, geometry …

The aim of this paper is to introduce a fast and efficient new two-grid method to solve the d-
dimensional (d= 1, 2, 3) Poisson elliptic equations. The finite difference equations at all …

Glaciological modelling is a computationally challenging task due to its high cost and
complexity associated with large spatial-and long time-scale simulations. In this paper, we …

Accurate projections of the evolution of ice sheets in a changing climate require a fine
mesh/grid resolution in ice sheet models to correctly capture fundamental physical …