The state-of-the-art methods for solving the DNN-parametrized version of eq:pde use penalty methods, Lagrangian methods, and augmented Lagrangian methods [1]}, [2]}, [3]}, [4]}, which often involve delicate tuning of multiple hyperparameters and could lead to infeasible solutions. By contrast, PyGRANSO stops iterations by rigorous check of constraint violation and stationarity.
