Parametric shape optimization with differentiable FEM simulation