Next: 1.5 Correspondence: operator-module
Up: 1 Solution of linear systems
Prev: 1.3 Solution by an iterative method
Contents
The numerical solution of large problems by domain decomposition techniques is very well adapted to current generation parallel computers. However, the efficiency of these techniques depend strongly on the algorithm chosen and its implementation.
The approach proposed in the MODULEF code divides the computational domain into non-structured sub-domains of arbitrary form, and reduces the initial problem to an interface problem. The corresponding operator (the Steklov-Poincaré operator at the continuous level, the Schur complement matrix at the discrete level) is inverted by an preconditioned conjugate gradient algorithm . This algorithm requires, at each step, the solution of a Dirichlet and a Neumann problem on each sub-domain.
For a detailed description of the algorithm implemented in MODULEF, consult [2].