Comparison of numerical methods in the contrast imaging problem in NMR

In this article, the contrast imaging problem in nuclear magnetic resonance is modeled as a Mayer problem in optimal control. A first synthesis of locally optimal solutions is given in the single-input case using geometric methods and the HamPath software, both based on Pontryagin´s maximum principle. We then compare these results using direct methods implemented with the Bocop toolbox and a moment-based approach, making a first step towards global optimality.