Constraint.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 Constraint
  :public dscomponents.ACoreFun {
    public:

        ComputeUnassembled = false;

        
    private:

        fsys;
    public:

        idx_p_elems_unass;

        Sigma;

        fDim_unass;

        
    public:

        Constraint(sys) {
            this = this@dscomponents.ACoreFun(sys);
        }


        function  configUpdated() {
            sys = this.fsys;
            mc = sys.Model.Config;
            if ~isempty(mc)
                this.fDim= sys.NumAlgebraicDofs;
                this.xDim= sys.NumAlgebraicDofs;
                this.JSparsityPattern= this.computeSparsityPattern;
                this.precomputeUnassembledData;
            end
        }


        function  setSystem(sys) {
            setSystem@dscomponents.ACoreFun(this, sys);
            if isa(sys," models.ReducedSecondOrderSystem ")
                this.fsys= sys.Model.FullModel.System;
            else
                this.fsys= sys;
            end
        }


        function  evaluateCoreFun() {
            error(" Do not call. evaluate is implemented directly. ");
        }


        function projected = project(V,W) {
            projected = this.clone;
            projected = project@dscomponents.ACoreFun(this, V, W, projected);
        }


        function copy = clone() {
            copy = muscle.ConstraintsFun(this.System);
            copy = clone@dscomponents.ACoreFun(this, copy);
        }
    private:

        function  precomputeUnassembledData() {
            sys = this.System;
            mc = sys.Model.Config;
            fem = mc.PressureFEM;
            geo = fem.Geometry;

            [I, ~] = find(fem.Sigma);

            n = numel(I);
            S = sparse(I,1:n,ones(n,1),geo.NumNodes,n);

            /*  Take out nodes with dirichlet BC on output side
 *             S(sys.idx_v_bc_local,:) = [];
             * Find corresponding unassembled dofs that would be ignored
             * (due to dirichlet velocity values, pressure dirichlet not
             * implemented)
 *             bc_unass = find(sum(S,1) == 0);
             * Remove them, too. The unassembled evaluation also removes the
             * corresponding entries of the unassembled vector.
 *             S(:,bc_unass) = []; */

            this.Sigma= S;

            this.fDim_unass= size(S,2);

            /*  Create boolean array that determines which unassembled dofs
             * belong to which element
             * dg dofs */
            hlp = repmat(1:geo.NumElements,geo.DofsPerElement,1);
            hlp = hlp(:);
            ass = false(geo.NumElements,length(hlp));
            for k = 1:geo.NumElements
                ass(k,:) = hlp == k;
            end
            this.idx_p_elems_unass= ass;
        }

    

};
}
}