next up previous contents
Next: A sample LMI Up: Examples Previous: A Lovász function

A sample truss problem

>> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
>> % Sample truss problem
>> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
>> 
>> load testing/truss/truss1     % from Nemirovskii
>> setpars                       % sets scalefac = 100
>> initvars
>> sdp

tau =   0.9990,     scalefac =      100

iter   p_step      d_step     p_infeas    d_infeas      X.Z         pobj        dobj
  0   0.000e+00   0.000e+00   7.803e+02   3.603e+02   1.300e+05   1.000e+02   0.000e+00
  1   1.000e+00   7.344e-01   6.356e-14   9.570e+01   1.391e+04   2.565e+02  -6.643e+01
  2   1.000e+00   1.000e+00   2.359e-12   8.789e-15   2.470e+02   2.471e+02   1.014e-01
  3   6.313e-01   1.000e+00   8.669e-13   1.194e-16   9.229e+01   9.189e+01  -4.013e-01
  4   8.211e-01   1.000e+00   1.635e-13   1.688e-16   1.668e+01   1.745e+01   7.700e-01
  5   5.064e-03   4.190e-01   1.635e-13   2.267e-15   1.004e+01   1.746e+01   7.424e+00
  6   1.000e+00   9.086e-01   2.853e-13   1.088e-15   1.254e+00   1.006e+01   8.809e+00
  7   9.886e-01   9.977e-01   2.754e-13   2.432e-15   1.919e-02   9.018e+00   8.999e+00
  8   9.990e-01   9.990e-01   1.535e-13   2.035e-15   1.948e-05   9.000e+00   9.000e+00
  9   9.990e-01   9.990e-01   9.392e-14   1.936e-15   1.948e-08   9.000e+00   9.000e+00
 10   9.990e-01   9.990e-01   3.879e-14   2.220e-15   1.971e-11   9.000e+00   9.000e+00
fsdp: stop since error reduced to desired value

sdp: elapsed time               =   1.30909 seconds
sdp: elapsed cpu time           =   1.26000 seconds
sdp: flops                      =   8.99130e+04
sdp: Number of iterations       =   10
sdp: final value of X.Z         =   1.971e-11
sdp: final primal infeasibility =   3.879e-14
sdp: final dual infeasibility   =   2.220e-15
sdp: primal objective value     =   8.9999963153051237e+00
sdp: dual objective value       =   8.9999963152853795e+00
>> primalcond(A,blk,X,1.0e-06);  % check if primal degenerate
primalcond =         Inf
>> dualcond(A,blk,Z,1.0e-06);    % check if dual degenerate
dualcond =   8.920e+00



Madhu Nayakkankuppam
Fri Mar 28 00:48:56 EST 1997