07-Jan-2022 20:03:53 fem3d_pack_test(): MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2 Test fem3d_pack(). BASIS_MN_TET4_TEST: Verify basis functions for element TET4. Number of nodes = 4 Tetrahedron Nodes: 1.000000 0.814724 0.905792 0.126987 2.000000 0.913376 0.632359 0.097540 3.000000 0.278498 0.546882 0.957507 4.000000 0.964889 0.157613 0.970593 The basis function values at basis nodes should form the identity matrix. 1.000000 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 1.000000 The basis function values at ANY point P should sum to 1: ------------P------------- -----------------PHI---------------- PHI_SUM 0.663787 0.685177 0.449899 0.401377 0.203537 0.335588 0.059498 1.000000 0.753096 0.500316 0.593261 0.136528 0.296432 0.256444 0.310596 1.000000 0.688824 0.496299 0.729955 0.264991 0.014431 0.343141 0.377436 1.000000 0.711749 0.607419 0.486947 0.263911 0.294630 0.288950 0.152509 1.000000 0.586586 0.698682 0.515726 0.418958 0.109416 0.451279 0.020346 1.000000 BASIS_MN_TET10_TEST: Verify basis functions for element TET10. Number of nodes = 10. Tetrahedron Nodes: 1.0000 0.2769 0.0462 0.0971 2.0000 0.8235 0.6948 0.3171 3.0000 0.9502 0.0344 0.4387 4.0000 0.3816 0.7655 0.7952 The basis function values at basis nodes should form the identity matrix. 1.0000 -0.0000 -0.0000 -0.0000 0.0000 0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 1.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 -0.0000 1.0000 -0.0000 -0.0000 0.0000 0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 1.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 -0.0000 0.0000 -0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 -0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 -0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 -0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 -0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 The basis function values at ANY point P should sum to 1: ------------P------------- ----------------------------------------------------PHI----------------------------------------- PHI_SUM 0.6362 0.4857 0.4992 -0.0833 -0.1235 -0.1250 -0.0983 0.1170 0.2791 0.1065 0.1545 0.4048 0.3683 1.0000 0.5536 0.4492 0.4008 -0.1213 -0.1173 -0.0880 -0.1231 0.3657 0.1423 0.1338 0.3294 0.3504 0.1282 1.0000 0.4310 0.3685 0.3928 -0.0397 -0.0876 -0.0692 -0.1061 0.2069 0.0376 0.1514 0.6341 0.1574 0.1152 1.0000 0.5631 0.2239 0.3015 -0.0409 -0.1093 -0.1235 -0.0836 0.2937 0.1791 0.5050 0.1931 0.0685 0.1178 1.0000 0.5271 0.3457 0.4213 -0.1089 -0.0887 -0.1241 -0.1163 0.1568 0.1056 0.3109 0.4296 0.1459 0.2893 1.0000 BASIS_BRICK8_TEST: Verify basis functions for element BRICK8. Number of nodes = 8 The basis function values at basis nodes should form the identity matrix. 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000 The basis function values at ANY point P should sum to 1: ------------P------------- PHI_SUM 0.8909 0.9593 0.5472 1.0000 0.1386 0.1493 0.2575 1.0000 0.8407 0.2543 0.8143 1.0000 0.2435 0.9293 0.3500 1.0000 0.1966 0.2511 0.6160 1.0000 BASIS_BRICK20_TEST(): Verify basis functions for element BRICK20. Number of nodes = 20 The basis function values at basis nodes should form the identity matrix. 1.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.000 1.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 1.000 -0.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 -0.000 1.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 -0.000 -0.000 1.000 -0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 -0.000 -0.000 -0.000 1.000 -0.000 -0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 1.000 -0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.000 0.000 0.000 -0.000 -0.000 -0.000 -0.000 -0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 0.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 -0.000 0.000 -0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.000 0.000 -0.000 -0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 0.000 -0.000 -0.000 -0.000 0.000 -0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 -0.000 0.000 -0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 -0.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000 -0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 -0.000 -0.000 -0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000 The basis function values at ANY point P should sum to 1: ------------P------------- PHI_SUM 0.4733 0.3517 0.8308 1.0000 0.5853 0.5497 0.9172 1.0000 0.2858 0.7572 0.7537 1.0000 0.3804 0.5678 0.0759 1.0000 0.0540 0.5308 0.7792 1.0000 BASIS_BRICK27_TEST: Verify basis functions for element BRICK27. Number of nodes = 27 The basis function values at basis nodes should form the identity matrix. 1.000 -0.000 0.000 -0.000 -0.000 0.000 -0.000 0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 0.000 0.000 -0.000 -0.000 0.000 -0.000 0.000 0.000 1.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000 -0.000 0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 0.000 -0.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 -0.000 0.000 0.000 1.000 0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000 -0.000 -0.000 1.000 -0.000 0.000 0.000 -0.000 0.000 0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000 -0.000 -0.000 -0.000 0.000 -0.000 0.000 0.000 0.000 -0.000 0.000 -0.000 0.000 -0.000 0.000 -0.000 1.000 -0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000 -0.000 0.000 -0.000 0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 0.000 0.000 -0.000 0.000 -0.000 0.000 0.000 -0.000 -0.000 0.000 1.000 -0.000 -0.000 -0.000 0.000 0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 0.000 0.000 0.000 -0.000 0.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 -0.000 1.000 0.000 0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 0.000 0.000 0.000 -0.000 -0.000 -0.000 -0.000 -0.000 0.000 -0.000 0.000 -0.000 0.000 -0.000 0.000 0.000 -0.000 0.000 -0.000 1.000 0.000 -0.000 -0.000 -0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 0.000 1.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 1.000 0.000 -0.000 0.000 0.000 -0.000 -0.000 0.000 -0.000 -0.000 -0.000 0.000 -0.000 0.000 0.000 0.000 -0.000 -0.000 0.000 -0.000 0.000 0.000 -0.000 0.000 -0.000 -0.000 -0.000 0.000 1.000 -0.000 0.000 -0.000 0.000 0.000 0.000 -0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 0.000 -0.000 0.000 -0.000 0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 1.000 -0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 0.000 -0.000 -0.000 -0.000 0.000 1.000 -0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 -0.000 0.000 0.000 0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000 -0.000 0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000 -0.000 -0.000 -0.000 0.000 0.000 -0.000 -0.000 1.000 0.000 0.000 0.000 -0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 -0.000 -0.000 0.000 -0.000 0.000 -0.000 0.000 -0.000 0.000 0.000 0.000 -0.000 -0.000 0.000 -0.000 0.000 -0.000 1.000 0.000 -0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 0.000 -0.000 -0.000 -0.000 0.000 0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 -0.000 0.000 0.000 0.000 -0.000 -0.000 0.000 1.000 -0.000 0.000 -0.000 -0.000 -0.000 0.000 -0.000 -0.000 0.000 -0.000 0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 0.000 -0.000 0.000 0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 1.000 0.000 -0.000 -0.000 0.000 -0.000 -0.000 -0.000 0.000 -0.000 0.000 -0.000 0.000 -0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 0.000 -0.000 0.000 -0.000 -0.000 -0.000 0.000 1.000 0.000 0.000 0.000 -0.000 -0.000 0.000 -0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000 -0.000 0.000 -0.000 0.000 0.000 -0.000 0.000 -0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000 1.000 -0.000 0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 0.000 -0.000 0.000 0.000 0.000 -0.000 -0.000 -0.000 1.000 0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 -0.000 -0.000 -0.000 -0.000 0.000 0.000 -0.000 -0.000 0.000 -0.000 -0.000 -0.000 0.000 0.000 0.000 0.000 -0.000 0.000 0.000 0.000 1.000 -0.000 0.000 -0.000 -0.000 0.000 -0.000 -0.000 0.000 -0.000 0.000 0.000 -0.000 -0.000 0.000 -0.000 -0.000 -0.000 0.000 0.000 -0.000 0.000 -0.000 0.000 0.000 0.000 0.000 -0.000 1.000 0.000 -0.000 -0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000 -0.000 0.000 0.000 0.000 -0.000 -0.000 -0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 -0.000 -0.000 0.000 1.000 0.000 0.000 0.000 -0.000 0.000 -0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 -0.000 0.000 -0.000 0.000 -0.000 0.000 0.000 -0.000 0.000 0.000 -0.000 -0.000 0.000 1.000 -0.000 -0.000 0.000 -0.000 0.000 0.000 -0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 0.000 -0.000 -0.000 0.000 0.000 -0.000 -0.000 -0.000 0.000 0.000 -0.000 0.000 1.000 The basis function values at ANY point P should sum to 1: ------------P------------- PHI_SUM 0.9340 0.1299 0.5688 1.0000 0.4694 0.0119 0.3371 1.0000 0.1622 0.7943 0.3112 1.0000 0.5285 0.1656 0.6020 1.0000 0.2630 0.6541 0.6892 1.0000 FEM3D_PACK_TEST03(): For an order 4 tetrahedron, PHYSICAL_TO_REFERENCE_TET4 maps a physical point to a reference point. REFERENCE_TO_PHYSICAL_TET4 maps a reference point to a physical point. ( R S T ) ==> ( X Y Z ) ==> ( R2 S2 T2 ) 0.0576 0.1582 0.4918 2.5512 2.8255 3.9185 0.0576 0.1582 0.4918 0.0326 0.6761 0.0688 3.5422 1.4614 2.8378 0.0326 0.6761 0.0688 0.0007 0.4547 0.1454 3.3079 1.8361 3.4892 0.0007 0.4547 0.1454 0.0066 0.9156 0.0434 3.8589 1.1713 2.1966 0.0066 0.9156 0.0434 0.0289 0.5068 0.1880 3.2611 1.8692 3.2338 0.0289 0.5068 0.1880 0.0779 0.2937 0.0327 3.1052 1.7717 3.9302 0.0779 0.2937 0.0327 0.2470 0.3573 0.3697 2.4935 2.3822 3.0644 0.2470 0.3573 0.3697 0.1686 0.1850 0.5973 2.2505 3.0096 3.5106 0.1686 0.1850 0.5973 0.1531 0.3377 0.2151 2.8164 2.0925 3.4655 0.1531 0.3377 0.2151 0.3308 0.1832 0.2686 2.2530 2.3541 3.5201 0.3308 0.1832 0.2686 fem3d_pack_test() Normal end of execution. 07-Jan-2022 20:03:54