pdfmpc/update_mv.m at main · mazenalamir/pdfmpc · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
%---------------------------------------------------------
% Author : Mazen Alamir
% CNRS, UNiversity of Grenoble-Alpes
% Gipsa Lab. Last modification March 2017
%---------------------------------------------------------

function [u_sol,p_sol,lesp]=update_mv(mv,subset,param)
coder.allowpcode('plain');
np=length(mv);
Nloop=max(2,fix((param.Nev-1)/(4*np*param.Niter)));
lesp=zeros(np,Nloop*np);
ind=1;
pmin=zeros(np,1);
pmax=zeros(np,1);
for i=1:np,
    lesp(i,1)=mv(i).p;
    pmin(i)=mv(i).pmin;
    pmax(i)=mv(i).pmax;
end
npr=length(subset);
ell=subset(1);
for i=1:Nloop-1,
    for sig=1:npr,
        param.p=lesp(:,ind);param.ell=ell;
        param.pmin=pmin;param.pmax=pmax;
        mv(ell).p=lesp(ell,ind);
        mv(ell)=update_sc(mv(ell),param,param.Niter);
        lesp(:,ind+1)=lesp(:,ind);lesp(ell,ind+1)=mv(ell).p;
        if (ell<npr),
            ell=subset(sig+1);
        else
            ell=subset(1);
        end
        ind=ind+1;
    end
end
lesp=lesp(:,1:ind-1);
p_sol=lesp(:,end);
u_sol=user_control_profile(p_sol,param.ode,param.uparam);
return