A widely used approach to image segmentation is to define corresponding segmentation energies and to compute shapes that are minimizers of these energies. In this work, we introduce a flexible and efficient numerical framework for minimization of such