clear Ahelp Astiff G Mass F ff t

m=par.m; n=par.n; n1=par.n1;
maska=par.maska; BI=par.BI; DELTA=par.DELTA; f=par.f; 

vol=100; compl = 100; par.cmp=compl;

t=sdpvar(m,1);
lambda = sdpvar(1,1);
setsdpvar(t,(vol/m).*ones(m,1));
setsdpvar(lambda,1.0e-5);

Ahelp=zeros(n,n);
for i=1:m
   Ahelp=Ahelp+t(i)*BI(i,:)'*BI(i,:);
end
Astiff=Ahelp(maska,maska);

ff=par.f; 

G=zeros(n1,n1);
for i=1:m
    G=G+t(i)*DELTA{i}'*DELTA{i};
end
Mass=G;

Mde = sparse(n1,n1);

%dead mass
  Mde(33,33) = 10; Mde(34,34) = 10;

F = set(100>t>0); 

Alarge = [compl -ff'; -ff Astiff];
F = F + set(Alarge>0)

lambda_bar = 1.0e-3;
F = F+ set(Astiff - lambda_bar*(Mass+Mde));

solvesdp(F,sum(t),sdpsettings('solver','penbmi','usex0',1,'savedebug',1,'showpro',1));

t = double(t);
K = double(Astiff);
M = double(Mass)+Mde;

%mypic1(par,t);
picvib1(par,t,K,M);