Solvemoment error "maximum size requested"

30 views
Skip to first unread message

researcher

unread,
Mar 26, 2020, 10:01:40 AM3/26/20
to YALMIP
I am trying to use solvemoment for the nonconvex polynomial program attached below (which tries to maximize the determinant of a non-symmetric matrix). However, it crashes after 10 minutes, with the error "Requested 46981200x47 (16.5GB) array exceeds maximum
array size preference" (see full output in the attached file)

Is there a way to reduce the size or simplify the moment relaxation problem?

Thanks!
problem.m
output.txt

Johan Löfberg

unread,
Mar 26, 2020, 1:27:36 PM3/26/20
to YALMIP
First, one typically never works explicitly with a det operator since it is such an awful mathematical object. For that reason, YALMIP does not construct an explicit representation of the det by default. What you want is 

obj=-det(A,'polynomial')

Then you get the degree 20 and 3260 term polynomial you want.

A degree 20 polynomial in 10 variables requires O(10^(20/2)) monomials in the basis, thus working with matrices of size 10^10 x 10^10. Ain't gonna happen.

To get anywhere here, you have to do a lot of research on reduced representations for this particular instance, or simply solve it using some other strategy than moments which is nice in theory but horrible in practice most often.
Reply all
Reply to author
Forward
0 new messages