dual of rank constrained SDPs












0












$begingroup$


The standard rank constrained SDP is as follows:
$$
min , Tr(CX),\ text{ s.t } Tr(A_iX)= b_i text{ for } i = 1,2,...,m ,,\ X >= 0,\ Rank(X) <= r
$$



I was interested in a lower bounding a particular rank constrained SDP. I was trying out the following : 1) Formulate above using the non-linear formulation of Burer-Monteiro . 2) Find a dual feasible solution and report obj value.



The non-linear version of the above would be



$$
min , Tr(CVV^T),\ text{ s.t } Tr(A_iVV^T)= b_i text{ for } i = 1,2,...,m \V in mathbb{R}^{n times r}.
$$



However, taking the Lagrangian Dual of this doesn't seem to be fruitful, at least for my SDP of interest.



Is there another way out? If not analytically, but say numerically ? When I use a nonlinear solver like 'fmincon' on MATLAB, I end up getting a large duality gap and a bad lower bound.










share|cite|improve this question









$endgroup$












  • $begingroup$
    The paper by Burer and Monteiro has some leads on solving the nonlinear version. Can you show the dual formulation that you derived?
    $endgroup$
    – LinAlg
    Jan 4 at 14:53
















0












$begingroup$


The standard rank constrained SDP is as follows:
$$
min , Tr(CX),\ text{ s.t } Tr(A_iX)= b_i text{ for } i = 1,2,...,m ,,\ X >= 0,\ Rank(X) <= r
$$



I was interested in a lower bounding a particular rank constrained SDP. I was trying out the following : 1) Formulate above using the non-linear formulation of Burer-Monteiro . 2) Find a dual feasible solution and report obj value.



The non-linear version of the above would be



$$
min , Tr(CVV^T),\ text{ s.t } Tr(A_iVV^T)= b_i text{ for } i = 1,2,...,m \V in mathbb{R}^{n times r}.
$$



However, taking the Lagrangian Dual of this doesn't seem to be fruitful, at least for my SDP of interest.



Is there another way out? If not analytically, but say numerically ? When I use a nonlinear solver like 'fmincon' on MATLAB, I end up getting a large duality gap and a bad lower bound.










share|cite|improve this question









$endgroup$












  • $begingroup$
    The paper by Burer and Monteiro has some leads on solving the nonlinear version. Can you show the dual formulation that you derived?
    $endgroup$
    – LinAlg
    Jan 4 at 14:53














0












0








0





$begingroup$


The standard rank constrained SDP is as follows:
$$
min , Tr(CX),\ text{ s.t } Tr(A_iX)= b_i text{ for } i = 1,2,...,m ,,\ X >= 0,\ Rank(X) <= r
$$



I was interested in a lower bounding a particular rank constrained SDP. I was trying out the following : 1) Formulate above using the non-linear formulation of Burer-Monteiro . 2) Find a dual feasible solution and report obj value.



The non-linear version of the above would be



$$
min , Tr(CVV^T),\ text{ s.t } Tr(A_iVV^T)= b_i text{ for } i = 1,2,...,m \V in mathbb{R}^{n times r}.
$$



However, taking the Lagrangian Dual of this doesn't seem to be fruitful, at least for my SDP of interest.



Is there another way out? If not analytically, but say numerically ? When I use a nonlinear solver like 'fmincon' on MATLAB, I end up getting a large duality gap and a bad lower bound.










share|cite|improve this question









$endgroup$




The standard rank constrained SDP is as follows:
$$
min , Tr(CX),\ text{ s.t } Tr(A_iX)= b_i text{ for } i = 1,2,...,m ,,\ X >= 0,\ Rank(X) <= r
$$



I was interested in a lower bounding a particular rank constrained SDP. I was trying out the following : 1) Formulate above using the non-linear formulation of Burer-Monteiro . 2) Find a dual feasible solution and report obj value.



The non-linear version of the above would be



$$
min , Tr(CVV^T),\ text{ s.t } Tr(A_iVV^T)= b_i text{ for } i = 1,2,...,m \V in mathbb{R}^{n times r}.
$$



However, taking the Lagrangian Dual of this doesn't seem to be fruitful, at least for my SDP of interest.



Is there another way out? If not analytically, but say numerically ? When I use a nonlinear solver like 'fmincon' on MATLAB, I end up getting a large duality gap and a bad lower bound.







optimization nonlinear-optimization duality-theorems semidefinite-programming






share|cite|improve this question













share|cite|improve this question











share|cite|improve this question




share|cite|improve this question










asked Jan 4 at 12:15









Maharshi RayMaharshi Ray

1




1












  • $begingroup$
    The paper by Burer and Monteiro has some leads on solving the nonlinear version. Can you show the dual formulation that you derived?
    $endgroup$
    – LinAlg
    Jan 4 at 14:53


















  • $begingroup$
    The paper by Burer and Monteiro has some leads on solving the nonlinear version. Can you show the dual formulation that you derived?
    $endgroup$
    – LinAlg
    Jan 4 at 14:53
















$begingroup$
The paper by Burer and Monteiro has some leads on solving the nonlinear version. Can you show the dual formulation that you derived?
$endgroup$
– LinAlg
Jan 4 at 14:53




$begingroup$
The paper by Burer and Monteiro has some leads on solving the nonlinear version. Can you show the dual formulation that you derived?
$endgroup$
– LinAlg
Jan 4 at 14:53










0






active

oldest

votes











Your Answer





StackExchange.ifUsing("editor", function () {
return StackExchange.using("mathjaxEditing", function () {
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
});
});
}, "mathjax-editing");

StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "69"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});

function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
noCode: true, onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});


}
});














draft saved

draft discarded


















StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3061593%2fdual-of-rank-constrained-sdps%23new-answer', 'question_page');
}
);

Post as a guest















Required, but never shown

























0






active

oldest

votes








0






active

oldest

votes









active

oldest

votes






active

oldest

votes
















draft saved

draft discarded




















































Thanks for contributing an answer to Mathematics Stack Exchange!


  • Please be sure to answer the question. Provide details and share your research!

But avoid



  • Asking for help, clarification, or responding to other answers.

  • Making statements based on opinion; back them up with references or personal experience.


Use MathJax to format equations. MathJax reference.


To learn more, see our tips on writing great answers.




draft saved


draft discarded














StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3061593%2fdual-of-rank-constrained-sdps%23new-answer', 'question_page');
}
);

Post as a guest















Required, but never shown





















































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown

































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown







Popular posts from this blog

Can a sorcerer learn a 5th-level spell early by creating spell slots using the Font of Magic feature?

Does disintegrating a polymorphed enemy still kill it after the 2018 errata?

A Topological Invariant for $pi_3(U(n))$