catch22 / quantbox Goto Github PK

% prepare Ψ ⊗ Ψ for a Haar random pure state Ψ
[U, ~, ~] = svd(rand(3) + 1i*rand(3));
sigma = U * diag([1 0 0]) * U';
sigma = sigma / trace(sigma);
rho = kron(sigma, sigma);
dims = [3 3];

% try to find 2:1 extension that is 1:0-PPT
[RHO_EXT, ~, ~, ~, INFO] = sym_extension(rho, dims, [2 1], [1 0], [1 1], 'sdpt3');

error:  total dimension of C should be > length(b)
error: called from
    validate at line 143 column 5
    sqlp at line 159 column 38
    solve_sdp>solve_sdp_with_sdpt3 at line 201 column 17
    solve_sdp at line 177 column 3
    sym_extension at line 120 column 14