  We propose new domain decomposition methods for systems of partial
differential equations in two and three dimensions. The algorithms are derived
with the help of the Smith factorization of the operator. This could also be
validated by numerical experiments.
