<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman, new york, times, serif;font-size:12pt"><div>Hi All:<br><br>Thanks to you (especially Stefan) I have made great progress. I still have a couple of issues to resolve. In my maximization problem I have 102 choice variables. I provide the derivative of the objective and the bounds and the hessian too. I seem to have a problem only with the jacobian. I am providing a part of the output. Please note that the only problem is that the value of my derivative and that computed by the derivatives checker is off by one -e.g., my jac_g(3,0) is the checkers jac_g(2,1)., my jag_g(2,1) is the checkers jac_g(3,1) etc. There seem to be no other problems I am using C-style indexing. Any insights would be welcome. Thank you.<br><br>Krish<br><br><br><br>List of user-set
options:<br><br> Name Value used<br> derivative_test = second-order yes<br> output_file = krish.txt
yes<br> print_user_options = yes yes<br><br>******************************************************************************<br>This program contains Ipopt, a library for large-scale nonlinear optimization.<br> Ipopt is released as open source code under the Common Public License (CPL).<br> For more information visit http://projects.coin-or.org/Ipopt<br>******************************************************************************<br><br>NOTE: You are using Ipopt by default with the MUMPS linear solver.<br> Other linear solvers might be more efficient (see Ipopt documentation).<br><br><br><br>Starting derivative checker.<br><br>* jac_g
[ 2, 0] = 0.0000000000000000e+000 v ~ 1.0000000000000000e+000 [1.000e+000]<br>* jac_g [ 3, 0] = 1.0000000000000000e+000 v ~ 0.0000000000000000e+000 [1.000e+000]<br>* jac_g [ 2, 1] = -1.2653525307998922e-002 v ~ 1.0000000000000000e+000 [1.013e+000]<br>* jac_g [ 3, 1] = 1.0000000000000000e+000 v ~ -1.2653525189448844e-002 [1.013e+000]<br>* jac_g [ 2, 2] = -1.3398148523156232e-002 v ~ 1.0000000000000000e+000 [1.013e+000]<br>* jac_g [ 3, 2] = 1.0000000000000000e+000 v ~ -1.3398148302617985e-002 [1.013e+000]<br>* jac_g [ 2, 3] = -1.4125931249430622e-002 v ~ 1.0000000000000000e+000 [1.014e+000]<br>* jac_g [
3, 3] = 1.0000000000000000e+000 v ~ -1.4125930980490864e-002 [1.014e+000]<br>* jac_g [ 2, 4] = -1.2426902681357370e-002 v ~ 1.0000000000000000e+000 [1.012e+000]<br>----------------deleted as repetitive--------------------------------------------------------------<br>* jac_g [ 3, 100] = 0.0000000000000000e+000 ~ 1.1953086350574571e-002 [1.195e-002]<br>* jac_g [ 2, 101] = 1.1310855712182842e-002 v ~ 0.0000000000000000e+000 [1.131e-002]<br>* jac_g [ 3, 101] = 0.0000000000000000e+000 ~ 1.1310855984464352e-002 [1.131e-002]<br><br>Derivative checker detected 202 error(s).<br><br>Too few degrees of freedom (n_x = 0, n_c = 1).<br> Trying fixed_variable_treatment = RELAX_BOUNDS<br><br>Number of nonzeros in equality constraint
Jacobian...: 102<br>Number of nonzeros in inequality constraint Jacobian.: 2499<br>Number of nonzeros in Lagrangian Hessian.............: 5253<br><br>Total number of variables............................: 102<br> variables with only lower bounds: 0<br> variables with lower and upper bounds: 102<br> variables with only upper bounds: 0<br>Total number of equality
constraints.................: 1<br>Total number of inequality constraints...............: 25<br> inequality constraints with only lower bounds: 0<br> inequality constraints with lower and upper bounds: 25<br> inequality constraints with only upper bounds: 0<br><br> <br></div></div><br>
</body></html>