CVX Not Running Quadratic Optimisation

1

Hello,

I’m trying to implement a simple quadratic optimisation my lecturer set me and there is an issue with my system which has completely stumped him so I would appreciate any assistance this forum can offer.

I am trying to implement this system:

image
image933×293 22.7 KB

This is my MATLAB (run with R2020b) code:
% Solving Least Squares Problem
clear all;
clc;

% Generate data
n = 10;
m = n;
P = randn(n);
P = P*P’;
q = randn(n,1);
r = randn;
G = randn(m,n);
h = randn(m,1);
A = randn(n,n);
b = randn(n,1);

% Create and solve problem
cvx_begin
variable x(n)
dual variables y1 y2
minimize(0.5x’Px + q’x +r)
subject to
y1 : Gx <= h;
y2 : Ax == b;
cvx_end

x

This is the error I’m getting:
Check for missing argument or incorrect argument data type in call to function ‘vec’.

Error in cvx/quad_form (line 43)
v = vec( v );

Error in * (line 261)
[ z2, success ] = quad_form( xx, P, Q, R );

Error in Quadratic_Problem (line 46)
minimize( 0.5*x’Px + q’*x + r );

I have installed CVX correctly as it runs for a least squares example program. However, it is highlighting that there is something wrong in the vec() function.

I would appreciate any guidance/support you can offer and if you need anymore information then please tell me and I’ll make sure to add it to the original post.

Josh

2

Your program runs without error for me, although the chances of your random inputs being feasible may not be very high,

Are you using CVX 2.2.? If not,switch to that. Otherwise, perhaps there is a conflict with a different version of vec in your MATLAB path.

Look at the output of
which -all vec

3

Hi Mark,

Many thanks for your suggestions. I have checked my CVX version and it’s 2.2. I did have multiple vec.m because I had YALMIP installed for a previous project. I have removed it from the path and the same error occurs despite restarting MATLAB and verifying the number of vec() functions:

which -all vec
D:\cvx-w64\cvx\functions\vec_\vec.m
D:\cvx-w64\cvx\functions@cvx\vec.m % cvx method

4

Nevermind Mark, it now runs without any issues.

Thank you very much, you have no idea how grateful I am!

5

What did the problem turn out to be? I.e., how did you get it to no longer be a problem?

6

I seemed to be the multiple vec()s you suggested. I’ve no clue why it didn’t initially work but after leaving it for 10 mins, it seemed to sort itself out. It’s not the handiest solution but I removed YALMIP from the path, left it and it now works.

7

Try putting YALMIP in the MATLAB path after CVX.

If you do so, I suggest removing geomean.m and logdet.m from the CVX installation; Because if you do not, those CVX files redirect to geo_mean and log_det. Which means that if YALMIP calls geomean and logdet, that would be converted into call to CVX’s geo_mean and log_det, and that would not end nicely. After removing geomean.m and logdet.m, you will still have geo_mean and log_det, which is what you actually need in CVX.

If you do all this, as far as I know, CVX and YALMIP should coexist and both work properly and completely.

And yes, i learned all this the hard way.

8

Hello Mark, I had a similar issue and tried removing the geomean as suggested but it still didn’t work for me. I also removed my YALMIP as well but that hasn’t also worked so far.

9

You only need to remove CVX’s geomean if you use YALMIP. If you are using CVX, use its command, which is geo_mean. I am not aware of there being any other geo_mean.

if in doubt, use
which -all name_of_function
to see all functions in order of appearance (precedence) in your MATLAB path which are named name_of_function.

10

I have a similar problem which only appears if the problem (XX) is sufficiently large

I did check

which -all vec.m
/home/tk/Matlab/cvx/functions/@cvx/vec.m % cvx method

But I have a vec.p file

which -all vec
/usr/local/MATLAB/R2024a/toolbox/shared/msblks_measurement/msblks_measurement/@eyeDiagramSI/vec.p % eyeDiagramSI method
/home/tk/Matlab/cvx/functions/@cvx/vec.m % cvx method

Here is the execution given error followed by one given no error

[u,slack,X0Input,Covxt]= SMPC(Phi,Gamma,eye(n),Q,C,D,R,LambdaX*(C’*C),T,…
G,g,F,f,100,pb,xh,P,ProbStateConsMP);

[u,slack,X0Input,Covxt]= SMPC(Phi,Gamma,eye(n),Q,C,D,R,LambdaX*(C’*C),T,…

                   G,g,F,f,100,pb,xh,P,ProbStateConsMP);

Warning: A non-empty cvx problem already exists in this scope.
It is being overwritten.

In cvxprob (line 28)
In cvx_begin (line 41)
In SMPC (line 121)
Incorrect number or types of inputs or outputs for function vec.
Error in cvx/quad_form (line 43)
v = vec( v );
Error in * (line 261)
[ z2, success ] = quad_form( xx, P, Q, R );
Error in SMPC (line 125)
V= X’SSX+U’TTU+ProbStateCons.SlackPenalty*s;

[u,slack,X0Input,Covxt]= SMPC(Phi,Gamma,eye(n),Q,C,D,R,LambdaX*(C’*C),T,…
G,g,F,f,10,pb,xh,P,ProbStateConsMP);

[u,slack,X0Input,Covxt]= SMPC(Phi,Gamma,eye(n),Q,C,D,R,LambdaX*(C’*C),T,…

                   G,g,F,f,10,pb,xh,P,ProbStateConsMP);

Warning: A non-empty cvx problem already exists in this scope.
It is being overwritten.

In cvxprob (line 28)
In cvx_begin (line 41)
In SMPC (line 121)

11

Either try moving CVX’s vec ahead of the vec.p in the MATLAB path, or change the quad_form code to v = v(:);

12

Thanks I’ll try that.

13

cvx_begin quiet
variable w(NumofAntenna) complex;
minimize( w’ * R_x * w );
subject to
norm((a0 * w’),2) <= 1;
cvx_ end
the above is my code , I am having error in vec function but my friend does not have the problem in executing the same code with cvx tool. the error is
Error in cvx/quad_form (line 43)
v = vec( v );
^^^^^^^^^^^^^
Error in * (line 261)
[ z2, success ] = quad_form( xx, P, Q, R );
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Error in mam_code (line 184)
minimize( w’ * R_x * w );
^^^^^^^^^^^^^^^^^^^^^^^^^
please can you help with this

14

There might be another vec in your MATLAB path ahead of the CVX version.

Look at
which -all vec

Move the CVX directories ahead of the other vec in your MATLAB path.
If you have both YALMIP and CVX installed, I recommend CVX be ahead of YA:MIP in the MATLAB path.

Actually, I have no new suggestions beyond those I gave previously in this thread to other users.