A parallel realisation of the Conjugate Gradient Method with Schur-Complement preconditioning, based on a domain decomposition approach, is described in detail. Special kinds of solvers for the resulting interiour and coupling systems are presented. A large range of numerical results is used to demonstrate the properties and behaviour of this solvers in practical situations