A constrained-optimization based half-quadratic algorithm for robustly fitting sets of linearly parametrized curves