EntireTA.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 EntireTA
  :public models.muscle.AMuscleConfig {
    public:

        EntireTA(varargin) {
            this = this@models.muscle.AMuscleConfig(varargin[:]);
            this.init;

            /* % Muscle fibre weights
 *             geo = s.geo27;
 *             types = [0 .2 .4 .6 .8 1];
 *             ftw = zeros(this.FEM.GaussPointsPerElem,length(types),geo.NumElements);
 *             % Test: Use only slow-twitch muscles
 *             ftw(:,1,:) = .4;
 *             ftw(:,2,:) = .05;
 *             ftw(:,3,:) = .05;
 *             ftw(:,4,:) = .1;
 *             ftw(:,5,:) = .2;
 *             ftw(:,6,:) = .2;
 *             this.FibreTypeWeights = ftw;
 *             p = models.motorunit.Pool;
 *             p.FibreTypes = types;
 *             this.Pool = p; */
        }
        function  configureModel(m) {
            configureModel@models.muscle.AMuscleConfig(this, m);
            /*  Overload this method to set model-specific quantities like
             * simulation time etc */
            m.T= 200;
            m.dt= 1;
            m.DefaultMu(1) = .1;
            m.DefaultMu(2) = 50;
/*              m.DefaultMu(4) = 4; */
            os = m.ODESolver;
            os.RelTol= .01;
            os.AbsTol= .08;
        }

    
    protected:


        function geo = getGeometry() {
            s = load(fullfile(fileparts(which(mfilename))," EntireTA.mat "));
            geo = s.geo27;
        }


        function displ_dir = setPositionDirichletBC(displ_dir) {
            geo = this.FEM.Geometry;
            /*  Front face
 *             displ_dir(:,geo.Elements(6,geo.MasterFaces(3,:))) = true;
             * Back faces */
            displ_dir(:,geo.Elements(8,geo.MasterFaces(4,:))) = true;
            displ_dir(:,geo.Elements(8,geo.MasterFaces(2,:))) = true;
        }
        function anull = seta0(anull) {
            anull(1,:,:) = 1;
        }
    public: /* ( Static ) */

        static function  test_EntireTA() {
            m = models.muscle.Model(models.muscle.examples.EntireTA);
            m.simulateAndPlot;
        }

    

};
}
}
}