Source imaging maps back boundary measurements to underlying generators within the domain; e.g., retrieving the parameters of the generating dipoles from electrical potential measurements on the scalp such as in electroencephalography EEG). Fitting such a parametric source model is non-linear in the positions of the sources and renewed interest in mathematical imaging has led to several promising approaches.