where G i is the column-reduced version of G and Θ is a p i ≤ p by N matrix depending on the G i and priors computed from the previous source estimate D^i−1. A similar approach was used in [31] where the source region was contracted several times but at each iteration, LORETA was used to estimate the source tomography.
