FuncSum.m

Go to the documentation of this file.

namespace general{
namespace functions{


/* (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 FuncSum
  :public general.functions.AFunGen {
    private:

        fungens;

        
    public:

        FuncSum(varargin) {
            if length(varargin) < 2
                error(" Must have at least two summands! ");
            end
            this.fungens= varargin;
        }


        function [fhandle , dfhandle ] = getFunction() {
            f = this.fungens;
            names = cell(1,length(f));
            for k = 1:length(f)
                mc = metaclass(f[k]);
                [~,name] = fileparts(which(mc.Name));
                names[k] = sprintf(" %s%d ",name,k);
                eval(sprintf(" %s = f{%d}.getFunction; ",names[k],k));
            end
            /*  Add (t) argument */
            names = cellfun(@(s)sprintf(" %s(t) ",s),names," UniformOutput ",false);
            /*  Join with + */
            sumstr = Utils.implode(names,"  +  ");
            /*  Eval to readable handle */
            eval(sprintf(" fhandle = @(t)%s; ",sumstr));
            dfhandle = [];
        }


        function str = getConfigStr() {
            str = this.fungens[1].getConfigStr;
            for k = 2:length(this.fungens)
                str = sprintf(" / %s ",this.fungens[k].getConfigStr);
            end
        }

    
    public: /* ( Static ) */

        static function res = test_FuncSum() {
            f1 = general.functions.PiecewiseLinear(-1:.3:1);
            f1h = f1.getFunction;
            w = rand(3,1)*5;
            f2 = general.functions.Sinus(w(1),w(2),w(3));
            f2h = f2.getFunction;
            w = rand(3,1)*5;
            f3 = general.functions.Sinus(w(1),w(2),w(3));
            f3h = f3.getFunction;

            fs = general.functions.FuncSum(f1,f2,f3);
            fsh = fs.getFunction;

            res = true;
            for x = -1:.1:1
                res = res && f1h(x)+f2h(x)+f3h(x) == fsh(x);
            end
            fs.plot;
        }

    

};
}
}


/* /*  Speed te*/
 Speed test
 * all(1:3) = 0;
 * for k=1:1000
 *     slen = min(len)*.99 + rand*(max(len)*1.01-min(len)*.99);
 *     t = tic;
 *     %sproc = spline(len,proc,slen);
 *     sproc = pchip(len,proc,slen);
 *     all(1) = all(1) + toc(t);
 *     l = rand*3;
 *     t = tic;
 *     fl(l);
 *     all(2) = all(2) + toc(t);
 *     t = tic;
 *     polyval(pol,slen);
 *     all(3) = all(3) + toc(t);
 * end
 * all / 1000