17 September 2021 9:31:15.774 AM FEM_BASIS_TEST: FORTRAN90 version. Test the FEM_BASIS library. TEST01 FEM_BASIS_1D evaluates an arbitrary basis function over an interval. I J X L(I,J)(X) 2 1 0.6667 1.00000 0 3 0.0000 -0.00000 1 2 0.3333 0.00000 2 1 0.6667 1.00000 3 0 1.0000 0.00000 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.00000 0 0 3 0.0000 0.0000 0.00000 1 0 2 0.3333 0.0000 1.00000 2 0 1 0.6667 0.0000 0.00000 3 0 0 1.0000 0.0000 -0.00000 0 1 2 0.0000 0.3333 0.00000 1 1 1 0.3333 0.3333 0.111022E-15 2 1 0 0.6667 0.3333 -0.166533E-15 0 2 1 0.0000 0.6667 0.00000 1 2 0 0.3333 0.6667 -0.166533E-15 0 3 0 0.0000 1.0000 -0.00000 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.00000 0 0 0 4 0.0000 0.0000 0.0000 -0.00000 1 0 0 3 0.2500 0.0000 0.0000 -0.00000 2 0 0 2 0.5000 0.0000 0.0000 -0.00000 3 0 0 1 0.7500 0.0000 0.0000 -0.00000 4 0 0 0 1.0000 0.0000 0.0000 -0.00000 0 1 0 3 0.0000 0.2500 0.0000 -0.00000 1 1 0 2 0.2500 0.2500 0.0000 -0.00000 2 1 0 1 0.5000 0.2500 0.0000 -0.00000 3 1 0 0 0.7500 0.2500 0.0000 -0.00000 0 2 0 2 0.0000 0.5000 0.0000 -0.00000 1 2 0 1 0.2500 0.5000 0.0000 -0.00000 2 2 0 0 0.5000 0.5000 0.0000 -0.00000 0 3 0 1 0.0000 0.7500 0.0000 -0.00000 1 3 0 0 0.2500 0.7500 0.0000 -0.00000 0 4 0 0 0.0000 1.0000 0.0000 -0.00000 0 0 1 3 0.0000 0.0000 0.2500 0.00000 1 0 1 2 0.2500 0.0000 0.2500 0.00000 2 0 1 1 0.5000 0.0000 0.2500 0.00000 3 0 1 0 0.7500 0.0000 0.2500 0.00000 0 1 1 2 0.0000 0.2500 0.2500 0.00000 1 1 1 1 0.2500 0.2500 0.2500 0.00000 2 1 1 0 0.5000 0.2500 0.2500 0.00000 0 2 1 1 0.0000 0.5000 0.2500 0.00000 1 2 1 0 0.2500 0.5000 0.2500 0.00000 0 3 1 0 0.0000 0.7500 0.2500 0.00000 0 0 2 2 0.0000 0.0000 0.5000 0.00000 1 0 2 1 0.2500 0.0000 0.5000 1.00000 2 0 2 0 0.5000 0.0000 0.5000 0.00000 0 1 2 1 0.0000 0.2500 0.5000 0.00000 1 1 2 0 0.2500 0.2500 0.5000 0.00000 0 2 2 0 0.0000 0.5000 0.5000 0.00000 0 0 3 1 0.0000 0.0000 0.7500 0.00000 1 0 3 0 0.2500 0.0000 0.7500 0.00000 0 1 3 0 0.0000 0.2500 0.7500 0.00000 0 0 4 0 0.0000 0.0000 1.0000 0.00000 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.00000 0 3 0.0000 -0.00000 1 2 0.3333 0.00000 2 1 0.6667 1.00000 3 0 1.0000 0.00000 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.00000 0 0 3 0.0000 0.0000 0.00000 1 0 2 0.3333 0.0000 1.00000 2 0 1 0.6667 0.0000 0.00000 3 0 0 1.0000 0.0000 -0.00000 0 1 2 0.0000 0.3333 0.00000 1 1 1 0.3333 0.3333 0.00000 2 1 0 0.6667 0.3333 -0.00000 0 2 1 0.0000 0.6667 0.00000 1 2 0 0.3333 0.6667 -0.00000 0 3 0 0.0000 1.0000 -0.00000 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.00000 0 0 0 4 0.0000 0.0000 0.0000 -0.00000 1 0 0 3 0.2500 0.0000 0.0000 -0.00000 2 0 0 2 0.5000 0.0000 0.0000 -0.00000 3 0 0 1 0.7500 0.0000 0.0000 -0.00000 4 0 0 0 1.0000 0.0000 0.0000 -0.00000 0 1 0 3 0.0000 0.2500 0.0000 -0.00000 1 1 0 2 0.2500 0.2500 0.0000 -0.00000 2 1 0 1 0.5000 0.2500 0.0000 -0.00000 3 1 0 0 0.7500 0.2500 0.0000 -0.00000 0 2 0 2 0.0000 0.5000 0.0000 -0.00000 1 2 0 1 0.2500 0.5000 0.0000 -0.00000 2 2 0 0 0.5000 0.5000 0.0000 -0.00000 0 3 0 1 0.0000 0.7500 0.0000 -0.00000 1 3 0 0 0.2500 0.7500 0.0000 -0.00000 0 4 0 0 0.0000 1.0000 0.0000 -0.00000 0 0 1 3 0.0000 0.0000 0.2500 0.00000 1 0 1 2 0.2500 0.0000 0.2500 0.00000 2 0 1 1 0.5000 0.0000 0.2500 0.00000 3 0 1 0 0.7500 0.0000 0.2500 0.00000 0 1 1 2 0.0000 0.2500 0.2500 0.00000 1 1 1 1 0.2500 0.2500 0.2500 0.00000 2 1 1 0 0.5000 0.2500 0.2500 0.00000 0 2 1 1 0.0000 0.5000 0.2500 0.00000 1 2 1 0 0.2500 0.5000 0.2500 0.00000 0 3 1 0 0.0000 0.7500 0.2500 0.00000 0 0 2 2 0.0000 0.0000 0.5000 0.00000 1 0 2 1 0.2500 0.0000 0.5000 1.00000 2 0 2 0 0.5000 0.0000 0.5000 0.00000 0 1 2 1 0.0000 0.2500 0.5000 0.00000 1 1 2 0 0.2500 0.2500 0.5000 0.00000 0 2 2 0 0.0000 0.5000 0.5000 0.00000 0 0 3 1 0.0000 0.0000 0.7500 0.00000 1 0 3 0 0.2500 0.0000 0.7500 0.00000 0 1 3 0 0.0000 0.2500 0.7500 0.00000 0 0 4 0 0.0000 0.0000 1.0000 0.00000 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.0000 0.0000 0.5000 1.00000 0 0 2 0 2 0.0000 0.0000 0.0000 -0.00000 0 0 2 1 1 0.0000 0.0000 0.5000 -0.00000 0 0 2 2 0 0.0000 0.0000 1.0000 -0.00000 0 1 1 0 2 0.0000 0.5000 0.0000 -0.00000 0 1 1 1 1 0.0000 0.5000 0.5000 -0.00000 0 1 1 2 0 0.0000 0.5000 1.0000 -0.00000 0 2 0 0 2 0.0000 1.0000 0.0000 -0.00000 0 2 0 1 1 0.0000 1.0000 0.5000 -0.00000 0 2 0 2 0 0.0000 1.0000 1.0000 -0.00000 1 0 1 0 2 0.5000 0.0000 0.0000 0.00000 1 0 1 1 1 0.5000 0.0000 0.5000 0.00000 1 0 1 2 0 0.5000 0.0000 1.0000 0.00000 1 1 0 0 2 0.5000 0.5000 0.0000 0.00000 1 1 0 1 1 0.5000 0.5000 0.5000 0.00000 1 1 0 2 0 0.5000 0.5000 1.0000 0.00000 2 0 0 0 2 1.0000 0.0000 0.0000 0.00000 2 0 0 1 1 1.0000 0.0000 0.5000 1.00000 2 0 0 2 0 1.0000 0.0000 1.0000 0.00000 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.6667 0.0000 1.0000 1.00000 0 0 3 0 1 0.0000 0.0000 0.0000 -0.00000 0 0 3 1 0 0.0000 0.0000 1.0000 -0.00000 0 1 2 0 1 0.0000 0.3333 0.0000 -0.00000 0 1 2 1 0 0.0000 0.3333 1.0000 -0.00000 0 2 1 0 1 0.0000 0.6667 0.0000 -0.00000 0 2 1 1 0 0.0000 0.6667 1.0000 -0.00000 0 3 0 0 1 0.0000 1.0000 0.0000 -0.00000 0 3 0 1 0 0.0000 1.0000 1.0000 -0.00000 1 0 2 0 1 0.3333 0.0000 0.0000 0.00000 1 0 2 1 0 0.3333 0.0000 1.0000 0.00000 1 1 1 0 1 0.3333 0.3333 0.0000 0.00000 1 1 1 1 0 0.3333 0.3333 1.0000 0.00000 1 2 0 0 1 0.3333 0.6667 0.0000 0.00000 1 2 0 1 0 0.3333 0.6667 1.0000 0.00000 2 0 1 0 1 0.6667 0.0000 0.0000 0.00000 2 0 1 1 0 0.6667 0.0000 1.0000 1.00000 2 1 0 0 1 0.6667 0.3333 0.0000 0.00000 2 1 0 1 0 0.6667 0.3333 1.0000 0.166533E-15 3 0 0 0 1 1.0000 0.0000 0.0000 0.00000 3 0 0 1 0 1.0000 0.0000 1.0000 0.00000 FEM_BASIS_TEST: Normal end of execution. 17 September 2021 9:31:15.774 AM