[Bonmin] understanding bonmin output
Giuseppe Aprea
giuseppe.aprea at gmail.com
Fri Jul 3 15:02:08 EDT 2009
Hi, I am solving a MINLP with bonmin. I am using B-Hyb
and the example file which is given with the software is my template.
The program ends catching TNLPSolver::UnsolvedError
and below you can see the output from the last two steps. It finally speaks
about "small primal infeasibility but I know that there is at least one
"trivial solution"; I mean there is an array which fulfills all
constraints. Can
anyone, plese, give me a suggestion to better understand this output?
By the way, I put:
print_solution yes
in options file but I can't see any solution array in the output, How can I
retrieve it?
Shall I use:
bb.bestSolution();
to retrive the pointer to the first element?
cheers,
giuseppe
**************************************************
*** Update HessianMatrix for Iteration 1015:
**************************************************
**************************************************
*** Summary of Iteration 1015 for original NLP:
**************************************************
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
1015r-7.4998033e+01 1.40e-06 1.14e-13 -9.6 4.78e-07 - 1.00e+00
1.00e+00s 19 FL
**************************************************
*** Beginning Iteration 1015 from the following point:
**************************************************
Primal infeasibility for restoration phase problem = 1.3960840661297880e-06
Dual infeasibility for restoration phase problem = 1.1368683772161603e-13
||curr_x||_inf = 9.9193648021086943e+08
||curr_s||_inf = 9.7255641900768816e+08
||curr_y_c||_inf = 2.8204645278437250e+00
||curr_y_d||_inf = 4.9928828305278242e-05
||curr_z_L||_inf = 1.0028204645278437e+03
||curr_z_U||_inf = 2.8204645278437268e+00
||curr_v_L||_inf = 1.0079063885684439e-17
||curr_v_U||_inf = 4.9928828305278520e-05
***Current NLP Values for Iteration (Restoration phase problem) 1015:
(scaled) (unscaled)
Objective...............: 1.3021940861321653e-10 1.3021940861321653e-10
Dual infeasibility......: 1.1368683772161603e-13 1.1368683772161603e-13
Constraint violation....: 1.3960840661297880e-06 1.3960840661297880e-06
Complementarity.........: 2.7660570686466480e-10 2.7660570686466480e-10
Overall NLP error.......: 1.3960840661297880e-06 1.3960840661297880e-06
**************************************************
*** Update Barrier Parameter for Iteration 1015:
**************************************************
Remaining in fixed mu mode.
Barrier Parameter: 2.766057e-10
**************************************************
*** Solving the Primal Dual System for Iteration 1015:
**************************************************
Number of doubles for MUMPS to hold factorization (INFO(9)) = 167090
Number of integers for MUMPS to hold factorization (INFO(10)) = 137037
MUMPS returned INFO(1) = -10 matrix is singular.
Factorization failed with retval = 1
Number of doubles for MUMPS to hold factorization (INFO(9)) = 166449
Number of integers for MUMPS to hold factorization (INFO(10)) = 136955
Factorization successful.
Number of trial factorizations performed: 2
Perturbation parameters: delta_x=0.000000e+00 delta_s=0.000000e+00
delta_c=4.078170e-11 delta_d=4.078170e-11
residual_ratio = 3.782836e-12
Factorization successful.
residual_ratio = 2.244350e-13
**************************************************
*** Finding Acceptable Trial Point for Iteration 1015:
**************************************************
--> Starting line search in iteration 1015 <--
The current filter has 6 entries.
Trying soft restoration phase step with step length 1.000000e+00
Checking acceptability for trial step size alpha_primal_test= 0.000000e+00:
New values of barrier function = 1.2585792238523389e-05
(reference 1.2585791362940986e-05):
New values of constraint violation = 1.9299189705205889e-05
(reference 2.1041560280738349e-05):
Checking sufficient reduction...
Succeeded...
Checking filter acceptability...
Failed...
Primal-dual error at current point: 6.2754482887496734e-09
Primal-dual error at trial point : 5.7558015386434678e-09
Trial step accepted.
orig_curr_theta = 1.89e-05, orig_trial_theta = 1.93e-05
Point does not provide sufficient reduction w.r.t the original theta
(orig_theta_max=1.887143e-06).
**************************************************
*** Update HessianMatrix for Iteration 1016:
**************************************************
**************************************************
*** Summary of Iteration 1016 for original NLP:
**************************************************
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls
1016r-7.4998033e+01 1.00e-06 1.14e-13 -9.6 5.75e-07 - 1.00e+00
1.00e+00s 0 FLA
**************************************************
*** Beginning Iteration 1016 from the following point:
**************************************************
Primal infeasibility for restoration phase problem = 9.9954729648743523e-07
Dual infeasibility for restoration phase problem = 1.1368683772161603e-13
||curr_x||_inf = 9.9193648021086931e+08
||curr_s||_inf = 9.7255641900768793e+08
||curr_y_c||_inf = 2.8204648379838653e+00
||curr_y_d||_inf = 4.9928828305278242e-05
||curr_z_L||_inf = 1.0028204648379839e+03
||curr_z_U||_inf = 2.8204648379838670e+00
||curr_v_L||_inf = 1.0079063885684386e-17
||curr_v_U||_inf = 4.9928828305278520e-05
***Current NLP Values for Iteration (Restoration phase problem) 1016:
(scaled) (unscaled)
Objective...............: 1.3021940856162432e-10 1.3021940856162432e-10
Dual infeasibility......: 1.1368683772161603e-13 1.1368683772161603e-13
Constraint violation....: 9.9954729648743523e-07 9.9954729648743523e-07
Complementarity.........: 2.7660570686410599e-10 2.7660570686410599e-10
Overall NLP error.......: 9.9954729648743523e-07 9.9954729648743523e-07
**************************************************
*** Update Barrier Parameter for Iteration 1016:
**************************************************
Remaining in fixed mu mode.
Barrier Parameter: 2.766057e-10
**************************************************
*** Solving the Primal Dual System for Iteration 1016:
**************************************************
Number of doubles for MUMPS to hold factorization (INFO(9)) = 167090
Number of integers for MUMPS to hold factorization (INFO(10)) = 137037
MUMPS returned INFO(1) = -10 matrix is singular.
Factorization failed with retval = 1
Number of doubles for MUMPS to hold factorization (INFO(9)) = 166449
Number of integers for MUMPS to hold factorization (INFO(10)) = 136955
Factorization successful.
Number of trial factorizations performed: 2
Perturbation parameters: delta_x=0.000000e+00 delta_s=0.000000e+00
delta_c=4.078170e-11 delta_d=4.078170e-11
residual_ratio = 2.425931e-11
Factorization successful.
residual_ratio = 7.912670e-13
**************************************************
*** Finding Acceptable Trial Point for Iteration 1016:
**************************************************
--> Starting line search in iteration 1016 <--
Storing current iterate as backup acceptable point.
The current filter has 6 entries.
Trying soft restoration phase step with step length 1.000000e+00
Checking acceptability for trial step size alpha_primal_test= 0.000000e+00:
New values of barrier function = 1.2585792491532923e-05
(reference 1.2585792238523389e-05):
New values of constraint violation = 1.9522375838284858e-05
(reference 1.9299189705205889e-05):
Checking sufficient reduction...
Failed...
Primal-dual error at current point: 5.7558015386434678e-09
Primal-dual error at trial point : 5.8223639139496865e-09
Trial step rejected.
Restoration phase converged to a point with small primal infeasibility.
Number of Iterations....: 1016
(scaled) (unscaled)
Objective...............: -7.4998032932067574e+01 -7.4998032932067574e+01
Dual infeasibility......: 9.9999999956378893e-01 9.9999999956378893e-01
Constraint violation....: 9.9954729648743523e-07 9.9954729648743523e-07
Complementarity.........: 2.7660570686410599e-10 2.7660570686410599e-10
Overall NLP error.......: 9.9999999956378893e-01 9.9999999956378893e-01
Number of objective function evaluations = 3186
Number of objective gradient evaluations = 860
Number of equality constraint evaluations = 3190
Number of inequality constraint evaluations = 3190
Number of equality constraint Jacobian evaluations = 1024
Number of inequality constraint Jacobian evaluations = 1024
Number of Lagrangian Hessian evaluations = 1017
Total CPU secs in IPOPT (w/o function evaluations) = 1656.757
Total CPU secs in NLP function evaluations = 264.214
EXIT: Restoration Failed!
Problem status: 3
Objective value: 1.79769e+308
More information about the Bonmin
mailing list