A robust alternating direction method of multipliers numerical scheme for solving geometric inverse problems in a shape optimization setting