No label defined (Q4937289)

The author considers stiff elliptic boundary value problems of the form \[ -\text{div}(\overline A(x)\nabla u)= f\quad\text{in }D,\quad u= 0\quad\text{on }\partial D,\tag{1} \] where \(D\) denotes a bounded domain in \(\mathbb{R}^2\). The \(2\times 2\) matrix \(\overline A\) can be represented as \(\overline A(x)= I+ \omega(x)A(x)\), with the identity matrix \(I\), a smooth, symmetric and positive definite matrix \(A\), and some parameter \(\omega\in [1,\infty]\) that is constant in certain isolated subdomains of \(D\). The author constructs an iterative solver of a special piecewise linear finite element approximation to (1) the convergence rate of which depends neither on the coefficient jumps nor on the discretization parameter. The numerical results presented confirm this theoretical result.