{"action":"opened","issue":{"url":"https://api.github.com/repos/sqlp/sedumi/issues/42","repository_url":"https://api.github.com/repos/sqlp/sedumi","labels_url":"https://api.github.com/repos/sqlp/sedumi/issues/42/labels{/name}","comments_url":"https://api.github.com/repos/sqlp/sedumi/issues/42/comments","events_url":"https://api.github.com/repos/sqlp/sedumi/issues/42/events","html_url":"https://github.com/sqlp/sedumi/issues/42","id":395090081,"node_id":"MDU6SXNzdWUzOTUwOTAwODE=","number":42,"title":"Incorrect solution for LP","user":{"login":"blegat","id":1048205,"node_id":"MDQ6VXNlcjEwNDgyMDU=","avatar_url":"https://avatars3.githubusercontent.com/u/1048205?v=4","gravatar_id":"","url":"https://api.github.com/users/blegat","html_url":"https://github.com/blegat","followers_url":"https://api.github.com/users/blegat/followers","following_url":"https://api.github.com/users/blegat/following{/other_user}","gists_url":"https://api.github.com/users/blegat/gists{/gist_id}","starred_url":"https://api.github.com/users/blegat/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/blegat/subscriptions","organizations_url":"https://api.github.com/users/blegat/orgs","repos_url":"https://api.github.com/users/blegat/repos","events_url":"https://api.github.com/users/blegat/events{/privacy}","received_events_url":"https://api.github.com/users/blegat/received_events","type":"User","site_admin":false},"labels":[],"state":"open","locked":false,"assignee":null,"assignees":[],"milestone":null,"comments":0,"created_at":"2019-01-01T14:57:50Z","updated_at":"2019-01-01T14:57:50Z","closed_at":null,"author_association":"NONE","body":"Consider the problem\r\n```\r\nmin -x2\r\n- x1 -  x2 = 0\r\n  x1 - 3x2 = 0\r\nx1 free\r\nx2 >= 0\r\n```\r\nwith dual\r\n```\r\nmax 0\r\n-y1 +  y2 = 0\r\n2y1 + 3y2 >= 1\r\n```\r\nThe only feasible solution of the primal is `(0, 0)` and for the dual, the feasible solution is the ray `(1, 1)` starting at `(1/5, 1/5)`.\r\nSeDuMi returns the correct `x` `(0, 0)` but an incorrect `ỳ` `(-1.2, 0.6)`.\r\n\r\nThe problem is the [13th test problem of MOI](https://github.com/JuliaOpt/MathOptInterface.jl/blob/999247c5e548c2235e49bec196809dc679c1948e/src/Test/contlinear.jl#L1276) and it is used to verify the correctness of the [Julia SeDuMi wrapper](https://github.com/blegat/SeDuMi.jl).\r\nHere is what I get when using `sedumi` directly:\r\n```julia\r\njulia> A = [-1.0 -2.0; 1.0 -3.0]\r\n2×2 Array{Float64,2}:\r\n -1.0  -2.0\r\n  1.0  -3.0\r\n\r\njulia> b = [0.0, 0.0]\r\n2-element Array{Float64,1}:\r\n 0.0\r\n 0.0\r\n\r\njulia> c = [0.0; -1.0]\r\n2-element Array{Float64,1}:\r\n  0.0\r\n -1.0\r\n\r\njulia> using SeDuMi\r\n\r\njulia> sedumi(A, b, c, SeDuMi.Cone(1, 1))\r\n>> SeDuMi 1.32 by AdvOL, 2005-2008 and Jos F. Sturm, 1998-2003.\r\nAlg = 2: xz-corrector, Adaptive Step-Differentiation, theta = 0.250, beta = 0.500\r\nPut 1 free variables in a quadratic cone\r\neqs m = 2, order n = 4, dim = 4, blocks = 2\r\nnnz(A) = 4 + 0, nnz(ADA) = 4, nnz(L) = 3\r\n it :     b*y       gap    delta  rate   t/tP*  t/tD*   feas cg cg  prec\r\n  0 :            1.25E+01 0.000\r\n  1 :   0.00E+00 1.44E-06 0.000 0.0000 1.0000 1.0000   1.00  1  1  1.0E-06\r\n  2 :   0.00E+00 5.75E-10 0.000 0.0004 0.9999 0.9999   1.00  1  2  3.5E-10\r\n\r\niter seconds digits       c*x               b*y\r\n  2      0.4   Inf -3.2329888893e-11  0.0000000000e+00\r\n|Ax-b| =   9.8e-10, [Ay-c]_+ =   1.6E-16, |x|=  1.4e+00, |y|=  1.3e+00\r\n\r\nDetailed timing (sec)\r\n   Pre          IPM          Post\r\n3.601E-01    2.107E-01    2.660E-02    \r\nMax-norms: ||b||=0, ||c|| = 1,\r\nCholesky |add|=0, |skip| = 1, ||L.L|| = 2.\r\n([-4.67511e-10, 3.23299e-11], [-1.2, 0.6], Dict{String,Any}(\"timing\"=>[0.360066 0.210701 0.0266042],\"dinf\"=>0.0,\"pinf\"=>0.0,\"cpusec\"=>1.66,\"wallsec\"=>0.597372,\"feasratio\"=>1.0,\"r0\"=>9.75092e-10,\"iter\"=>2.0,\"numerr\"=>0.0))\r\n```\r\n\r\nAny idea what's the issue ?"}}
{"actor":{"display_login":"blegat"}}