Re: glpk problem,please check code

[Ben Abbott]

> On Feb 12, 2017, at 10:54 AM, saiko7 <address@hidden> wrote:
> 
> function [alpha]=mclp(B,c,d);
> 
> a=B(:,1:29);
> r=size(a,1);
> c=size(a,2);
> Q=gallery("lehmer",r);Q2=[Q zeros(r,r);zeros(r,r) zeros(r,r)];
> alpha=zeros(r,1);
> e=ones(r,1);e2=[e;zeros(r,1)];
> u=zeros(r,1);u2=[u;zeros(r,1)];
> v=zeros(r,1);v2=[v;zeros(r,1)];
> y=B(:,30);
> p=sum(y==0);q=sum(y==1);
> 
> E=[eye(q,q) zeros(q,p);zeros(p,q) -eye(p,p)];E2=[E zeros(r,r);zeros(r,r)
> zeros(r,r)];
> alpha=zeros(r,1);
> b=a*a'+e*e';
> H=[b zeros(r,r);zeros(r,r) Q];
> alpha2=[alpha;u];alpha2(1,1)=0.000000001;
> 
> [alpha2 obj
> inf]=qp(alpha2,H,zeros(2*r,1),[],[],[],[],(-Q2*u2)-(d*e2),E2,-c*e2);
> 
> ihave checked the dimensions of the matrices invoved ,but keep getting the
> message 
> mclp(a,3,1)
> error: glpk: A must be a real valued 8 by 12 matrix
> error: called from
>    glpk at line 506 column 5
>    qp at line 370 column 29
>    mclp at line 20 column 16
> 
> where a is a 4*30 matrix

I’m seeing a possible problem with glpk as well. What version of Octave are you 
using and what is the value for “a” when you call mclp(a,3,1)?

Ben