07-Jan-2022 20:05:02 fem_basis_test(): MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2 Test fem_basis(). FEM_BASIS_TEST01 FEM_BASIS_1D evaluates an arbitrary basis function over an interval. I J X L(I,J)(X) 2 1 0.6667 1 0 3 0.0000 -0 1 2 0.3333 0 2 1 0.6667 1 3 0 1.0000 -0 FEM_BASIS_TEST02 FEM_BASIS_2D evaluates an arbitrary triangular basis function. I J K X Y L(I,J,K)(X,Y) 1 0 2 0.3333 0.0000 1 0 0 3 0.0000 0.0000 0 1 0 2 0.3333 0.0000 1 2 0 1 0.6667 0.0000 -0 3 0 0 1.0000 0.0000 0 0 1 2 0.0000 0.3333 0 1 1 1 0.3333 0.3333 -0 2 1 0 0.6667 0.3333 0 0 2 1 0.0000 0.6667 -0 1 2 0 0.3333 0.6667 0 0 3 0 0.0000 1.0000 0 FEM_BASIS_TEST03 FEM_BASIS_3D evaluates an arbitrary tetrahedral basis function. I J K L X Y Z L(I,J,K,L)(X,Y,Z) 1 0 2 1 0.2500 0.0000 0.5000 1 0 0 0 4 0.0000 0.0000 0.0000 -0 1 0 0 3 0.2500 0.0000 0.0000 -0 2 0 0 2 0.5000 0.0000 0.0000 -0 3 0 0 1 0.7500 0.0000 0.0000 -0 4 0 0 0 1.0000 0.0000 0.0000 0 0 1 0 3 0.0000 0.2500 0.0000 -0 1 1 0 2 0.2500 0.2500 0.0000 -0 2 1 0 1 0.5000 0.2500 0.0000 -0 3 1 0 0 0.7500 0.2500 0.0000 0 0 2 0 2 0.0000 0.5000 0.0000 -0 1 2 0 1 0.2500 0.5000 0.0000 -0 2 2 0 0 0.5000 0.5000 0.0000 0 0 3 0 1 0.0000 0.7500 0.0000 -0 1 3 0 0 0.2500 0.7500 0.0000 0 0 4 0 0 0.0000 1.0000 0.0000 0 0 0 1 3 0.0000 0.0000 0.2500 0 1 0 1 2 0.2500 0.0000 0.2500 0 2 0 1 1 0.5000 0.0000 0.2500 0 3 0 1 0 0.7500 0.0000 0.2500 -0 0 1 1 2 0.0000 0.2500 0.2500 0 1 1 1 1 0.2500 0.2500 0.2500 0 2 1 1 0 0.5000 0.2500 0.2500 -0 0 2 1 1 0.0000 0.5000 0.2500 0 1 2 1 0 0.2500 0.5000 0.2500 -0 0 3 1 0 0.0000 0.7500 0.2500 -0 0 0 2 2 0.0000 0.0000 0.5000 0 1 0 2 1 0.2500 0.0000 0.5000 1 2 0 2 0 0.5000 0.0000 0.5000 -0 0 1 2 1 0.0000 0.2500 0.5000 0 1 1 2 0 0.2500 0.2500 0.5000 -0 0 2 2 0 0.0000 0.5000 0.5000 -0 0 0 3 1 0.0000 0.0000 0.7500 0 1 0 3 0 0.2500 0.0000 0.7500 -0 0 1 3 0 0.0000 0.2500 0.7500 -0 0 0 4 0 0.0000 0.0000 1.0000 -0 FEM_BASIS_TEST04 FEM_BASIS_MD evaluates an arbitrary basis function over an M-dimensional simplex. I J X L(I,J)(X) 2 1 0.6667 1 0 3 0.0000 -0 1 2 0.3333 0 2 1 0.6667 1 3 0 1.0000 0 FEM_BASIS_TEST05 FEM_BASIS_MD evaluates an arbitrary basis function over an M-dimensional simplex. I J K X Y L(I,J,K)(X,Y) 1 0 2 0.3333 0.0000 1 0 0 3 0.0000 0.0000 0 1 0 2 0.3333 0.0000 1 2 0 1 0.6667 0.0000 0 3 0 0 1.0000 0.0000 -0 0 1 2 0.0000 0.3333 0 1 1 1 0.3333 0.3333 0 2 1 0 0.6667 0.3333 -0 0 2 1 0.0000 0.6667 0 1 2 0 0.3333 0.6667 -0 0 3 0 0.0000 1.0000 -0 FEM_BASIS_TEST06 FEM_BASIS_MD evaluates an arbitrary basis function over an M-dimensional simplex. I J K L X Y Z L(I,J,K,L)(X,Y,Z) 1 0 2 1 0.2500 0.0000 0.5000 1 0 0 0 4 0.0000 0.0000 0.0000 -0 1 0 0 3 0.2500 0.0000 0.0000 -0 2 0 0 2 0.5000 0.0000 0.0000 -0 3 0 0 1 0.7500 0.0000 0.0000 -0 4 0 0 0 1.0000 0.0000 0.0000 -0 0 1 0 3 0.0000 0.2500 0.0000 -0 1 1 0 2 0.2500 0.2500 0.0000 -0 2 1 0 1 0.5000 0.2500 0.0000 -0 3 1 0 0 0.7500 0.2500 0.0000 -0 0 2 0 2 0.0000 0.5000 0.0000 -0 1 2 0 1 0.2500 0.5000 0.0000 -0 2 2 0 0 0.5000 0.5000 0.0000 -0 0 3 0 1 0.0000 0.7500 0.0000 -0 1 3 0 0 0.2500 0.7500 0.0000 -0 0 4 0 0 0.0000 1.0000 0.0000 -0 0 0 1 3 0.0000 0.0000 0.2500 0 1 0 1 2 0.2500 0.0000 0.2500 0 2 0 1 1 0.5000 0.0000 0.2500 0 3 0 1 0 0.7500 0.0000 0.2500 0 0 1 1 2 0.0000 0.2500 0.2500 0 1 1 1 1 0.2500 0.2500 0.2500 0 2 1 1 0 0.5000 0.2500 0.2500 0 0 2 1 1 0.0000 0.5000 0.2500 0 1 2 1 0 0.2500 0.5000 0.2500 0 0 3 1 0 0.0000 0.7500 0.2500 0 0 0 2 2 0.0000 0.0000 0.5000 0 1 0 2 1 0.2500 0.0000 0.5000 1 2 0 2 0 0.5000 0.0000 0.5000 0 0 1 2 1 0.0000 0.2500 0.5000 0 1 1 2 0 0.2500 0.2500 0.5000 0 0 2 2 0 0.0000 0.5000 0.5000 0 0 0 3 1 0.0000 0.0000 0.7500 0 1 0 3 0 0.2500 0.0000 0.7500 0 0 1 3 0 0.0000 0.2500 0.7500 0 0 0 4 0 0.0000 0.0000 1.0000 0 FEM_BASIS_TEST07 FEM_BASIS_PRISM_TRIANGLE evaluates an arbitrary basis function over a right triangular prism. Here, we generate basis functions which can be up to degree 2 in X and Y, and up to degree 2 in Z. Choose a node N1, define the basis function associated with that node, and then evaluate it at all other nodes. I1 I2 I3 J1 J2 X Y Z B(X,Y,Z) 2 0 0 1 1 1.000000 0.000000 0.500000 1 0 0 2 0 2 0.000000 0.000000 0.000000 -0 0 0 2 1 1 0.000000 0.000000 0.500000 -0 0 0 2 2 0 0.000000 0.000000 1.000000 0 0 1 1 0 2 0.000000 0.500000 0.000000 -0 0 1 1 1 1 0.000000 0.500000 0.500000 -0 0 1 1 2 0 0.000000 0.500000 1.000000 0 0 2 0 0 2 0.000000 1.000000 0.000000 -0 0 2 0 1 1 0.000000 1.000000 0.500000 -0 0 2 0 2 0 0.000000 1.000000 1.000000 0 1 0 1 0 2 0.500000 0.000000 0.000000 0 1 0 1 1 1 0.500000 0.000000 0.500000 0 1 0 1 2 0 0.500000 0.000000 1.000000 -0 1 1 0 0 2 0.500000 0.500000 0.000000 0 1 1 0 1 1 0.500000 0.500000 0.500000 0 1 1 0 2 0 0.500000 0.500000 1.000000 -0 2 0 0 0 2 1.000000 0.000000 0.000000 0 2 0 0 1 1 1.000000 0.000000 0.500000 1 2 0 0 2 0 1.000000 0.000000 1.000000 -0 FEM_BASIS_TEST08 FEM_BASIS_PRISM_TRIANGLE evaluates an arbitrary basis function over a right triangular prism. Here, we generate basis functions which can be up to degree 3 in X and Y, and up to degree 1 in Z. Choose a node N1, define the basis function associated with that node, and then evaluate it at all other nodes. I1 I2 I3 J1 J2 X Y Z B(X,Y,Z) 2 0 1 1 0 0.666667 0.000000 1.000000 1 0 0 3 0 1 0.000000 0.000000 0.000000 -0 0 0 3 1 0 0.000000 0.000000 1.000000 -0 0 1 2 0 1 0.000000 0.333333 0.000000 -0 0 1 2 1 0 0.000000 0.333333 1.000000 -0 0 2 1 0 1 0.000000 0.666667 0.000000 -0 0 2 1 1 0 0.000000 0.666667 1.000000 -0 0 3 0 0 1 0.000000 1.000000 0.000000 0 0 3 0 1 0 0.000000 1.000000 1.000000 0 1 0 2 0 1 0.333333 0.000000 0.000000 0 1 0 2 1 0 0.333333 0.000000 1.000000 0 1 1 1 0 1 0.333333 0.333333 0.000000 0 1 1 1 1 0 0.333333 0.333333 1.000000 0 1 2 0 0 1 0.333333 0.666667 0.000000 -0 1 2 0 1 0 0.333333 0.666667 1.000000 -0 2 0 1 0 1 0.666667 0.000000 0.000000 0 2 0 1 1 0 0.666667 0.000000 1.000000 1 2 1 0 0 1 0.666667 0.333333 0.000000 -0 2 1 0 1 0 0.666667 0.333333 1.000000 -0 3 0 0 0 1 1.000000 0.000000 0.000000 -0 3 0 0 1 0 1.000000 0.000000 1.000000 -0 fem_basis_test(): Normal end of execution. 07-Jan-2022 20:05:02