undefexa.pro~orig

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

donlp2, v3, 05/29/98, copyright P. Spellucci
Thu Feb 24 16:56:06 2000
undefexa



================================================================================
             1-th iteration step
   scf=  1.0000e+00 psist=  0.0000e+00   psi=  0.0000e+00  upsi=  0.0000e+00
  fxst= -2.0000e+00    fx= -2.0000e+00
  x=
   1.0000e+00   1.0000e+00
 valid permutation of x

  1   2 

  multipliers: first estimate
  u =
    2  -1.2000e+00


  del=  1.00000e-01  b2n0=  8.94427e-01   b2n=  8.94427e-01   gfn=  2.82843e+00


 values of restrictions
 (   2    0.0000e+00    2.2361e+00)  
 

   diag[r]=
   1.0000e+00

  multipliers: second estimate
  u =
    2  -1.2000e+00
  condition number of r     1.000000000000000e+00
  condition number of a     1.000000000000000e+00


  current scaling,  scf =   1.0000e+00 clow =            1 eta =   0.0000e+00


  scalres=
   1.0000e+00   2.6000e+00   1.0000e+00   1.0000e+00


  phase=  0  scf0=  1.0000e+00
  d =
  -1.0400e+00  -8.0000e-02


 start unimin

    phi= -2.0000e+00   dphi= -2.2400e+00    psi=  0.0000e+00 tau0/2=  5.0000e-01
     fx= -2.0000e+00  dscal=  1.0000e+00    scf=  1.0000e+00   upsi=  0.0000e+00
    sig=  1.2500e-01     fx= -2.4093e+00    psi=  0.0000e+00   upsi=  0.0000e+00


 end unimin

 sig=  1.2500e-01  num. f-evaluations 3
 list of inactive hit constraints
   1     1  
BFGS-update as in Pantoja and Mayne
  tk  =   5.000000e-01
 xsik =   3.790695e+01



================================================================================
             2-th iteration step
   scf=  1.0000e+00 psist=  0.0000e+00   psi=  0.0000e+00  upsi=  0.0000e+00
  fxst= -2.0000e+00    fx= -2.4093e+00
  x=
   8.7000e-01   9.9000e-01
 valid permutation of x

  1   2 


  del=  6.69016e-02  b2n0=  5.62347e+01   b2n=  5.61338e+00   gfn=  5.61338e+00
  condition number of r     -1.000000000000000e+00
  condition number of a     1.300102158716682e+01


  current scaling,  scf =   1.0000e+00 clow =            1 eta =   1.1600e-02


  scalres=
   1.0000e+00   2.6000e+00   1.0000e+00   1.0000e+00


  phase=  0  scf0=  1.0000e+00
  d =
  -9.1516e+00  -1.4880e+00


 start unimin

    phi= -2.4093e+00   dphi= -4.4102e+01    psi=  0.0000e+00 tau0/2=  5.0000e-01
     fx= -2.4093e+00  dscal=  1.3946e-02    scf=  1.0000e+00   upsi=  0.0000e+00
    sig=  8.3163e-03     fx= -2.8951e+00    psi=  0.0000e+00   upsi=  0.0000e+00


 end unimin

 sig=  8.3163e-03  num. f-evaluations 4
 list of inactive hit constraints
   1     1     1  
BFGS-update as in Pantoja and Mayne
  tk  =   5.000000e-01
 xsik =   9.366103e+01



================================================================================
             3-th iteration step
   scf=  1.0000e+00 psist=  0.0000e+00   psi=  0.0000e+00  upsi=  0.0000e+00
  fxst= -2.0000e+00    fx= -2.8951e+00
  x=
   7.9389e-01   9.7762e-01
 valid permutation of x

  1   2 

  multipliers: first estimate
  u =
    1   7.5895e+00    2   9.1663e-14


  del=  7.02021e-01  b2n0=  1.85109e+01   b2n=  9.05340e-14   gfn=  9.55798e+00


 values of restrictions
 (   1    2.6129e-02    1.2594e+00)  (   2    2.5036e-01    2.1961e+00)  
 

   diag[r]=
  -1.0000e+00   1.2997e-02

  multipliers: second estimate
  u =
    1   7.5895e+00    2   9.1663e-14
  condition number of r     7.694124972521753e+01
  condition number of a     4.662701264246899e+00


  current scaling,  scf =   1.0000e+00 clow =            1 eta =   1.1600e-02


  scalres=
   1.5379e+01   2.6000e+00   1.0000e+00   1.0000e+00


  phase=  0  scf0=  1.0000e+00
  d =
  -2.2355e-01   2.4238e-01


 start unimin

    phi= -2.8951e+00   dphi= -1.9830e-01    psi=  0.0000e+00 tau0/2=  5.0000e-01
     fx= -2.8951e+00  dscal=  1.0000e+00    scf=  1.0000e+00   upsi=  0.0000e+00
    sig=  1.0000e+00     fx= -3.7674e+00    psi=  9.8604e-01   upsi=  1.1293e-01
    sig=  5.0000e-01     fx= -3.1197e+00    psi=  7.4094e-03   upsi=  4.8179e-04


 end unimin

 sig=  5.0000e-01  num. f-evaluations 2
 list of inactive hit constraints
none
BFGS-update as in Pantoja and Mayne
  tk  =   1.087231e-01
 xsik =   2.013281e+01



================================================================================
             4-th iteration step
   scf=  1.0000e+00 psist=  0.0000e+00   psi=  7.4094e-03  upsi=  4.8179e-04
  fxst= -2.0000e+00    fx= -3.1197e+00
  x=
   6.8212e-01   1.0988e+00
 valid permutation of x

  2   1 

  multipliers: first estimate
  u =
    1   9.3832e+00    2  -3.5434e-16


  del=  1.00000e-01  b2n0=  1.11022e-16   b2n=  3.70864e-16   gfn=  1.21355e+01


 values of restrictions
 (   1   -4.8179e-04    1.2933e+00)  (   2    1.1049e-01    2.4145e+00)  
 

   diag[r]=
   1.0000e+00   4.7626e-03

  multipliers: second estimate
  u =
    1   9.3832e+00    2  -3.5434e-16
  condition number of r     2.099675724276054e+02
  condition number of a     5.969988975168789e+05


  current scaling,  scf =   1.0000e+00 clow =            1 eta =   1.1600e-02


  scalres=
   1.8966e+01   2.6000e+00   1.0000e+00   1.0000e+00


  phase=  2  scf0=  1.0000e+00
  d =
   4.4109e-02  -7.0349e-02


  dd=
   0.0000e+00   0.0000e+00


 start unimin

    phi= -3.1106e+00   dphi= -4.6171e-03    psi=  9.1379e-03 tau0/2=  5.0000e-01
     fx= -3.1197e+00  dscal=  1.0000e+00    scf=  1.0000e+00   upsi=  4.8179e-04
    sig=  1.0000e+00     fx= -3.1446e+00    psi=  5.8854e-02   upsi=  3.1031e-03
    sig=  5.0000e-01     fx= -3.1247e+00    psi=  1.9282e-02   upsi=  1.0167e-03
    sig=  1.2500e-01     fx= -3.1196e+00    psi=  8.9152e-03   upsi=  4.7005e-04


 end unimin

 sig=  1.2500e-01  num. f-evaluations 3
 list of inactive hit constraints
none
BFGS-update as in Pantoja and Mayne
  tk  =   6.894646e-03
 xsik =   0.000000e+00



================================================================================
             5-th iteration step
   scf=  1.0000e+00 psist=  0.0000e+00   psi=  8.9152e-03  upsi=  4.7005e-04
  fxst= -2.0000e+00    fx= -3.1196e+00
  x=
   6.8763e-01   1.0900e+00
 valid permutation of x

  2   1 

  multipliers: first estimate
  u =
    1   9.3823e+00    2  -4.6621e-16


  del=  5.17942e-02  b2n0=  3.45968e-12   b2n=  1.51463e-15   gfn=  1.20918e+01


 values of restrictions
 (   1   -4.7005e-04    1.2888e+00)  (   2    1.2423e-01    2.3985e+00)  
 

   diag[r]=
  -1.0000e+00  -9.9995e-01

  multipliers: second estimate
  u =
    1   9.3823e+00    2  -4.6621e-16
  condition number of r     1.000050119836917e+00
  condition number of a     6.822601719592174e+08


  current scaling,  scf =   1.0000e+00 clow =            1 eta =   1.1600e-02


  scalres=
   1.8966e+01   2.6000e+00   1.0000e+00   1.0000e+00


  phase=  2  scf0=  1.0000e+00
  d =
   5.1192e-02  -8.0466e-02


  dd=
   0.0000e+00   0.0000e+00


 start unimin

    phi= -3.1107e+00   dphi= -4.5051e-03    psi=  8.9152e-03 tau0/2=  5.0000e-01
     fx= -3.1196e+00  dscal=  1.0000e+00    scf=  1.0000e+00   upsi=  4.7005e-04
    sig=  1.0000e+00     fx= -3.1544e+00    psi=  7.8127e-02   upsi=  4.1192e-03
    sig=  5.0000e-01     fx= -3.1271e+00    psi=  2.3989e-02   upsi=  1.2648e-03
    sig=  1.2500e-01     fx= -3.1197e+00    psi=  9.0216e-03   upsi=  4.7566e-04
    sig=  6.2500e-02     fx= -3.1195e+00    psi=  8.6632e-03   upsi=  4.5676e-04


 end unimin

 sig=  6.2500e-02  num. f-evaluations 4
 list of inactive hit constraints
none
restart with scaled unit matrix
restart with scaled unit matrix



================================================================================
             6-th iteration step
   scf=  1.0000e+00 psist=  0.0000e+00   psi=  8.6632e-03  upsi=  4.5676e-04
  fxst= -2.0000e+00    fx= -3.1195e+00
  x=
   6.9083e-01   1.0850e+00
 valid permutation of x