関数 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]