>> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% >> % 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