EstimatorParamSweep.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 [Y , E ] = EstimatorParamSweep(ea,mu,unused1,param,rowvec<double> paramvals) {


/*  Validity checks */
if ~isa(ea," EstimatorAnalyzer ")
    error(" The first argument must be a EstimatorAnalyzer instance. ");
end
if ischar(param)
    pidx = ea.ReducedModel.getParamIndexFromName(param);
    pname = param;
elseif isposintscalar(param)
    pidx = param;
    pname = ea.ReducedModel.System.Params(pidx).Name;
else
    error(" param must be either a char array (parameter name) or positive parameter index ");
end
npar = length(paramvals);
if npar == 0
    error(" Parameter values to sweep must be given. ");
end

/*  % Iterate over parameter values
 * mu(pidx) = paramvals(1);
 * [errs, ctimes, relerrs] = ea.getErrorEstimates(mu, inputidx, true);
 * fprintf('Estimator analyzer parameter sweep 1/%d with %s=%f\n',npar,pname,mu(pidx));
 * for idx = 2:npar
 *     mu(pidx) = paramvals(idx);
 *     [errs(:,:,idx), ctimes(idx,:), relerrs(:,:,idx)] = ea.getErrorEstimates(mu, inputidx, true);
 * 
 *     fprintf('Estimator analyzer parameter sweep %d/%d with %s=%f\n',idx,npar,pname,mu(pidx));
 * end */
load tmp;

/* % Prepare plot */
[T, MU] = meshgrid(ea.Model.Times,paramvals);
tit = sprintf(" Error estimator analysis for parameter sweep of '%s' from %1.3f to %1.3f, base \\mu = [%s] ",...
    pname,paramvals(1),paramvals(end),num2str(mu^t));

/*  %% Upper bound
 * if rmodel.ErrorEstimator.Enabled
 *     obj = mesh(T,MU,Y+E,'EdgeColor','none','FaceColor','red');
 *     alpha(obj,.3);
 * %     mesh(T,MU,Y+E,'EdgeColor','none','FaceColor','red');
 *     title(tit); axis tight; xlabel('t'); ylabel(pname);
 * end */

/* % y plot */
p = figure;
ax = gca(p);
rotate3d(ax," on ");
view(ax,9,16);
inv = figure(" Visible "," off ");
iax = gca(inv);
hold on;
/*  cm = jet; */
cd = cell.empty;
for idx = 1:size(errs,1)/* #ok */

    Z = squeeze(errs(idx,:,:))^t;
    /* surf(ax,T,MU,Z,'EdgeColor','none','FaceColor',cm(idx*floor(size(cm,1)/size(errs,1)),:)); */
    h = surf(iax,T,MU,Z," EdgeColor "," none ");
    cd[idx] = get(h," CData ");
    set(h," Parent ",ax);
end
title(tit);
emin = min(errs(:));
/* axis(ax,[0 ea.Model.T 1 npar emin*.9 max(emin,1e4)]); */
xlabel(" Time t "); ylabel(sprintf(" Parameter %s value ",pname)); zlabel(" Error estimates ");
zlim(ax,[emin*.9 max(emin,1e4)]);

cm = jet;
colormap(ax,cm);
for idx = 1:size(errs,1)
    ncd = cd[idx];
    /* ncd = ncd-min(ncd(:)); */
    ncd = size(cm,1)*max(ncd(:))./ncd;
    set(h," CData ",ncd," CDataMapping "," direct ");
end
close(inv);

/*  %% Lower bound
 * if rmodel.ErrorEstimator.Enabled
 *     figure;
 *     colormap jet;
 *     obj = surf(T,MU,Y-E,'EdgeColor','none','FaceColor','red');
 *     alpha(obj,.3);
 * %     mesh(T,MU,Y-E,'EdgeColor','none','FaceColor','red');
 *     hold on;
 * end */

}