Cube20Node.m

Go to the documentation of this file.

namespace fem{
namespace geometry{


/* (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 Cube20Node
  :public fem.geometry.BaseGeometry {
    public:


        Cube20Node(pts,cubes) {
            if nargin < 2
                this = fem.geometry.RegularHex20Grid;
                return;
            elseif size(unique(pts" , "rows^t),1) ~= size(pts,2);
                error(" Please provide unique points! ");
            end
            this.Nodes= pts;
            this.Elements= cubes;
            this.DofsPerElement= 20;
            /* EdgeIndices = [1 3 6 8 13 15 18 20]; */

            /* % Compute edges */
            e = int16.empty(0,2);
            for i=1:size(cubes,1)
                hlp = cubes(i,[1 2 1 4 1 9 2 3 3 5 3 10 4 6 5 8 6 7 ...
                    6 11 7 8 8 12 9 13 10 15 11 18 12 20 13 14 13 16 ...
                    14 15 15 17 16 18 17 20 18 19 19 20]);
                e(end+1:end+24,:) = reshape(hlp" ,2,[]) ";
            end
            e = unique(e," rows "," stable ");
            this.Edges= e;

            /* % Set Face indices */
            this.MasterFaces= [  1 4 6 9 11 13 16 18
                            3 5 8 10 12 15 17 20
                            1:3 9 10 13:15
                            6:8 11 12 18:20
                            1:8
                            13:20];
             this.PatchFacesIdx= [ 1 4 6 11 18 16 13 9
                                    3 5 8 12 20 17 15 10
                                    1 2 3 10 15 14 13 9
                                    6 7 8 12 20 19 18 11
                                    1 2 3 5 8 7 6 4
                                    13 14 15 17 20 19 18 16];
            this.PatchesPerFace= 1;
            this.Faces= this.computeFaces;

            /*  Sanity checks */
            this.ReverseAxesIndices= [3 2 1 5 4 8 7 6 10 9 12 11 15 14 13 17 16 20 19 18
                                       6:8 4 5 1:3 11 12 9 10 18:20 16 17 13:15
                                       13:20 9:12 1:8];
            this.OrientationCheckIndices(:,:,1) = [1:3; 6:8; 13:15; 18:20];
            this.OrientationCheckIndices(:,:,2) = [1 4 6; 3 5 8; 13 16 18; 15 17 20];
            this.OrientationCheckIndices(:,:,3) = [1 9 13; 3 10 15; 6 11 18; 8 12 20];
            this.checkOrientation;
        }


        function cube8 = toCube8Node() {
            elems20 = this.Elements;
            nodes20 = this.Nodes;
            elems8 = elems20(:,[1 3 6 8 13 15 18 20]);
            usednodes = unique(elems8(:)," stable ");
            nodes8 = nodes20(:,usednodes);
            invidx(usednodes) = 1:length(usednodes);
            elems8 = invidx(elems8);
            cube8 = fem.geometry.Cube8Node(nodes8,elems8);
        }
};
}
}