Currently, emg3d only returns one misfit and one gradient (per direction). TODO: Implement the possibility to return the misfit and the gradient for each source-receiver pair.