Interferometric synthetic-aperture radar (InSAR) interferograms contain valuable information about the fault systems hidden beneath the surface of the Earth. In a new approach, we aim to fit InSAR ground deformation data using a volumetric distribution of multiple seismic point sources whose parameters are found by a genetic algorithm. The resulting source distribution could provide another useful tool in solving the difficult problem of accurately mapping earthquake faults. To test the algorithm, we first apply it to synthetic data, followed by applications to an ALOS-2 InSAR interferogram. We report first results and discuss advantages and disadvantages of this approach.