  A parallelizable iterative procedure based on domain decomposition is
presented and analyzed for weak Galerkin finite element methods for second
order elliptic equations. The convergence analysis is established for the
decomposition of the domain into individual elements associated to the weak
Galerkin methods or into larger subdomains. A series of numerical tests are
illustrated to verify the theory developed in this paper.
