2 March 2010 9:03:54.010 AM FEM3D_PACK_PRB FORTRAN90 version Test routines in the FEM3D_PACK library. BASIS_MN_TET4_TEST: Verify basis functions for element TET4. Number of nodes = 4 Tetrahedron Nodes: 1 0.218 0.956 0.830 2 0.562 0.415 0.066 3 0.258 0.110 0.044 4 0.634 0.062 0.450 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.3673 0.4015 0.2141 0.2053 0.3860 0.4078 0.0009 1.0000 0.3141 0.7484 0.5735 0.6617 0.2586 0.0697 0.0100 1.0000 0.3804 0.6472 0.4245 0.4693 0.4593 0.0672 0.0041 1.0000 0.5016 0.4042 0.2682 0.1589 0.5570 0.0694 0.2147 1.0000 0.3687 0.4304 0.4194 0.3511 0.1140 0.2953 0.2396 1.0000 BASIS_MN_TET10_TEST: Verify basis functions for element TET10. Number of nodes = 10. Tetrahedron Nodes: 1 0.218418 0.956318 0.829509 2 0.561695 0.415307 0.661187E-01 3 0.257578 0.109957 0.438290E-01 4 0.633966 0.617272E-01 0.449539 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.3673 0.4015 0.2141 -0.1210 -0.0880 -0.0752 -0.0009 0.3169 0.6296 0.3348 0.0008 0.0015 0.0015 1.0000 0.3141 0.7484 0.5735 0.2139 -0.1249 -0.0600 -0.0098 0.6844 0.0721 0.1845 0.0266 0.0104 0.0028 1.0000 0.3804 0.6472 0.4245 -0.0288 -0.0374 -0.0582 -0.0041 0.8623 0.1236 0.1262 0.0077 0.0075 0.0011 1.0000 0.5016 0.4042 0.2682 -0.1084 0.0636 -0.0598 -0.1225 0.3541 0.1546 0.0441 0.1364 0.4783 0.0596 1.0000 0.3687 0.4304 0.4194 -0.1046 -0.0880 -0.1209 -0.1248 0.1601 0.1346 0.4147 0.3366 0.1093 0.2830 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.2184 0.9563 0.8295 1.0000 0.5617 0.4153 0.0661 1.0000 0.2576 0.1100 0.0438 1.0000 0.6340 0.0617 0.4495 1.0000 0.4013 0.7547 0.7973 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.2184 0.9563 0.8295 1.0000 0.5617 0.4153 0.0661 1.0000 0.2576 0.1100 0.0438 1.0000 0.6340 0.0617 0.4495 1.0000 0.4013 0.7547 0.7973 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.2184 0.9563 0.8295 1.0000 0.5617 0.4153 0.0661 1.0000 0.2576 0.1100 0.0438 1.0000 0.6340 0.0617 0.4495 1.0000 0.4013 0.7547 0.7973 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.6530 0.0192 0.0802 1.6329 2.1413 3.5562 0.6530 0.0192 0.0802 0.1227 0.3794 0.1895 2.9445 1.9995 3.4268 0.1227 0.3794 0.1895 0.4363 0.0636 0.3885 1.8024 2.7135 3.5481 0.4363 0.0636 0.3885 0.1183 0.0365 0.0293 2.7706 2.0222 4.6247 0.1183 0.0365 0.0293 0.0138 0.1341 0.3020 2.8045 2.4698 4.2680 0.0138 0.1341 0.3020 0.0208 0.0237 0.2865 2.6957 2.5493 4.6008 0.0208 0.0237 0.2865 0.2890 0.0197 0.4669 1.9748 2.9142 3.8961 0.2890 0.0197 0.4669 0.0792 0.5366 0.1496 3.2285 1.7626 3.0820 0.0792 0.5366 0.1496 0.0966 0.5111 0.0596 3.2582 1.6081 3.2139 0.0966 0.5111 0.0596 0.3663 0.0599 0.2031 2.1241 2.3463 3.8845 0.3663 0.0599 0.2031 FEM3D_PACK_PRB Normal end of execution. 2 March 2010 9:03:54.018 AM