lyapunov2_shift.m

 関数 solvesdp のオプション shift を利用した,lyapunov2_eps.m と等価な記述.

実行結果

>> lyapunov2_shift

SeDuMi 1.3 by AdvOL, 2005-2008 and Jos F. Sturm, 1998-2003.
Alg = 2: xz-corrector, theta = 0.250, beta = 0.500
Put 1 free variables in a quadratic cone
eqs m = 3, order n = 7, dim = 11, blocks = 4
nnz(A) = 11 + 0, nnz(ADA) = 9, nnz(L) = 6
 it :     b*y       gap    delta  rate   t/tP*  t/tD*   feas cg cg  prec
  0 :            2.53E+001 0.000
  1 :  0.00E+000 5.36E+000 0.000 0.2117 0.9000 0.9000   0.79  1  1  4.5E+000
  2 :  0.00E+000 9.86E-001 0.000 0.1840 0.9000 0.9000   0.33  1  1  1.1E+000
  3 :  0.00E+000 2.12E-002 0.000 0.0215 0.9900 0.9900   1.03  1  1  6.8E-002
  4 :  0.00E+000 1.27E-005 0.000 0.0006 0.9999 0.9999   1.00  1  1  1.4E-005
  5 :  0.00E+000 9.87E-007 0.437 0.0776 0.9900 0.9900   0.30  1  1  2.1E-006
  6 :  0.00E+000 7.55E-008 0.000 0.0765 0.9900 0.9900  -0.52  1  1  1.1E-006
  7 :  0.00E+000 4.55E-011 0.319 0.0006 0.9999 0.9999  -0.97  1  1  6.7E-006
  8 :  0.00E+000 2.99E-013 0.000 0.0066 0.9990 0.9990  -1.06  2  2  1.0E-006
  9 :  0.00E+000 7.34E-014 0.000 0.2457 0.9000 0.9000  -1.00  3  3  1.1E-006
 10 :  0.00E+000 2.85E-014 0.000 0.3884 0.9000 0.9000  -1.00  3  3  1.1E-006
Run into numerical problems.

iter seconds digits       c*x               b*y
 10      0.3   6.9 -2.6498073166e+001  0.0000000000e+000
|Ax-b| =  7.4e-007, [Ay-c]_+ =  1.1E-007, |x|= 2.4e+008, |y|= 9.0e-001

Detailed timing (sec)
   Pre          IPM          Post
7.001E-003    2.080E-001    3.007E-003    
Max-norms: ||b||=0, ||c|| = 1,
Cholesky |add|=0, |skip| = 0, ||L.L|| = 1.08759.
sol = 
    yalmiptime: 0.1220
    solvertime: 0.2190
          info: 'Infeasible problem (SeDuMi-1.3)'
       problem: 1
        dimacs: [7.7443e-007 0 0 5.6295e-008 -0.9636 -0.9636]
    
前のページ (kawata_ebihara.zip) に戻る