24 December 2011 01:38:00 PM FEM3D_PACK_PRB: C++ version Test the FEM3D_PACK library. BASIS_MN_TET4_TEST: Verify basis functions for element TET4. Number of nodes = 4 Tetrahedron Nodes: 7.82637e-06 0.131538 0.45865 0.532767 0.0470446 0.678865 0.934693 0.383502 The basis function values at basis nodes should form the identity matrix. 1.0000 -0.0000 -0.0000 0.0000 -0.0000 1.0000 -0.0000 0.0000 -0.0000 -0.0000 1.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.3544 0.2534 0.6536 0.5736 0.0239 0.0369 0.3656 1.0000 0.0744 0.3350 0.7121 0.5944 0.0068 0.3396 0.0592 1.0000 0.4621 0.4488 0.5142 0.1591 0.2618 0.2245 0.3546 1.0000 0.4045 0.3328 0.5460 0.3993 0.2487 0.0434 0.3086 1.0000 0.3860 0.4798 0.5313 0.1491 0.2513 0.3263 0.2732 1.0000 BASIS_MN_TET10_TEST: Verify basis functions for element TET10. Number of nodes = 10. Tetrahedron Nodes: 0 0.2184 0.9563 0.8295 1 0.5617 0.4153 0.0661 2 0.2576 0.1100 0.0438 3 0.6340 0.0617 0.4495 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.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.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 1.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 1.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 1.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_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.0000-0.0000-0.0000-0.0000-0.0000-0.0000-0.0000-0.0000 0.0000 0.0000 0.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 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 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 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 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 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 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 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 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 -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-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-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-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-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-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-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-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 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 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_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.0000-0.0000 0.0000-0.0000-0.0000 0.0000-0.0000 0.0000 0.0000-0.0000-0.0000 0.0000 0.0000-0.0000 0.0000-0.0000-0.0000 0.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 0.0000 0.0000-0.0000-0.0000 0.0000-0.0000-0.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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.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-0.0000 0.0000-0.0000-0.0000-0.0000 0.0000-0.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 0.0000 0.0000-0.0000-0.0000 0.0000 0.0000-0.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-0.0000 0.0000 0.0000 0.0000-0.0000 0.0000-0.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-0.0000-0.0000-0.0000 0.0000 0.0000 0.0000 0.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 0.0000-0.0000-0.0000-0.0000-0.0000-0.0000 0.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 0.0000 0.0000 0.0000 0.0000-0.0000-0.0000-0.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-0.0000 0.0000 0.0000-0.0000 0.0000 0.0000-0.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 0.0000-0.0000 0.0000 0.0000 0.0000-0.0000 -0.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-0.0000 0.0000 0.0000-0.0000 0.0000 -0.0000 0.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-0.0000 0.0000 0.0000 0.0000 -0.0000-0.0000 0.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 0.0000-0.0000-0.0000 -0.0000-0.0000-0.0000-0.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 0.0000 0.0000 -0.0000 0.0000 0.0000-0.0000 0.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-0.0000 -0.0000 0.0000-0.0000 0.0000-0.0000 0.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 -0.0000-0.0000 0.0000 0.0000-0.0000-0.0000 0.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 0.0000 0.0000-0.0000-0.0000 0.0000 0.0000-0.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-0.0000 0.0000-0.0000 0.0000-0.0000 0.0000-0.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-0.0000-0.0000 0.0000-0.0000 0.0000-0.0000 0.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-0.0000 0.0000-0.0000 0.0000-0.0000-0.0000 0.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-0.0000 0.0000 0.0000-0.0000-0.0000 0.0000-0.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 0.0000-0.0000-0.0000 0.0000-0.0000-0.0000-0.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 0.0000 0.0000 0.0000-0.0000-0.0000-0.0000 0.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-0.0000 0.0000 0.0000-0.0000-0.0000 0.0000-0.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-0.0000-0.0000 0.0000 0.0000-0.0000 0.0000-0.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_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. 24 December 2011 01:38:00 PM