starting from a previouly calculated generation

[ZeusMarti]

Hi there,

I have started using MG-GPO to optimize the dynamical aperture and beam lifetime. First of all, thanks for sharing this type of content, it is very nice, I really apreciate it.

I was trying to start from a previously calculated generation 2 of a diferent run. I did:

load('generation_2.mat');
gbest = MGGPO(evaluate,predict,Npop,Ngen,Nobj,Nvar,da);

It produced ar error because f0 iwas not defined I did the following modification in MGGPO (line 43):

elseif nargin > 6
    %{
    da = varargin{1}; 
    %     da.Xmat=[da.Xmat xt];
    %     da.fa = [da.fa; ynew];
    %     da.nf = nf;
    %     da.invKmat = inv(Kmat+sigy^2*eye(nf))
    gbest = da.gbest;
    v0 = da.v0;
    %}
    % zeus modification:
    vrange1 = ones(Nvar,1)*[0,1,1e-6]*1;
    l_limit = vrange1(:,1);
    u_limit = vrange1(:,2);
    gen = varargin{1};
    da=gen.da;
    gbest = da.gbest;
    v0 = da.v0;
    f0=gen.f0;
end

Then I run:

gen=importdata('generation_2.mat');
gbest = MGGPO(evaluate,predict,Npop,Ngen,Nobj,Nvar,gen);

This time the code calculated the evaluate function (it takes a long time to evaluate the dynamical aperture and lifetime) but it produced an error because pbest was not defined... Maybe there are more things missing. Could you take a look at it?

Also I know I'm not understanding what the field da contains, could you
briefly describe it?

Thanks a lot

[hshang315]

Xiaobiao: Can you answer his question? Thanks Hairong

…

On Tue, Oct 26, 2021 at 3:41 AM Zeus Martí Diaz ***@***.***> wrote: Hi there, I have started using MG-GPO to optimize the dynamical aperture and beam lifetime. First of all, thanks for sharing this type of content, it is very nice, I really apreciate it. I was trying to start from a previously calculated generation 2 of a diferent run. I did: load('generation_2.mat'); gbest = MGGPO(evaluate,predict,Npop,Ngen,Nobj,Nvar,da); It produced ar error because *f0* iwas not defined I did the following modification in MGGPO (line 43): elseif nargin > 6 %{ da = varargin{1}; % da.Xmat=[da.Xmat xt]; % da.fa = [da.fa; ynew]; % da.nf = nf; % da.invKmat = inv(Kmat+sigy^2*eye(nf)) gbest = da.gbest; v0 = da.v0; %} % zeus modification: vrange1 = ones(Nvar,1)*[0,1,1e-6]*1; l_limit = vrange1(:,1); u_limit = vrange1(:,2); gen = varargin{1}; da=gen.da; gbest = da.gbest; v0 = da.v0; f0=gen.f0; end Then I run: gen=importdata('generation_2.mat'); gbest = MGGPO(evaluate,predict,Npop,Ngen,Nobj,Nvar,gen); This time the code calculated the evaluate function (it takes a long time to evaluate the dynamical aperture and lifetime) but it produced an error because *pbest* was not defined... Maybe there are more things missing. Could you take a look at it? Also I know I'm not understanding what the field *da* contains, could you briefly describe it? Thanks a lot — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub <#1>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AE5RPD2J43FELNTMFC6RB4DUIZZVVANCNFSM5GXGSRBQ> . Triage notifications on the go with GitHub Mobile for iOS <https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675> or Android <https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub>.

[hshang315]

Hi, Zeus, It’s good to hear that you are using MG-GPO for DA/MA optimization. Thanks for your interest. The code indeed wasn’t ready to start from a saved generation, even though it is easy to make it so. I just make some changes to MGGPO.m and the test code ‘run.m’. I think it works now. I tested it with the option of not using GPML (should be the same for the other option). I uploaded the files to GitHub. They are also attached to this message. Please don’t hesitate to send any more questions. Best regards, Xiaobiao From: Lucy Shang ***@***.***> Sent: Tuesday, October 26, 2021 6:17 AM To: SPEAR3-ML/MG-GPO ***@***.***>; Huang, Xiaobiao ***@***.***> Subject: Re: [SPEAR3-ML/MG-GPO] starting from a previouly calculated generation (Issue #1) Xiaobiao: Can you answer his question? Thanks Hairong On Tue, Oct 26, 2021 at 3:41 AM Zeus Martí Diaz ***@***.******@***.***>> wrote: Hi there, I have started using MG-GPO to optimize the dynamical aperture and beam lifetime. First of all, thanks for sharing this type of content, it is very nice, I really apreciate it. I was trying to start from a previously calculated generation 2 of a diferent run. I did: load('generation_2.mat'); gbest = MGGPO(evaluate,predict,Npop,Ngen,Nobj,Nvar,da); It produced ar error because f0 iwas not defined I did the following modification in MGGPO (line 43): elseif nargin > 6 %{ da = varargin{1}; % da.Xmat=[da.Xmat xt]; % da.fa = [da.fa; ynew]; % da.nf<http://da.nf> = nf; % da.invKmat = inv(Kmat+sigy^2*eye(nf)) gbest = da.gbest; v0 = da.v0; %} % zeus modification: vrange1 = ones(Nvar,1)*[0,1,1e-6]*1; l_limit = vrange1(:,1); u_limit = vrange1(:,2); gen = varargin{1}; da=gen.da; gbest = da.gbest; v0 = da.v0; f0=gen.f0; end Then I run: gen=importdata('generation_2.mat'); gbest = MGGPO(evaluate,predict,Npop,Ngen,Nobj,Nvar,gen); This time the code calculated the evaluate function (it takes a long time to evaluate the dynamical aperture and lifetime) but it produced an error because pbest was not defined... Maybe there are more things missing. Could you take a look at it? Also I know I'm not understanding what the field da contains, could you briefly describe it? Thanks a lot — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub<#1>, or unsubscribe<https://github.com/notifications/unsubscribe-auth/AE5RPD2J43FELNTMFC6RB4DUIZZVVANCNFSM5GXGSRBQ>. Triage notifications on the go with GitHub Mobile for iOS<https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675> or Android<https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub>.