himmelbl.pro~orig

sqp程序包。用sqp算法实现非线性约束的优化求解

donlp2, v3, 05/29/98, copyright P. Spellucci
Thu Feb 24 16:53:27 2000
himmelblau6

     n=        45    nh=        16    ng=        45

  epsx= 1.000e-05 sigsm= 1.490e-08

startvalue
  1.0000000e-01   1.0000000e-01   1.0000000e-01   1.0000000e-01   1.0000000e-01 
  1.0000000e-01   1.0000000e-01   1.0000000e-01   1.0000000e-01   1.0000000e-01 
  1.0000000e-01   1.0000000e-01   1.0000000e-01   1.0000000e-01   1.0000000e-01 
  1.0000000e-01   1.0000000e-01   1.0000000e-01   1.0000000e-01   1.0000000e-01 
  1.0000000e-01   1.0000000e-01   1.0000000e-01   1.0000000e-01   1.0000000e-01 
  1.0000000e-01   1.0000000e-01   1.0000000e-01   1.0000000e-01   1.0000000e-01 
  1.0000000e-01   1.0000000e-01   1.0000000e-01   1.0000000e-01   1.0000000e-01 
  1.0000000e-01   1.0000000e-01   1.0000000e-01   1.0000000e-01   1.0000000e-01 
  1.0000000e-01   1.0000000e-01   1.0000000e-01   1.0000000e-01   1.0000000e-01 

  eps=  2.22e-16  tol= 1.98e-323 del0=  1.00e-05 delm=  3.67e-11 tau0=  1.00e+04
  tau=  1.00e-01   sd=  1.00e-01   sw=  3.67e-11  rho=  1.00e-06 rho1=  1.00e-10
 scfm=  1.00e+04  c1d=  1.00e-02 epdi=  1.00e-08
  nre=        45 anal=         1

 termination reason:
 KT-conditions satisfied, no further correction computed
 evaluations of f                         2285
 evaluations of grad f                     467
 evaluations of constraints              36560
 evaluations of grads of constraints        16
 final scaling of objective           4.860436e-02
 norm of grad(f)                      1.307726e+02
 lagrangian violation                 6.025803e-04
 feasibility violation                7.770066e-17
 dual feasibility violation           0.000000e+00
 optimizer runtime sec's              1.413000e+01


 optimal value of f =  -1.91038283278016e+03

 optimal solution  x =
  1.08366286278979e-06  2.55014516130919e-01  3.70480160720013e+00
  2.53736400510651e-01  6.52946016337137e-01  1.24636755497584e-03
  3.93083840470973e-04  1.00000000000000e-06  1.00000000000000e-06
  7.59876867033251e-02  8.75229861646824e-02  2.57231846609046e-02
  4.41431636273205e+01  2.17046709006796e-02  1.74754724187125e-06
  2.69066107911506e-05  1.55000000000000e-02  1.00000000000000e-06
  7.36832904773706e-05  3.83089593972778e-05  1.00000000000000e-06
  1.00000000000000e-06  4.06231329667491e-03  6.07013835317554e-04
  2.25668153390954e-02  2.60974738014518e+00  1.33477095750609e-03
  1.00000000000000e-06  1.49600740425417e-06  2.11275000000000e-02
  2.26850000000000e-03  1.00000000000000e-06  1.00000000000000e-06
  1.00000000000000e-06  1.00000000000000e-06  2.91250186801307e-05
  5.16083805723171e-03  3.89003692408816e-03  1.00000000000000e-06
  1.00000000000000e-06  8.00000000000000e-06  1.30265664455865e-02
  2.53358125076617e-03  2.97402507611457e-05  2.25974923885426e-06

  multipliers are relativ to scf=1
  nr.    constraint      normgrad (or 1)   multiplier
    1   0.0000000e+00    5.7445626e+00   -1.5173461e+01
    2   5.5511151e-17    3.3166248e+00   -1.0494741e+01
    3   0.0000000e+00    1.7320508e+00   -1.1648677e+01
    4   0.0000000e+00    3.8729833e+00   -2.0816420e+01
    5   0.0000000e+00    3.3166248e+00   -1.8593342e+01
    6   1.3877788e-17    1.4142136e+00   -6.4850370e+00
    7   0.0000000e+00    1.4142136e+00   -6.1423999e+00
    8   0.0000000e+00    1.4142136e+00   -4.7703003e+00
    9   0.0000000e+00    1.0000000e+00   -8.0747705e+00
   10   0.0000000e+00    1.0000000e+00   -4.8361983e+00
   11   0.0000000e+00    2.2360680e+00   -9.3275833e+00
   12  -6.9388939e-18    3.3166248e+00   -1.0065494e-01
   13  -9.3830075e-19    5.7445626e+00   -5.6498230e-01
   14   0.0000000e+00    5.7445626e+00   -2.3543459e+00
   15  -4.3368087e-19    4.2426407e+00   -1.7772919e-01
   16  -8.4703295e-22    4.2426407e+00   -7.3175340e-02
   17   8.3662863e-08    1.0000000e+00    0.0000000e+00
   18   2.5501352e-01    1.0000000e+00    0.0000000e+00
   19   3.7048006e+00    1.0000000e+00    0.0000000e+00
   20   2.5373540e-01    1.0000000e+00    0.0000000e+00
   21   6.5294502e-01    1.0000000e+00    0.0000000e+00
   22   1.2453676e-03    1.0000000e+00    0.0000000e+00
   23   3.9208384e-04    1.0000000e+00    0.0000000e+00
   24   0.0000000e+00    1.0000000e+00    3.2943638e+00
   25   0.0000000e+00    1.0000000e+00    8.6997644e-01
   26   7.5986687e-02    1.0000000e+00    0.0000000e+00
   27   8.7521986e-02    1.0000000e+00    0.0000000e+00
   28   2.5722185e-02    1.0000000e+00    0.0000000e+00
   29   4.4143163e+01    1.0000000e+00    0.0000000e+00
   30   2.1703671e-02    1.0000000e+00    0.0000000e+00
   31   7.4754724e-07    1.0000000e+00    0.0000000e+00
   32   2.5906611e-05    1.0000000e+00    0.0000000e+00
   33   1.5499000e-02    1.0000000e+00    0.0000000e+00
   34   0.0000000e+00    1.0000000e+00    1.0828932e+01
   35   7.2683290e-05    1.0000000e+00    0.0000000e+00
   36   3.7308959e-05    1.0000000e+00    0.0000000e+00
   37   0.0000000e+00    1.0000000e+00    6.0218906e+00
   38   0.0000000e+00    1.0000000e+00    3.7988131e+00
   39   4.0613133e-03    1.0000000e+00    0.0000000e+00
   40   6.0601384e-04    1.0000000e+00    0.0000000e+00
   41   2.2565815e-02    1.0000000e+00    0.0000000e+00
   42   2.6097464e+00    1.0000000e+00    0.0000000e+00
   43   1.3337710e-03    1.0000000e+00    0.0000000e+00
   44   0.0000000e+00    1.0000000e+00    2.2699737e+00
   45   4.9600740e-07    1.0000000e+00    0.0000000e+00
   46   2.1126500e-02    1.0000000e+00    0.0000000e+00
   47   2.2675000e-03    1.0000000e+00    0.0000000e+00
   48   0.0000000e+00    1.0000000e+00    3.7542759e+00
   49   0.0000000e+00    1.0000000e+00    1.6152081e+01
   50   0.0000000e+00    1.0000000e+00    2.9063250e+01
   51   0.0000000e+00    1.0000000e+00    3.8323859e+01
   52   2.8125019e-05    1.0000000e+00    0.0000000e+00
   53   5.1598381e-03    1.0000000e+00    0.0000000e+00
   54   3.8890369e-03    1.0000000e+00    0.0000000e+00
   55   0.0000000e+00    1.0000000e+00    4.0783706e+00
   56   0.0000000e+00    1.0000000e+00    5.1967232e-02
   57   7.0000000e-06    1.0000000e+00    0.0000000e+00
   58   1.3025566e-02    1.0000000e+00    0.0000000e+00
   59   2.5325813e-03    1.0000000e+00    0.0000000e+00
   60   2.8740251e-05    1.0000000e+00    0.0000000e+00
   61   1.2597492e-06    1.0000000e+00    0.0000000e+00

 evaluations of restrictions and their gradients
 (  2285,     1) (  2285,     1) (  2285,     1) (  2285,     1) (  2285,     1)
 (  2285,     1) (  2285,     1) (  2285,     1) (  2285,     1) (  2285,     1)
 (  2285,     1) (  2285,     1) (  2285,     1) (  2285,     1) (  2285,     1)
 (  2285,     1) (     0,     0) (     0,     0) (     0,     0) (     0,     0)
 (     0,     0) (     0,     0) (     0,     0) (     0,     0) (     0,     0)
 (     0,     0) (     0,     0) (     0,     0) (     0,     0) (     0,     0)
 (     0,     0) (     0,     0) (     0,     0) (     0,     0) (     0,     0)
 (     0,     0) (     0,     0) (     0,     0) (     0,     0) (     0,     0)
 (     0,     0) (     0,     0) (     0,     0) (     0,     0) (     0,     0)
 (     0,     0) (     0,     0) (     0,     0) (     0,     0) (     0,     0)
 (     0,     0) (     0,     0) (     0,     0) (     0,     0) (     0,     0)
 (     0,     0) (     0,     0) (     0,     0) (     0,     0) (     0,     0)
 (     0,     0)
last estimate of condition of active gradients  2.301e+02
last estimate of condition of approx. hessian   1.190e+09
iterative steps total             466
# of restarts                       1
# of full regular updates         465
# of updates                      465
# of full regularized SQP-steps    12