The Hessian matrix of every function in this model is constant, so it will be calculated once in the entire optimal process based on interior point method, which speeds up each iteration.