Long.m

Go to the documentation of this file.

namespace models{
namespace muscle{
namespace examples{


/* (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 Long
  :public models.muscle.AMuscleConfig {
    public:

        Long(varargin) {
            this = this@models.muscle.AMuscleConfig(varargin[:]);
            this.addOption(" Devi ",1);
            this.init;
        }


        function  configureModel(m) {
            configureModel@models.muscle.AMuscleConfig(this, m);
            m.T= 50;
            m.dt= 1;
            m.DefaultMu(1) = .1;
            m.DefaultMu(2) = 20;
        }

    
    protected:


        function geo = getGeometry() {
            geo = fem.geometry.RegularHex20Grid(-10:10:10,-40:10:40, [0 10], this.Options.Devi);
        }
        function displ_dir = setPositionDirichletBC(displ_dir) {
            geo = this.FEM.Geometry;
            for k = [8 16]
                displ_dir(:,geo.Elements(k,[6:8 11 12 18:20])) = true;
            end
        }
        function anull = seta0(anull) {
            fe = this.FEM;
            if fe.GaussPointsPerElem ~= 27
                warning(" a0 designed for 27 gauss points! ");
            end
            x = linspace(0,1,8*3);
            basea0 = [sin(x*pi); cos(x*pi)]; /* ; zeros(size(x)) */

            front = fe.GaussPoints(2,:) < 0;
            mid = fe.GaussPoints(2,:) == 0;
            back = fe.GaussPoints(2,:) > 0;
            /*  Direction is x */
            for m = 1:8
                off = (m-1)*3;
                anull([1 3],front,[m m+8]) = basea0(2,off+1);
                anull([1 3],mid,[m m+8]) = basea0(2,off+2);
                anull([1 3],back,[m m+8]) = basea0(2,off+3);

                anull(2,front,[m m+8]) = basea0(1,off+1);
                anull(2,mid,[m m+8]) = basea0(1,off+2);
                anull(2,back,[m m+8]) = basea0(1,off+3);
            end
        }

    
    public: /* ( Static ) */

        static function  test_Long() {
            m = models.muscle.Model(models.muscle.examples.Long);
            m.simulateAndPlot;
        }

    

};
}
}
}