Then an efficient penalty-based GRG (generalized reduced gradient) approach is used to solve the transformed problem, in which the adjoin equation method is applied to evaluate the gradient of the penalty term associated with the stability constraints.