<div dir="ltr"><div><div>All,<br><br></div>I have some questions regarding the KKT system that is being factorized by whichever linear solver is interfaced with IPOPT (Mumps, MA27, MA57). The particular linear solver is not in question, it is the KKT matrix that IPOPT is sending to the linear solver. Specifically, it says in Zavala's book, "Computational Strategies for the Optimal Operation of Large-Scale Chemical Processes" that IPOPT factorizes the augmented matrix, where the augmented matrix includes the barrier Hessian. I want to know where the following quantities are computed in IPOPT:<br>
<br></div><div>(1) The actual Hessian, H_i<br></div><div>(2) The barrier Hessian, W_i, where W_i = H_i + Sigma_i<br></div><div>(3) Sigma_i, where Sigma_i = inverse(X_i)*V_i where X_i is a diagonal matrix of the current estimate of the decision vector and V_i is a diagonal matrix of the Lagrange multipliers of the simple bound constraints x_min <= x <= x_max<br>
<br></div><div>The reason I ask this question is that in some problems we are solving IPOPT appears to be sending to MA57 augmented matrices whose elements have extremely wide disparities in magnitude. For example, in some problems which IPOPT seems to solve perfectly fine we are seeing many orders of magnitude difference in elements of the augmented matrix and we cannot understand why the elements are so widely separated in magnitude. <br>
<br></div><div>In addition to the above requests, I am grateful if somebody can point me to the exact location in IPOPT where the augmented matrix is constructed and sent to the linear solver. <br></div><div><br></div><div>
Thank you for any help you can provide.<br><br></div><div>Regards,<br><br>Anil Rao<br><br></div><div><br><br></div><div><div><div>-- <br>Anil V. Rao<br>
Associate Professor<br>
Department of Mechanical and Aerospace Engineering<br>
University of Florida<br>
Gainesville, FL 32611-6250<br>Tel: (352) 672-1529<br>E-mail: <a href="mailto:anilvrao@gmail.com" target="_blank">anilvrao@gmail.com</a><br>
</div></div></div></div>