SubMeshModelConfig.m

Go to the documentation of this file.

namespace models{
namespace muscle{


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


        
    public:

        FullConfig;


        SubToFullNodeIdx;

        SubToFullElemIdx;

        SubToFullFaceIdx;

        
    public:

        SubMeshModelConfig(full,elems,varargin) {
            geo = full.FEM.Geometry;
            [subgeo, nodeidx, faceidx] = geo.getSubMesh(elems, varargin[:]);
            this = this@models.muscle.AMuscleConfig(subgeo);

            /*  Store original config */
            this.FullConfig= full;
            /*  Store index of elems/nodes in original geometry */
            this.SubToFullNodeIdx= nodeidx;
            this.SubToFullElemIdx= elems;
            this.SubToFullFaceIdx= faceidx;

            /*  Copy other quantities */
            this.NeumannCoordinateSystem= full.NeumannCoordinateSystem;
            this.a0CoordinateSystem= full.a0CoordinateSystem;
            this.FibreTypeWeights= full.FibreTypeWeights;
            this.Pool= full.Pool;
        }


        function  configureModel(models.BaseFullModel model) {
            this.FullConfig.configureModel(model);
        }


        function  prepareSimulation(colvec<double> mu,integer inputidx) {
            this.FullConfig.prepareSimulation(mu, inputidx);
        }


        function P = getBoundaryPressure(elemidx,faceidx) {
            fullelemidx = this.SubToFullElemIdx(elemidx);
            fullfaceidx = this.SubToFullFaceIdx(faceidx);
            P = this.FullConfig.getBoundaryPressure(fullelemidx, fullfaceidx);
        }


        function u = getInputs() {
            u = this.FullConfig.getInputs;
        }


        function x0 = getX0(x0) {
            x0 = this.FullConfig.getX0(x0);
        }
    protected:

        function anull = seta0(unused1) {
            fullanull = this.FullConfig.geta0;
            anull = fullanull(:,:,this.SubToFullElemIdx);
        }


        function [velo_dir , velo_dir_val ] = setVelocityDirichletBC(unused1,unused2) {
            [~, velo_dir, velo_dir_val] = this.FullConfig.getBC;
            velo_dir = velo_dir(:,this.SubToFullNodeIdx);
            velo_dir_val = velo_dir_val(:,this.SubToFullNodeIdx);
        }



        function displ_dir = setPositionDirichletBC(unused1) {
            displ_dir = this.FullConfig.getBC;
            displ_dir = displ_dir(:,this.SubToFullNodeIdx);
        }

    

};
}
}