In this code we implement an extension of the augmented Lagrangian (AL) method proposed in: Vermeylen, C., Van Barel, M.