recoverSolution -- factor a rank one positive semidefinite matrix

• Usage:

  sol = recoverSolution(mon,X)

  sol = recoverSolution(sdpsol)

• Inputs:
  • mon, a matrix, of monomials
  • X, a matrix, the moment matrix
  • sdpsol, an instance of the type SDPResult,
• Outputs:
  • sol, a list, the optimal solution

i1 : R = RR[x,y];

i2 : mon = matrix {{1},{x},{y}};

             3      1
o2 : Matrix R  <-- R

i3 : X = matrix(RR, {{1,0,1},{0,0,0},{1,0,1}} );

                3         3
o3 : Matrix RR    <-- RR
              53        53

i4 : sol = recoverSolution(mon,X)

o4 = {x => 0, y => 1}

o4 : List