CubicToLinear.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 CubicToLinear
  :public general.functions.AFunGen {
    public:

        M;

        lam0;

        
    public:

        CubicToLinear(lam0,M) {
            if nargin < 2
                M = 100;
                if nargin < 1
                    lam0 = 2;
                end
            end
            this.lam0= lam0;
            this.M= M;
        }


        function [fhandle , dfhandle ] = getFunction() {
            m = this.M;
            l = this.lam0;
            a = m/(l-1)^2; 
            fl = a*(l.^3/3-l.^2+l-1/3);
            fhandle = @(t)(t>=1 & t<l).*a.*(t.^3/3-t.^2+t-1/3) + (t>=l).*(fl + (t-l)*m);
            dfhandle = @(t)(t>=1 & t<l).*a.*(t.^2-2*t+1) + (t>=l)*m;
        }


        function str = getConfigStr() {
            str = sprintf(" lam0: %g, M: %g ",this.lam0,this.M);
        }


        function pm = plot(varargin) {
            varargin = [[" R ",[1 1.2*this.lam0]] varargin];
            
            i.addParamValue('R',0);
            i.parse(varargin[:]);
            res = i.Results;
            pm = plot@general.functions.AFunGen(this, varargin[:]);
            if (res.R(2) > this.lam0)
                mc = metaclass(this);
                for k = 1:length(pm.Figures)
                    if ishandle(pm.Figures(k)) && strcmpi(get(pm.Figures(k)," Tag "),...
                            strrep(mc.Name," . "," _ "))
                        f = this.getFunction;
                        ax = get(pm.Figures(k)," Children ");
                        hold(ax(1)," on ");
                        plot(ax(1),this.lam0,f(this.lam0)," rx "," MarkerSize ",16);
                        break;
                    end
                end
            end
        }

    
    public: /* ( Static ) */

        static function qfungen = fromLinearizedMarkert(mfungen) {
            mfun = mfungen.getFunction;
            m = mfungen.max_modulus;
            t0 = 2*mfungen.t0;
            ft0 = mfun(t0);
            l0 = 2*(t0-ft0/m)-1;
            qfungen = general.functions.CubicToLinear(l0,m);
        }


        static function res = test_CubicToLinear() {
            f = general.functions.CubicToLinear;
            f.plot;
            f = general.functions.CubicToLinear(4,5e6);
            f.plot;
            res = true;
        }

    

};
}
}