PointerCoreFun.m

Go to the documentation of this file.

namespace dscomponents{


/* (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 PointerCoreFun
  :public dscomponents.ACoreFun {
    private:

        target;

        
    public:

        PointerCoreFun(function_handle funPtr,xdim,logical timedep) {

            /*  Create ACoreFun without associated system (cannot be used
             * here) */
            this = this@dscomponents.ACoreFun([]);
            if nargin < 3
                /*  Assume worst case: set time dependency to true! */
                timedep = true;
            end

            /*  Creates a new wrapper for a core function handle. */
            if ~isa(funPtr," function_handle ")
                error(" Argument funPtr must be a function handle. ");
            elseif nargin(funPtr) ~= 3
                error(" funPtr nargin must equal three (= x,t,mu). ");
            end
            this.CustomProjection= false;
            this.TimeDependent= timedep;
            this.target= funPtr;
            this.xDim= xdim;
            this.fDim= size(funPtr(rand(xdim,1),0,rand(100,1)),1);
        }
        function fx = evaluateCoreFun(colvec<double> x,double t) {
            fx = this.target(x, t, this.mu);
        }
        function proj = project(V,W) {

            /*  no need to call clone here as the only property gets changed
             * anyways. */
            newfun = @(z,t,mu)W^t * this.target(V*z,t,mu);
            proj = dscomponents.PointerCoreFun(newfun, size(V,2));
        }
    public: /* ( Sealed ) */

        function copy = clone() {
            copy = dscomponents.PointerCoreFun(this.target);
        }

    
    public: /* ( Static ) */

        static function  test_PointerCoreFun() {
            m = models.BaseFullModel;
            fun = @(x,t,mu)x*3.*t;
            m.System.x0= dscomponents.ConstInitialValue(1);
            m.System.f= dscomponents.PointerCoreFun(fun,1);
            m.simulate;
            clear m;
        }

    

};
}