LinearSplitOfOne.m

Go to the documentation of this file.

/* (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 LinearSplitOfOne
  :public handle {
    public:

        char VarName = "x";
        integer Splits = 1;
        double IntervalLength = 1;
    public:


        LinearSplitOfOne(char varname,integer splits,double len) {
            if nargin > 0
                this.VarName= varname;
                if nargin > 1
                    this.Splits= splits;
                    if nargin > 2
                        this.IntervalLength= len;
                    end
                end
            end
        }
        function str = getFunStr(idx) {
            if idx > this.Splits || idx < 0
                error(" Max %d functions (starting at 0) for %d splits ",this.Splits+1,this.Splits);
            end
            ilen = this.IntervalLength / this.Splits;
            if idx == 0
                str = sprintf(" (1-@@varname@@/%g).*(@@varname@@<%g) ",ilen,ilen);
            elseif idx == this.Splits
                str = sprintf(" ((@@varname@@-%g)/%g).*(@@varname@@>=%g) ",(idx-1)*ilen,ilen,(idx-1)*ilen);
            else
                /*  zero to one */
                zto = sprintf(" (@@varname@@>=%g).*(@@varname@@<%g).*((@@varname@@-%g)/%g) ",...
                    (idx-1)*ilen,idx*ilen,(idx-1)*ilen,ilen);
                /*  one to zero */
                otz = sprintf(" (@@varname@@>=%g).*(@@varname@@<%g).*(1-(@@varname@@-%g)/%g) ",...
                    idx*ilen,(idx+1)*ilen,idx*ilen,ilen);
                str = [zto " + " otz];
            end
            str = strrep(str," @@varname@@ ",this.VarName);
        }
        function cellfuns = getAllFunctions() {
            funs = cell(this.Splits+1,1);
            for idx = 1:this.Splits+1
                funs[idx] = this.getFunStr(idx-1);
            end
        }
        function handle = getAllFunctionsParsed() {
            funs = this.getAllFunctions;
            str = [" @( " this.VarName " )[ "];
            for idx = 1:this.Splits+1
                str = [str funs[idx] " ; "];/* #ok */

            end
            str = [str " ]; "];
            handle = eval(str);
        }
#if 0 //mtoc++: 'set.Splits'
function  Splits(value) {
            if value < 1
                error(" Minimum of one split necessary. ");
            end
            this.Splits= value;
        }

#endif

    
    public: /* ( Static ) */

        static function res = test_LinearSplitOfOne() {
            num = 20;
            s = randi(50,1,num);
            l = rand(1,num);
            res = true;
            for k=1:num
                lso = LinearSplitOfOne(" x ",s(k),l(k));
                f = lso.getAllFunctionsParsed;

                pos = linspace(0,l(k),200);
                for i=1:length(pos)
                    if i < 200
                        res = res && sum(f(pos(i))) == 1;
                    else
                        fprintf(" Splits: %d, len: %g, diff at end: %g\n ",s(k),l(k),sum(f(pos(i)))-1);
                    end
                end
            end
        }

    
};