grain_validation_minimization.m

function [ all,fval] = grain_validation_minimization( B,G,E,D,all0)
%all0=horzcat(G.r(division,:),G.x,D(1).pos,D(1).tilt);%,D(2).pos,D(2).tilt);
% lall=horzcat(all0(1:3)-0.66,all0(4:6)-10,D(1).lpostol,-D(1).tilttol);%,D(2).lpostol,-D(2).tilttol);
% uall=horzcat(all0(1:3)+0.66,all0(4:6)+10,D(1).upostol,D(1).tilttol);%,D(2).upostol,D(2).tilttol);
% options = optimoptions(@fmincon,'Display','off','FunValCheck','off','UseParallel',true);
% [all,fval,exitflag,output]=fmincon(@fun,all0,[],[],[],[],lall,uall,[],options);
%all0=horzcat(G.r(division,:),G.x,D(1).pos(1:2),D(2).pos,D(1).tilt,D(2).tilt);
lall=horzcat(all0(1:3)-0.66,all0(4:6)-10,D(1).lpostol(1:2),D(2).lpostol,-D(1).tilttol,-D(2).tilttol);
uall=horzcat(all0(1:3)+0.66,all0(4:6)+10,D(1).upostol(1:2),D(2).upostol,D(1).tilttol,D(2).tilttol);
options = optimoptions(@fmincon,'Display','off','FunValCheck','off','UseParallel',true);
[all,fval]=fmincon(@fun,all0,[],[],[],[],lall,uall,[],options);
all= [all(1:8),0,all(9:17)];
function [cost]=fun(allx)
    r=allx(1:3);
    x=allx(4:6);
    Dpos1=allx(7:8);
    Dpos1(3)=0;
    Dpos2=allx(9:11);
    Dtilt1=allx(12:14);
    Dtilt2=allx(15:17);
    [L] =laue_continuous_source(B,G,r);% Get the reflections for all divisions and omegas
    [spositions1,spositions2] =generate_diff_spots(B,D,L,x,Dpos1,Dpos2,Dtilt1,Dtilt2);
    [temp1,cost,ref1,ref2] = compare_G_and_E(B,E,spositions1,spositions2);
end
end