DEIMEstimatorAnalyzer.m

Go to the documentation of this file.

/* (Autoinserted by mtoc++)
 * This source code has been filtered by the mtoc++ executable,
 * which generates code that can be processed by the doxygen documentation tool.
 *
 * On the other hand, it can neither be interpreted by MATLAB, nor can it be compiled with a C++ compiler.
 * Except for the comments, the function bodies of your M-file functions are untouched.
 * Consequently, the FILTER_SOURCE_FILES doxygen switch (default in our Doxyfile.template) will produce
 * attached source files that are highly readable by humans.
 *
 * Additionally, links in the doxygen generated documentation to the source code of functions and class members refer to
 * the correct locations in the source code browser.
 * However, the line numbers most likely do not correspond to the line numbers in the original MATLAB source files.
 */

function varargout = DEIMEstimatorAnalyzer(varargin) {


/*  Edit the above text to modify the response to help DEIMEstimatorAnalyzer */

/*  Last Modified by GUIDE v2.5 12-Oct-2012 13:53:00 */

/*  Begin initialization code - DO NOT EDIT */
gui_Singleton = 0;
gui_State = struct(" gui_Name ",       mfilename, ...
                   " gui_Singleton ",  gui_Singleton, ...
                   " gui_OpeningFcn ", @DEIMEstimatorAnalyzer_OpeningFcn, ...
                   " gui_OutputFcn ",  @DEIMEstimatorAnalyzer_OutputFcn, ...
                   " gui_LayoutFcn ",  [] , ...
                   " gui_Callback ",   []);
if nargin && ischar(varargin[1])
    gui_State.gui_Callback= str2func(varargin[1]);
end

if nargout
    [varargout[1:nargout]] = gui_mainfcn(gui_State, varargin[:]);
else
    gui_mainfcn(gui_State, varargin[:]);
end
/*  End initialization code - DO NOT EDIT */


/*  --- Executes just before DEIMEstimatorAnalyzer is made visible. */
}
function  DEIMEstimatorAnalyzer>DEIMEstimatorAnalyzer_OpeningFcn(hObject,eventdata,handles,varargin) {

/*  Choose default command line output for DEIMEstimatorAnalyzer */
handles.output= hObject;

/*  Update handles structure */
guidata(hObject, handles);

/*  Parse inputs */

ip.addRequired('Model',@(r)isa(r," models.ReducedModel "));
ip.addParamValue('StartIdx',4,@isposintscalar);
ip.parse(varargin[:]);
res = ip.Results;
r = res.Model;
setappdata(handles.main," startidx ",res.StartIdx);

mu = r.System.mu;
if isempty(mu)
    mu = r.ParamSamples(:,1);
end
setappdata(handles.main," mu ",mu);
in = [];
if r.System.InputCount > 0
    in = 1;
end
setappdata(handles.main," inputidx ",in);

/*  Care condition: dont simulate very expensive by default! */
e = r.ErrorEstimator;
if ~isempty(e) && e.JacSimTransMaxSize > 0 && e.JacSimTransSize == 0
    e.JacSimTransSize= 1;
end
setappdata(handles.main," r ",r);
createModelParamSliders(handles,r);
modelToGUI(handles, r);

extra = ;
sr = r.FullModel.SpaceReducer;
if sr.IncludeTrajectoryFxiData
    extra = " InclFxiData:yes ";
end
if sr.IncludeFiniteDifferences
    extra = [extra " InclFinDiff:yes "];
end
n = sprintf(" Model '%s', %d-dim, reduced %d-dim, %s ",...
    r.Name,r.FullModel.Dimension,size(r.V,2),extra);
set(hObject," Name ",n);
set(handles.rbgAlpha," SelectionChangeFcn ",@(s,e)rbgAlpha_SelectionChanged(handles,s,e));
set(handles.rbgBeta," SelectionChangeFcn ",@(s,e)rbgBeta_SelectionChanged(handles,s,e));
set(handles.rbg_Param," SelectionChangeFcn ",@(s,e)rbgParam_SelectionChanged(handles,s,e));

reSimulate(handles);

}

function  DEIMEstimatorAnalyzer>createModelParamSliders(h,r) {
if r.System.ParamCount > 0
    /*  Setup training sample slider */
    set(h.slTrainParam," Min ",1," Max ",size(r.ParamSamples,2)," Value ",1);
    /*  Create model parameter slides */
    parent = h.rbg_Param;
    set(parent," Units "," pixels ");
    pos = get(parent," Position ");
    set(parent," Units "," normalized ");
    dist = 22;/* px */

    pcnt = 0;
    pdata = struct;
    pdata.sl= zeros(r.System.ParamCount,1);
    pdata.lbl= zeros(r.System.ParamCount,1);
    mu = getappdata(h.main," mu ");
    for pidx = 1:r.System.ParamCount
        p = r.System.Params(pidx);
        if p.HasRange
            /*  Top location */
            top = pos(4)-pcnt*dist-40-2*dist;

            /* % Create name labels */
            label = uicontrol(" Tag ",sprintf(" lblP%d ",pidx)," Style "," text ",...
                    " Parent ",parent," HorizontalAlignment "," left ");
            set(label," String ",[p.Name " : "]," Units "," pixels ",...
                " Position ",[10 top 60 14]);
            set(label," Units "," normalized ");

            /* % Labels for min/max values */
            label = uicontrol(" Tag ",[" runtime_lbll " num2str(pidx)]," Style "," text ",...
                    " Parent ",parent," HorizontalAlignment "," right ");
            set(label," String ",sprintf(" %g ",p.MinVal)," Units "," pixels ",...
                " Position ",[70 top 40 14]);
            set(label," Units "," normalized ");
            label = uicontrol(" Tag ",[" runtime_lblr " num2str(pidx)]," Style "," text ",...
                    " Parent ",parent," HorizontalAlignment "," left ");
            set(label," String ",sprintf(" %g ",p.MaxVal)," Units "," pixels ",...
                " Position ",[290 top 40 14]);
            set(label," Units "," normalized ");

            /* % Create slider */
            ctrl = uicontrol(" Tag ",sprintf(" slP%d ",pidx)," Parent ",...
                parent," Style "," slider "," UserData ",pidx);
            /*  Position */
            set(ctrl," Units "," pixels "," Position ",...
                [110 top 170 16]);
            set(ctrl," Units "," normalized ");
            /*  Range etc */
            set(ctrl," Min ",p.MinVal," Max ",p.MaxVal," Value ",mu(pidx));
            set(ctrl," SliderStep ",[0.01 0.1]);
            /*  Set callback & string */
            set(ctrl," Callback ",@(~,~)(updateUserParam(h)));
            pdata.sl(pidx) = ctrl;

            /* % Create value labels */
            label = uicontrol(" Tag ",sprintf(" lblPV%d ",pidx)," Style "," text ",...
                    " Parent ",parent," HorizontalAlignment "," left ");
            set(label," String ",sprintf(" %g ",mu(pidx))," Units "," pixels ",...
                " Position ",[340 top 40 14]);
            set(label," Units "," normalized ");
            pdata.lbl(pidx) = label;

            /*  increase position counter */
            pcnt = pcnt + 1;
        end
    end
    setappdata(h.main," pdata ",pdata);
end

}

function  DEIMEstimatorAnalyzer>updateUserParam(h) {
r = getappdata(h.main," r ");
pdata = getappdata(h.main," pdata ");
pc = r.System.ParamCount;
if pc > 0
    set(h.rbCustomParam," Value ",1);
    mu = zeros(pc,1);
    for pidx = 1:pc
        if r.FullModel.System.Params(pidx).HasRange
            mu(pidx) = get(pdata.sl(pidx)," Value ");
        else
            mu(pidx) = r.FullModel.System.Params(pidx).MinVal;
        end
        s = sprintf(" %g ",r.System.Params(pidx).Name,mu(pidx));
        set(pdata.lbl(pidx)," String ",s);
    end
    setappdata(h.main," mu ",mu);
end
reSimulate(h);

}

function  DEIMEstimatorAnalyzer>modelToGUI(h,r) {
    
    set(h.slM," Min ",1," Max ",r.System.f.MaxOrder," Value ", r.System.f.Order(1));
    set(h.lblM," String ",sprintf(" %d ",r.System.f.Order(1)));

    /*  m' estimates */
    e = r.ErrorEstimator;
    if ~isempty(e) && e.UseTrueDEIMErr
        set(h.rbTrueDEIM," Value ",1);
    else
        set(h.rbMD," Value ",1);
    end
    mv = r.System.f.MaxOrder-r.System.f.Order(1);
    if mv ~= 0, en=" on "; else en=" off "; end
    set(h.slMD," Min ",0," Max ",max(mv,1)," Enable ",en," Value ", r.System.f.Order(2));
    set(h.lblMD," String ",sprintf(" %d ",r.System.f.Order(2)));

    /*  mj JacMDEIM */
    if ~isempty(e)
        set(h.slJM," Min ",1," Max ",e.JacMatDEIMMaxOrder," Value ", e.JacMatDEIMOrder);
        set(h.lblJM," String ",sprintf(" %d ",e.JacMatDEIMOrder));
        /*  simtrans k */
        set(h.slK," Min ",0," Max ",e.JacSimTransMaxSize," Value ", e.JacSimTransSize);
        set(h.lblK," String ",sprintf(" %d ",e.JacSimTransSize));
        if e.UseTrueLogLipConst
            set(h.rbtrueloclip," Value ",1);
        elseif e.UseJacobianLogLipConst
            set(h.rbjacloclip," Value ",1);
        elseif e.UseFullJacobian
            set(h.rbfulljac," Value ",1);
        else
            set(h.rbjmst," Value ",1);
        end
    end

}

function  DEIMEstimatorAnalyzer>reSimulate(h) {
drawnow;
fprintf(" Simulating...  ");
r = getappdata(h.main," r ");
mu = getappdata(h.main," mu ");
in = getappdata(h.main," inputidx ");
s = struct;
[~, s.y, s.tf, s.x] = r.FullModel.simulate(mu,in);
[s.t, s.yr, s.tr, s.xr] = r.simulate(mu,in);
setappdata(h.main," s ",s);
rePlot(h);

}

function  DEIMEstimatorAnalyzer>rePlot(h) {
fprintf(" plotting...  ");
r = getappdata(h.main," r ");
s = getappdata(h.main," s ");
esterr = zeros(size(r.Times));
if get(h.rbOutput," Value ") /*  Output error plot */

    Cn = normest(r.FullModel.System.C.C);
    err = Cn*1e-7*Norm.L2(s.x-r.V*s.xr);
    /* err2 = Norm.L2(r.FullModel.System.computeOutput(s.x-r.V*s.xr));
     *trajnorm = Norm.L2(s.y-s.yr); */
    trajnorm = Norm.L2(s.y);
    if ~isempty(r.ErrorEstimator)
        esterr = r.ErrorEstimator.OutputError;
    end
else /*  State space error plot */

    err = Norm.L2(s.x-r.V*s.xr);
    trajnorm = Norm.L2(s.x);
    if ~isempty(r.ErrorEstimator)
        esterr = r.ErrorEstimator.StateError;
    end
end
doPlot(h.axerr,r.Times,err," b ");
doPlot(h.axrel,r.Times,err./trajnorm," b ");
if ~isempty(r.ErrorEstimator) && r.ErrorEstimator.Enabled
    hold(h.axerr," on ");
    doPlot(h.axerr,r.Times,esterr," r ");
    hold(h.axerr," off ");
    /*  Relative error */
    hold(h.axrel," on ");
    doPlot(h.axrel,r.Times,esterr./trajnorm," r ");
    hold(h.axrel," off ");
    /*  Effectivities */
    doPlot(h.axeff,r.Times,esterr./err," g ");
    eest = esterr(end);
else
    cla(h.axeff);
    eest = -1;
end
str = sprintf(" Times: Full %gs, Red.: %gs, Speedup: %g\n ",s.tf,s.tr,s.tf/s.tr);
str = [str sprintf(" Errors(T=%g): True: %g, Est.: %g,\nRel: %g, Est.Rel.: %g, Eff.:%g ",...
    r.Times(end),err(end),eest,err(end)/trajnorm(end),eest/trajnorm(end),eest/err(end))];
set(h.lblRes," String ",str);
if ~isempty(r.ErrorEstimator)
    doPlot(h.axhlp,r.Times,r.ErrorEstimator.LastAlpha," m ");
    doPlot(h.axhlp2,r.Times,r.ErrorEstimator.LastBeta," c ");
    axis(h.axerr," tight ");
    if r.ErrorEstimator.Enabled
        axis(h.axeff," tight ");
        axis(h.axhlp," tight ");
        axis(h.axhlp2," tight ");
    end
end
fprintf(" done.\n ");

}

function  DEIMEstimatorAnalyzer>doPlot(h,colvec<double> x,matrix<double> y,varargin) {
startidx = getappdata(get(h," Parent ")," startidx ");
LogPlot.cleverPlot(h,x(startidx:end),y(startidx:end),varargin[:]);

}

function  DEIMEstimatorAnalyzer>getDEIMErrorsOnTraj(h) {
r = getappdata(h.main," r ");
s = getappdata(h.main," s ");
mu = getappdata(h.main," mu ");
[ef,er,fxno] = testing.DEIM.getApproxErrorFullRed(r, s.xr, s.t, mu, r.V);
pm = PlotManager;
pm.LeaveOpen= true;
ax = pm.nextPlot(" deimerr "," True and projected approximation error on current trajectory "," time "," error ");
if max(ef)/min(ef) < 100 || max(er)/min(er) < 100
    plotfun = @plot;
else
    plotfun = @semilogy;
end
plotfun(ax,r.Times,ef," b ",r.Times,er," r ");
legend(" Full "," Projected ");
pm.done;

/*  --- Outputs from this function are returned to the command line. */
}

function varargout = DEIMEstimatorAnalyzer>DEIMEstimatorAnalyzer_OutputFcn(hObject,eventdata,handles) {

/*  Get default command line output from handles structure */
varargout[1] = handles.output;

}

function  DEIMEstimatorAnalyzer>rbgAlpha_SelectionChanged(h,source,unused1) {
rb = get(source," SelectedObject ");
r = getappdata(h.main," r ");
e = r.ErrorEstimator;
switch get(rb," UserData ")
    case 1
        e.UseTrueDEIMErr= true;
        e.Enabled= true;
    case 2
        e.UseTrueDEIMErr= false;
end
reSimulate(h);

}

function  DEIMEstimatorAnalyzer>rbgBeta_SelectionChanged(h,source,eventdata) {
rb = get(source," SelectedObject ");
r = getappdata(h.main," r ");
e = r.ErrorEstimator;
e.UseTrueLogLipConst= false;
e.UseJacobianLogLipConst= false;
e.UseFullJacobian= false;
e.UseJacobianNorm= false;
switch get(rb," UserData ")
    case 1
        e.UseTrueLogLipConst= true;
    case 2
        e.UseJacobianLogLipConst= true;
    case 3
        e.UseFullJacobian= true;
    case 4
        e.UseJacobianNorm= true;
end
reSimulate(h);

}

function  DEIMEstimatorAnalyzer>rbgParam_SelectionChanged(h,source,unused1) {
rb = get(source," SelectedObject ");
switch get(rb," UserData ")
    case 1
        slTrainParam_Callback(h.slTrainParam, [], h);
    case 2
        updateUserParam(h);
end

/*  --- Executes on slider movement. */
}

function  DEIMEstimatorAnalyzer>slMD_Callback(hObject,eventdata,handles) {
r = getappdata(handles.main," r ");
md = round(get(hObject," Value "));
if r.System.f.Order(2) ~= md
    fprintf(" Setting m "" =%d...  ",md);
    r.System.f.Order= [r.System.f.Order(1) md];
    r.ErrorEstimator.Enabled= md > 0;
    set(handles.lblMD," String ",sprintf(" %d ",md));
    set(handles.rbMD," Value ",1);
    reSimulate(handles);
end

/*  --- Executes during object creation, after setting all properties. */
}

function  DEIMEstimatorAnalyzer>slMD_CreateFcn(hObject,eventdata,handles) {
if isequal(get(hObject," BackgroundColor "), get(0," defaultUicontrolBackgroundColor "))
    set(hObject," BackgroundColor ",[.9 .9 .9]);
end


/*  --- Executes on slider movement. */
}

function  DEIMEstimatorAnalyzer>slJM_Callback(hObject,eventdata,handles) {
r = getappdata(handles.main," r ");
jm = round(get(hObject," Value "));
if r.ErrorEstimator.JacMatDEIMOrder ~= jm
    fprintf(" Setting m_j=%d...  ",jm);
    r.ErrorEstimator.JacMatDEIMOrder= jm;
    set(handles.lblJM," String ",sprintf(" %d ",jm));
    reSimulate(handles);
end

/*  --- Executes during object creation, after setting all properties. */
}

function  DEIMEstimatorAnalyzer>slJM_CreateFcn(hObject,eventdata,handles) {
if isequal(get(hObject," BackgroundColor "), get(0," defaultUicontrolBackgroundColor "))
    set(hObject," BackgroundColor ",[.9 .9 .9]);
end

/*  --- Executes on slider movement. */
}

function  DEIMEstimatorAnalyzer>slK_Callback(hObject,eventdata,handles) {
r = getappdata(handles.main," r ");
k = round(get(hObject," Value "));
if r.ErrorEstimator.JacSimTransSize ~= k
    fprintf(" Setting k=%d...  ",k);
    r.ErrorEstimator.JacSimTransSize= k;
    set(handles.lblK," String ",sprintf(" %d ",k));
    reSimulate(handles);
end

/*  --- Executes during object creation, after setting all properties. */
}

function  DEIMEstimatorAnalyzer>slK_CreateFcn(hObject,eventdata,handles) {
if isequal(get(hObject," BackgroundColor "), get(0," defaultUicontrolBackgroundColor "))
    set(hObject," BackgroundColor ",[.9 .9 .9]);
end


/*  --- Executes on slider movement. */
}

function  DEIMEstimatorAnalyzer>slM_Callback(hObject,eventdata,handles) {
r = getappdata(handles.main," r ");
m = round(get(hObject," Value "));
if r.System.f.Order(1) ~= m
    fprintf(" Setting m=%d...  ",m);
    r.System.f.Order= [m min(r.System.f.Order(2),r.System.f.MaxOrder-m)];
    set(handles.lblM," String ",sprintf(" %d ",m));

    /*  Update range for m' values */
    mv = r.System.f.MaxOrder-m;
    if mv ~= 0, en=" on "; else en=" off "; end
    set(handles.slMD," Min ",0," Max ",max(mv,1)," Enable ",en," Value ", r.System.f.Order(2));
    set(handles.lblMD," String ",sprintf(" %d ",r.System.f.Order(2)));
    reSimulate(handles);
end

/*  --- Executes during object creation, after setting all properties. */
}

function  DEIMEstimatorAnalyzer>slM_CreateFcn(hObject,eventdata,handles) {
if isequal(get(hObject," BackgroundColor "), get(0," defaultUicontrolBackgroundColor "))
    set(hObject," BackgroundColor ",[.9 .9 .9]);
end


/*  --- Executes on button press in btnDEIMErr. */
}

function  DEIMEstimatorAnalyzer>btnDEIMErr_Callback(hObject,eventdata,handles) {
getDEIMErrorsOnTraj(handles);


/*  --- Executes on button press in btnRedSummary. */
}

function  DEIMEstimatorAnalyzer>btnRedSummary_Callback(hObject,eventdata,h) {
r = getappdata(h.main," r ");
ma = ModelAnalyzer(r);
pm = ma.plotReductionOverview;
pm.setFigureNames(sprintf(" Reduction summary of '%s' ",get(h.main," Name ")));


/*  --- Executes on slider movement. */
}

function  DEIMEstimatorAnalyzer>slTrainParam_Callback(hObject,unused1,h) {

/*  Hints: get(hObject,'Value') returns position of slider
 *        get(hObject,'Min') and get(hObject,'Max') to determine range of slider */
r = getappdata(h.main," r ");
pidx = round(get(hObject," Value "));
mu = r.ParamSamples(:,pidx);
if ~isequal(getappdata(h.main," mu "),mu)
    setappdata(h.main," mu ",mu);
    set(h.rbTrainParam," Value ",1);
    pdata = getappdata(h.main," pdata ");
    for pidx = 1:length(mu)
        set(pdata.sl(pidx)," Value ",mu(pidx));
        set(pdata.lbl(pidx)," String ",sprintf(" %g ",mu(pidx)));
    end
    drawnow;
    reSimulate(h);
end

/*  --- Executes during object creation, after setting all properties. */
}

function  DEIMEstimatorAnalyzer>slTrainParam_CreateFcn(hObject,eventdata,handles) {

/*  Hint: slider controls usually have a light gray background. */
if isequal(get(hObject," BackgroundColor "), get(0," defaultUicontrolBackgroundColor "))
    set(hObject," BackgroundColor ",[.9 .9 .9]);
end


/*  --- Executes on button press in btnVisResult. */
}

function  DEIMEstimatorAnalyzer>btnVisResult_Callback(hObject,eventdata,h) {
r = getappdata(h.main," r ");
m = r.FullModel;
s = getappdata(h.main," s ");
pm = [];
if get(h.rbOutput," Value ")
    vr = s.yr;
    v = s.y;
    pm = PlotManager(false,1,2);
    pm.LeaveOpen= true;
    pfun = @(t,y)r.plot(t,y,pm);
    pfunm = @(t,y)m.plot(t,y,pm);
else
    vr = r.V*s.xr;
    v = s.x;
    pfun = @r.plotState;
    pfunm = @r.plotState;
end
if get(h.chkAbsErr," Value ")
    pfun(s.t,vr);
    set(gcf," Name "," Results of reduced simulation ");
    pfunm(s.t,abs(v-vr));
    set(gcf," Name "," Absolute errors ");
else
    pfunm(s.t,v);
    set(gcf," Name "," Results of full simulation ");
    pfun(s.t,vr);
    set(gcf," Name "," Results of reduced simulation ");
end
if ~isempty(pm)
    pm.done;
end


/*  --- Executes on button press in pushbutton5. */
}

function  DEIMEstimatorAnalyzer>pushbutton5_Callback(hObject,eventdata,h) {
r = getappdata(h.main," r ");
s = getappdata(h.main," s ");
v = Norm.L2(s.x-r.V*(r.V^t*s.x));
pm = PlotManager;
pm.LeaveOpen= true;
h = pm.nextPlot(," Trajectory subspace projection L2(state) error "," time "," error ");
LogPlot.cleverPlot(h,s.t,v," r "," LineWidth ",2);
pm.done;


/*  --- Executes on button press in btnApproxLogNorm. */
}

function  DEIMEstimatorAnalyzer>btnApproxLogNorm_Callback(hObject,eventdata,h) {
r = getappdata(h.main," r ");
/* m = r.FullModel; */
s = getappdata(h.main," s ");
x = s.xr;
use = 20;
sel = round(linspace(1,size(x,2),min(use,size(x,2))));
x = x(:,sel);
mui = r.getRandomParam(400,1);
res = testing.LogNorm.getApproxLogNormsAtPos(r, x, r.scaledTimes(sel),mui);
pm = PlotManager;
pm.LeaveOpen= true;
testing.LogNorm.getApproxLogNormsAtPos_plots(res, pm);
pm.done;


/*  --- Executes on button press in rbState. */
}

function  DEIMEstimatorAnalyzer>rbState_Callback(hObject,eventdata,h) {
if get(hObject," Value ")
    set(h.rbOutput," Value ",0);
end
rePlot(h);

/*  --- Executes on button press in rbOutput. */
}

function  DEIMEstimatorAnalyzer>rbOutput_Callback(hObject,eventdata,h) {
if get(hObject," Value ")
    set(h.rbState," Value ",0);
end
rePlot(h);

}