DEIM.m

Go to the documentation of this file.

namespace approx{


/* (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 DEIM
  :public approx.BaseApprox,
   public general.DEIM {
    public:

        DEIM(sys) {
            this = this@general.DEIM;
            this = this@approx.BaseApprox(sys);

            this.CustomProjection= true;
            this.TimeDependent= true;
        }


        function  approximateSystemFunction(models.BaseFullModel model) {
            this.computeDEIM(model.System.f, model.Data.ApproxTrainData.fxi);
        }


        function fx = evaluate(colvec<double> x,double t) {
            fx = evaluate@general.DEIM(this, x, t);
        }


        function fx = evaluateMulti(colvec<double> x,double t,colvec<double> mu) {
            fx = evaluateMulti@general.DEIM(this, x, t, mu);
        }


        function  prepareSimulation(colvec<double> mu) {
            prepareSimulation@approx.BaseApprox(this, mu);
            /*  Forward the parameter setting to the inner ACompEvalCoreFun
             * (if already set) */
            if ~isempty(this.f)
                this.f.prepareSimulation(mu);
            end
        }


        function J = getStateJacobian(colvec<double> x,double t) {
            J = getStateJacobian@general.DEIM(this, x, t);
            /*  Finite difference one:
             *J2 = getStateJacobian@dscomponents.ACoreFun(this, x, t, mu); */
        }


        function fx = evaluateCoreFun() {
            error(" Do not call me! ");
        }
        function projected = project(V,W) {
            projected = this.clone;
            projected = project@approx.BaseApprox(this, V, W, projected);
            projected = project@general.DEIM(this, V, W, projected);
        }


        function copy = clone() {
            copy = approx.DEIM(this.System);
            copy = clone@approx.BaseApprox(this, copy);
            copy = clone@general.DEIM(this, copy);
        }

    
    protected: /* ( Static ) */

        static function obj = loadobj(obj,varargin) {
            obj = loadobj@general.DEIM(obj, varargin[:]);
            obj = loadobj@approx.BaseApprox(obj, varargin[:]);
        }

    
    public: /* ( Static ) */

        static function res = test_DEIM() {
            m = models.pcd.PCDModel(1);
            m.EnableTrajectoryCaching= false;
            m.Approx= approx.DEIM(m.System);
            m.Approx.MaxOrder= 40;
            m.System.Params(1).Desired = 10;
            m.SpaceReducer= spacereduction.PODGreedy;
            m.offlineGenerations;
            m.Approx.Order= [20 4];
            r = m.buildReducedModel;
            mu = m.getRandomParam;
            r.simulate(mu);
            r.System.f.Order= [15 2];
            r.simulate(mu);
            r.System.f.Order= 40;
            r.simulate(mu);
            res = true;
        }


};
}