Tue Oct 19 11:34:32 2021 fem_basis_test(): Python version: 3.6.9 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. Tue Oct 19 11:34:32 2021