Cube8Node.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 Cube8Node
  :public fem.geometry.BaseGeometry {
    public:


        Cube8Node(pts,cubes) {
            if nargin < 2
                this = fem.geometry.RegularHex8Grid;
                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= 8;

            /* % Compute edges */
            e = int16.empty(0,2);
            for i=1:size(cubes,1)
                hlp = cubes(i,[1 2 1 3 1 5 3 4 2 4 4 8 3 7 ...
                    8 7 5 7 6 2 6 5 6 8]);
                e(end+1:end+12,:) = reshape(hlp" ,2,[]) ";
            end
            e = unique(e," rows "," stable ");
            this.Edges= e;

            /* % Set Face indices */
            this.MasterFaces= [1 3 5 7
                            2 4 6 8
                            1 2 5 6
                            3 4 7 8
                            1 2 3 4
                            5 6 7 8];
            this.PatchFacesIdx= [1 3 7 5
                            2 4 8 6
                            1 2 6 5
                            3 4 8 7
                            1 2 4 3
                            5 6 8 7];
            this.PatchesPerFace= 1;
            this.Faces= this.computeFaces;

            /*  Sanity checks */
            this.ReverseAxesIndices= [2 1 4 3 6 5 8 7
                                       3 4 1 2 7 8 5 6
                                       5:8 1:4];
            this.OrientationCheckIndices(:,:,1) = [1 2; 3 4; 5 6; 7 8];
            this.OrientationCheckIndices(:,:,2) = [1 3; 2 4; 5 7; 6 8];
            this.OrientationCheckIndices(:,:,3) = [1 5; 2 6; 3 7; 4 8];
            this.checkOrientation;
        }


        function cube20 = toCube20Node() {
            elems8 = this.Elements;
            nodes8 = this.Nodes;
            nc = size(elems8,1);
            /*  Transformation matrix for corners to corner+edges locations */
            i = [1 2  2  3  4  4  5  5 6  7  7 8  9  9 10 10 11 11 12 12 13 14 14 15 16 16 17 17 18 19 19 20];  
            j = [1 1  2  2  1  3  2  4 3  3  4 4  1  5 2  6  3  7  4  8  5  5  6  6  5  7  6  8  7  7  8  8];
            s = [1 .5 .5 1 .5 .5 .5 .5 1 .5 .5 1 .5 .5 .5 .5 .5 .5 .5 .5 1  .5 .5 1  .5 .5 .5 .5 1  .5 .5 1];
            T = sparse(i,j,s,20,8);
            nodes20 = zeros(3,nc*20);

            /*  Iterate all cubes and collect nodes */
            for cidx = 1:nc
                pos = 20*(cidx-1)+1:20*cidx;
                cube = elems8(cidx,:);
                nodes20(:,pos) = nodes8(:,cube)*T^t;
            end
            [nodes20, ~, elems20] = unique(nodes20" , "rows" , "stable^t);
            cube20 = fem.geometry.Cube20Node(nodes20" ,reshape(elems20,20,[]) ");
        }
        function cube27 = toCube27Node() {
            elems8 = this.Elements;
            nodes8 = this.Nodes;
            nc = size(elems8,1);
            h = .5;
            v = .25;
            /*  Transformation matrix for corners to corner+edges locations
             * Nodes 1-13 */
            i = [1 2  2  3  4  4  5 5 5 5 6 6 7 8 8 9 10 10 11 11 11 11 12 12 13 13 13 13];
            j = [1 1  2  2  1  3  1 2 3 4 2 4 3 3 4 4 1  5  1  2  5  6  2  6  1  3  5  7]; 
            s = [1 h  h  1  h  h  v v v v h h 1 h h 1 .5 .5 v  v  v  v  .5 .5 v  v  v  v];
            /*  Node 14 */
            i = [i ones(1,8)*14];
            j = [j 1:8];
            s = [s ones(1,8)/8];
            /*  Nodes 15-27 */
            i = [i 15 15 15 15 16 16 17 17 17 17 18 18 19 20 20 21 22 22 23 23 23 23 24 24 25 26 26 27];
            j = [j 2  4  6  8  3  7  3  4  7  8  4  8  5  5  6  6  5  7  5  6  7  8  6  8  7  7  8  8]; 
            s = [s v  v  v  v  h  h  v  v  v  v  h  h  1  h  h  1  .5 .5 v  v  v  v  .5 .5 1  h  h  1];
            T = sparse(i,j,s,27,8);
            nodes27 = zeros(3,nc*27);

            /*  Iterate all cubes and collect nodes */
            for cidx = 1:nc
                pos = 27*(cidx-1)+1:27*cidx;
                cube = elems8(cidx,:);
                nodes27(:,pos) = nodes8(:,cube)*T^t;
            end
            [nodes27, ~, elems27] = unique(nodes27" , "rows" , "stable^t);
            cube27 = fem.geometry.Cube27Node(nodes27" ,reshape(elems27,27,[]) ");
        }
};
}
}