A $\texttt{SPIRED}$ code for the reconstruction of spin distribution

In Nuclear Magnetic Resonance (NMR), it is of crucial importance to have an accurate knowledge of the sample probability distribution corresponding to inhomogeneities of the magnetic fields. An accurate identification of the sample distribution requires a set of experimental data that is sufficiently rich to extract all fundamental information. These data depend strongly on the control fields (and their number) used experimentally. In this work, we present and analyze a greedy reconstruction algorithm, and provide the corresponding $\texttt{SPIRED}$ code, for the computation of a set of control functions allowing the generation of data that are appropriate for the accurate reconstruction of a sample distribution. In particular, the focus is on NMR and the Bloch system with inhomogeneities in the magnetic fields in all spatial directions. Numerical examples illustrate this general study.