We present a Sobolev gradient type preconditioning for iterative methods used in solving second order semilinear elliptic systems; the n-tuple of independent Laplacians acts as a preconditioning operator in Sobolev spaces. The theoretical iteration is done at the continuous level, providing a linearization approach that reduces the original problem to a system of linear Poisson equations. The method obtained preserves linear convergence when a polynomial growth of the lower order reaction type terms is involved. For the proof of linear convergence for systems with mixed boundary conditions, we use suitable energy spaces. We use Sobolev embedding estimates in the construction of the exact algorithm. The numerical implementation has focus on a direct and elementary realization, for which a detailed discussion and some examples are given.
|Number of pages||26|
|Journal||Electronic Journal of Differential Equations|
|Publication status||Published - máj. 21 2004|
ASJC Scopus subject areas