DefaultEstimator.m

Go to the documentation of this file.

namespace error{


/* (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.
 */

class DefaultEstimator
  :public error.BaseEstimator {
    public:

        RealOutputError;
    public:

        DefaultEstimator(models.ReducedModel rmodel) {

            this.ExtraODEDims= 0;
            if nargin == 1
                this.Enabled= true;
                this.offlineComputations(rmodel.FullModel);
                this = this.prepareForReducedModel(rmodel);
            else
                this.Enabled= false;
            end
        }
        function copy = clone() {
            copy = error.DefaultEstimator;
            copy = clone@error.BaseEstimator(this, copy);
            /*  No local properties there. */
        }
        function eint = evalODEPart() {
            eint = [];
        }
        function e0 = getE0(colvec<double> mu) {
            e0 = [];
        }
        function ct = postProcess(colvec<double> x,double t,integer inputidx) {

            m = this.ReducedModel.FullModel;

            /*  Compute full solution */
            [~, yf, ct, xf] = m.simulate(this.mu, inputidx);
            st = tic;
            xr = x(1:end-this.ExtraODEDims,:);
            if ~isempty(this.ReducedModel.V)
                diff = xf-this.ReducedModel.V*xr;
            else
                diff = xf-xr;
            end
            this.StateError= Norm.LG(diff,this.ReducedModel.G);

            /*  The reduced model's output C contains scaling! */
            yr = this.ReducedModel.System.computeOutput(xr);
            this.RealOutputError= Norm.L2(yf-yr);
            ct = ct + toc(st);
            ct = ct + postProcess@error.BaseEstimator(this, x, t, inputidx);
        }
        function ct = prepareConstants(colvec<double> mu,integer inputidx) {

            /*  Nothing to do here */
            ct = prepareConstants@error.BaseEstimator(this, mu, inputidx);
        }
    public: /* ( Static ) */

        static function errmsg = validModelForEstimator(varargin) {
            errmsg = [];
        }


};
}