Strong Underrelaxation in Kaczmarz's Method for Inconsistent Systems
We present a parallel preconditioning method for the iterative solution of the time-harmonic elastic wave equation which makes use of higher-order spectral elements to reduce pollution error. In particular, the method leverages perfectly matched layer boundary conditions to efficiently approximate the Schur complement matrices of a block LDLT factorization. Both sequential and parallel versions of the algorithm are discussed and results for large-scale problems from exploration geophysics are presented....
Convection-diffusion problems posed on the unit square and with solutions displaying exponential layers are solved using a sparse grid Galerkin finite element method with Shishkin meshes. Writing for the maximum number of mesh intervals in each coordinate direction, our “combination” method simply adds or subtracts solutions that have been computed by the Galerkin FEM on , and meshes. It is shown that the combination FEM yields (up to a factor ) the same order of accuracy in the associated...
For the symmetric Pareto Eigenvalue Complementarity Problem (EiCP), by reformulating it as a constrained optimization problem on a differentiable Rayleigh quotient function, we present a class of descent methods and prove their convergence. The main features include: using nonlinear complementarity functions (NCP functions) and Rayleigh quotient gradient as the descent direction, and determining the step size with exact linear search. In addition, these algorithms are further extended to solve the...