MarkertLawOriginal.m

Go to the documentation of this file.

namespace models{
namespace muscle{
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 MarkertLawOriginal
  :public general.functions.AFunGen {
    public:

        b;

        d;

        
    public:

        MarkertLawOriginal(b,d) {
            if nargin < 2
                d = 16.6;
                if nargin < 1
                    b = 7990;
                end
            end
            this.b= b;
            this.d= d;
            this.xLabel= " \lambda [-] ";
            this.yLabel= " Pressure [MPa] ";
        }


        function [fhandle , dfhandle , fbdhandle , dfbdhandle ] = getFunction() {
            b1 = this.b;
            d1 = this.d;            

            /* % Regular function
 *             fhandle = @(t)b1./t.^2.*(t.^d1-1);
 *             dfhandle = @(t)b1./t.^3.*((d1-2).*t.^d1+2);
 *             fbdhandle = @(t,b1,d1)b1./t.^2.*(t.^d1-1);
 *             dfbdhandle = @(t,b1,d1)b1./t.^3.*((d1-2).*t.^d1+2); */

            /* % Variant with support [1, \infty] */
            fhandle = @(t)max(0,b1./t.^2.*(t.^d1-1));
            dfhandle = @(t)(t>1).*(b1./t.^3.*((d1-2).*t.^d1+2));
            fbdhandle = @(t,b1,d1)max(0,b1./t.^2.*(t.^d1-1));
            dfbdhandle = @(t,b1,d1)max(0,b1./t.^3.*((d1-2).*t.^d1+2));

            /* % Variant without dlam/dC term (1/2\lambda)
 *             fhandle = @(t)b1./t.*(t.^d1-1);
 *             dfhandle = @(t)b1./t.^2.*((d1-2).*t.^d1+2);
 *             fbdhandle = @(t,b1,d1)b1./t.*(t.^d1-1);
 *             dfbdhandle = @(t,b1,d1)b1./t.^2.*((d1-2).*t.^d1+2); */
        }


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


        function pm = plot(range,varargin) {
            if nargin < 2
                range = [.9 1.2];
            end
            pm = plot@general.functions.AFunGen(this, " R ", range, varargin[:]);
        }

    
    public: /* ( Static ) */

        static function  createScratchPlot() {
            range = [.9 1.4];
            so = general.functions.MarkertLawOriginal(7990,16.6);
            pm = so.plot(range," b ");
            ax = get(gcf," Children ");
            ax = flipud(ax);
            hold(ax(1)," on ");
            hold(ax(2)," on ");
            ylim(ax(1),[0 16500]);
            ylim(ax(2),[0 4.1e5]);
            s = general.functions.MarkertLaw(7990,16.6,1e9);
            s.plot(range,ax," r ");
            s = general.functions.MarkertLaw(1000*7990,4.5,1e9);
            s.plot(range,ax," g ");
            legend(ax(1)," Original b=7990,d=16.6 "," Modified b=7990,d=16.6 "," Modified b=7990000,d=4.5 ");/* ,'Location','NorthWest'); */

            pm.savePlots(" /home/dwirtz/software/muscle/latex/img "," Format "," eps ");
        }

    

};
}
}
}