GLEstimator.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 GLEstimator
  :public error.BaseCompLemmaEstimator {
    private:

        cf;

        f;

        
    public:

        function copy = clone() {
            copy = error.GLEstimator;
            copy = clone@error.BaseCompLemmaEstimator(this, copy);
            copy.cf= this.cf;
            copy.f= this.f;
        }
        function b = getBeta(colvec<double> x,double t) {
            if ~isempty(this.cf)
                b = this.cf;
            else
                b = this.f.getGlobalLipschitz(t, this.mu);
            end
        }


        function ct = prepareConstants(colvec<double> mu,integer inputidx) {

            ct = prepareConstants@error.BaseCompLemmaEstimator(this, mu, inputidx);
            /*  Standard case: the approx function is a kernel expansion. it
             * can also be that the system's core function is already a
             * kernel expansion */
            st = tic;
            fm = this.ReducedModel.FullModel;
            if ~isempty(fm.Approx)
                /*  Get full d x N coeff matrix of approx function */
                this.f= fm.Approx;
            else
                /*  Get full d x N coeff matrix of core function */
                this.f= fm.System.f;
            end
            /*  Only precompute the global lipschitz constant if it is (t,\mu)-independent */
            if isa(this.f," kernels.KernelExpansion ") || (isa(this.f," kernels.ParamTimeKernelExpansion ") && ...
                    isa(this.f.ParamKernel," kernels.NoKernel ") && isa(this.f.TimeKernel," kernels.NoKernel "))
                this.cf= this.f.getGlobalLipschitz(0, mu);
            else
                this.cf= [];
            end
            ct = ct + toc(st);
        }
        function ct = postProcess(colvec<double> x,double t,integer inputidx) {

            st = tic;
            this.StateError(1,:) = x(end,:);
            ct = postProcess@error.BaseCompLemmaEstimator(this, x, t, inputidx) + toc(st);
        }
    public: /* ( Static ) */

        static function errmsg = validModelForEstimator(models.BaseFullModel model) {
            errmsg = validModelForEstimator@error.BaseCompLemmaEstimator(model);
            if isempty(errmsg) && ~isa(model.System.f.Expansion," dscomponents.IGlobalLipschitz ")
                errmsg = " The full model "" s core function must implement the dscomponents.IGlobalLipschitz interface for this error estimator. ";
            end
        }
};
}