17 September 2021 9:31:15.510 AM FEM3D_PACK_TEST FORTRAN90 version Test the FEM3D_PACK library. BASIS_MN_TET4_TEST: Verify basis functions for element TET4. Number of nodes = 4 Tetrahedron Nodes: 1 0.576 0.433 0.125 2 0.947 0.677 0.358 3 0.301 0.306 0.558 4 0.360 0.054 0.928 The basis function values at basis nodes should form the identity matrix. 1.000 -0.000 -0.000 -0.000 -0.000 1.000 -0.000 -0.000 -0.000 -0.000 1.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---------------- PHI_SUM 0.7458 0.5168 0.4359 0.1154 0.6260 0.1152 0.1434 1.0000 0.5013 0.2923 0.5240 0.3767 0.1126 0.1008 0.4099 1.0000 0.5322 0.2982 0.6355 0.1070 0.2686 0.1440 0.4804 1.0000 0.5876 0.3863 0.5228 0.1564 0.3518 0.2138 0.2781 1.0000 0.4799 0.3557 0.3977 0.4536 0.0716 0.3384 0.1363 1.0000 BASIS_MN_TET10_TEST: Verify basis functions for element TET10. Number of nodes = 10. Tetrahedron Nodes: 1 0.453203 0.473228 0.771334 2 0.502755 0.374764 0.518830 3 0.213829 0.177931E-01 0.557857 4 0.747917 0.796065 0.853130 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 1.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 1.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 1.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 1.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 1.000 The basis function values at ANY point P should sum to 1: ------------P------------- ----------------------------------------------------PHI----------------------------------------- PHI_SUM 0.5754 0.5178 0.6593 -0.0532 0.0208 -0.0457 -0.0969 0.1260 0.1059 0.0123 0.0893 0.7665 0.0750 1.0000 0.3680 0.2291 0.5735 -0.0874 -0.0822 -0.0307 -0.0225 0.1789 0.7406 0.2109 0.0107 0.0375 0.0442 1.0000 0.5027 0.4278 0.6561 -0.0857 -0.1034 -0.1246 -0.1200 0.1554 0.3342 0.1036 0.1318 0.4251 0.2835 1.0000 0.4199 0.3190 0.6272 -0.1180 -0.1195 -0.0961 -0.0994 0.2307 0.4476 0.2825 0.1044 0.1653 0.2025 1.0000 0.3906 0.2886 0.6648 -0.0949 -0.0053 0.1149 -0.1241 0.0027 0.0127 0.3038 0.1380 0.0058 0.6464 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.3463 0.4476 0.4709 1.0000 0.0949 0.0444 0.1850 1.0000 0.6552 0.6317 0.8475 1.0000 0.7773 0.1043 0.8428 1.0000 0.7427 0.2847 0.7630 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.6932 0.8338 0.5977 1.0000 0.3579 0.6460 0.0800 1.0000 0.4006 0.1150 0.5823 1.0000 0.8993 0.9823 0.6632 1.0000 0.8420 0.7959 0.2797 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. -0. 0. -0. -0. 0. -0. 0. 0. -0. -0. 0. 0. -0. 0. -0. -0. 0. 0. -0. 0. 0. -0. -0. 0. -0. 0. 0. 1. -0. -0. -0. -0. 0. 0. -0. 0. 0. 0. 0. 0. -0. -0. 0. -0. -0. -0. -0. -0. 0. 0. 0. 0. -0. 0. 0. 1. 0. -0. -0. -0. -0. -0. -0. -0. -0. 0. 0. 0. 0. 0. 0. 0. 0. 0. -0. -0. -0. -0. -0. 0. 0. -0. -0. 1. -0. 0. 0. -0. 0. 0. 0. -0. 0. -0. -0. 0. -0. -0. -0. 0. -0. 0. 0. 0. -0. 0. -0. 0. -0. 0. -0. 1. -0. 0. -0. 0. -0. -0. 0. -0. 0. -0. 0. 0. -0. -0. 0. 0. -0. 0. 0. -0. 0. -0. 0. 0. -0. -0. 0. 1. -0. -0. -0. 0. 0. 0. -0. -0. 0. 0. -0. 0. 0. 0. -0. 0. -0. -0. -0. -0. 0. 0. 0. 0. 0. 0. 0. 1. 0. -0. -0. -0. -0. -0. -0. -0. -0. -0. -0. -0. -0. 0. 0. 0. 0. 0. 0. -0. 0. -0. -0. 0. 0. -0. -0. 1. 0. 0. 0. -0. -0. 0. 0. -0. 0. 0. 0. -0. -0. -0. -0. -0. 0. -0. 0. -0. 0. -0. 0. 0. -0. 0. -0. 1. 0. -0. -0. -0. 0. -0. 0. -0. -0. 0. 0. 0. 0. 0. -0. -0. -0. 0. -0. -0. 0. 0. 0. 0. -0. -0. 0. 1. -0. 0. -0. -0. 0. 0. -0. -0. 0. -0. -0. 0. 0. -0. 0. 0. -0. -0. 0. 0. -0. 0. -0. -0. 0. 0. -0. 1. 0. -0. 0. 0. -0. -0. 0. -0. -0. -0. 0. -0. 0. 0. 0. -0. -0. 0. -0. 0. 0. -0. 0. -0. -0. -0. 0. 1. -0. 0. -0. 0. 0. 0. -0. -0. 0. -0. -0. 0. 0. -0. 0. -0. 0. -0. 0. 0. -0. 0. -0. -0. 0. 0. -0. 1. -0. 0. -0. 0. -0. -0. 0. -0. 0. -0. -0. 0. 0. 0. -0. -0. 0. 0. 0. 0. -0. -0. 0. -0. -0. -0. 0. 1. -0. -0. -0. 0. 0. 0. 0. -0. 0. 0. 0. -0. -0. -0. -0. -0. -0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 1. 0. -0. -0. -0. -0. -0. -0. -0. -0. -0. 0. 0. -0. 0. 0. -0. 0. -0. -0. 0. -0. -0. -0. 0. 0. -0. -0. 1. 0. 0. 0. -0. 0. 0. 0. 0. -0. -0. -0. -0. 0. -0. 0. -0. 0. -0. 0. 0. 0. -0. -0. 0. -0. 0. -0. 1. 0. -0. -0. 0. -0. -0. 0. 0. 0. -0. -0. -0. 0. 0. -0. -0. 0. 0. 0. 0. -0. 0. 0. 0. -0. -0. 0. 1. -0. 0. -0. -0. -0. 0. -0. -0. 0. -0. 0. 0. -0. -0. 0. 0. -0. 0. -0. 0. 0. 0. -0. -0. 0. 0. -0. 1. 0. -0. -0. 0. -0. -0. -0. 0. -0. 0. -0. 0. -0. 0. -0. 0. -0. -0. 0. 0. 0. -0. 0. -0. -0. -0. 0. 1. 0. 0. 0. -0. -0. 0. -0. 0. -0. 0. -0. -0. 0. -0. 0. -0. 0. 0. -0. 0. -0. 0. -0. 0. -0. -0. 0. 1. -0. 0. 0. -0. -0. 0. 0. -0. 0. -0. -0. 0. -0. 0. -0. -0. 0. 0. -0. 0. -0. 0. 0. 0. -0. -0. -0. 1. 0. -0. -0. 0. 0. 0. 0. -0. -0. -0. -0. 0. 0. -0. -0. 0. -0. -0. -0. 0. 0. 0. 0. -0. 0. 0. 0. 1. -0. 0. -0. -0. 0. -0. -0. 0. -0. 0. 0. -0. -0. 0. -0. -0. -0. 0. 0. -0. 0. -0. 0. 0. 0. 0. -0. 1. 0. -0. -0. 0. -0. 0. -0. -0. 0. -0. 0. 0. 0. -0. -0. -0. 0. -0. 0. -0. -0. 0. 0. -0. -0. -0. 0. 1. 0. 0. 0. -0. 0. -0. 0. -0. 0. -0. -0. 0. 0. -0. -0. 0. -0. 0. -0. 0. 0. -0. 0. 0. -0. -0. 0. 1. -0. -0. 0. -0. 0. 0. -0. 0. -0. 0. -0. -0. 0. 0. -0. 0. -0. -0. 0. 0. -0. -0. -0. 0. 0. -0. 0. 1. The basis function values at ANY point P should sum to 1: ------------P------------- PHI_SUM 0.3975 0.2390 0.6530 1.0000 0.9486 0.4470 0.5309 1.0000 0.5347 0.9256 0.7222 1.0000 0.5897 0.5575 0.3522 1.0000 0.2945 0.5184 0.9290 1.0000 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.0507 0.5036 0.4119 2.9903 2.3203 2.9761 0.0507 0.5036 0.4119 0.4460 0.2522 0.0181 2.3421 1.7840 3.3332 0.4460 0.2522 0.0181 0.1401 0.4094 0.2672 2.8620 2.1251 3.2245 0.1401 0.4094 0.2672 0.3168 0.0536 0.0006 2.4194 1.9476 4.2051 0.3168 0.0536 0.0006 0.3668 0.2152 0.1254 2.3561 2.0356 3.4953 0.3668 0.2152 0.1254 0.1528 0.0586 0.1164 2.6366 2.1742 4.4024 0.1528 0.0586 0.1164 0.0582 0.2176 0.4678 2.6334 2.7180 3.7631 0.0582 0.2176 0.4678 0.5990 0.0814 0.0388 1.8445 1.9963 3.5190 0.5990 0.0814 0.0388 0.1184 0.2056 0.6617 2.3070 3.1179 3.4846 0.1184 0.2056 0.6617 0.3838 0.4023 0.1611 2.4736 1.9199 2.8642 0.3838 0.4023 0.1611 FEM3D_PACK_TEST Normal end of execution. 17 September 2021 9:31:15.511 AM