The Cumulative Reconstructor is an accurate, extremely fast reconstruction algorithm for Shack-Hartmann wavefront sensor data. But it has shown an unacceptable high noise propagation for large apertures. Therefore, in this paper we describe a domain decomposition approach to deal with this drawback. We show that this adaptation of the algorithm gives the same reconstruction quality as the original algorithm and leads to a significant improvement with respect to noise propagation. The method is combined with an integral control and compared to the classical matrix vector multiplication algorithm on an end-to-end simulation of a single conjugate adaptive optics system. The reconstruction time is 20n (number of subapertures), and the method is parallelizable.