05 July 2013 10:13:07 AM FEM2D_PACK_PRB: C++ version Test the routines in the FEM2D_PACK library. TEST01 BANDWIDTH_MESH computes the geometric bandwidth: of a finite element mesh. NELEMX = 2 NELEMY = 6 ELEMENT_ORDER = 6 ELEMENT_NUM = 24 GRID_PRINT: Element -> Node table. Element order = 6 Number of elements = 24 # 0 1 2 3 4 5 0 1 3 11 2 7 6 1 13 11 3 12 7 8 2 3 5 13 4 9 8 3 15 13 5 14 9 10 4 11 13 21 12 17 16 5 23 21 13 22 17 18 6 13 15 23 14 19 18 7 25 23 15 24 19 20 8 21 23 31 22 27 26 9 33 31 23 32 27 28 10 23 25 33 24 29 28 11 35 33 25 34 29 30 12 31 33 41 32 37 36 13 43 41 33 42 37 38 14 33 35 43 34 39 38 15 45 43 35 44 39 40 16 41 43 51 42 47 46 17 53 51 43 52 47 48 18 43 45 53 44 49 48 19 55 53 45 54 49 50 20 51 53 61 52 57 56 21 63 61 53 62 57 58 22 53 55 63 54 59 58 23 65 63 55 64 59 60 Lower bandwidth ML = 10 Upper bandwidth MU = 10 Total bandwidth M = 21 TEST02 For the Navier Stokes variables associated with a mesh of T6 elements, NS_T6_VAR_COUNT counts variables, NS_T6_VAR_SET sets them, BANDWIDTH_VAR computes the variable bandwidth. NELEMX = 2 NELEMY = 6 ELEMENT_ORDER = 6 ELEMENT_NUM = 24 NODE_NUM = 65 GRID_PRINT: Element -> Node table. Element order = 6 Number of elements = 24 # 0 1 2 3 4 5 0 1 3 11 2 7 6 1 13 11 3 12 7 8 2 3 5 13 4 9 8 3 15 13 5 14 9 10 4 11 13 21 12 17 16 5 23 21 13 22 17 18 6 13 15 23 14 19 18 7 25 23 15 24 19 20 8 21 23 31 22 27 26 9 33 31 23 32 27 28 10 23 25 33 24 29 28 11 35 33 25 34 29 30 12 31 33 41 32 37 36 13 43 41 33 42 37 38 14 33 35 43 34 39 38 15 45 43 35 44 39 40 16 41 43 51 42 47 46 17 53 51 43 52 47 48 18 43 45 53 44 49 48 19 55 53 45 54 49 50 20 51 53 61 52 57 56 21 63 61 53 62 57 58 22 53 55 63 54 59 58 23 65 63 55 64 59 60 Number of variables VAR_NUM = 151 VAR_NODE pointer vector: 1 1 2 4 3 6 4 9 5 11 6 14 7 16 8 18 9 20 10 22 11 24 12 27 13 29 14 32 15 34 16 37 17 39 18 41 19 43 20 45 21 47 22 50 23 52 24 55 25 57 26 60 27 62 28 64 29 66 30 68 31 70 32 73 33 75 34 78 35 80 36 83 37 85 38 87 39 89 40 91 41 93 42 96 43 98 44 101 45 103 46 106 47 108 48 110 49 112 50 114 51 116 52 119 53 121 54 124 55 126 56 129 57 131 58 133 59 135 60 137 61 139 62 142 63 144 64 147 65 149 66 152 Node U_Var V_Var P_Var 0 1 2 3 1 4 5 2 6 7 8 3 9 10 4 11 12 13 5 14 15 6 16 17 7 18 19 8 20 21 9 22 23 10 24 25 26 11 27 28 12 29 30 31 13 32 33 14 34 35 36 15 37 38 16 39 40 17 41 42 18 43 44 19 45 46 20 47 48 49 21 50 51 22 52 53 54 23 55 56 24 57 58 59 25 60 61 26 62 63 27 64 65 28 66 67 29 68 69 30 70 71 72 31 73 74 32 75 76 77 33 78 79 34 80 81 82 35 83 84 36 85 86 37 87 88 38 89 90 39 91 92 40 93 94 95 41 96 97 42 98 99 100 43 101 102 44 103 104 105 45 106 107 46 108 109 47 110 111 48 112 113 49 114 115 50 116 117 118 51 119 120 52 121 122 123 53 124 125 54 126 127 128 55 129 130 56 131 132 57 133 134 58 135 136 59 137 138 60 139 140 141 61 142 143 62 144 145 146 63 147 148 64 149 150 151 Lower bandwidth ML = 25 Upper bandwidth MU = 25 Total bandwidth M = 51 TEST03 BASIS_11_T3_TEST - Test the T3 basis functions. BASIS_11_T4_TEST - Test the T4 basis functions. BASIS_11_T6_TEST - Test the T6 basis functions. BASIS_11_T3_TEST: Verify basis functions for element T3. Number of nodes = 3 Physical Nodes: 2 0 4 3 0 4 The basis function values at basis nodes should form the identity matrix. 1 0 0 -0 1 0 0 0 1 The X and Y derivatives should sum to 0. dPhidX sum, dPhidY sum: 0 0 0 0 0 0 BASIS_11_T4_TEST: Verify basis functions for element T4. Number of nodes = 4 Physical Nodes: 2 0 4 3 0 4 2 2.33333 The basis function values at basis nodes should form the identity matrix. 1 0 0 0 0 1 0 -0 0 0 1 0 -1.11022e-16 -5.55112e-17 -5.55112e-17 1 The X and Y derivatives should sum to 0. dPhidX sum, dPhidY sum: 0 0 0 0 0 0 2.08167e-17 1.73472e-17 BASIS_11_T6_TEST: Verify basis functions for element T6. Number of nodes = 6 Physical Nodes: 0 2 0 1 4 3 2 0 4 3 3 1.5 4 2 3.5 5 1 2 The basis function values at basis nodes should form the identity matrix. 1 0 0 -0 0 -0 0 1 -0 -0 0 0 0 0 1 0 -0 -0 0 -0 -0 1 -0 0 -0 0 0 0 1 -0 -0 -0 -0 -0 0 1 The X and Y derivatives should sum to 0. dPhidX sum, dPhidY sum: 0 0 0 0 0 0 0 0 -2.77556e-17 0 0 0 TEST04n Test the computation of basis functions by evaluating them at the nodes that define the basis functions. BASIS_MN_Q4_TEST - for the Q4 element. BASIS_MN_T3_TEST - for the T3 element. BASIS_MN_T4_TEST - for the T4 element. BASIS_MN_T6_TEST - for the T6 element. BASIS_MN_Q4_TEST: Verify basis functions for element Q4. Number of nodes = 4 Physical Nodes: 3 1 5 1 5 4 3 4 The basis function values at basis nodes should form the identity matrix. 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 The X and Y derivatives should sum to 0. dPhidX sum, dPhidY sum: 0 0 0 0 0 0 0 0 BASIS_MN_T3_TEST: Verify basis functions for element T3. Number of nodes = 3 Physical Nodes: 2 0 4 3 0 4 The basis function values at basis nodes should form the identity matrix. 1 0 0 -0 1 0 0 0 1 The X and Y derivatives should sum to 0. dPhidX sum, dPhidY sum: 0 0 0 0 0 0 BASIS_MN_T4_TEST: Verify basis functions for element T4. Number of nodes = 4 Physical Nodes: 2 0 4 2 0 4 2 2 The basis function values at basis nodes should form the identity matrix. 1 0 0 0 0 1 0 -0 0 0 1 0 0 0 0 1 The X and Y derivatives should sum to 0. dPhidX sum, dPhidY sum: 0 0 0 0 0 0 0 0 BASIS_MN_T6_TEST: Verify basis functions for element T6. Number of nodes = 6 Physical Nodes: 0 2 0 1 4 3 2 0 4 3 3 1.5 4 2 3.5 5 1 2 The basis function values at basis nodes should form the identity matrix. 1 0 0 -0 0 -0 -0 1 -0 -0 0 -0 -0 -0 1 -0 -0 -0 -0 0 -0 1 0 -0 -0 -0 0 -0 1 0 0 -0 -0 0 -0 1 The X and Y derivatives should sum to 0. dPhidX sum, dPhidY sum: 0 0 0 0 -1.11022e-16 0 0 0 -2.77556e-17 0 0 0 TEST05 DERIVATIVE_AVERAGE_T3 averages the spatial derivatives of a finite element function at the nodes. C X Y dCdX(computed) dCdY(computed) dCdX(exact) dCdY(exact) -0.974455 -0.707107 0.707107 -0.573115 -0.737595 0.380122 -0.918725 -0.87098 -0.532032 0.846724 -0.649263 -0.680676 0.617845 -0.859062 -0.613237 -0.330279 0.943883 -0.719895 -0.571736 0.842763 -0.612216 -0.222061 -0.111964 0.993712 -0.0677977 -0.539862 0.981281 -0.222056 0.222061 0.111964 0.993712 0.12476 1.1947 0.981281 0.222056 0.613237 0.330279 0.943883 1.26831 0.622382 0.842763 0.612216 0.87098 0.532032 0.846724 -0.236644 -2.75476 0.617845 0.859062 0.974455 0.707107 0.707107 -0.370105 0.923911 0.380122 0.918725 -1.65528 -0.989949 0.989949 -0.104417 -0.15575 0.537757 -1.65519 -1.63039 -0.744845 1.18541 -0.0344645 0.266935 1.0331 -1.60708 -1.22505 -0.462391 1.32144 0.318631 0.00563349 1.56282 -1.17896 -0.458247 -0.15675 1.3912 0.119043 0.489082 1.9117 -0.434357 0.458247 0.15675 1.3912 6.86787 14.9048 1.9117 0.434357 1.22505 0.462391 1.32144 6.16566 11.7814 1.56282 1.17896 1.63039 0.744845 1.18541 -1.97411 -2.55127 1.0331 1.60708 1.65528 0.989949 0.989949 -1.09411 0.176054 0.537757 1.65519 -2.50452 -1.27279 1.27279 0.314138 0.475937 0.475653 -2.43339 -2.71761 -0.957658 1.5241 -0.105249 0.209391 1.33668 -2.49296 -2.17685 -0.594502 1.69899 0.0054792 -0.314391 2.39131 -1.9032 -0.84061 -0.201536 1.78868 -0.206591 -0.22215 3.13463 -0.716097 0.84061 0.201536 1.78868 6.64042 15.1651 3.13463 0.716097 2.17685 0.594502 1.69899 -0.0302642 -0.784298 2.39131 1.9032 2.71761 0.957658 1.5241 -0.787231 -1.7835 1.33668 2.49296 2.50452 1.27279 1.27279 -0.263208 1.29966 0.475653 2.43339 -3.41961 -1.55563 1.55563 -0.433069 0.00324528 0.0366892 -3.11091 -4.11657 -1.17047 1.86279 -0.359672 0.260016 1.35232 -3.43102 -3.52901 -0.726614 2.07654 0.470846 -0.453528 3.22293 -2.75907 -1.40922 -0.246322 2.18617 -0.0045631 0.942909 4.63507 -1.06614 1.40922 0.246322 2.18617 -5.39071 -13.0176 4.63507 1.06614 3.52901 0.726614 2.07654 -5.70315 -17.1507 3.22293 2.75907 4.11657 1.17047 1.86279 -1.94389 -2.37511 1.35232 3.43102 3.41961 1.55563 1.55563 -1.58897 0.428659 0.0366892 3.11091 -4.22401 -1.83848 1.83848 0.145085 0.898265 -0.893997 -3.54601 -5.74404 -1.38328 2.20148 -0.205636 0.463511 0.90347 -4.32579 -5.31617 -0.858726 2.4541 0.166724 -0.81529 3.93518 -3.71555 -2.2029 -0.291108 2.58365 -0.536682 -0.226744 6.39441 -1.48309 2.2029 0.291108 2.58365 -5.70408 -12.6417 6.39441 1.48309 5.31617 0.858726 2.4541 0.765323 -0.659221 3.93518 3.71555 5.74404 1.38328 2.20148 -0.0339366 0.282309 0.90347 4.32579 4.22401 1.83848 1.83848 -0.598544 1.97957 -0.893997 3.54601 -4.68738 -2.12132 2.12132 -1.22451 -0.0928882 -2.3541 -3.61579 -7.45009 -1.5961 2.54017 -1.45253 -0.876968 -0.16323 -5.07872 -7.54362 -0.990837 2.83165 0.26328 -1.76023 4.39391 -4.73727 -3.25894 -0.335893 2.98114 0.996952 4.36828 8.39053 -1.96524 3.25894 0.335893 2.98114 3.8615 5.37918 8.39053 1.96524 7.54362 0.990837 2.83165 3.53784 0.664369 4.39391 4.73727 7.45009 1.5961 2.54017 -0.397988 1.00846 -0.16323 5.07872 4.68738 2.12132 2.12132 -4.22515 0.976767 -2.3541 3.61579 TEST07 ELEMENTS_EPS creates an Encapsulated PostScript file containing an image of a mesh. TEST08 ELEMENTS_EPS creates an Encapsulated PostScript file containing an image of a mesh. TEST09 ELEMENTS_EPS creates an Encapsulated PostScript file containing an image of a T4 mesh. TEST10 ELEMENTS_EPS creates an Encapsulated PostScript file containing an image of a mesh. TEST105 GRID_NODES_01 computes a regular grid in the unit square. NUM_X = 5 NUM_Y = 3 NODE_NUM = 15 0 0 0 1 0.25 0 2 0.5 0 3 0.75 0 4 1 0 5 0 0.5 6 0.25 0.5 7 0.5 0.5 8 0.75 0.5 9 1 0.5 10 0 1 11 0.25 1 12 0.5 1 13 0.75 1 14 1 1 TEST11 GRID_TEST tests the grid routines. GRID_TEST: Test the grid routine for element Q4 GRID_PRINT: Element -> Node table. Element order = 4 Number of elements = 6 # 0 1 2 3 0 1 2 6 5 1 2 3 7 6 2 3 4 8 7 3 5 6 10 9 4 6 7 11 10 5 7 8 12 11 Grid width is 5 GRID_TEST: Test the grid routine for element Q8 GRID_PRINT: Element -> Node table. Element order = 8 Number of elements = 6 # 0 1 2 3 4 5 6 7 0 1 3 14 12 2 9 13 8 1 3 5 16 14 4 10 15 9 2 5 7 18 16 6 11 17 10 3 12 14 25 23 13 20 24 19 4 14 16 27 25 15 21 26 20 5 16 18 29 27 17 22 28 21 Grid width is 13 GRID_TEST: Test the grid routine for element Q9 GRID_PRINT: Element -> Node table. Element order = 9 Number of elements = 6 # 0 1 2 3 4 5 6 7 8 0 1 3 17 15 2 10 16 8 9 1 3 5 19 17 4 12 18 10 11 2 5 7 21 19 6 14 20 12 13 3 15 17 31 29 16 24 30 22 23 4 17 19 33 31 18 26 32 24 25 5 19 21 35 33 20 28 34 26 27 Grid width is 16 GRID_TEST: Test the grid routine for element Q12 GRID_PRINT: Element -> Node table. Element order = 12 Number of elements = 6 # 0 1 2 3 4 5 6 7 8 9 10 11 0 1 2 3 4 11 12 15 16 19 20 21 22 1 4 5 6 7 12 13 16 17 22 23 24 25 2 7 8 9 10 13 14 17 18 25 26 27 28 3 19 20 21 22 29 30 33 34 37 38 39 40 4 22 23 24 25 30 31 34 35 40 41 42 43 5 25 26 27 28 31 32 35 36 43 44 45 46 Grid width is 21 GRID_TEST: Test the grid routine for element Q16 GRID_PRINT: Element -> Node table. Element order = 16 Number of elements = 6 # 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 0 1 2 3 4 11 12 13 14 21 22 23 24 31 32 33 34 1 4 5 6 7 14 15 16 17 24 25 26 27 34 35 36 37 2 7 8 9 10 17 18 19 20 27 28 29 30 37 38 39 40 3 31 32 33 34 41 42 43 44 51 52 53 54 61 62 63 64 4 34 35 36 37 44 45 46 47 54 55 56 57 64 65 66 67 5 37 38 39 40 47 48 49 50 57 58 59 60 67 68 69 70 Grid width is 33 GRID_TEST: Test the grid routine for element QL GRID_PRINT: Element -> Node table. Element order = 6 Number of elements = 6 # 0 1 2 3 4 5 0 1 2 3 8 9 10 1 3 4 5 10 11 12 2 5 6 7 12 13 14 3 8 9 10 15 16 17 4 10 11 12 17 18 19 5 12 13 14 19 20 21 Grid width is 9 GRID_TEST: Test the grid routine for element T3 GRID_PRINT: Element -> Node table. Element order = 3 Number of elements = 12 # 0 1 2 0 1 2 5 1 6 5 2 2 2 3 6 3 7 6 3 4 3 4 7 5 8 7 4 6 5 6 9 7 10 9 6 8 6 7 10 9 11 10 7 10 7 8 11 11 12 11 8 Grid width is 4 GRID_TEST: Test the grid routine for element T6 GRID_PRINT: Element -> Node table. Element order = 6 Number of elements = 12 # 0 1 2 3 4 5 0 1 3 15 2 9 8 1 17 15 3 16 9 10 2 3 5 17 4 11 10 3 19 17 5 18 11 12 4 5 7 19 6 13 12 5 21 19 7 20 13 14 6 15 17 29 16 23 22 7 31 29 17 30 23 24 8 17 19 31 18 25 24 9 33 31 19 32 25 26 10 19 21 33 20 27 26 11 35 33 21 34 27 28 Grid width is 14 GRID_TEST: Test the grid routine for element T10 GRID_PRINT: Element -> Node table. Element order = 10 Number of elements = 12 # 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 13 22 31 21 10 11 1 34 33 32 31 22 13 4 14 24 23 2 4 5 6 7 16 25 34 24 13 14 3 37 36 35 34 25 16 7 17 27 26 4 7 8 9 10 19 28 37 27 16 17 5 40 39 38 37 28 19 10 20 30 29 6 31 32 33 34 43 52 61 51 40 41 7 64 63 62 61 52 43 34 44 54 53 8 34 35 36 37 46 55 64 54 43 44 9 67 66 65 64 55 46 37 47 57 56 10 37 38 39 40 49 58 67 57 46 47 11 70 69 68 67 58 49 40 50 60 59 Grid width is 30 TEST12 INTERP_TEST tests the interpolating power of the element. INTERP_TEST for element "Q4". Number of nodes = 4 Interpolate R^0 * S^0 (R,S,U): 0 0 1 (R,S,U): 1 0 1 (R,S,U): 1 1 1 (R,S,U): 0 1 1 (R,S): 0.218418 0.956318 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0 0 0 (R,S): 0.829509 0.561695 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0 0 0 (R,S): 0.415307 0.0661187 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0 0 0 (R,S): 0.257578 0.109957 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0 0 0 (R,S): 0.043829 0.633966 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0 0 0 Interpolate R^0 * S^1 (R,S,U): 0 0 0 (R,S,U): 1 0 0 (R,S,U): 1 1 1 (R,S,U): 0 1 1 (R,S): 0.0617272 0.449539 (U ,U* ,Error): 0.449539 0.449539 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 (R,S): 0.401306 0.754673 (U ,U* ,Error): 0.754673 0.754673 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 (R,S): 0.797287 0.00183837 (U ,U* ,Error): 0.00183837 0.00183837 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 (R,S): 0.897504 0.350752 (U ,U* ,Error): 0.350752 0.350752 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 (R,S): 0.0945448 0.0136169 (U ,U* ,Error): 0.0136169 0.0136169 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 Interpolate R^1 * S^0 (R,S,U): 0 0 0 (R,S,U): 1 0 1 (R,S,U): 1 1 1 (R,S,U): 0 1 0 (R,S): 0.859097 0.840847 (U ,U* ,Error): 0.859097 0.859097 0 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 0 0 (R,S): 0.123104 0.00751236 (U ,U* ,Error): 0.123104 0.123104 0 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 0 0 (R,S): 0.260303 0.912484 (U ,U* ,Error): 0.260303 0.260303 0 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 0 0 (R,S): 0.113664 0.351629 (U ,U* ,Error): 0.113664 0.113664 1.38778e-17 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 0 0 (R,S): 0.822887 0.267132 (U ,U* ,Error): 0.822887 0.822887 1.11022e-16 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 0 0 Interpolate R^1 * S^1 (R,S,U): 0 0 0 (R,S,U): 1 0 0 (R,S,U): 1 1 1 (R,S,U): 0 1 0 (R,S): 0.692066 0.561662 (U ,U* ,Error): 0.388708 0.388708 0 (Ur,Ur*,Error): 0.561662 0.561662 0 (Us,Us*,Error): 0.692066 0.692066 0 (R,S): 0.861216 0.453794 (U ,U* ,Error): 0.390814 0.390814 0 (Ur,Ur*,Error): 0.453794 0.453794 0 (Us,Us*,Error): 0.861216 0.861216 0 (R,S): 0.911977 0.597917 (U ,U* ,Error): 0.545286 0.545286 0 (Ur,Ur*,Error): 0.597917 0.597917 0 (Us,Us*,Error): 0.911977 0.911977 0 (R,S): 0.188955 0.761492 (U ,U* ,Error): 0.143887 0.143887 0 (Ur,Ur*,Error): 0.761492 0.761492 0 (Us,Us*,Error): 0.188955 0.188955 0 (R,S): 0.396988 0.185314 (U ,U* ,Error): 0.0735676 0.0735676 0 (Ur,Ur*,Error): 0.185314 0.185314 0 (Us,Us*,Error): 0.396988 0.396988 0 INTERP_TEST for element "Q8". Number of nodes = 8 Interpolate R^0 * S^0 (R,S,U): 0 0 1 (R,S,U): 1 0 1 (R,S,U): 1 1 1 (R,S,U): 0 1 1 (R,S,U): 0.5 0 1 (R,S,U): 1 0.5 1 (R,S,U): 0.5 1 1 (R,S,U): 0 0.5 1 (R,S): 0.218418 0.956318 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 5.55112e-17 5.55112e-17 (Us,Us*,Error): 0 -4.44089e-16 4.44089e-16 (R,S): 0.829509 0.561695 (U ,U* ,Error): 1 1 1.11022e-16 (Ur,Ur*,Error): 0 -1.11022e-16 1.11022e-16 (Us,Us*,Error): 0 0 0 (R,S): 0.415307 0.0661187 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 1.11022e-16 1.11022e-16 (Us,Us*,Error): 0 4.44089e-16 4.44089e-16 (R,S): 0.257578 0.109957 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 -2.77556e-16 2.77556e-16 (Us,Us*,Error): 0 0 0 (R,S): 0.043829 0.633966 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 3.33067e-16 3.33067e-16 (Us,Us*,Error): 0 0 0 Interpolate R^0 * S^1 (R,S,U): 0 0 0 (R,S,U): 1 0 0 (R,S,U): 1 1 1 (R,S,U): 0 1 1 (R,S,U): 0.5 0 0 (R,S,U): 1 0.5 0.5 (R,S,U): 0.5 1 1 (R,S,U): 0 0.5 0.5 (R,S): 0.0617272 0.449539 (U ,U* ,Error): 0.449539 0.449539 5.55112e-17 (Ur,Ur*,Error): 0 1.11022e-16 1.11022e-16 (Us,Us*,Error): 1 1 0 (R,S): 0.401306 0.754673 (U ,U* ,Error): 0.754673 0.754673 1.11022e-16 (Ur,Ur*,Error): 0 1.11022e-16 1.11022e-16 (Us,Us*,Error): 1 1 0 (R,S): 0.797287 0.00183837 (U ,U* ,Error): 0.00183837 0.00183837 2.1684e-19 (Ur,Ur*,Error): 0 8.67362e-19 8.67362e-19 (Us,Us*,Error): 1 1 2.22045e-16 (R,S): 0.897504 0.350752 (U ,U* ,Error): 0.350752 0.350752 0 (Ur,Ur*,Error): 0 2.22045e-16 2.22045e-16 (Us,Us*,Error): 1 1 2.22045e-16 (R,S): 0.0945448 0.0136169 (U ,U* ,Error): 0.0136169 0.0136169 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 Interpolate R^1 * S^0 (R,S,U): 0 0 0 (R,S,U): 1 0 1 (R,S,U): 1 1 1 (R,S,U): 0 1 0 (R,S,U): 0.5 0 0.5 (R,S,U): 1 0.5 1 (R,S,U): 0.5 1 0.5 (R,S,U): 0 0.5 0 (R,S): 0.859097 0.840847 (U ,U* ,Error): 0.859097 0.859097 1.11022e-16 (Ur,Ur*,Error): 1 1 2.22045e-16 (Us,Us*,Error): 0 -2.77556e-16 2.77556e-16 (R,S): 0.123104 0.00751236 (U ,U* ,Error): 0.123104 0.123104 1.38778e-17 (Ur,Ur*,Error): 1 1 1.11022e-16 (Us,Us*,Error): 0 2.77556e-17 2.77556e-17 (R,S): 0.260303 0.912484 (U ,U* ,Error): 0.260303 0.260303 0 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 -1.11022e-16 1.11022e-16 (R,S): 0.113664 0.351629 (U ,U* ,Error): 0.113664 0.113664 1.38778e-17 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 2.77556e-17 2.77556e-17 (R,S): 0.822887 0.267132 (U ,U* ,Error): 0.822887 0.822887 0 (Ur,Ur*,Error): 1 1 2.22045e-16 (Us,Us*,Error): 0 3.33067e-16 3.33067e-16 Interpolate R^0 * S^2 (R,S,U): 0 0 0 (R,S,U): 1 0 0 (R,S,U): 1 1 1 (R,S,U): 0 1 1 (R,S,U): 0.5 0 0 (R,S,U): 1 0.5 0.25 (R,S,U): 0.5 1 1 (R,S,U): 0 0.5 0.25 (R,S): 0.692066 0.561662 (U ,U* ,Error): 0.315465 0.315465 5.55112e-17 (Ur,Ur*,Error): 0 2.77556e-17 2.77556e-17 (Us,Us*,Error): 1.12332 1.12332 2.22045e-16 (R,S): 0.861216 0.453794 (U ,U* ,Error): 0.205929 0.205929 5.55112e-17 (Ur,Ur*,Error): 0 -1.94289e-16 1.94289e-16 (Us,Us*,Error): 0.907588 0.907588 2.22045e-16 (R,S): 0.911977 0.597917 (U ,U* ,Error): 0.357505 0.357505 1.66533e-16 (Ur,Ur*,Error): 0 -2.77556e-17 2.77556e-17 (Us,Us*,Error): 1.19583 1.19583 4.44089e-16 (R,S): 0.188955 0.761492 (U ,U* ,Error): 0.57987 0.57987 1.11022e-16 (Ur,Ur*,Error): 0 -3.60822e-16 3.60822e-16 (Us,Us*,Error): 1.52298 1.52298 2.22045e-16 (R,S): 0.396988 0.185314 (U ,U* ,Error): 0.0343413 0.0343413 3.46945e-17 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0.370628 0.370628 5.55112e-17 Interpolate R^1 * S^1 (R,S,U): 0 0 0 (R,S,U): 1 0 0 (R,S,U): 1 1 1 (R,S,U): 0 1 0 (R,S,U): 0.5 0 0 (R,S,U): 1 0.5 0.5 (R,S,U): 0.5 1 0.5 (R,S,U): 0 0.5 0 (R,S): 0.574366 0.367027 (U ,U* ,Error): 0.210808 0.210808 0 (Ur,Ur*,Error): 0.367027 0.367027 5.55112e-17 (Us,Us*,Error): 0.574366 0.574366 0 (R,S): 0.617205 0.361529 (U ,U* ,Error): 0.223137 0.223137 0 (Ur,Ur*,Error): 0.361529 0.361529 0 (Us,Us*,Error): 0.617205 0.617205 1.11022e-16 (R,S): 0.21293 0.714471 (U ,U* ,Error): 0.152132 0.152132 0 (Ur,Ur*,Error): 0.714471 0.714471 2.22045e-16 (Us,Us*,Error): 0.21293 0.21293 0 (R,S): 0.117707 0.299329 (U ,U* ,Error): 0.0352331 0.0352331 6.93889e-18 (Ur,Ur*,Error): 0.299329 0.299329 0 (Us,Us*,Error): 0.117707 0.117707 1.38778e-17 (R,S): 0.825003 0.82466 (U ,U* ,Error): 0.680347 0.680347 2.22045e-16 (Ur,Ur*,Error): 0.82466 0.82466 1.11022e-16 (Us,Us*,Error): 0.825003 0.825003 0 Interpolate R^2 * S^0 (R,S,U): 0 0 0 (R,S,U): 1 0 1 (R,S,U): 1 1 1 (R,S,U): 0 1 0 (R,S,U): 0.5 0 0.25 (R,S,U): 1 0.5 1 (R,S,U): 0.5 1 0.25 (R,S,U): 0 0.5 0 (R,S): 0.0618618 0.710781 (U ,U* ,Error): 0.00382688 0.00382688 1.73472e-18 (Ur,Ur*,Error): 0.123724 0.123724 1.11022e-16 (Us,Us*,Error): 0 0 0 (R,S): 0.0882833 0.777994 (U ,U* ,Error): 0.00779395 0.00779395 5.20417e-18 (Ur,Ur*,Error): 0.176567 0.176567 1.66533e-16 (Us,Us*,Error): 0 0 0 (R,S): 0.745303 0.308675 (U ,U* ,Error): 0.555477 0.555477 1.11022e-16 (Ur,Ur*,Error): 1.49061 1.49061 0 (Us,Us*,Error): 0 2.77556e-17 2.77556e-17 (R,S): 0.899373 0.763537 (U ,U* ,Error): 0.808872 0.808872 2.22045e-16 (Ur,Ur*,Error): 1.79875 1.79875 2.22045e-16 (Us,Us*,Error): 0 -8.32667e-17 8.32667e-17 (R,S): 0.761731 0.40697 (U ,U* ,Error): 0.580234 0.580234 1.11022e-16 (Ur,Ur*,Error): 1.52346 1.52346 2.22045e-16 (Us,Us*,Error): 0 1.38778e-16 1.38778e-16 Interpolate R^1 * S^2 (R,S,U): 0 0 0 (R,S,U): 1 0 0 (R,S,U): 1 1 1 (R,S,U): 0 1 0 (R,S,U): 0.5 0 0 (R,S,U): 1 0.5 0.25 (R,S,U): 0.5 1 0.5 (R,S,U): 0 0.5 0 (R,S): 0.938749 0.562088 (U ,U* ,Error): 0.296592 0.296592 0 (Ur,Ur*,Error): 0.315943 0.315943 2.77556e-16 (Us,Us*,Error): 1.05532 1.05532 4.44089e-16 (R,S): 0.01782 0.501103 (U ,U* ,Error): 0.00447469 0.00447469 3.46945e-18 (Ur,Ur*,Error): 0.251104 0.251104 1.11022e-16 (Us,Us*,Error): 0.0178593 0.0178593 0 (R,S): 0.0419093 0.368851 (U ,U* ,Error): 0.00570179 0.00570179 6.07153e-18 (Ur,Ur*,Error): 0.136051 0.136051 8.32667e-17 (Us,Us*,Error): 0.0309165 0.0309165 3.46945e-18 (R,S): 0.271724 0.858573 (U ,U* ,Error): 0.2003 0.2003 5.55112e-17 (Ur,Ur*,Error): 0.737147 0.737147 2.22045e-16 (Us,Us*,Error): 0.466589 0.466589 0 (R,S): 0.0290366 0.0174423 (U ,U* ,Error): 8.83388e-06 8.83388e-06 1.38913e-19 (Ur,Ur*,Error): 0.000304233 0.000304233 7.75205e-18 (Us,Us*,Error): 0.00101293 0.00101293 6.93889e-18 Interpolate R^2 * S^1 (R,S,U): 0 0 0 (R,S,U): 1 0 0 (R,S,U): 1 1 1 (R,S,U): 0 1 0 (R,S,U): 0.5 0 0 (R,S,U): 1 0.5 0.5 (R,S,U): 0.5 1 0.25 (R,S,U): 0 0.5 0 (R,S): 0.152384 0.114319 (U ,U* ,Error): 0.00265457 0.00265457 1.73472e-18 (Ur,Ur*,Error): 0.0348406 0.0348406 4.16334e-17 (Us,Us*,Error): 0.0232208 0.0232208 1.38778e-17 (R,S): 0.353907 0.119308 (U ,U* ,Error): 0.0149433 0.0149433 1.38778e-17 (Ur,Ur*,Error): 0.0844478 0.0844478 1.38778e-17 (Us,Us*,Error): 0.12525 0.12525 1.11022e-16 (R,S): 0.206653 0.212924 (U ,U* ,Error): 0.00909299 0.00909299 6.93889e-18 (Ur,Ur*,Error): 0.0880026 0.0880026 5.55112e-17 (Us,Us*,Error): 0.0427054 0.0427054 1.38778e-17 (R,S): 0.612948 0.809519 (U ,U* ,Error): 0.30414 0.30414 0 (Ur,Ur*,Error): 0.992385 0.992385 1.11022e-16 (Us,Us*,Error): 0.375705 0.375705 1.11022e-16 (R,S): 0.58709 0.215492 (U ,U* ,Error): 0.0742744 0.0742744 0 (Ur,Ur*,Error): 0.253026 0.253026 0 (Us,Us*,Error): 0.344674 0.344674 1.66533e-16 INTERP_TEST for element "Q9". Number of nodes = 9 Interpolate R^0 * S^0 (R,S,U): 0 0 1 (R,S,U): 1 0 1 (R,S,U): 1 1 1 (R,S,U): 0 1 1 (R,S,U): 0.5 0 1 (R,S,U): 1 0.5 1 (R,S,U): 0.5 1 1 (R,S,U): 0 0.5 1 (R,S,U): 0.5 0.5 1 (R,S): 0.218418 0.956318 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 -1.66533e-16 1.66533e-16 (Us,Us*,Error): 0 4.44089e-16 4.44089e-16 (R,S): 0.829509 0.561695 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 4.44089e-16 4.44089e-16 (Us,Us*,Error): 0 0 0 (R,S): 0.415307 0.0661187 (U ,U* ,Error): 1 1 2.22045e-16 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0 0 0 (R,S): 0.257578 0.109957 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 2.22045e-16 2.22045e-16 (Us,Us*,Error): 0 0 0 (R,S): 0.043829 0.633966 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 8.88178e-16 8.88178e-16 (Us,Us*,Error): 0 0 0 Interpolate R^0 * S^1 (R,S,U): 0 0 0 (R,S,U): 1 0 0 (R,S,U): 1 1 1 (R,S,U): 0 1 1 (R,S,U): 0.5 0 0 (R,S,U): 1 0.5 0.5 (R,S,U): 0.5 1 1 (R,S,U): 0 0.5 0.5 (R,S,U): 0.5 0.5 0.5 (R,S): 0.0617272 0.449539 (U ,U* ,Error): 0.449539 0.449539 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 (R,S): 0.401306 0.754673 (U ,U* ,Error): 0.754673 0.754673 0 (Ur,Ur*,Error): 0 -5.55112e-17 5.55112e-17 (Us,Us*,Error): 1 1 2.22045e-16 (R,S): 0.797287 0.00183837 (U ,U* ,Error): 0.00183837 0.00183837 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 (R,S): 0.897504 0.350752 (U ,U* ,Error): 0.350752 0.350752 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 (R,S): 0.0945448 0.0136169 (U ,U* ,Error): 0.0136169 0.0136169 3.46945e-18 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 Interpolate R^1 * S^0 (R,S,U): 0 0 0 (R,S,U): 1 0 1 (R,S,U): 1 1 1 (R,S,U): 0 1 0 (R,S,U): 0.5 0 0.5 (R,S,U): 1 0.5 1 (R,S,U): 0.5 1 0.5 (R,S,U): 0 0.5 0 (R,S,U): 0.5 0.5 0.5 (R,S): 0.859097 0.840847 (U ,U* ,Error): 0.859097 0.859097 0 (Ur,Ur*,Error): 1 1 2.22045e-16 (Us,Us*,Error): 0 -3.33067e-16 3.33067e-16 (R,S): 0.123104 0.00751236 (U ,U* ,Error): 0.123104 0.123104 1.38778e-17 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 0 0 (R,S): 0.260303 0.912484 (U ,U* ,Error): 0.260303 0.260303 0 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 0 0 (R,S): 0.113664 0.351629 (U ,U* ,Error): 0.113664 0.113664 1.38778e-17 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 2.77556e-17 2.77556e-17 (R,S): 0.822887 0.267132 (U ,U* ,Error): 0.822887 0.822887 1.11022e-16 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 1.11022e-16 1.11022e-16 Interpolate R^0 * S^2 (R,S,U): 0 0 0 (R,S,U): 1 0 0 (R,S,U): 1 1 1 (R,S,U): 0 1 1 (R,S,U): 0.5 0 0 (R,S,U): 1 0.5 0.25 (R,S,U): 0.5 1 1 (R,S,U): 0 0.5 0.25 (R,S,U): 0.5 0.5 0.25 (R,S): 0.692066 0.561662 (U ,U* ,Error): 0.315465 0.315465 5.55112e-17 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1.12332 1.12332 0 (R,S): 0.861216 0.453794 (U ,U* ,Error): 0.205929 0.205929 2.77556e-17 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0.907588 0.907588 1.11022e-16 (R,S): 0.911977 0.597917 (U ,U* ,Error): 0.357505 0.357505 0 (Ur,Ur*,Error): 0 -2.22045e-16 2.22045e-16 (Us,Us*,Error): 1.19583 1.19583 2.22045e-16 (R,S): 0.188955 0.761492 (U ,U* ,Error): 0.57987 0.57987 1.11022e-16 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1.52298 1.52298 0 (R,S): 0.396988 0.185314 (U ,U* ,Error): 0.0343413 0.0343413 6.93889e-18 (Ur,Ur*,Error): 0 -2.77556e-17 2.77556e-17 (Us,Us*,Error): 0.370628 0.370628 0 Interpolate R^1 * S^1 (R,S,U): 0 0 0 (R,S,U): 1 0 0 (R,S,U): 1 1 1 (R,S,U): 0 1 0 (R,S,U): 0.5 0 0 (R,S,U): 1 0.5 0.5 (R,S,U): 0.5 1 0.5 (R,S,U): 0 0.5 0 (R,S,U): 0.5 0.5 0.25 (R,S): 0.574366 0.367027 (U ,U* ,Error): 0.210808 0.210808 0 (Ur,Ur*,Error): 0.367027 0.367027 1.11022e-16 (Us,Us*,Error): 0.574366 0.574366 1.11022e-16 (R,S): 0.617205 0.361529 (U ,U* ,Error): 0.223137 0.223137 2.77556e-17 (Ur,Ur*,Error): 0.361529 0.361529 0 (Us,Us*,Error): 0.617205 0.617205 0 (R,S): 0.21293 0.714471 (U ,U* ,Error): 0.152132 0.152132 0 (Ur,Ur*,Error): 0.714471 0.714471 1.11022e-16 (Us,Us*,Error): 0.21293 0.21293 2.77556e-17 (R,S): 0.117707 0.299329 (U ,U* ,Error): 0.0352331 0.0352331 1.38778e-17 (Ur,Ur*,Error): 0.299329 0.299329 5.55112e-17 (Us,Us*,Error): 0.117707 0.117707 1.38778e-17 (R,S): 0.825003 0.82466 (U ,U* ,Error): 0.680347 0.680347 0 (Ur,Ur*,Error): 0.82466 0.82466 1.11022e-16 (Us,Us*,Error): 0.825003 0.825003 2.22045e-16 Interpolate R^2 * S^0 (R,S,U): 0 0 0 (R,S,U): 1 0 1 (R,S,U): 1 1 1 (R,S,U): 0 1 0 (R,S,U): 0.5 0 0.25 (R,S,U): 1 0.5 1 (R,S,U): 0.5 1 0.25 (R,S,U): 0 0.5 0 (R,S,U): 0.5 0.5 0.25 (R,S): 0.0618618 0.710781 (U ,U* ,Error): 0.00382688 0.00382688 8.67362e-18 (Ur,Ur*,Error): 0.123724 0.123724 0 (Us,Us*,Error): 0 0 0 (R,S): 0.0882833 0.777994 (U ,U* ,Error): 0.00779395 0.00779395 1.21431e-17 (Ur,Ur*,Error): 0.176567 0.176567 1.11022e-16 (Us,Us*,Error): 0 -2.77556e-17 2.77556e-17 (R,S): 0.745303 0.308675 (U ,U* ,Error): 0.555477 0.555477 0 (Ur,Ur*,Error): 1.49061 1.49061 0 (Us,Us*,Error): 0 -5.55112e-17 5.55112e-17 (R,S): 0.899373 0.763537 (U ,U* ,Error): 0.808872 0.808872 1.11022e-16 (Ur,Ur*,Error): 1.79875 1.79875 4.44089e-16 (Us,Us*,Error): 0 3.05311e-16 3.05311e-16 (R,S): 0.761731 0.40697 (U ,U* ,Error): 0.580234 0.580234 0 (Ur,Ur*,Error): 1.52346 1.52346 2.22045e-16 (Us,Us*,Error): 0 2.77556e-17 2.77556e-17 Interpolate R^1 * S^2 (R,S,U): 0 0 0 (R,S,U): 1 0 0 (R,S,U): 1 1 1 (R,S,U): 0 1 0 (R,S,U): 0.5 0 0 (R,S,U): 1 0.5 0.25 (R,S,U): 0.5 1 0.5 (R,S,U): 0 0.5 0 (R,S,U): 0.5 0.5 0.125 (R,S): 0.938749 0.562088 (U ,U* ,Error): 0.296592 0.296592 0 (Ur,Ur*,Error): 0.315943 0.315943 0 (Us,Us*,Error): 1.05532 1.05532 0 (R,S): 0.01782 0.501103 (U ,U* ,Error): 0.00447469 0.00447469 0 (Ur,Ur*,Error): 0.251104 0.251104 0 (Us,Us*,Error): 0.0178593 0.0178593 3.46945e-18 (R,S): 0.0419093 0.368851 (U ,U* ,Error): 0.00570179 0.00570179 2.60209e-18 (Ur,Ur*,Error): 0.136051 0.136051 2.77556e-17 (Us,Us*,Error): 0.0309165 0.0309165 6.93889e-18 (R,S): 0.271724 0.858573 (U ,U* ,Error): 0.2003 0.2003 5.55112e-17 (Ur,Ur*,Error): 0.737147 0.737147 0 (Us,Us*,Error): 0.466589 0.466589 1.11022e-16 (R,S): 0.0290366 0.0174423 (U ,U* ,Error): 8.83388e-06 8.83388e-06 3.04932e-20 (Ur,Ur*,Error): 0.000304233 0.000304233 8.13152e-19 (Us,Us*,Error): 0.00101293 0.00101293 6.93889e-18 Interpolate R^2 * S^1 (R,S,U): 0 0 0 (R,S,U): 1 0 0 (R,S,U): 1 1 1 (R,S,U): 0 1 0 (R,S,U): 0.5 0 0 (R,S,U): 1 0.5 0.5 (R,S,U): 0.5 1 0.25 (R,S,U): 0 0.5 0 (R,S,U): 0.5 0.5 0.125 (R,S): 0.152384 0.114319 (U ,U* ,Error): 0.00265457 0.00265457 1.73472e-18 (Ur,Ur*,Error): 0.0348406 0.0348406 1.38778e-17 (Us,Us*,Error): 0.0232208 0.0232208 1.38778e-17 (R,S): 0.353907 0.119308 (U ,U* ,Error): 0.0149433 0.0149433 0 (Ur,Ur*,Error): 0.0844478 0.0844478 0 (Us,Us*,Error): 0.12525 0.12525 0 (R,S): 0.206653 0.212924 (U ,U* ,Error): 0.00909299 0.00909299 0 (Ur,Ur*,Error): 0.0880026 0.0880026 1.38778e-17 (Us,Us*,Error): 0.0427054 0.0427054 1.38778e-17 (R,S): 0.612948 0.809519 (U ,U* ,Error): 0.30414 0.30414 5.55112e-17 (Ur,Ur*,Error): 0.992385 0.992385 1.11022e-16 (Us,Us*,Error): 0.375705 0.375705 0 (R,S): 0.58709 0.215492 (U ,U* ,Error): 0.0742744 0.0742744 0 (Ur,Ur*,Error): 0.253026 0.253026 0 (Us,Us*,Error): 0.344674 0.344674 5.55112e-17 Interpolate R^2 * S^2 (R,S,U): 0 0 0 (R,S,U): 1 0 0 (R,S,U): 1 1 1 (R,S,U): 0 1 0 (R,S,U): 0.5 0 0 (R,S,U): 1 0.5 0.25 (R,S,U): 0.5 1 0.25 (R,S,U): 0 0.5 0 (R,S,U): 0.5 0.5 0.0625 (R,S): 0.768056 0.723297 (U ,U* ,Error): 0.308617 0.308617 0 (Ur,Ur*,Error): 0.80363 0.80363 0 (Us,Us*,Error): 0.853361 0.853361 0 (R,S): 0.448019 0.855176 (U ,U* ,Error): 0.146793 0.146793 2.77556e-17 (Ur,Ur*,Error): 0.655296 0.655296 0 (Us,Us*,Error): 0.343304 0.343304 5.55112e-17 (R,S): 0.945017 0.909057 (U ,U* ,Error): 0.738009 0.738009 1.11022e-16 (Ur,Ur*,Error): 1.5619 1.5619 2.22045e-16 (Us,Us*,Error): 1.62368 1.62368 4.44089e-16 (R,S): 0.519726 0.0301946 (U ,U* ,Error): 0.000246268 0.000246268 0 (Ur,Ur*,Error): 0.000947684 0.000947684 7.58942e-19 (Us,Us*,Error): 0.016312 0.016312 3.1225e-17 (R,S): 0.481067 0.292313 (U ,U* ,Error): 0.0197746 0.0197746 3.46945e-18 (Ur,Ur*,Error): 0.0822113 0.0822113 0 (Us,Us*,Error): 0.135297 0.135297 0 INTERP_TEST for element "Q12". Number of nodes = 12 Interpolate R^0 * S^0 (R,S,U): 0 0 1 (R,S,U): 0.333333 0 1 (R,S,U): 0.666667 0 1 (R,S,U): 1 0 1 (R,S,U): 0 0.333333 1 (R,S,U): 1 0.333333 1 (R,S,U): 0 0.666667 1 (R,S,U): 1 0.666667 1 (R,S,U): 0 1 1 (R,S,U): 0.333333 1 1 (R,S,U): 0.666667 1 1 (R,S,U): 1 1 1 (R,S): 0.218418 0.956318 (U ,U* ,Error): 1 1 4.44089e-16 (Ur,Ur*,Error): 0 1.05471e-15 1.05471e-15 (Us,Us*,Error): 0 1.77636e-15 1.77636e-15 (R,S): 0.829509 0.561695 (U ,U* ,Error): 1 1 4.44089e-16 (Ur,Ur*,Error): 0 1.44329e-15 1.44329e-15 (Us,Us*,Error): 0 2.72005e-15 2.72005e-15 (R,S): 0.415307 0.0661187 (U ,U* ,Error): 1 1 1.11022e-16 (Ur,Ur*,Error): 0 9.71445e-17 9.71445e-17 (Us,Us*,Error): 0 4.996e-16 4.996e-16 (R,S): 0.257578 0.109957 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 9.71445e-17 9.71445e-17 (Us,Us*,Error): 0 -6.66134e-16 6.66134e-16 (R,S): 0.043829 0.633966 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 -9.99201e-16 9.99201e-16 (Us,Us*,Error): 0 1.56472e-15 1.56472e-15 Interpolate R^0 * S^1 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0.333333 (R,S,U): 1 0.333333 0.333333 (R,S,U): 0 0.666667 0.666667 (R,S,U): 1 0.666667 0.666667 (R,S,U): 0 1 1 (R,S,U): 0.333333 1 1 (R,S,U): 0.666667 1 1 (R,S,U): 1 1 1 (R,S): 0.0617272 0.449539 (U ,U* ,Error): 0.449539 0.449539 1.66533e-16 (Ur,Ur*,Error): 0 -2.22045e-16 2.22045e-16 (Us,Us*,Error): 1 1 2.22045e-16 (R,S): 0.401306 0.754673 (U ,U* ,Error): 0.754673 0.754673 2.22045e-16 (Ur,Ur*,Error): 0 9.99201e-16 9.99201e-16 (Us,Us*,Error): 1 1 1.77636e-15 (R,S): 0.797287 0.00183837 (U ,U* ,Error): 0.00183837 0.00183837 4.33681e-19 (Ur,Ur*,Error): 0 9.54098e-18 9.54098e-18 (Us,Us*,Error): 1 1 4.44089e-16 (R,S): 0.897504 0.350752 (U ,U* ,Error): 0.350752 0.350752 0 (Ur,Ur*,Error): 0 6.66134e-16 6.66134e-16 (Us,Us*,Error): 1 1 2.22045e-16 (R,S): 0.0945448 0.0136169 (U ,U* ,Error): 0.0136169 0.0136169 0 (Ur,Ur*,Error): 0 -5.20417e-18 5.20417e-18 (Us,Us*,Error): 1 1 2.22045e-16 Interpolate R^1 * S^0 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0.333333 (R,S,U): 0.666667 0 0.666667 (R,S,U): 1 0 1 (R,S,U): 0 0.333333 0 (R,S,U): 1 0.333333 1 (R,S,U): 0 0.666667 0 (R,S,U): 1 0.666667 1 (R,S,U): 0 1 0 (R,S,U): 0.333333 1 0.333333 (R,S,U): 0.666667 1 0.666667 (R,S,U): 1 1 1 (R,S): 0.859097 0.840847 (U ,U* ,Error): 0.859097 0.859097 2.22045e-16 (Ur,Ur*,Error): 1 1 1.33227e-15 (Us,Us*,Error): 0 8.88178e-16 8.88178e-16 (R,S): 0.123104 0.00751236 (U ,U* ,Error): 0.123104 0.123104 0 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 9.02056e-17 9.02056e-17 (R,S): 0.260303 0.912484 (U ,U* ,Error): 0.260303 0.260303 0 (Ur,Ur*,Error): 1 1 6.66134e-16 (Us,Us*,Error): 0 5.55112e-16 5.55112e-16 (R,S): 0.113664 0.351629 (U ,U* ,Error): 0.113664 0.113664 0 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 8.32667e-17 8.32667e-17 (R,S): 0.822887 0.267132 (U ,U* ,Error): 0.822887 0.822887 1.11022e-16 (Ur,Ur*,Error): 1 1 3.10862e-15 (Us,Us*,Error): 0 3.33067e-16 3.33067e-16 Interpolate R^0 * S^2 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0.111111 (R,S,U): 1 0.333333 0.111111 (R,S,U): 0 0.666667 0.444444 (R,S,U): 1 0.666667 0.444444 (R,S,U): 0 1 1 (R,S,U): 0.333333 1 1 (R,S,U): 0.666667 1 1 (R,S,U): 1 1 1 (R,S): 0.692066 0.561662 (U ,U* ,Error): 0.315465 0.315465 1.11022e-16 (Ur,Ur*,Error): 0 2.27596e-15 2.27596e-15 (Us,Us*,Error): 1.12332 1.12332 8.88178e-16 (R,S): 0.861216 0.453794 (U ,U* ,Error): 0.205929 0.205929 0 (Ur,Ur*,Error): 0 2.44249e-15 2.44249e-15 (Us,Us*,Error): 0.907588 0.907588 5.55112e-16 (R,S): 0.911977 0.597917 (U ,U* ,Error): 0.357505 0.357505 2.22045e-16 (Ur,Ur*,Error): 0 6.66134e-16 6.66134e-16 (Us,Us*,Error): 1.19583 1.19583 1.33227e-15 (R,S): 0.188955 0.761492 (U ,U* ,Error): 0.57987 0.57987 0 (Ur,Ur*,Error): 0 3.33067e-16 3.33067e-16 (Us,Us*,Error): 1.52298 1.52298 1.33227e-15 (R,S): 0.396988 0.185314 (U ,U* ,Error): 0.0343413 0.0343413 3.46945e-17 (Ur,Ur*,Error): 0 8.32667e-17 8.32667e-17 (Us,Us*,Error): 0.370628 0.370628 1.11022e-16 Interpolate R^1 * S^1 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0 (R,S,U): 1 0.333333 0.333333 (R,S,U): 0 0.666667 0 (R,S,U): 1 0.666667 0.666667 (R,S,U): 0 1 0 (R,S,U): 0.333333 1 0.333333 (R,S,U): 0.666667 1 0.666667 (R,S,U): 1 1 1 (R,S): 0.574366 0.367027 (U ,U* ,Error): 0.210808 0.210808 2.77556e-17 (Ur,Ur*,Error): 0.367027 0.367027 5.55112e-17 (Us,Us*,Error): 0.574366 0.574366 2.22045e-16 (R,S): 0.617205 0.361529 (U ,U* ,Error): 0.223137 0.223137 8.32667e-17 (Ur,Ur*,Error): 0.361529 0.361529 7.21645e-16 (Us,Us*,Error): 0.617205 0.617205 2.22045e-16 (R,S): 0.21293 0.714471 (U ,U* ,Error): 0.152132 0.152132 0 (Ur,Ur*,Error): 0.714471 0.714471 5.55112e-16 (Us,Us*,Error): 0.21293 0.21293 0 (R,S): 0.117707 0.299329 (U ,U* ,Error): 0.0352331 0.0352331 6.93889e-18 (Ur,Ur*,Error): 0.299329 0.299329 2.22045e-16 (Us,Us*,Error): 0.117707 0.117707 8.32667e-17 (R,S): 0.825003 0.82466 (U ,U* ,Error): 0.680347 0.680347 1.11022e-16 (Ur,Ur*,Error): 0.82466 0.82466 2.88658e-15 (Us,Us*,Error): 0.825003 0.825003 2.22045e-16 Interpolate R^2 * S^0 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0.111111 (R,S,U): 0.666667 0 0.444444 (R,S,U): 1 0 1 (R,S,U): 0 0.333333 0 (R,S,U): 1 0.333333 1 (R,S,U): 0 0.666667 0 (R,S,U): 1 0.666667 1 (R,S,U): 0 1 0 (R,S,U): 0.333333 1 0.111111 (R,S,U): 0.666667 1 0.444444 (R,S,U): 1 1 1 (R,S): 0.0618618 0.710781 (U ,U* ,Error): 0.00382688 0.00382688 6.93889e-18 (Ur,Ur*,Error): 0.123724 0.123724 0 (Us,Us*,Error): 0 1.249e-16 1.249e-16 (R,S): 0.0882833 0.777994 (U ,U* ,Error): 0.00779395 0.00779395 1.38778e-17 (Ur,Ur*,Error): 0.176567 0.176567 3.33067e-16 (Us,Us*,Error): 0 1.38778e-16 1.38778e-16 (R,S): 0.745303 0.308675 (U ,U* ,Error): 0.555477 0.555477 0 (Ur,Ur*,Error): 1.49061 1.49061 6.66134e-16 (Us,Us*,Error): 0 -2.22045e-16 2.22045e-16 (R,S): 0.899373 0.763537 (U ,U* ,Error): 0.808872 0.808872 1.11022e-16 (Ur,Ur*,Error): 1.79875 1.79875 1.77636e-15 (Us,Us*,Error): 0 2.66454e-15 2.66454e-15 (R,S): 0.761731 0.40697 (U ,U* ,Error): 0.580234 0.580234 1.11022e-16 (Ur,Ur*,Error): 1.52346 1.52346 2.22045e-16 (Us,Us*,Error): 0 0 0 Interpolate R^0 * S^3 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0.037037 (R,S,U): 1 0.333333 0.037037 (R,S,U): 0 0.666667 0.296296 (R,S,U): 1 0.666667 0.296296 (R,S,U): 0 1 1 (R,S,U): 0.333333 1 1 (R,S,U): 0.666667 1 1 (R,S,U): 1 1 1 (R,S): 0.938749 0.562088 (U ,U* ,Error): 0.177588 0.177588 2.77556e-17 (Ur,Ur*,Error): 0 -1.11022e-15 1.11022e-15 (Us,Us*,Error): 0.94783 0.94783 1.11022e-16 (R,S): 0.01782 0.501103 (U ,U* ,Error): 0.125829 0.125829 2.77556e-17 (Ur,Ur*,Error): 0 -3.05311e-16 3.05311e-16 (Us,Us*,Error): 0.753313 0.753313 3.33067e-16 (R,S): 0.0419093 0.368851 (U ,U* ,Error): 0.0501824 0.0501824 3.46945e-17 (Ur,Ur*,Error): 0 8.32667e-17 8.32667e-17 (Us,Us*,Error): 0.408152 0.408152 2.77556e-16 (R,S): 0.271724 0.858573 (U ,U* ,Error): 0.632894 0.632894 0 (Ur,Ur*,Error): 0 2.22045e-16 2.22045e-16 (Us,Us*,Error): 2.21144 2.21144 0 (R,S): 0.0290366 0.0174423 (U ,U* ,Error): 5.30652e-06 5.30652e-06 8.96161e-19 (Ur,Ur*,Error): 0 1.73472e-17 1.73472e-17 (Us,Us*,Error): 0.000912699 0.000912699 7.87131e-17 Interpolate R^1 * S^2 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0 (R,S,U): 1 0.333333 0.111111 (R,S,U): 0 0.666667 0 (R,S,U): 1 0.666667 0.444444 (R,S,U): 0 1 0 (R,S,U): 0.333333 1 0.333333 (R,S,U): 0.666667 1 0.666667 (R,S,U): 1 1 1 (R,S): 0.152384 0.114319 (U ,U* ,Error): 0.00199147 0.00199147 7.80626e-18 (Ur,Ur*,Error): 0.0130688 0.0130688 3.98986e-17 (Us,Us*,Error): 0.0348406 0.0348406 5.55112e-17 (R,S): 0.353907 0.119308 (U ,U* ,Error): 0.00503764 0.00503764 2.60209e-18 (Ur,Ur*,Error): 0.0142344 0.0142344 7.45931e-17 (Us,Us*,Error): 0.0844478 0.0844478 2.77556e-17 (R,S): 0.206653 0.212924 (U ,U* ,Error): 0.00936894 0.00936894 1.04083e-17 (Ur,Ur*,Error): 0.0453366 0.0453366 3.46945e-17 (Us,Us*,Error): 0.0880026 0.0880026 1.249e-16 (R,S): 0.612948 0.809519 (U ,U* ,Error): 0.401677 0.401677 1.11022e-16 (Ur,Ur*,Error): 0.655321 0.655321 1.33227e-15 (Us,Us*,Error): 0.992385 0.992385 6.66134e-16 (R,S): 0.58709 0.215492 (U ,U* ,Error): 0.0272625 0.0272625 1.04083e-17 (Ur,Ur*,Error): 0.0464366 0.0464366 2.70617e-16 (Us,Us*,Error): 0.253026 0.253026 1.11022e-16 Interpolate R^2 * S^1 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0 (R,S,U): 1 0.333333 0.333333 (R,S,U): 0 0.666667 0 (R,S,U): 1 0.666667 0.666667 (R,S,U): 0 1 0 (R,S,U): 0.333333 1 0.111111 (R,S,U): 0.666667 1 0.444444 (R,S,U): 1 1 1 (R,S): 0.768056 0.723297 (U ,U* ,Error): 0.42668 0.42668 5.55112e-17 (Ur,Ur*,Error): 1.11107 1.11107 2.22045e-16 (Us,Us*,Error): 0.589911 0.589911 6.66134e-16 (R,S): 0.448019 0.855176 (U ,U* ,Error): 0.171652 0.171652 2.77556e-17 (Ur,Ur*,Error): 0.76627 0.76627 2.22045e-16 (Us,Us*,Error): 0.200721 0.200721 8.60423e-16 (R,S): 0.945017 0.909057 (U ,U* ,Error): 0.811841 0.811841 0 (Ur,Ur*,Error): 1.71815 1.71815 1.33227e-15 (Us,Us*,Error): 0.893058 0.893058 0 (R,S): 0.519726 0.0301946 (U ,U* ,Error): 0.00815602 0.00815602 1.56125e-17 (Ur,Ur*,Error): 0.0313858 0.0313858 2.77556e-17 (Us,Us*,Error): 0.270115 0.270115 1.66533e-16 (R,S): 0.481067 0.292313 (U ,U* ,Error): 0.0676486 0.0676486 2.77556e-17 (Ur,Ur*,Error): 0.281244 0.281244 5.55112e-17 (Us,Us*,Error): 0.231425 0.231425 8.32667e-17 Interpolate R^3 * S^0 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0.037037 (R,S,U): 0.666667 0 0.296296 (R,S,U): 1 0 1 (R,S,U): 0 0.333333 0 (R,S,U): 1 0.333333 1 (R,S,U): 0 0.666667 0 (R,S,U): 1 0.666667 1 (R,S,U): 0 1 0 (R,S,U): 0.333333 1 0.037037 (R,S,U): 0.666667 1 0.296296 (R,S,U): 1 1 1 (R,S): 0.90264 0.667842 (U ,U* ,Error): 0.735434 0.735434 1.11022e-16 (Ur,Ur*,Error): 2.44428 2.44428 4.44089e-16 (Us,Us*,Error): 0 5.55112e-16 5.55112e-16 (R,S): 0.412278 0.156948 (U ,U* ,Error): 0.0700762 0.0700762 5.55112e-17 (Ur,Ur*,Error): 0.50992 0.50992 1.11022e-16 (Us,Us*,Error): 0 5.55112e-17 5.55112e-17 (R,S): 0.833282 0.964404 (U ,U* ,Error): 0.578596 0.578596 0 (Ur,Ur*,Error): 2.08307 2.08307 8.88178e-16 (Us,Us*,Error): 0 1.77636e-15 1.77636e-15 (R,S): 0.74079 0.456099 (U ,U* ,Error): 0.406523 0.406523 2.77556e-16 (Ur,Ur*,Error): 1.64631 1.64631 4.44089e-16 (Us,Us*,Error): 0 -6.66134e-16 6.66134e-16 (R,S): 0.653561 0.406827 (U ,U* ,Error): 0.279164 0.279164 1.66533e-16 (Ur,Ur*,Error): 1.28143 1.28143 2.22045e-16 (Us,Us*,Error): 0 1.77636e-15 1.77636e-15 Interpolate R^1 * S^3 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0 (R,S,U): 1 0.333333 0.037037 (R,S,U): 0 0.666667 0 (R,S,U): 1 0.666667 0.296296 (R,S,U): 0 1 0 (R,S,U): 0.333333 1 0.333333 (R,S,U): 0.666667 1 0.666667 (R,S,U): 1 1 1 (R,S): 0.540539 0.832281 (U ,U* ,Error): 0.311628 0.311628 5.55112e-17 (Ur,Ur*,Error): 0.576514 0.576514 6.66134e-16 (Us,Us*,Error): 1.12328 1.12328 0 (R,S): 0.145756 0.717128 (U ,U* ,Error): 0.0537547 0.0537547 4.16334e-17 (Ur,Ur*,Error): 0.3688 0.3688 4.44089e-16 (Us,Us*,Error): 0.224875 0.224875 5.55112e-17 (R,S): 0.775651 0.362262 (U ,U* ,Error): 0.0368751 0.0368751 5.55112e-17 (Ur,Ur*,Error): 0.0475409 0.0475409 2.498e-16 (Us,Us*,Error): 0.305374 0.305374 3.33067e-16 (R,S): 0.531111 0.379977 (U ,U* ,Error): 0.0291378 0.0291378 3.46945e-17 (Ur,Ur*,Error): 0.0548619 0.0548619 1.52656e-16 (Us,Us*,Error): 0.230049 0.230049 3.05311e-16 (R,S): 0.269285 0.877418 (U ,U* ,Error): 0.1819 0.1819 5.55112e-17 (Ur,Ur*,Error): 0.675491 0.675491 8.88178e-16 (Us,Us*,Error): 0.621937 0.621937 0 Interpolate R^3 * S^1 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0 (R,S,U): 1 0.333333 0.333333 (R,S,U): 0 0.666667 0 (R,S,U): 1 0.666667 0.666667 (R,S,U): 0 1 0 (R,S,U): 0.333333 1 0.037037 (R,S,U): 0.666667 1 0.296296 (R,S,U): 1 1 1 (R,S): 0.761285 0.913675 (U ,U* ,Error): 0.403119 0.403119 5.55112e-17 (Ur,Ur*,Error): 1.58857 1.58857 4.44089e-16 (Us,Us*,Error): 0.441206 0.441206 3.88578e-16 (R,S): 0.135794 0.291195 (U ,U* ,Error): 0.000729167 0.000729167 9.75782e-18 (Ur,Ur*,Error): 0.0161089 0.0161089 2.77556e-17 (Us,Us*,Error): 0.00250405 0.00250405 1.29671e-16 (R,S): 0.118211 0.771781 (U ,U* ,Error): 0.00127487 0.00127487 2.1684e-19 (Ur,Ur*,Error): 0.0323542 0.0323542 2.98372e-16 (Us,Us*,Error): 0.00165186 0.00165186 9.75782e-18 (R,S): 0.320633 0.877413 (U ,U* ,Error): 0.028922 0.028922 6.93889e-18 (Ur,Ur*,Error): 0.270609 0.270609 1.66533e-16 (Us,Us*,Error): 0.0329628 0.0329628 2.08167e-17 (R,S): 0.679476 0.949176 (U ,U* ,Error): 0.297761 0.297761 5.55112e-17 (Ur,Ur*,Error): 1.31467 1.31467 0 (Us,Us*,Error): 0.313705 0.313705 0 INTERP_TEST for element "Q16". Number of nodes = 16 Interpolate R^0 * S^0 (R,S,U): 0 0 1 (R,S,U): 0.333333 0 1 (R,S,U): 0.666667 0 1 (R,S,U): 1 0 1 (R,S,U): 0 0.333333 1 (R,S,U): 0.333333 0.333333 1 (R,S,U): 0.666667 0.333333 1 (R,S,U): 1 0.333333 1 (R,S,U): 0 0.666667 1 (R,S,U): 0.333333 0.666667 1 (R,S,U): 0.666667 0.666667 1 (R,S,U): 1 0.666667 1 (R,S,U): 0 1 1 (R,S,U): 0.333333 1 1 (R,S,U): 0.666667 1 1 (R,S,U): 1 1 1 (R,S): 0.218418 0.956318 (U ,U* ,Error): 1 1 3.33067e-16 (Ur,Ur*,Error): 0 -9.99201e-16 9.99201e-16 (Us,Us*,Error): 0 6.18949e-15 6.18949e-15 (R,S): 0.829509 0.561695 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 -5.4956e-15 5.4956e-15 (Us,Us*,Error): 0 -5.34295e-16 5.34295e-16 (R,S): 0.415307 0.0661187 (U ,U* ,Error): 1 1 2.22045e-16 (Ur,Ur*,Error): 0 -4.51028e-16 4.51028e-16 (Us,Us*,Error): 0 1.02349e-15 1.02349e-15 (R,S): 0.257578 0.109957 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 -2.91434e-16 2.91434e-16 (Us,Us*,Error): 0 -1.96544e-15 1.96544e-15 (R,S): 0.043829 0.633966 (U ,U* ,Error): 1 1 2.22045e-16 (Ur,Ur*,Error): 0 -2.08167e-15 2.08167e-15 (Us,Us*,Error): 0 -1.85615e-15 1.85615e-15 Interpolate R^0 * S^1 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0.333333 (R,S,U): 0.333333 0.333333 0.333333 (R,S,U): 0.666667 0.333333 0.333333 (R,S,U): 1 0.333333 0.333333 (R,S,U): 0 0.666667 0.666667 (R,S,U): 0.333333 0.666667 0.666667 (R,S,U): 0.666667 0.666667 0.666667 (R,S,U): 1 0.666667 0.666667 (R,S,U): 0 1 1 (R,S,U): 0.333333 1 1 (R,S,U): 0.666667 1 1 (R,S,U): 1 1 1 (R,S): 0.0617272 0.449539 (U ,U* ,Error): 0.449539 0.449539 5.55112e-17 (Ur,Ur*,Error): 0 -5.6552e-16 5.6552e-16 (Us,Us*,Error): 1 1 8.88178e-16 (R,S): 0.401306 0.754673 (U ,U* ,Error): 0.754673 0.754673 0 (Ur,Ur*,Error): 0 1.02002e-15 1.02002e-15 (Us,Us*,Error): 1 1 0 (R,S): 0.797287 0.00183837 (U ,U* ,Error): 0.00183837 0.00183837 1.51788e-18 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 2.22045e-16 (R,S): 0.897504 0.350752 (U ,U* ,Error): 0.350752 0.350752 1.11022e-16 (Ur,Ur*,Error): 0 -1.70003e-15 1.70003e-15 (Us,Us*,Error): 1 1 3.33067e-16 (R,S): 0.0945448 0.0136169 (U ,U* ,Error): 0.0136169 0.0136169 0 (Ur,Ur*,Error): 0 3.51282e-17 3.51282e-17 (Us,Us*,Error): 1 1 2.22045e-16 Interpolate R^1 * S^0 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0.333333 (R,S,U): 0.666667 0 0.666667 (R,S,U): 1 0 1 (R,S,U): 0 0.333333 0 (R,S,U): 0.333333 0.333333 0.333333 (R,S,U): 0.666667 0.333333 0.666667 (R,S,U): 1 0.333333 1 (R,S,U): 0 0.666667 0 (R,S,U): 0.333333 0.666667 0.333333 (R,S,U): 0.666667 0.666667 0.666667 (R,S,U): 1 0.666667 1 (R,S,U): 0 1 0 (R,S,U): 0.333333 1 0.333333 (R,S,U): 0.666667 1 0.666667 (R,S,U): 1 1 1 (R,S): 0.859097 0.840847 (U ,U* ,Error): 0.859097 0.859097 1.11022e-16 (Ur,Ur*,Error): 1 1 2.22045e-15 (Us,Us*,Error): 0 2.22045e-16 2.22045e-16 (R,S): 0.123104 0.00751236 (U ,U* ,Error): 0.123104 0.123104 5.55112e-17 (Ur,Ur*,Error): 1 1 3.33067e-16 (Us,Us*,Error): 0 2.01228e-16 2.01228e-16 (R,S): 0.260303 0.912484 (U ,U* ,Error): 0.260303 0.260303 0 (Ur,Ur*,Error): 1 1 2.22045e-16 (Us,Us*,Error): 0 -1.41553e-15 1.41553e-15 (R,S): 0.113664 0.351629 (U ,U* ,Error): 0.113664 0.113664 5.55112e-17 (Ur,Ur*,Error): 1 1 1.11022e-16 (Us,Us*,Error): 0 -3.57353e-16 3.57353e-16 (R,S): 0.822887 0.267132 (U ,U* ,Error): 0.822887 0.822887 2.22045e-16 (Ur,Ur*,Error): 1 1 1.11022e-16 (Us,Us*,Error): 0 -8.32667e-17 8.32667e-17 Interpolate R^0 * S^2 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0.111111 (R,S,U): 0.333333 0.333333 0.111111 (R,S,U): 0.666667 0.333333 0.111111 (R,S,U): 1 0.333333 0.111111 (R,S,U): 0 0.666667 0.444444 (R,S,U): 0.333333 0.666667 0.444444 (R,S,U): 0.666667 0.666667 0.444444 (R,S,U): 1 0.666667 0.444444 (R,S,U): 0 1 1 (R,S,U): 0.333333 1 1 (R,S,U): 0.666667 1 1 (R,S,U): 1 1 1 (R,S): 0.692066 0.561662 (U ,U* ,Error): 0.315465 0.315465 5.55112e-17 (Ur,Ur*,Error): 0 -9.02056e-16 9.02056e-16 (Us,Us*,Error): 1.12332 1.12332 4.44089e-16 (R,S): 0.861216 0.453794 (U ,U* ,Error): 0.205929 0.205929 0 (Ur,Ur*,Error): 0 -3.05311e-16 3.05311e-16 (Us,Us*,Error): 0.907588 0.907588 1.11022e-16 (R,S): 0.911977 0.597917 (U ,U* ,Error): 0.357505 0.357505 5.55112e-17 (Ur,Ur*,Error): 0 2.35922e-15 2.35922e-15 (Us,Us*,Error): 1.19583 1.19583 4.44089e-16 (R,S): 0.188955 0.761492 (U ,U* ,Error): 0.57987 0.57987 0 (Ur,Ur*,Error): 0 -3.85109e-16 3.85109e-16 (Us,Us*,Error): 1.52298 1.52298 6.66134e-16 (R,S): 0.396988 0.185314 (U ,U* ,Error): 0.0343413 0.0343413 2.77556e-17 (Ur,Ur*,Error): 0 -8.32667e-17 8.32667e-17 (Us,Us*,Error): 0.370628 0.370628 1.66533e-16 Interpolate R^1 * S^1 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0 (R,S,U): 0.333333 0.333333 0.111111 (R,S,U): 0.666667 0.333333 0.222222 (R,S,U): 1 0.333333 0.333333 (R,S,U): 0 0.666667 0 (R,S,U): 0.333333 0.666667 0.222222 (R,S,U): 0.666667 0.666667 0.444444 (R,S,U): 1 0.666667 0.666667 (R,S,U): 0 1 0 (R,S,U): 0.333333 1 0.333333 (R,S,U): 0.666667 1 0.666667 (R,S,U): 1 1 1 (R,S): 0.574366 0.367027 (U ,U* ,Error): 0.210808 0.210808 5.55112e-17 (Ur,Ur*,Error): 0.367027 0.367027 0 (Us,Us*,Error): 0.574366 0.574366 2.22045e-16 (R,S): 0.617205 0.361529 (U ,U* ,Error): 0.223137 0.223137 0 (Ur,Ur*,Error): 0.361529 0.361529 2.22045e-16 (Us,Us*,Error): 0.617205 0.617205 4.44089e-16 (R,S): 0.21293 0.714471 (U ,U* ,Error): 0.152132 0.152132 0 (Ur,Ur*,Error): 0.714471 0.714471 1.11022e-16 (Us,Us*,Error): 0.21293 0.21293 2.498e-16 (R,S): 0.117707 0.299329 (U ,U* ,Error): 0.0352331 0.0352331 1.38778e-17 (Ur,Ur*,Error): 0.299329 0.299329 0 (Us,Us*,Error): 0.117707 0.117707 1.38778e-17 (R,S): 0.825003 0.82466 (U ,U* ,Error): 0.680347 0.680347 1.11022e-16 (Ur,Ur*,Error): 0.82466 0.82466 1.66533e-15 (Us,Us*,Error): 0.825003 0.825003 5.55112e-16 Interpolate R^2 * S^0 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0.111111 (R,S,U): 0.666667 0 0.444444 (R,S,U): 1 0 1 (R,S,U): 0 0.333333 0 (R,S,U): 0.333333 0.333333 0.111111 (R,S,U): 0.666667 0.333333 0.444444 (R,S,U): 1 0.333333 1 (R,S,U): 0 0.666667 0 (R,S,U): 0.333333 0.666667 0.111111 (R,S,U): 0.666667 0.666667 0.444444 (R,S,U): 1 0.666667 1 (R,S,U): 0 1 0 (R,S,U): 0.333333 1 0.111111 (R,S,U): 0.666667 1 0.444444 (R,S,U): 1 1 1 (R,S): 0.0618618 0.710781 (U ,U* ,Error): 0.00382688 0.00382688 1.38778e-17 (Ur,Ur*,Error): 0.123724 0.123724 1.38778e-16 (Us,Us*,Error): 0 0 0 (R,S): 0.0882833 0.777994 (U ,U* ,Error): 0.00779395 0.00779395 1.38778e-17 (Ur,Ur*,Error): 0.176567 0.176567 1.94289e-16 (Us,Us*,Error): 0 -2.77556e-17 2.77556e-17 (R,S): 0.745303 0.308675 (U ,U* ,Error): 0.555477 0.555477 1.11022e-16 (Ur,Ur*,Error): 1.49061 1.49061 2.22045e-16 (Us,Us*,Error): 0 -6.10623e-16 6.10623e-16 (R,S): 0.899373 0.763537 (U ,U* ,Error): 0.808872 0.808872 1.11022e-16 (Ur,Ur*,Error): 1.79875 1.79875 0 (Us,Us*,Error): 0 1.77636e-15 1.77636e-15 (R,S): 0.761731 0.40697 (U ,U* ,Error): 0.580234 0.580234 0 (Ur,Ur*,Error): 1.52346 1.52346 4.44089e-16 (Us,Us*,Error): 0 1.04083e-16 1.04083e-16 Interpolate R^0 * S^3 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0.037037 (R,S,U): 0.333333 0.333333 0.037037 (R,S,U): 0.666667 0.333333 0.037037 (R,S,U): 1 0.333333 0.037037 (R,S,U): 0 0.666667 0.296296 (R,S,U): 0.333333 0.666667 0.296296 (R,S,U): 0.666667 0.666667 0.296296 (R,S,U): 1 0.666667 0.296296 (R,S,U): 0 1 1 (R,S,U): 0.333333 1 1 (R,S,U): 0.666667 1 1 (R,S,U): 1 1 1 (R,S): 0.938749 0.562088 (U ,U* ,Error): 0.177588 0.177588 5.55112e-17 (Ur,Ur*,Error): 0 2.77556e-16 2.77556e-16 (Us,Us*,Error): 0.94783 0.94783 2.22045e-16 (R,S): 0.01782 0.501103 (U ,U* ,Error): 0.125829 0.125829 0 (Ur,Ur*,Error): 0 -2.91434e-16 2.91434e-16 (Us,Us*,Error): 0.753313 0.753313 5.55112e-16 (R,S): 0.0419093 0.368851 (U ,U* ,Error): 0.0501824 0.0501824 6.93889e-18 (Ur,Ur*,Error): 0 2.08167e-17 2.08167e-17 (Us,Us*,Error): 0.408152 0.408152 5.55112e-17 (R,S): 0.271724 0.858573 (U ,U* ,Error): 0.632894 0.632894 1.11022e-16 (Ur,Ur*,Error): 0 -9.4369e-16 9.4369e-16 (Us,Us*,Error): 2.21144 2.21144 4.44089e-16 (R,S): 0.0290366 0.0174423 (U ,U* ,Error): 5.30652e-06 5.30652e-06 5.66326e-18 (Ur,Ur*,Error): 0 -1.56125e-17 1.56125e-17 (Us,Us*,Error): 0.000912699 0.000912699 6.70037e-17 Interpolate R^1 * S^2 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0 (R,S,U): 0.333333 0.333333 0.037037 (R,S,U): 0.666667 0.333333 0.0740741 (R,S,U): 1 0.333333 0.111111 (R,S,U): 0 0.666667 0 (R,S,U): 0.333333 0.666667 0.148148 (R,S,U): 0.666667 0.666667 0.296296 (R,S,U): 1 0.666667 0.444444 (R,S,U): 0 1 0 (R,S,U): 0.333333 1 0.333333 (R,S,U): 0.666667 1 0.666667 (R,S,U): 1 1 1 (R,S): 0.152384 0.114319 (U ,U* ,Error): 0.00199147 0.00199147 1.30104e-17 (Ur,Ur*,Error): 0.0130688 0.0130688 2.08167e-17 (Us,Us*,Error): 0.0348406 0.0348406 4.16334e-17 (R,S): 0.353907 0.119308 (U ,U* ,Error): 0.00503764 0.00503764 1.30104e-17 (Ur,Ur*,Error): 0.0142344 0.0142344 5.20417e-18 (Us,Us*,Error): 0.0844478 0.0844478 2.77556e-17 (R,S): 0.206653 0.212924 (U ,U* ,Error): 0.00936894 0.00936894 6.93889e-18 (Ur,Ur*,Error): 0.0453366 0.0453366 1.38778e-17 (Us,Us*,Error): 0.0880026 0.0880026 6.93889e-17 (R,S): 0.612948 0.809519 (U ,U* ,Error): 0.401677 0.401677 5.55112e-17 (Ur,Ur*,Error): 0.655321 0.655321 6.66134e-16 (Us,Us*,Error): 0.992385 0.992385 3.33067e-16 (R,S): 0.58709 0.215492 (U ,U* ,Error): 0.0272625 0.0272625 3.46945e-18 (Ur,Ur*,Error): 0.0464366 0.0464366 4.16334e-17 (Us,Us*,Error): 0.253026 0.253026 5.55112e-17 Interpolate R^2 * S^1 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0 (R,S,U): 0.333333 0.333333 0.037037 (R,S,U): 0.666667 0.333333 0.148148 (R,S,U): 1 0.333333 0.333333 (R,S,U): 0 0.666667 0 (R,S,U): 0.333333 0.666667 0.0740741 (R,S,U): 0.666667 0.666667 0.296296 (R,S,U): 1 0.666667 0.666667 (R,S,U): 0 1 0 (R,S,U): 0.333333 1 0.111111 (R,S,U): 0.666667 1 0.444444 (R,S,U): 1 1 1 (R,S): 0.768056 0.723297 (U ,U* ,Error): 0.42668 0.42668 0 (Ur,Ur*,Error): 1.11107 1.11107 2.22045e-16 (Us,Us*,Error): 0.589911 0.589911 1.11022e-16 (R,S): 0.448019 0.855176 (U ,U* ,Error): 0.171652 0.171652 2.77556e-17 (Ur,Ur*,Error): 0.76627 0.76627 5.55112e-16 (Us,Us*,Error): 0.200721 0.200721 2.77556e-17 (R,S): 0.945017 0.909057 (U ,U* ,Error): 0.811841 0.811841 3.33067e-16 (Ur,Ur*,Error): 1.71815 1.71815 6.66134e-16 (Us,Us*,Error): 0.893058 0.893058 1.9984e-15 (R,S): 0.519726 0.0301946 (U ,U* ,Error): 0.00815602 0.00815602 1.73472e-18 (Ur,Ur*,Error): 0.0313858 0.0313858 2.08167e-17 (Us,Us*,Error): 0.270115 0.270115 1.11022e-16 (R,S): 0.481067 0.292313 (U ,U* ,Error): 0.0676486 0.0676486 1.38778e-17 (Ur,Ur*,Error): 0.281244 0.281244 0 (Us,Us*,Error): 0.231425 0.231425 0 Interpolate R^3 * S^0 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0.037037 (R,S,U): 0.666667 0 0.296296 (R,S,U): 1 0 1 (R,S,U): 0 0.333333 0 (R,S,U): 0.333333 0.333333 0.037037 (R,S,U): 0.666667 0.333333 0.296296 (R,S,U): 1 0.333333 1 (R,S,U): 0 0.666667 0 (R,S,U): 0.333333 0.666667 0.037037 (R,S,U): 0.666667 0.666667 0.296296 (R,S,U): 1 0.666667 1 (R,S,U): 0 1 0 (R,S,U): 0.333333 1 0.037037 (R,S,U): 0.666667 1 0.296296 (R,S,U): 1 1 1 (R,S): 0.90264 0.667842 (U ,U* ,Error): 0.735434 0.735434 0 (Ur,Ur*,Error): 2.44428 2.44428 4.44089e-16 (Us,Us*,Error): 0 5.55112e-16 5.55112e-16 (R,S): 0.412278 0.156948 (U ,U* ,Error): 0.0700762 0.0700762 1.38778e-17 (Ur,Ur*,Error): 0.50992 0.50992 1.11022e-16 (Us,Us*,Error): 0 -2.47198e-17 2.47198e-17 (R,S): 0.833282 0.964404 (U ,U* ,Error): 0.578596 0.578596 2.22045e-16 (Ur,Ur*,Error): 2.08307 2.08307 8.88178e-16 (Us,Us*,Error): 0 -1.77636e-15 1.77636e-15 (R,S): 0.74079 0.456099 (U ,U* ,Error): 0.406523 0.406523 0 (Ur,Ur*,Error): 1.64631 1.64631 2.22045e-16 (Us,Us*,Error): 0 -5.51642e-16 5.51642e-16 (R,S): 0.653561 0.406827 (U ,U* ,Error): 0.279164 0.279164 5.55112e-17 (Ur,Ur*,Error): 1.28143 1.28143 6.66134e-16 (Us,Us*,Error): 0 -8.42208e-16 8.42208e-16 Interpolate R^1 * S^3 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0 (R,S,U): 0.333333 0.333333 0.0123457 (R,S,U): 0.666667 0.333333 0.0246914 (R,S,U): 1 0.333333 0.037037 (R,S,U): 0 0.666667 0 (R,S,U): 0.333333 0.666667 0.0987654 (R,S,U): 0.666667 0.666667 0.197531 (R,S,U): 1 0.666667 0.296296 (R,S,U): 0 1 0 (R,S,U): 0.333333 1 0.333333 (R,S,U): 0.666667 1 0.666667 (R,S,U): 1 1 1 (R,S): 0.540539 0.832281 (U ,U* ,Error): 0.311628 0.311628 5.55112e-17 (Ur,Ur*,Error): 0.576514 0.576514 6.66134e-16 (Us,Us*,Error): 1.12328 1.12328 2.22045e-16 (R,S): 0.145756 0.717128 (U ,U* ,Error): 0.0537547 0.0537547 2.08167e-17 (Ur,Ur*,Error): 0.3688 0.3688 0 (Us,Us*,Error): 0.224875 0.224875 1.11022e-16 (R,S): 0.775651 0.362262 (U ,U* ,Error): 0.0368751 0.0368751 0 (Ur,Ur*,Error): 0.0475409 0.0475409 1.11022e-16 (Us,Us*,Error): 0.305374 0.305374 4.44089e-16 (R,S): 0.531111 0.379977 (U ,U* ,Error): 0.0291378 0.0291378 6.93889e-18 (Ur,Ur*,Error): 0.0548619 0.0548619 6.245e-17 (Us,Us*,Error): 0.230049 0.230049 1.94289e-16 (R,S): 0.269285 0.877418 (U ,U* ,Error): 0.1819 0.1819 0 (Ur,Ur*,Error): 0.675491 0.675491 1.11022e-16 (Us,Us*,Error): 0.621937 0.621937 1.11022e-16 Interpolate R^2 * S^2 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0 (R,S,U): 0.333333 0.333333 0.0123457 (R,S,U): 0.666667 0.333333 0.0493827 (R,S,U): 1 0.333333 0.111111 (R,S,U): 0 0.666667 0 (R,S,U): 0.333333 0.666667 0.0493827 (R,S,U): 0.666667 0.666667 0.197531 (R,S,U): 1 0.666667 0.444444 (R,S,U): 0 1 0 (R,S,U): 0.333333 1 0.111111 (R,S,U): 0.666667 1 0.444444 (R,S,U): 1 1 1 (R,S): 0.761285 0.913675 (U ,U* ,Error): 0.483813 0.483813 0 (Ur,Ur*,Error): 1.27104 1.27104 2.22045e-16 (Us,Us*,Error): 1.05905 1.05905 4.44089e-16 (R,S): 0.135794 0.291195 (U ,U* ,Error): 0.00156362 0.00156362 1.95156e-18 (Ur,Ur*,Error): 0.0230292 0.0230292 3.46945e-18 (Us,Us*,Error): 0.0107393 0.0107393 1.73472e-17 (R,S): 0.118211 0.771781 (U ,U* ,Error): 0.00832345 0.00832345 1.38778e-17 (Ur,Ur*,Error): 0.140824 0.140824 2.77556e-17 (Us,Us*,Error): 0.0215695 0.0215695 4.16334e-17 (R,S): 0.320633 0.877413 (U ,U* ,Error): 0.0791451 0.0791451 1.38778e-17 (Ur,Ur*,Error): 0.493681 0.493681 2.77556e-16 (Us,Us*,Error): 0.180406 0.180406 8.32667e-17 (R,S): 0.679476 0.949176 (U ,U* ,Error): 0.41595 0.41595 5.55112e-17 (Ur,Ur*,Error): 1.22433 1.22433 6.66134e-16 (Us,Us*,Error): 0.876445 0.876445 5.55112e-16 Interpolate R^3 * S^1 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0 (R,S,U): 0.333333 0.333333 0.0123457 (R,S,U): 0.666667 0.333333 0.0987654 (R,S,U): 1 0.333333 0.333333 (R,S,U): 0 0.666667 0 (R,S,U): 0.333333 0.666667 0.0246914 (R,S,U): 0.666667 0.666667 0.197531 (R,S,U): 1 0.666667 0.666667 (R,S,U): 0 1 0 (R,S,U): 0.333333 1 0.037037 (R,S,U): 0.666667 1 0.296296 (R,S,U): 1 1 1 (R,S): 0.793361 0.0213662 (U ,U* ,Error): 0.0106694 0.0106694 0 (Ur,Ur*,Error): 0.0403451 0.0403451 2.77556e-17 (Us,Us*,Error): 0.499359 0.499359 2.22045e-16 (R,S): 0.102278 0.991406 (U ,U* ,Error): 0.00106072 0.00106072 2.25514e-17 (Ur,Ur*,Error): 0.0311129 0.0311129 1.73472e-16 (Us,Us*,Error): 0.00106992 0.00106992 3.29597e-17 (R,S): 0.56326 0.704528 (U ,U* ,Error): 0.125899 0.125899 0 (Ur,Ur*,Error): 0.670558 0.670558 3.33067e-16 (Us,Us*,Error): 0.178701 0.178701 2.77556e-17 (R,S): 0.994146 0.604334 (U ,U* ,Error): 0.593782 0.593782 0 (Ur,Ur*,Error): 1.79184 1.79184 1.11022e-15 (Us,Us*,Error): 0.982539 0.982539 8.88178e-16 (R,S): 0.0461452 0.562029 (U ,U* ,Error): 5.52253e-05 5.52253e-05 2.98156e-19 (Ur,Ur*,Error): 0.00359032 0.00359032 7.32921e-17 (Us,Us*,Error): 9.82605e-05 9.82605e-05 1.87025e-17 Interpolate R^2 * S^3 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0 (R,S,U): 0.333333 0.333333 0.00411523 (R,S,U): 0.666667 0.333333 0.0164609 (R,S,U): 1 0.333333 0.037037 (R,S,U): 0 0.666667 0 (R,S,U): 0.333333 0.666667 0.0329218 (R,S,U): 0.666667 0.666667 0.131687 (R,S,U): 1 0.666667 0.296296 (R,S,U): 0 1 0 (R,S,U): 0.333333 1 0.111111 (R,S,U): 0.666667 1 0.444444 (R,S,U): 1 1 1 (R,S): 0.0215154 0.609135 (U ,U* ,Error): 0.000104626 0.000104626 4.74338e-19 (Ur,Ur*,Error): 0.00972567 0.00972567 2.77556e-17 (Us,Us*,Error): 0.000515284 0.000515284 6.50521e-19 (R,S): 0.728789 0.748407 (U ,U* ,Error): 0.222646 0.222646 5.55112e-17 (Ur,Ur*,Error): 0.611004 0.611004 2.22045e-16 (Us,Us*,Error): 0.892482 0.892482 4.44089e-16 (R,S): 0.469634 0.135194 (U ,U* ,Error): 0.000544991 0.000544991 7.80626e-18 (Ur,Ur*,Error): 0.00232092 0.00232092 5.63785e-18 (Us,Us*,Error): 0.0120935 0.0120935 1.04083e-17 (R,S): 0.202348 0.86012 (U ,U* ,Error): 0.026054 0.026054 6.93889e-18 (Ur,Ur*,Error): 0.257517 0.257517 1.66533e-16 (Us,Us*,Error): 0.0908734 0.0908734 8.32667e-17 (R,S): 0.0406244 0.773973 (U ,U* ,Error): 0.000765156 0.000765156 4.44523e-18 (Ur,Ur*,Error): 0.0376698 0.0376698 1.38778e-17 (Us,Us*,Error): 0.00296583 0.00296583 8.67362e-19 Interpolate R^3 * S^2 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0 (R,S,U): 0.333333 0.333333 0.00411523 (R,S,U): 0.666667 0.333333 0.0329218 (R,S,U): 1 0.333333 0.111111 (R,S,U): 0 0.666667 0 (R,S,U): 0.333333 0.666667 0.0164609 (R,S,U): 0.666667 0.666667 0.131687 (R,S,U): 1 0.666667 0.444444 (R,S,U): 0 1 0 (R,S,U): 0.333333 1 0.037037 (R,S,U): 0.666667 1 0.296296 (R,S,U): 1 1 1 (R,S): 0.155941 0.895727 (U ,U* ,Error): 0.0030425 0.0030425 2.29851e-17 (Ur,Ur*,Error): 0.0585318 0.0585318 2.08167e-17 (Us,Us*,Error): 0.00679336 0.00679336 2.77556e-17 (R,S): 0.487127 0.148682 (U ,U* ,Error): 0.0025553 0.0025553 2.60209e-18 (Ur,Ur*,Error): 0.015737 0.015737 6.93889e-18 (Us,Us*,Error): 0.0343728 0.0343728 6.93889e-18 (R,S): 0.893345 0.445207 (U ,U* ,Error): 0.141313 0.141313 0 (Ur,Ur*,Error): 0.474552 0.474552 0 (Us,Us*,Error): 0.634818 0.634818 3.33067e-16 (R,S): 0.599347 0.228689 (U ,U* ,Error): 0.0112597 0.0112597 1.73472e-18 (Ur,Ur*,Error): 0.0563596 0.0563596 3.46945e-17 (Us,Us*,Error): 0.0984714 0.0984714 4.16334e-17 (R,S): 0.5723 0.639618 (U ,U* ,Error): 0.0766853 0.0766853 1.38778e-17 (Ur,Ur*,Error): 0.401985 0.401985 1.66533e-16 (Us,Us*,Error): 0.239785 0.239785 5.55112e-17 Interpolate R^3 * S^3 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0 (R,S,U): 0.333333 0.333333 0.00137174 (R,S,U): 0.666667 0.333333 0.0109739 (R,S,U): 1 0.333333 0.037037 (R,S,U): 0 0.666667 0 (R,S,U): 0.333333 0.666667 0.0109739 (R,S,U): 0.666667 0.666667 0.0877915 (R,S,U): 1 0.666667 0.296296 (R,S,U): 0 1 0 (R,S,U): 0.333333 1 0.037037 (R,S,U): 0.666667 1 0.296296 (R,S,U): 1 1 1 (R,S): 0.0659845 0.00211991 (U ,U* ,Error): 2.73704e-12 2.73704e-12 5.23167e-20 (Ur,Ur*,Error): 1.2444e-10 1.2444e-10 2.58908e-19 (Us,Us*,Error): 3.87333e-09 3.87333e-09 2.06749e-17 (R,S): 0.629405 0.412842 (U ,U* ,Error): 0.0175445 0.0175445 0 (Ur,Ur*,Error): 0.0836242 0.0836242 2.77556e-17 (Us,Us*,Error): 0.127491 0.127491 2.77556e-17 (R,S): 0.628397 0.469879 (U ,U* ,Error): 0.0257431 0.0257431 0 (Ur,Ur*,Error): 0.122899 0.122899 2.77556e-17 (Us,Us*,Error): 0.16436 0.16436 0 (R,S): 0.258509 0.754493 (U ,U* ,Error): 0.00741978 0.00741978 8.67362e-19 (Ur,Ur*,Error): 0.0861068 0.0861068 1.38778e-17 (Us,Us*,Error): 0.0295024 0.0295024 3.46945e-18 (R,S): 0.768513 0.403978 (U ,U* ,Error): 0.0299246 0.0299246 0 (Ur,Ur*,Error): 0.116815 0.116815 4.16334e-17 (Us,Us*,Error): 0.222224 0.222224 1.94289e-16 INTERP_TEST for element "QL". Number of nodes = 6 Interpolate R^0 * S^0 (R,S,U): 0 0 1 (R,S,U): 0.5 0 1 (R,S,U): 1 0 1 (R,S,U): 0 1 1 (R,S,U): 0.5 1 1 (R,S,U): 1 1 1 (R,S): 0.218418 0.956318 (U ,U* ,Error): 1 1 1.11022e-16 (Ur,Ur*,Error): 0 -6.93889e-17 6.93889e-17 (Us,Us*,Error): 0 0 0 (R,S): 0.829509 0.561695 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0 0 0 (R,S): 0.415307 0.0661187 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 1.17961e-16 1.17961e-16 (Us,Us*,Error): 0 -1.80411e-16 1.80411e-16 (R,S): 0.257578 0.109957 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 1.58727e-16 1.58727e-16 (Us,Us*,Error): 0 -1.94289e-16 1.94289e-16 (R,S): 0.043829 0.633966 (U ,U* ,Error): 1 1 2.22045e-16 (Ur,Ur*,Error): 0 3.33067e-16 3.33067e-16 (Us,Us*,Error): 0 4.16334e-17 4.16334e-17 Interpolate R^0 * S^1 (R,S,U): 0 0 0 (R,S,U): 0.5 0 0 (R,S,U): 1 0 0 (R,S,U): 0 1 1 (R,S,U): 0.5 1 1 (R,S,U): 1 1 1 (R,S): 0.0617272 0.449539 (U ,U* ,Error): 0.449539 0.449539 0 (Ur,Ur*,Error): 0 1.66533e-16 1.66533e-16 (Us,Us*,Error): 1 1 2.22045e-16 (R,S): 0.401306 0.754673 (U ,U* ,Error): 0.754673 0.754673 1.11022e-16 (Ur,Ur*,Error): 0 -5.55112e-17 5.55112e-17 (Us,Us*,Error): 1 1 0 (R,S): 0.797287 0.00183837 (U ,U* ,Error): 0.00183837 0.00183837 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 (R,S): 0.897504 0.350752 (U ,U* ,Error): 0.350752 0.350752 0 (Ur,Ur*,Error): 0 1.11022e-16 1.11022e-16 (Us,Us*,Error): 1 1 0 (R,S): 0.0945448 0.0136169 (U ,U* ,Error): 0.0136169 0.0136169 0 (Ur,Ur*,Error): 0 -5.20417e-18 5.20417e-18 (Us,Us*,Error): 1 1 0 Interpolate R^1 * S^0 (R,S,U): 0 0 0 (R,S,U): 0.5 0 0.5 (R,S,U): 1 0 1 (R,S,U): 0 1 0 (R,S,U): 0.5 1 0.5 (R,S,U): 1 1 1 (R,S): 0.859097 0.840847 (U ,U* ,Error): 0.859097 0.859097 0 (Ur,Ur*,Error): 1 1 2.22045e-16 (Us,Us*,Error): 0 0 0 (R,S): 0.123104 0.00751236 (U ,U* ,Error): 0.123104 0.123104 0 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 0 0 (R,S): 0.260303 0.912484 (U ,U* ,Error): 0.260303 0.260303 0 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 -1.38778e-17 1.38778e-17 (R,S): 0.113664 0.351629 (U ,U* ,Error): 0.113664 0.113664 2.77556e-17 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 0 0 (R,S): 0.822887 0.267132 (U ,U* ,Error): 0.822887 0.822887 1.11022e-16 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 0 0 Interpolate R^1 * S^1 (R,S,U): 0 0 0 (R,S,U): 0.5 0 0 (R,S,U): 1 0 0 (R,S,U): 0 1 0 (R,S,U): 0.5 1 0.5 (R,S,U): 1 1 1 (R,S): 0.692066 0.561662 (U ,U* ,Error): 0.388708 0.388708 0 (Ur,Ur*,Error): 0.561662 0.561662 0 (Us,Us*,Error): 0.692066 0.692066 0 (R,S): 0.861216 0.453794 (U ,U* ,Error): 0.390814 0.390814 0 (Ur,Ur*,Error): 0.453794 0.453794 5.55112e-17 (Us,Us*,Error): 0.861216 0.861216 0 (R,S): 0.911977 0.597917 (U ,U* ,Error): 0.545286 0.545286 0 (Ur,Ur*,Error): 0.597917 0.597917 1.11022e-16 (Us,Us*,Error): 0.911977 0.911977 0 (R,S): 0.188955 0.761492 (U ,U* ,Error): 0.143887 0.143887 0 (Ur,Ur*,Error): 0.761492 0.761492 1.11022e-16 (Us,Us*,Error): 0.188955 0.188955 2.77556e-17 (R,S): 0.396988 0.185314 (U ,U* ,Error): 0.0735676 0.0735676 0 (Ur,Ur*,Error): 0.185314 0.185314 0 (Us,Us*,Error): 0.396988 0.396988 0 Interpolate R^2 * S^0 (R,S,U): 0 0 0 (R,S,U): 0.5 0 0.25 (R,S,U): 1 0 1 (R,S,U): 0 1 0 (R,S,U): 0.5 1 0.25 (R,S,U): 1 1 1 (R,S): 0.574366 0.367027 (U ,U* ,Error): 0.329896 0.329896 0 (Ur,Ur*,Error): 1.14873 1.14873 0 (Us,Us*,Error): 0 1.38778e-17 1.38778e-17 (R,S): 0.617205 0.361529 (U ,U* ,Error): 0.380942 0.380942 0 (Ur,Ur*,Error): 1.23441 1.23441 0 (Us,Us*,Error): 0 0 0 (R,S): 0.21293 0.714471 (U ,U* ,Error): 0.0453392 0.0453392 6.93889e-18 (Ur,Ur*,Error): 0.42586 0.42586 5.55112e-17 (Us,Us*,Error): 0 0 0 (R,S): 0.117707 0.299329 (U ,U* ,Error): 0.0138549 0.0138549 1.73472e-18 (Ur,Ur*,Error): 0.235414 0.235414 5.55112e-17 (Us,Us*,Error): 0 0 0 (R,S): 0.825003 0.82466 (U ,U* ,Error): 0.68063 0.68063 0 (Ur,Ur*,Error): 1.65001 1.65001 0 (Us,Us*,Error): 0 0 0 Interpolate R^2 * S^1 (R,S,U): 0 0 0 (R,S,U): 0.5 0 0 (R,S,U): 1 0 0 (R,S,U): 0 1 0 (R,S,U): 0.5 1 0.25 (R,S,U): 1 1 1 (R,S): 0.0618618 0.710781 (U ,U* ,Error): 0.00272007 0.00272007 6.50521e-18 (Ur,Ur*,Error): 0.0879403 0.0879403 1.38778e-17 (Us,Us*,Error): 0.00382688 0.00382688 1.73472e-18 (R,S): 0.0882833 0.777994 (U ,U* ,Error): 0.00606364 0.00606364 1.30104e-17 (Ur,Ur*,Error): 0.137368 0.137368 2.77556e-17 (Us,Us*,Error): 0.00779395 0.00779395 1.21431e-17 (R,S): 0.745303 0.308675 (U ,U* ,Error): 0.171462 0.171462 0 (Ur,Ur*,Error): 0.460113 0.460113 0 (Us,Us*,Error): 0.555477 0.555477 0 (R,S): 0.899373 0.763537 (U ,U* ,Error): 0.617603 0.617603 0 (Ur,Ur*,Error): 1.37341 1.37341 0 (Us,Us*,Error): 0.808872 0.808872 0 (R,S): 0.761731 0.40697 (U ,U* ,Error): 0.236137 0.236137 2.77556e-17 (Ur,Ur*,Error): 0.620002 0.620002 0 (Us,Us*,Error): 0.580234 0.580234 0 INTERP_TEST for element "T3". Number of nodes = 3 Interpolate R^0 * S^0 (R,S,U): 0 0 1 (R,S,U): 1 0 1 (R,S,U): 0 1 1 (R,S): 0.781582 0.0436824 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0 0 0 (R,S): 0.170491 0.438305 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0 0 0 (R,S): 0.415307 0.0661187 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0 0 0 (R,S): 0.257578 0.109957 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0 0 0 (R,S): 0.043829 0.633966 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0 0 0 Interpolate R^0 * S^1 (R,S,U): 0 0 0 (R,S,U): 1 0 0 (R,S,U): 0 1 1 (R,S): 0.0617272 0.449539 (U ,U* ,Error): 0.449539 0.449539 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 (R,S): 0.598694 0.245327 (U ,U* ,Error): 0.245327 0.245327 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 (R,S): 0.797287 0.00183837 (U ,U* ,Error): 0.00183837 0.00183837 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 (R,S): 0.102496 0.649248 (U ,U* ,Error): 0.649248 0.649248 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 (R,S): 0.0945448 0.0136169 (U ,U* ,Error): 0.0136169 0.0136169 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 Interpolate R^1 * S^0 (R,S,U): 0 0 0 (R,S,U): 1 0 1 (R,S,U): 0 1 0 (R,S): 0.140903 0.159153 (U ,U* ,Error): 0.140903 0.140903 0 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 0 0 (R,S): 0.123104 0.00751236 (U ,U* ,Error): 0.123104 0.123104 0 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 0 0 (R,S): 0.739697 0.0875163 (U ,U* ,Error): 0.739697 0.739697 0 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 0 0 (R,S): 0.113664 0.351629 (U ,U* ,Error): 0.113664 0.113664 0 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 0 0 (R,S): 0.177113 0.732868 (U ,U* ,Error): 0.177113 0.177113 0 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 0 0 INTERP_TEST - Warning! Skipping test for element "T4". INTERP_TEST for element "T6". Number of nodes = 6 Interpolate R^0 * S^0 (R,S,U): 0 0 1 (R,S,U): 1 0 1 (R,S,U): 0 1 1 (R,S,U): 0.5 0 1 (R,S,U): 0.5 0.5 1 (R,S,U): 0 0.5 1 (R,S): 0.781582 0.0436824 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0 0 0 (R,S): 0.170491 0.438305 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0 0 0 (R,S): 0.415307 0.0661187 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 -5.55112e-17 5.55112e-17 (Us,Us*,Error): 0 2.22045e-16 2.22045e-16 (R,S): 0.257578 0.109957 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0 0 0 (R,S): 0.043829 0.633966 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0 0 0 Interpolate R^0 * S^1 (R,S,U): 0 0 0 (R,S,U): 1 0 0 (R,S,U): 0 1 1 (R,S,U): 0.5 0 0 (R,S,U): 0.5 0.5 0.5 (R,S,U): 0 0.5 0.5 (R,S): 0.0617272 0.449539 (U ,U* ,Error): 0.449539 0.449539 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 (R,S): 0.598694 0.245327 (U ,U* ,Error): 0.245327 0.245327 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 (R,S): 0.797287 0.00183837 (U ,U* ,Error): 0.00183837 0.00183837 2.1684e-19 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 (R,S): 0.102496 0.649248 (U ,U* ,Error): 0.649248 0.649248 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 (R,S): 0.0945448 0.0136169 (U ,U* ,Error): 0.0136169 0.0136169 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 0 Interpolate R^1 * S^0 (R,S,U): 0 0 0 (R,S,U): 1 0 1 (R,S,U): 0 1 0 (R,S,U): 0.5 0 0.5 (R,S,U): 0.5 0.5 0.5 (R,S,U): 0 0.5 0 (R,S): 0.140903 0.159153 (U ,U* ,Error): 0.140903 0.140903 0 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 0 0 (R,S): 0.123104 0.00751236 (U ,U* ,Error): 0.123104 0.123104 0 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 0 0 (R,S): 0.739697 0.0875163 (U ,U* ,Error): 0.739697 0.739697 0 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 0 0 (R,S): 0.113664 0.351629 (U ,U* ,Error): 0.113664 0.113664 0 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 0 0 (R,S): 0.177113 0.732868 (U ,U* ,Error): 0.177113 0.177113 2.77556e-17 (Ur,Ur*,Error): 1 1 0 (Us,Us*,Error): 0 0 0 Interpolate R^0 * S^2 (R,S,U): 0 0 0 (R,S,U): 1 0 0 (R,S,U): 0 1 1 (R,S,U): 0.5 0 0 (R,S,U): 0.5 0.5 0.25 (R,S,U): 0 0.5 0.25 (R,S): 0.307934 0.438338 (U ,U* ,Error): 0.19214 0.19214 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0.876675 0.876675 1.11022e-16 (R,S): 0.138784 0.546206 (U ,U* ,Error): 0.298341 0.298341 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1.09241 1.09241 0 (R,S): 0.088023 0.402083 (U ,U* ,Error): 0.161671 0.161671 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0.804166 0.804166 0 (R,S): 0.188955 0.761492 (U ,U* ,Error): 0.57987 0.57987 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1.52298 1.52298 2.22045e-16 (R,S): 0.396988 0.185314 (U ,U* ,Error): 0.0343413 0.0343413 6.93889e-18 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0.370628 0.370628 0 Interpolate R^1 * S^1 (R,S,U): 0 0 0 (R,S,U): 1 0 0 (R,S,U): 0 1 0 (R,S,U): 0.5 0 0 (R,S,U): 0.5 0.5 0.25 (R,S,U): 0 0.5 0 (R,S): 0.574366 0.367027 (U ,U* ,Error): 0.210808 0.210808 0 (Ur,Ur*,Error): 0.367027 0.367027 0 (Us,Us*,Error): 0.574366 0.574366 0 (R,S): 0.617205 0.361529 (U ,U* ,Error): 0.223137 0.223137 0 (Ur,Ur*,Error): 0.361529 0.361529 0 (Us,Us*,Error): 0.617205 0.617205 0 (R,S): 0.21293 0.714471 (U ,U* ,Error): 0.152132 0.152132 0 (Ur,Ur*,Error): 0.714471 0.714471 0 (Us,Us*,Error): 0.21293 0.21293 0 (R,S): 0.117707 0.299329 (U ,U* ,Error): 0.0352331 0.0352331 0 (Ur,Ur*,Error): 0.299329 0.299329 0 (Us,Us*,Error): 0.117707 0.117707 0 (R,S): 0.174997 0.17534 (U ,U* ,Error): 0.030684 0.030684 0 (Ur,Ur*,Error): 0.17534 0.17534 0 (Us,Us*,Error): 0.174997 0.174997 0 Interpolate R^2 * S^0 (R,S,U): 0 0 0 (R,S,U): 1 0 1 (R,S,U): 0 1 0 (R,S,U): 0.5 0 0.25 (R,S,U): 0.5 0.5 0.25 (R,S,U): 0 0.5 0 (R,S): 0.0618618 0.710781 (U ,U* ,Error): 0.00382688 0.00382688 1.73472e-18 (Ur,Ur*,Error): 0.123724 0.123724 0 (Us,Us*,Error): 0 0 0 (R,S): 0.0882833 0.777994 (U ,U* ,Error): 0.00779395 0.00779395 5.20417e-18 (Ur,Ur*,Error): 0.176567 0.176567 5.55112e-17 (Us,Us*,Error): 0 0 0 (R,S): 0.254697 0.691325 (U ,U* ,Error): 0.0648705 0.0648705 1.38778e-17 (Ur,Ur*,Error): 0.509394 0.509394 0 (Us,Us*,Error): 0 0 0 (R,S): 0.100627 0.236463 (U ,U* ,Error): 0.0101258 0.0101258 1.73472e-18 (Ur,Ur*,Error): 0.201254 0.201254 0 (Us,Us*,Error): 0 0 0 (R,S): 0.238269 0.59303 (U ,U* ,Error): 0.0567723 0.0567723 6.93889e-18 (Ur,Ur*,Error): 0.476539 0.476539 0 (Us,Us*,Error): 0 0 0 INTERP_TEST for element "T10". Number of nodes = 10 Interpolate R^0 * S^0 (R,S,U): 0 0 1 (R,S,U): 0.333333 0 1 (R,S,U): 0.666667 0 1 (R,S,U): 1 0 1 (R,S,U): 0 0.333333 1 (R,S,U): 0.333333 0.333333 1 (R,S,U): 0.666667 0.333333 1 (R,S,U): 0 0.666667 1 (R,S,U): 0.333333 0.666667 1 (R,S,U): 0 1 1 (R,S): 0.781582 0.0436824 (U ,U* ,Error): 1 1 2.22045e-16 (Ur,Ur*,Error): 0 8.88178e-16 8.88178e-16 (Us,Us*,Error): 0 4.44089e-16 4.44089e-16 (R,S): 0.170491 0.438305 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 7.77156e-16 7.77156e-16 (Us,Us*,Error): 0 2.60902e-15 2.60902e-15 (R,S): 0.415307 0.0661187 (U ,U* ,Error): 1 1 0 (Ur,Ur*,Error): 0 1.27676e-15 1.27676e-15 (Us,Us*,Error): 0 8.88178e-16 8.88178e-16 (R,S): 0.257578 0.109957 (U ,U* ,Error): 1 1 2.22045e-16 (Ur,Ur*,Error): 0 1.66533e-15 1.66533e-15 (Us,Us*,Error): 0 -3.05311e-16 3.05311e-16 (R,S): 0.043829 0.633966 (U ,U* ,Error): 1 1 2.22045e-16 (Ur,Ur*,Error): 0 -4.44089e-16 4.44089e-16 (Us,Us*,Error): 0 0 0 Interpolate R^0 * S^1 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0.333333 (R,S,U): 0.333333 0.333333 0.333333 (R,S,U): 0.666667 0.333333 0.333333 (R,S,U): 0 0.666667 0.666667 (R,S,U): 0.333333 0.666667 0.666667 (R,S,U): 0 1 1 (R,S): 0.0617272 0.449539 (U ,U* ,Error): 0.449539 0.449539 0 (Ur,Ur*,Error): 0 3.33067e-16 3.33067e-16 (Us,Us*,Error): 1 1 1.11022e-16 (R,S): 0.598694 0.245327 (U ,U* ,Error): 0.245327 0.245327 8.32667e-17 (Ur,Ur*,Error): 0 2.22045e-16 2.22045e-16 (Us,Us*,Error): 1 1 2.22045e-16 (R,S): 0.797287 0.00183837 (U ,U* ,Error): 0.00183837 0.00183837 1.0842e-18 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1 1 4.44089e-16 (R,S): 0.102496 0.649248 (U ,U* ,Error): 0.649248 0.649248 1.11022e-16 (Ur,Ur*,Error): 0 2.22045e-16 2.22045e-16 (Us,Us*,Error): 1 1 2.66454e-15 (R,S): 0.0945448 0.0136169 (U ,U* ,Error): 0.0136169 0.0136169 8.67362e-18 (Ur,Ur*,Error): 0 6.93889e-18 6.93889e-18 (Us,Us*,Error): 1 1 4.44089e-16 Interpolate R^1 * S^0 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0.333333 (R,S,U): 0.666667 0 0.666667 (R,S,U): 1 0 1 (R,S,U): 0 0.333333 0 (R,S,U): 0.333333 0.333333 0.333333 (R,S,U): 0.666667 0.333333 0.666667 (R,S,U): 0 0.666667 0 (R,S,U): 0.333333 0.666667 0.333333 (R,S,U): 0 1 0 (R,S): 0.140903 0.159153 (U ,U* ,Error): 0.140903 0.140903 0 (Ur,Ur*,Error): 1 1 1.11022e-16 (Us,Us*,Error): 0 2.25514e-17 2.25514e-17 (R,S): 0.123104 0.00751236 (U ,U* ,Error): 0.123104 0.123104 0 (Ur,Ur*,Error): 1 1 2.22045e-16 (Us,Us*,Error): 0 5.55112e-17 5.55112e-17 (R,S): 0.739697 0.0875163 (U ,U* ,Error): 0.739697 0.739697 0 (Ur,Ur*,Error): 1 1 8.88178e-16 (Us,Us*,Error): 0 1.11022e-16 1.11022e-16 (R,S): 0.113664 0.351629 (U ,U* ,Error): 0.113664 0.113664 1.38778e-17 (Ur,Ur*,Error): 1 1 2.22045e-16 (Us,Us*,Error): 0 0 0 (R,S): 0.177113 0.732868 (U ,U* ,Error): 0.177113 0.177113 2.77556e-17 (Ur,Ur*,Error): 1 1 1.11022e-16 (Us,Us*,Error): 0 2.22045e-16 2.22045e-16 Interpolate R^0 * S^2 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0.111111 (R,S,U): 0.333333 0.333333 0.111111 (R,S,U): 0.666667 0.333333 0.111111 (R,S,U): 0 0.666667 0.444444 (R,S,U): 0.333333 0.666667 0.444444 (R,S,U): 0 1 1 (R,S): 0.307934 0.438338 (U ,U* ,Error): 0.19214 0.19214 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0.876675 0.876675 3.33067e-16 (R,S): 0.138784 0.546206 (U ,U* ,Error): 0.298341 0.298341 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1.09241 1.09241 4.44089e-16 (R,S): 0.088023 0.402083 (U ,U* ,Error): 0.161671 0.161671 2.77556e-17 (Ur,Ur*,Error): 0 1.38778e-16 1.38778e-16 (Us,Us*,Error): 0.804166 0.804166 2.22045e-16 (R,S): 0.188955 0.761492 (U ,U* ,Error): 0.57987 0.57987 1.11022e-16 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 1.52298 1.52298 2.22045e-16 (R,S): 0.396988 0.185314 (U ,U* ,Error): 0.0343413 0.0343413 6.93889e-18 (Ur,Ur*,Error): 0 5.55112e-17 5.55112e-17 (Us,Us*,Error): 0.370628 0.370628 2.22045e-16 Interpolate R^1 * S^1 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0 (R,S,U): 0.333333 0.333333 0.111111 (R,S,U): 0.666667 0.333333 0.222222 (R,S,U): 0 0.666667 0 (R,S,U): 0.333333 0.666667 0.222222 (R,S,U): 0 1 0 (R,S): 0.574366 0.367027 (U ,U* ,Error): 0.210808 0.210808 0 (Ur,Ur*,Error): 0.367027 0.367027 1.66533e-16 (Us,Us*,Error): 0.574366 0.574366 0 (R,S): 0.617205 0.361529 (U ,U* ,Error): 0.223137 0.223137 2.77556e-17 (Ur,Ur*,Error): 0.361529 0.361529 5.55112e-17 (Us,Us*,Error): 0.617205 0.617205 0 (R,S): 0.21293 0.714471 (U ,U* ,Error): 0.152132 0.152132 0 (Ur,Ur*,Error): 0.714471 0.714471 1.11022e-16 (Us,Us*,Error): 0.21293 0.21293 8.32667e-17 (R,S): 0.117707 0.299329 (U ,U* ,Error): 0.0352331 0.0352331 6.93889e-18 (Ur,Ur*,Error): 0.299329 0.299329 5.55112e-17 (Us,Us*,Error): 0.117707 0.117707 0 (R,S): 0.174997 0.17534 (U ,U* ,Error): 0.030684 0.030684 1.04083e-17 (Ur,Ur*,Error): 0.17534 0.17534 0 (Us,Us*,Error): 0.174997 0.174997 2.77556e-17 Interpolate R^2 * S^0 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0.111111 (R,S,U): 0.666667 0 0.444444 (R,S,U): 1 0 1 (R,S,U): 0 0.333333 0 (R,S,U): 0.333333 0.333333 0.111111 (R,S,U): 0.666667 0.333333 0.444444 (R,S,U): 0 0.666667 0 (R,S,U): 0.333333 0.666667 0.111111 (R,S,U): 0 1 0 (R,S): 0.0618618 0.710781 (U ,U* ,Error): 0.00382688 0.00382688 1.9082e-17 (Ur,Ur*,Error): 0.123724 0.123724 5.55112e-17 (Us,Us*,Error): 0 0 0 (R,S): 0.0882833 0.777994 (U ,U* ,Error): 0.00779395 0.00779395 1.9082e-17 (Ur,Ur*,Error): 0.176567 0.176567 5.55112e-17 (Us,Us*,Error): 0 2.77556e-17 2.77556e-17 (R,S): 0.254697 0.691325 (U ,U* ,Error): 0.0648705 0.0648705 2.77556e-17 (Ur,Ur*,Error): 0.509394 0.509394 1.11022e-16 (Us,Us*,Error): 0 0 0 (R,S): 0.100627 0.236463 (U ,U* ,Error): 0.0101258 0.0101258 3.46945e-18 (Ur,Ur*,Error): 0.201254 0.201254 5.55112e-17 (Us,Us*,Error): 0 3.1225e-17 3.1225e-17 (R,S): 0.238269 0.59303 (U ,U* ,Error): 0.0567723 0.0567723 6.93889e-18 (Ur,Ur*,Error): 0.476539 0.476539 2.77556e-16 (Us,Us*,Error): 0 0 0 Interpolate R^0 * S^3 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0.037037 (R,S,U): 0.333333 0.333333 0.037037 (R,S,U): 0.666667 0.333333 0.037037 (R,S,U): 0 0.666667 0.296296 (R,S,U): 0.333333 0.666667 0.296296 (R,S,U): 0 1 1 (R,S): 0.0612505 0.437912 (U ,U* ,Error): 0.0839769 0.0839769 0 (Ur,Ur*,Error): 0 2.77556e-17 2.77556e-17 (Us,Us*,Error): 0.5753 0.5753 2.22045e-16 (R,S): 0.01782 0.501103 (U ,U* ,Error): 0.125829 0.125829 0 (Ur,Ur*,Error): 0 0 0 (Us,Us*,Error): 0.753313 0.753313 1.11022e-16 (R,S): 0.0419093 0.368851 (U ,U* ,Error): 0.0501824 0.0501824 6.93889e-18 (Ur,Ur*,Error): 0 -2.77556e-17 2.77556e-17 (Us,Us*,Error): 0.408152 0.408152 2.22045e-16 (R,S): 0.728276 0.141427 (U ,U* ,Error): 0.00282879 0.00282879 1.0842e-17 (Ur,Ur*,Error): 0 -1.38778e-17 1.38778e-17 (Us,Us*,Error): 0.0600052 0.0600052 1.38778e-16 (R,S): 0.0290366 0.0174423 (U ,U* ,Error): 5.30652e-06 5.30652e-06 1.21803e-18 (Ur,Ur*,Error): 0 3.46945e-18 3.46945e-18 (Us,Us*,Error): 0.000912699 0.000912699 1.36393e-16 Interpolate R^1 * S^2 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0 (R,S,U): 0.333333 0.333333 0.037037 (R,S,U): 0.666667 0.333333 0.0740741 (R,S,U): 0 0.666667 0 (R,S,U): 0.333333 0.666667 0.148148 (R,S,U): 0 1 0 (R,S): 0.152384 0.114319 (U ,U* ,Error): 0.00199147 0.00199147 0 (Ur,Ur*,Error): 0.0130688 0.0130688 8.67362e-18 (Us,Us*,Error): 0.0348406 0.0348406 2.77556e-17 (R,S): 0.353907 0.119308 (U ,U* ,Error): 0.00503764 0.00503764 8.67362e-19 (Ur,Ur*,Error): 0.0142344 0.0142344 5.20417e-18 (Us,Us*,Error): 0.0844478 0.0844478 0 (R,S): 0.206653 0.212924 (U ,U* ,Error): 0.00936894 0.00936894 5.20417e-18 (Ur,Ur*,Error): 0.0453366 0.0453366 2.08167e-17 (Us,Us*,Error): 0.0880026 0.0880026 1.38778e-17 (R,S): 0.387052 0.190481 (U ,U* ,Error): 0.0140434 0.0140434 3.46945e-18 (Ur,Ur*,Error): 0.036283 0.036283 1.38778e-17 (Us,Us*,Error): 0.147452 0.147452 2.77556e-17 (R,S): 0.58709 0.215492 (U ,U* ,Error): 0.0272625 0.0272625 1.73472e-17 (Ur,Ur*,Error): 0.0464366 0.0464366 2.08167e-17 (Us,Us*,Error): 0.253026 0.253026 0 Interpolate R^2 * S^1 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0 (R,S,U): 0.666667 0 0 (R,S,U): 1 0 0 (R,S,U): 0 0.333333 0 (R,S,U): 0.333333 0.333333 0.037037 (R,S,U): 0.666667 0.333333 0.148148 (R,S,U): 0 0.666667 0 (R,S,U): 0.333333 0.666667 0.0740741 (R,S,U): 0 1 0 (R,S): 0.231944 0.276703 (U ,U* ,Error): 0.014886 0.014886 0 (Ur,Ur*,Error): 0.128359 0.128359 2.77556e-17 (Us,Us*,Error): 0.0537979 0.0537979 1.38778e-17 (R,S): 0.551981 0.144824 (U ,U* ,Error): 0.0441254 0.0441254 6.93889e-18 (Ur,Ur*,Error): 0.15988 0.15988 2.77556e-17 (Us,Us*,Error): 0.304683 0.304683 0 (R,S): 0.0549825 0.0909431 (U ,U* ,Error): 0.000274928 0.000274928 1.0842e-19 (Ur,Ur*,Error): 0.0100006 0.0100006 0 (Us,Us*,Error): 0.00302308 0.00302308 3.90313e-18 (R,S): 0.519726 0.0301946 (U ,U* ,Error): 0.00815602 0.00815602 1.73472e-18 (Ur,Ur*,Error): 0.0313858 0.0313858 0 (Us,Us*,Error): 0.270115 0.270115 5.55112e-17 (R,S): 0.481067 0.292313 (U ,U* ,Error): 0.0676486 0.0676486 0 (Ur,Ur*,Error): 0.281244 0.281244 1.11022e-16 (Us,Us*,Error): 0.231425 0.231425 2.77556e-17 Interpolate R^3 * S^0 (R,S,U): 0 0 0 (R,S,U): 0.333333 0 0.037037 (R,S,U): 0.666667 0 0.296296 (R,S,U): 1 0 1 (R,S,U): 0 0.333333 0 (R,S,U): 0.333333 0.333333 0.037037 (R,S,U): 0.666667 0.333333 0.296296 (R,S,U): 0 0.666667 0 (R,S,U): 0.333333 0.666667 0.037037 (R,S,U): 0 1 0 (R,S): 0.0973602 0.332158 (U ,U* ,Error): 0.000922877 0.000922877 1.19262e-18 (Ur,Ur*,Error): 0.028437 0.028437 1.17961e-16 (Us,Us*,Error): 0 1.04083e-17 1.04083e-17 (R,S): 0.412278 0.156948 (U ,U* ,Error): 0.0700762 0.0700762 0 (Ur,Ur*,Error): 0.50992 0.50992 2.22045e-16 (Us,Us*,Error): 0 2.60209e-17 2.60209e-17 (R,S): 0.166718 0.0355958 (U ,U* ,Error): 0.00463394 0.00463394 3.46945e-18 (Ur,Ur*,Error): 0.083385 0.083385 0 (Us,Us*,Error): 0 -6.93889e-18 6.93889e-18 (R,S): 0.25921 0.543901 (U ,U* ,Error): 0.0174163 0.0174163 6.93889e-18 (Ur,Ur*,Error): 0.20157 0.20157 3.33067e-16 (Us,Us*,Error): 0 0 0 (R,S): 0.346439 0.593173 (U ,U* ,Error): 0.0415795 0.0415795 6.93889e-18 (Ur,Ur*,Error): 0.360059 0.360059 1.11022e-16 (Us,Us*,Error): 0 2.77556e-17 2.77556e-17 TEST13 MAP_TEST tests the map routines. MAP_TEST: The interpolation matrix for element Q4 The interpolation matrix: Col: 1 2 3 4 Row --- 1 1 0 -0 0 2 -1 0 0 1 3 -1 1 0 0 4 1 -1 1 -1 MAP_TEST: The interpolation matrix for element Q8 The interpolation matrix: Col: 1 2 3 4 5 Row --- 1 1 0 0 0 0 2 -3 0 0 -1 0 3 -3 -1 0 0 4 4 2 0 0 2 0 5 5 -1 -3 -1 -4 6 2 2 0 0 -4 7 -2 2 2 -2 0 8 -2 -2 2 2 4 Col: 6 7 8 Row --- 1 0 0 -0 2 0 0 4 3 0 0 -0 4 0 0 -4 5 4 4 -4 6 0 0 -0 7 -4 0 4 8 0 -4 0 MAP_TEST: The interpolation matrix for element Q9 The interpolation matrix: Col: 1 2 3 4 5 Row --- 1 1 0 0 0 0 2 -3 0 0 -1 0 3 -3 -1 0 0 4 4 2 0 0 2 0 5 9 3 1 3 -12 6 2 2 0 0 -4 7 -6 -2 -2 -6 8 8 -6 -6 -2 -2 12 9 4 4 4 4 -8 Col: 6 7 8 9 Row --- 1 0 0 0 -0 2 0 0 4 -0 3 0 0 0 -0 4 0 0 -4 -0 5 -4 -4 -12 16 6 0 0 0 -0 7 4 8 12 -16 8 8 4 8 -16 9 -8 -8 -8 16 MAP_TEST: The interpolation matrix for element Q12 The interpolation matrix: Col: 1 2 3 4 5 Row --- 1 1 0 -0 0 0 2 -5.5 0 -0 0 9 3 -5.5 9 -4.5 1 0 4 9 0 -0 0 -22.5 5 10 -9 4.5 -5.5 -9 6 9 -22.5 18 -4.5 0 7 -4.5 0 -0 0 13.5 8 -9 0 -0 9 22.5 9 -9 22.5 -18 4.5 0 10 -4.5 13.5 -13.5 4.5 0 11 4.5 0 -0 -4.5 -13.5 12 4.5 -13.5 13.5 -4.5 0 Col: 6 7 8 9 10 Row --- 1 0 0 0 0 0 2 0 -4.5 0 1 0 3 0 0 0 0 0 4 0 18 0 -4.5 0 5 9 4.5 -4.5 -5.5 9 6 -0 0 0 0 0 7 0 -13.5 0 4.5 0 8 -22.5 -18 18 4.5 0 9 0 0 0 9 -22.5 10 0 0 0 0 0 11 13.5 13.5 -13.5 -4.5 0 12 0 0 0 -4.5 13.5 Col: 11 12 Row --- 1 0 0 2 0 0 3 0 0 4 0 0 5 -4.5 1 6 0 0 7 0 0 8 0 -4.5 9 18 -4.5 10 0 0 11 0 4.5 12 -13.5 4.5 MAP_TEST: The interpolation matrix for element Q16 The interpolation matrix: Col: 1 2 3 4 5 Row --- 1 1 0 0 0 0 2 -5.5 0 0 0 9 3 -5.5 9 -4.5 1 1.76942e-15 4 9 0 0 0 -22.5 5 30.25 -49.5 24.75 -5.5 -49.5 6 9 -22.5 18 -4.5 4.19803e-15 7 -4.5 0 0 0 13.5 8 -49.5 81 -40.5 9 123.75 9 -49.5 123.75 -99 24.75 81 10 -4.5 13.5 -13.5 4.5 2.24126e-14 11 24.75 -40.5 20.25 -4.5 -74.25 12 81 -202.5 162 -40.5 -202.5 13 24.75 -74.25 74.25 -24.75 -40.5 14 -40.5 101.25 -81 20.25 121.5 15 -40.5 121.5 -121.5 40.5 101.25 16 20.25 -60.75 60.75 -20.25 -60.75 Col: 6 7 8 9 10 Row --- 1 0 0 0 0 0 2 0 0 0 -4.5 0 3 8.88178e-16 -8.43076e-15 5.77316e-15 -2.04003e-15 3.3723e-15 4 0 0 0 18 0 5 81 -40.5 9 24.75 -40.5 6 -1.77636e-14 2.86646e-14 -1.5099e-14 3.64292e-15 1.68615e-15 7 0 0 0 -13.5 0 8 -202.5 101.25 -22.5 -99 162 9 -202.5 162 -40.5 -40.5 101.25 10 -5.15143e-14 3.70953e-14 -7.99361e-15 -2.19824e-14 5.39568e-14 11 121.5 -60.75 13.5 74.25 -121.5 12 506.25 -405 101.25 162 -405 13 121.5 -121.5 40.5 20.25 -60.75 14 -303.75 243 -60.75 -121.5 303.75 15 -303.75 303.75 -101.25 -81 243 16 182.25 -182.25 60.75 60.75 -182.25 Col: 11 12 13 14 15 Row --- 1 0 0 0 0 0 2 0 0 1 0 0 3 1.77636e-15 -3.10862e-15 6.12936e-16 -1.77636e-15 5.6205e-16 4 0 0 -4.5 0 0 5 20.25 -4.5 -5.5 9 -4.5 6 -1.42109e-14 8.88178e-15 -3.21271e-15 3.55271e-15 5.6205e-16 7 0 0 4.5 0 0 8 -81 18 24.75 -40.5 20.25 9 -81 20.25 9 -22.5 18 10 -4.08562e-14 8.88178e-15 6.57961e-15 -1.77636e-14 1.46133e-14 11 60.75 -13.5 -24.75 40.5 -20.25 12 324 -81 -40.5 101.25 -81 13 60.75 -20.25 -4.5 13.5 -13.5 14 -243 60.75 40.5 -101.25 81 15 -243 81 20.25 -60.75 60.75 16 182.25 -60.75 -20.25 60.75 -60.75 Col: 16 Row --- 1 -0 2 0 3 6.01371e-16 4 -0 5 1 6 -9.02056e-16 7 0 8 -4.5 9 -4.5 10 -3.42936e-15 11 4.5 12 20.25 13 4.5 14 -20.25 15 -20.25 16 20.25 MAP_TEST: The interpolation matrix for element QL The interpolation matrix: Col: 1 2 3 4 5 Row --- 1 1 -0 0 0 0 2 -1 -0 0 1 0 3 -3 4 -1 0 0 4 3 -4 1 -3 4 5 2 -4 2 0 0 6 -2 4 -2 2 -4 Col: 6 Row --- 1 0 2 0 3 -0 4 -1 5 0 6 2 MAP_TEST: The interpolation matrix for element T3 The interpolation matrix: Col: 1 2 3 Row --- 1 1 -0 0 2 -1 -0 1 3 -1 1 0 MAP_TEST: The interpolation matrix for element T6 The interpolation matrix: Col: 1 2 3 4 5 Row --- 1 1 -0 0 0 0 2 -3 -0 -1 0 0 3 -3 -1 0 4 0 4 2 0 2 0 0 5 4 0 0 -4 4 6 2 2 0 -4 0 Col: 6 Row --- 1 0 2 4 3 0 4 -4 5 -4 6 0 MAP_TEST: The interpolation matrix for element T10 The interpolation matrix: Col: 1 2 3 4 5 Row --- 1 1 -0 0 0 0 2 -5.5 -0 0 0 9 3 -5.5 9 -4.5 1 0 4 9 -0 0 0 -22.5 5 18 -22.5 4.5 -8.88178e-16 -22.5 6 9 -22.5 18 -4.5 0 7 -4.5 -0 0 0 13.5 8 -13.5 13.5 0 0 27 9 -13.5 27 -13.5 0 13.5 10 -4.5 13.5 -13.5 4.5 0 Col: 6 7 8 9 10 Row --- 1 0 0 0 0 0 2 0 0 -4.5 0 1 3 0 0 0 -0 0 4 0 0 18 0 -4.5 5 27 -4.5 4.5 -4.5 0 6 0 0 0 0 0 7 0 0 -13.5 0 4.5 8 -27 0 -13.5 13.5 0 9 -27 13.5 0 0 0 10 0 0 0 0 0 TEST135 MASS_MATRIX_T3 computes the mass matrix for a finite element system using T3 elements (linear triangles). The T3 mass matrix: Col: 1 2 3 4 5 Row --- 1 0.0208333 0.0104167 0 0.0104167 0 2 0.0104167 0.0625 0.0104167 0.0208333 0.0208333 3 0 0.0104167 0.0416667 0 0.0208333 4 0.0104167 0.0208333 0 0.0625 0.0208333 5 0 0.0208333 0.0208333 0.0208333 0.125 6 0 0 0.0104167 0 0.0208333 7 0 0 0 0.0104167 0.0208333 8 0 0 0 0 0.0208333 9 0 0 0 0 0 Col: 6 7 8 9 Row --- 1 0 0 0 0 2 0 0 0 0 3 0.0104167 0 0 0 4 0 0.0104167 0 0 5 0.0208333 0.0208333 0.0208333 0 6 0.0625 0 0.0208333 0.0104167 7 0 0.0416667 0.0104167 0 8 0.0208333 0.0104167 0.0625 0.0104167 9 0.0104167 0 0.0104167 0.0208333 TEST14 MASS_MATRIX_T6 computes the mass matrix for a finite element system using T6 elements (quadratic triangles). The T6 mass matrix: Col: 1 2 3 4 5 Row --- 1 0.0166667 8.68446e-17 -0.00277778 8.84709e-17 -0.0111111 2 8.7712e-17 0.0888889 6.33174e-17 0.0444444 0.0444444 3 -0.00277778 6.28837e-17 0.0333333 -0.0111111 1.54824e-16 4 8.84709e-17 0.0444444 -0.0111111 0.0888889 0.0444444 5 -0.0111111 0.0444444 1.57426e-16 0.0444444 0.177778 6 0 0 6.59195e-17 0 0.0444444 7 -0.00277778 -0.0111111 -0.00555556 6.54858e-17 1.54445e-16 8 0 0 -0.0111111 0 0.0444444 9 0 0 -0.00277778 0 -0.0111111 Col: 6 7 8 9 Row --- 1 0 -0.00277778 0 0 2 0 -0.0111111 0 0 3 6.54858e-17 -0.00555556 -0.0111111 -0.00277778 4 0 6.59195e-17 0 0 5 0.0444444 1.55746e-16 0.0444444 -0.0111111 6 0.0888889 -0.0111111 0.0444444 8.84709e-17 7 -0.0111111 0.0333333 6.28837e-17 -0.00277778 8 0.0444444 6.33174e-17 0.0888889 8.7712e-17 9 8.84709e-17 -0.00277778 8.68446e-17 0.0166667 TEST15 For an order 3 triangle, PHYSICAL_TO_REFERENCE_T3 maps a physical point to a reference point. REFERENCE_TO_PHYSICAL_T3 maps a reference point to a physical point. XSI ETA ==> X Y ==> XSI2 ETA2 0.781582 0.0436824 2.60685 1.17473 0.781582 0.0436824 0.170491 0.438305 1.77929 2.75322 0.170491 0.438305 0.415307 0.0661187 1.89673 1.26447 0.415307 0.0661187 0.257578 0.109957 1.62511 1.43983 0.257578 0.109957 0.043829 0.633966 1.72162 3.53586 0.043829 0.633966 0.0617272 0.449539 1.57299 2.79816 0.0617272 0.449539 0.598694 0.245327 2.44271 1.98131 0.598694 0.245327 0.797287 0.00183837 2.59641 1.00735 0.797287 0.00183837 0.102496 0.649248 1.85424 3.59699 0.102496 0.649248 0.0945448 0.0136169 1.20271 1.05447 0.0945448 0.0136169 TEST16 For an order 6 triangle, REFERENCE_TO_PHYSICAL_T6 maps a reference point to a physical point. XSI ETA ==> X Y 0 0 0 0 1 0 2 0 0 1 0 4 0.5 0 1 0 0.5 0.5 1 1 0 0.5 0 2 0.25 0.75 0.5 2.25 0.75 0.25 1.5 0.25 0.4 0.1 0.8 0.24 0.3 0.2 0.6 0.56 0.2 0.3 0.4 0.96 0.1 0.4 0.2 1.44 0.1 0.1 0.2 0.36 0.2 0.2 0.4 0.64 0.3 0.3 0.6 0.84 0.4 0.4 0.8 0.96 TEST18 SHAPE_TEST tests the shape routines. SHAPE_TEST: Verify shape functions of type Q4 Element order = 4 Basis function values at basis nodes should form the identity matrix. 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 The R and S derivatives should sum to 0. dTdR sum, dTdS sum: 0 0 0 0 0 0 0 0 SHAPE_TEST: Verify shape functions of type Q8 Element order = 8 Basis function values at basis nodes should form the identity matrix. 1 -0 -0 -0 0 0 0 0 0 1 -0 0 -0 0 -0 -0 -0 0 1 0 0 -0 -0 0 0 0 -0 1 -0 -0 0 -0 0 -0 -0 -0 1 0 0 0 0 -0 0 0 -0 1 -0 -0 0 0 0 -0 -0 -0 1 -0 0 -0 -0 -0 0 0 0 1 The R and S derivatives should sum to 0. dTdR sum, dTdS sum: 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 SHAPE_TEST: Verify shape functions of type Q9 Element order = 9 Basis function values at basis nodes should form the identity matrix. 1 -0 0 -0 0 -0 -0 0 0 0 1 -0 -0 -0 0 0 0 -0 0 0 1 0 -0 -0 -0 -0 0 0 -0 -0 1 0 0 0 -0 -0 -0 0 -0 0 1 0 -0 -0 0 -0 -0 0 0 0 1 -0 0 -0 -0 0 0 -0 0 -0 1 0 -0 -0 0 -0 0 -0 -0 0 1 0 0 -0 0 -0 -0 0 0 -0 1 The R and S derivatives should sum to 0. dTdR sum, dTdS sum: 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 SHAPE_TEST: Verify shape functions of type Q12 Element order = 12 Basis function values at basis nodes should form the identity matrix. 1 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 1 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 1 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 1 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 1 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 1 0 0 -0 0 -0 -0 0 0 -0 -0 0 -0 1 The R and S derivatives should sum to 0. dTdR sum, dTdS sum: 0 0 1.33227e-15 4.44089e-16 2.77556e-15 0 1.77636e-15 0 0 1.33227e-15 2.44249e-15 1.33227e-15 0 2.77556e-15 2.22045e-15 2.77556e-15 0 1.77636e-15 1.33227e-15 2.22045e-15 2.77556e-15 1.77636e-15 1.77636e-15 1.77636e-15 SHAPE_TEST: Verify shape functions of type Q16 Element order = 16 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 1 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 1 -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 1 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 1 -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 1 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 1 -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 1 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 1 The R and S derivatives should sum to 0. dTdR sum, dTdS sum: 0 0 8.88178e-16 -1.55431e-15 2.66454e-15 -1.55431e-15 -2.66454e-15 -6.66134e-16 0 8.88178e-16 4.44089e-16 7.77156e-16 2.44249e-15 7.77156e-16 -4.44089e-15 1.22125e-15 0 2.22045e-15 4.44089e-16 2.22045e-15 2.44249e-15 2.22045e-15 -4.44089e-15 2.22045e-15 -6.66134e-16 -5.32907e-15 7.77156e-16 -2.66454e-15 2.44249e-15 -2.66454e-15 -3.55271e-15 -3.55271e-15 SHAPE_TEST: Verify shape functions of type QL Element order = 6 Basis function values at basis nodes should form the identity matrix. 1 0 -0 0 0 -0 -0 1 0 -0 0 0 0 -0 1 0 -0 0 -0 -0 0 1 0 -0 0 -0 -0 -0 1 0 -0 0 -0 0 -0 1 The R and S derivatives should sum to 0. dTdR sum, dTdS sum: 0 0 0 0 0 0 0 0 0 0 0 0 SHAPE_TEST: Verify shape functions of type T3 Element order = 3 Basis function values at basis nodes should form the identity matrix. 1 0 0 0 1 0 0 0 1 The R and S derivatives should sum to 0. dTdR sum, dTdS sum: 0 0 0 0 0 0 SHAPE_TEST: Verify shape functions of type T6 Element order = 6 Basis function values at basis nodes should form the identity matrix. 1 -0 -0 0 0 0 -0 1 -0 0 0 0 -0 -0 1 0 0 0 0 0 -0 1 0 0 -0 0 0 0 1 0 0 -0 0 0 0 1 The R and S derivatives should sum to 0. dTdR sum, dTdS sum: 0 0 0 0 0 0 0 0 0 0 0 0 SHAPE_TEST: Verify shape functions of type T10 Element order = 10 Basis function values at basis nodes should form the identity matrix. 1 0 -0 0 0 0 -0 -0 -0 0 0 1 -0 0 0 0 0 -0 -0 0 -0 0 1 -0 0 0 0 -0 -0 0 0 -0 0 1 -0 0 0 -0 -0 0 0 0 -0 0 1 0 -0 0 0 0 -0 0 -0 0 0 1 0 0 0 0 5.55112e-17 -1.66533e-16 1.66533e-16 -0 -8.32667e-17 3.33067e-16 1 0 0 0 -0 0 -0 0 0 0 -0 1 0 -0 1.11022e-16 -1.66533e-16 -0 0 -3.33067e-16 6.66134e-16 0 3.33067e-16 1 -0 0 -0 -0 0 -0 0 -0 0 0 1 The R and S derivatives should sum to 0. dTdR sum, dTdS sum: 8.88178e-16 8.88178e-16 1.33227e-15 8.88178e-16 9.99201e-16 8.88178e-16 -1.77636e-15 8.88178e-16 8.88178e-16 1.33227e-15 8.88178e-16 8.88178e-16 0 4.44089e-15 0 9.99201e-16 1.77636e-15 1.22125e-15 1.77636e-15 -1.77636e-15 TEST19 SPHERE_GRID_Q4_ELEMENT sets up a grid of Q4 quadrilaterals on a sphere. SPHERE_GRID_Q4_ELEMENT_NUM returns the number of elements in the grid SPHERE_GRID_Q4_NODE_NUM returns the number of nodes in the grid. SPHERE_GRID_Q4_NODE_XYZ returns the coordinates of nodes in the grid. Expected number of nodes = 58 Expected number of elements = 64 The elements and their nodes: 1 1 1 3 2 2 1 1 4 3 3 1 1 5 4 4 1 1 6 5 5 1 1 7 6 6 1 1 8 7 7 1 1 9 8 8 1 1 2 9 9 2 3 11 10 10 3 4 12 11 11 4 5 13 12 12 5 6 14 13 13 6 7 15 14 14 7 8 16 15 15 8 9 17 16 16 9 2 10 17 17 10 11 19 18 18 11 12 20 19 19 12 13 21 20 20 13 14 22 21 21 14 15 23 22 22 15 16 24 23 23 16 17 25 24 24 17 10 18 25 25 18 19 27 26 26 19 20 28 27 27 20 21 29 28 28 21 22 30 29 29 22 23 31 30 30 23 24 32 31 31 24 25 33 32 32 25 18 26 33 33 26 27 35 34 34 27 28 36 35 35 28 29 37 36 36 29 30 38 37 37 30 31 39 38 38 31 32 40 39 39 32 33 41 40 40 33 26 34 41 41 34 35 43 42 42 35 36 44 43 43 36 37 45 44 44 37 38 46 45 45 38 39 47 46 46 39 40 48 47 47 40 41 49 48 48 41 34 42 49 49 42 43 51 50 50 43 44 52 51 51 44 45 53 52 52 45 46 54 53 53 46 47 55 54 54 47 48 56 55 55 48 49 57 56 56 49 42 50 57 57 50 51 58 58 58 51 52 58 58 59 52 53 58 58 60 53 54 58 58 61 54 55 58 58 62 55 56 58 58 63 56 57 58 58 64 57 50 58 58 The node coordinates: 1 0 0 -1 2 0.382683 0 -0.92388 3 0.270598 0.270598 -0.92388 4 2.34326e-17 0.382683 -0.92388 5 -0.270598 0.270598 -0.92388 6 -0.382683 4.68652e-17 -0.92388 7 -0.270598 -0.270598 -0.92388 8 -7.02978e-17 -0.382683 -0.92388 9 0.270598 -0.270598 -0.92388 10 0.707107 0 -0.707107 11 0.5 0.5 -0.707107 12 4.32978e-17 0.707107 -0.707107 13 -0.5 0.5 -0.707107 14 -0.707107 8.65956e-17 -0.707107 15 -0.5 -0.5 -0.707107 16 -1.29893e-16 -0.707107 -0.707107 17 0.5 -0.5 -0.707107 18 0.92388 0 -0.382683 19 0.653281 0.653281 -0.382683 20 5.65713e-17 0.92388 -0.382683 21 -0.653281 0.653281 -0.382683 22 -0.92388 1.13143e-16 -0.382683 23 -0.653281 -0.653281 -0.382683 24 -1.69714e-16 -0.92388 -0.382683 25 0.653281 -0.653281 -0.382683 26 1 0 6.12323e-17 27 0.707107 0.707107 6.12323e-17 28 6.12323e-17 1 6.12323e-17 29 -0.707107 0.707107 6.12323e-17 30 -1 1.22465e-16 6.12323e-17 31 -0.707107 -0.707107 6.12323e-17 32 -1.83697e-16 -1 6.12323e-17 33 0.707107 -0.707107 6.12323e-17 34 0.92388 0 0.382683 35 0.653281 0.653281 0.382683 36 5.65713e-17 0.92388 0.382683 37 -0.653281 0.653281 0.382683 38 -0.92388 1.13143e-16 0.382683 39 -0.653281 -0.653281 0.382683 40 -1.69714e-16 -0.92388 0.382683 41 0.653281 -0.653281 0.382683 42 0.707107 0 0.707107 43 0.5 0.5 0.707107 44 4.32978e-17 0.707107 0.707107 45 -0.5 0.5 0.707107 46 -0.707107 8.65956e-17 0.707107 47 -0.5 -0.5 0.707107 48 -1.29893e-16 -0.707107 0.707107 49 0.5 -0.5 0.707107 50 0.382683 0 0.92388 51 0.270598 0.270598 0.92388 52 2.34326e-17 0.382683 0.92388 53 -0.270598 0.270598 0.92388 54 -0.382683 4.68652e-17 0.92388 55 -0.270598 -0.270598 0.92388 56 -7.02978e-17 -0.382683 0.92388 57 0.270598 -0.270598 0.92388 58 0 0 1 TEST20 SPHERE_GRID_Q9_ELEMENT sets up a grid of Q9 quadrilaterals on a sphere. SPHERE_GRID_Q9_ELEMENT_NUM returns the number of elements in the grid SPHERE_GRID_Q9_NODE_NUM returns the number of nodes in the grid. SPHERE_GRID_Q9_NODE_XYZ returns the coordinates of nodes in the grid. Expected number of nodes = 44 Expected number of elements = 12 The elements and their nodes: 1 1 1 10 8 1 4 9 2 3 2 1 1 12 10 1 6 11 4 5 3 1 1 8 12 1 2 13 6 7 4 8 10 22 20 9 16 21 14 15 5 10 12 24 22 11 18 23 16 17 6 12 8 20 24 13 14 25 18 19 7 20 22 34 32 21 28 33 26 27 8 22 24 36 34 23 30 35 28 29 9 24 20 32 36 25 26 37 30 31 10 32 34 44 44 33 40 44 38 39 11 34 36 44 44 35 42 44 40 41 12 36 32 44 44 37 38 44 42 43 The node coordinates: 1 0 0 -1 2 0.382683 0 -0.92388 3 0.191342 0.331414 -0.92388 4 -0.191342 0.331414 -0.92388 5 -0.382683 4.68652e-17 -0.92388 6 -0.191342 -0.331414 -0.92388 7 0.191342 -0.331414 -0.92388 8 0.707107 0 -0.707107 9 0.353553 0.612372 -0.707107 10 -0.353553 0.612372 -0.707107 11 -0.707107 8.65956e-17 -0.707107 12 -0.353553 -0.612372 -0.707107 13 0.353553 -0.612372 -0.707107 14 0.92388 0 -0.382683 15 0.46194 0.800103 -0.382683 16 -0.46194 0.800103 -0.382683 17 -0.92388 1.13143e-16 -0.382683 18 -0.46194 -0.800103 -0.382683 19 0.46194 -0.800103 -0.382683 20 1 0 6.12323e-17 21 0.5 0.866025 6.12323e-17 22 -0.5 0.866025 6.12323e-17 23 -1 1.22465e-16 6.12323e-17 24 -0.5 -0.866025 6.12323e-17 25 0.5 -0.866025 6.12323e-17 26 0.92388 0 0.382683 27 0.46194 0.800103 0.382683 28 -0.46194 0.800103 0.382683 29 -0.92388 1.13143e-16 0.382683 30 -0.46194 -0.800103 0.382683 31 0.46194 -0.800103 0.382683 32 0.707107 0 0.707107 33 0.353553 0.612372 0.707107 34 -0.353553 0.612372 0.707107 35 -0.707107 8.65956e-17 0.707107 36 -0.353553 -0.612372 0.707107 37 0.353553 -0.612372 0.707107 38 0.382683 0 0.92388 39 0.191342 0.331414 0.92388 40 -0.191342 0.331414 0.92388 41 -0.382683 4.68652e-17 0.92388 42 -0.191342 -0.331414 0.92388 43 0.191342 -0.331414 0.92388 44 0 0 1 TEST21 SPHERE_GRID_Q16_ELEMENT sets up a grid of Q16 quadrilaterals on a sphere. SPHERE_GRID_Q16_ELEMENT_NUM returns the number of elements in the grid SPHERE_GRID_Q16_NODE_NUM returns the number of nodes in the grid. SPHERE_GRID_Q16_NODE_XYZ returns the coordinates of nodes in the grid. Expected number of nodes = 32 Expected number of elements = 4 The elements and their nodes, listed in a way that suggests their geometry: 4 32 32 32 32 29 30 31 26 23 24 25 20 17 18 19 14 3 32 32 32 32 26 27 28 29 20 21 22 23 14 15 16 17 2 17 18 19 14 11 12 13 8 5 6 7 2 1 1 1 1 1 14 15 16 17 8 9 10 11 2 3 4 5 1 1 1 1 The node coordinates: 1 0 0 -1 2 0.5 0 -0.866025 3 0.25 0.433013 -0.866025 4 -0.25 0.433013 -0.866025 5 -0.5 6.12323e-17 -0.866025 6 -0.25 -0.433013 -0.866025 7 0.25 -0.433013 -0.866025 8 0.866025 0 -0.5 9 0.433013 0.75 -0.5 10 -0.433013 0.75 -0.5 11 -0.866025 1.06058e-16 -0.5 12 -0.433013 -0.75 -0.5 13 0.433013 -0.75 -0.5 14 1 0 6.12323e-17 15 0.5 0.866025 6.12323e-17 16 -0.5 0.866025 6.12323e-17 17 -1 1.22465e-16 6.12323e-17 18 -0.5 -0.866025 6.12323e-17 19 0.5 -0.866025 6.12323e-17 20 0.866025 0 0.5 21 0.433013 0.75 0.5 22 -0.433013 0.75 0.5 23 -0.866025 1.06058e-16 0.5 24 -0.433013 -0.75 0.5 25 0.433013 -0.75 0.5 26 0.5 0 0.866025 27 0.25 0.433013 0.866025 28 -0.25 0.433013 0.866025 29 -0.5 6.12323e-17 0.866025 30 -0.25 -0.433013 0.866025 31 0.25 -0.433013 0.866025 32 0 0 1 TEST22 SPHERE_GRID_T3_ELEMENT sets up a grid of T3 triangles on a sphere. SPHERE_GRID_T3_ELEMENT_NUM returns the number of elements in the grid SPHERE_GRID_T3_NODE_NUM returns the number of nodes in the grid. SPHERE_GRID_T3_NODE_XYZ returns the coordinates of nodes in the grid. Expected number of nodes = 58 Expected number of elements = 112 The elements and their nodes: 1 3 2 1 2 4 3 1 3 5 4 1 4 6 5 1 5 7 6 1 6 8 7 1 7 9 8 1 8 2 9 1 9 2 3 10 10 11 10 3 11 3 4 11 12 12 11 4 13 4 5 12 14 13 12 5 15 5 6 13 16 14 13 6 17 6 7 14 18 15 14 7 19 7 8 15 20 16 15 8 21 8 9 16 22 17 16 9 23 9 2 17 24 10 17 2 25 10 11 18 26 19 18 11 27 11 12 19 28 20 19 12 29 12 13 20 30 21 20 13 31 13 14 21 32 22 21 14 33 14 15 22 34 23 22 15 35 15 16 23 36 24 23 16 37 16 17 24 38 25 24 17 39 17 10 25 40 18 25 10 41 18 19 26 42 27 26 19 43 19 20 27 44 28 27 20 45 20 21 28 46 29 28 21 47 21 22 29 48 30 29 22 49 22 23 30 50 31 30 23 51 23 24 31 52 32 31 24 53 24 25 32 54 33 32 25 55 25 18 33 56 26 33 18 57 26 27 34 58 35 34 27 59 27 28 35 60 36 35 28 61 28 29 36 62 37 36 29 63 29 30 37 64 38 37 30 65 30 31 38 66 39 38 31 67 31 32 39 68 40 39 32 69 32 33 40 70 41 40 33 71 33 26 41 72 34 41 26 73 34 35 42 74 43 42 35 75 35 36 43 76 44 43 36 77 36 37 44 78 45 44 37 79 37 38 45 80 46 45 38 81 38 39 46 82 47 46 39 83 39 40 47 84 48 47 40 85 40 41 48 86 49 48 41 87 41 34 49 88 42 49 34 89 42 43 50 90 51 50 43 91 43 44 51 92 52 51 44 93 44 45 52 94 53 52 45 95 45 46 53 96 54 53 46 97 46 47 54 98 55 54 47 99 47 48 55 100 56 55 48 101 48 49 56 102 57 56 49 103 49 42 57 104 50 57 42 105 50 51 58 106 51 52 58 107 52 53 58 108 53 54 58 109 54 55 58 110 55 56 58 111 56 57 58 112 57 50 58 The node coordinates: 1 0 0 -1 2 0.382683 0 -0.92388 3 0.270598 0.270598 -0.92388 4 2.34326e-17 0.382683 -0.92388 5 -0.270598 0.270598 -0.92388 6 -0.382683 4.68652e-17 -0.92388 7 -0.270598 -0.270598 -0.92388 8 -7.02978e-17 -0.382683 -0.92388 9 0.270598 -0.270598 -0.92388 10 0.707107 0 -0.707107 11 0.5 0.5 -0.707107 12 4.32978e-17 0.707107 -0.707107 13 -0.5 0.5 -0.707107 14 -0.707107 8.65956e-17 -0.707107 15 -0.5 -0.5 -0.707107 16 -1.29893e-16 -0.707107 -0.707107 17 0.5 -0.5 -0.707107 18 0.92388 0 -0.382683 19 0.653281 0.653281 -0.382683 20 5.65713e-17 0.92388 -0.382683 21 -0.653281 0.653281 -0.382683 22 -0.92388 1.13143e-16 -0.382683 23 -0.653281 -0.653281 -0.382683 24 -1.69714e-16 -0.92388 -0.382683 25 0.653281 -0.653281 -0.382683 26 1 0 6.12323e-17 27 0.707107 0.707107 6.12323e-17 28 6.12323e-17 1 6.12323e-17 29 -0.707107 0.707107 6.12323e-17 30 -1 1.22465e-16 6.12323e-17 31 -0.707107 -0.707107 6.12323e-17 32 -1.83697e-16 -1 6.12323e-17 33 0.707107 -0.707107 6.12323e-17 34 0.92388 0 0.382683 35 0.653281 0.653281 0.382683 36 5.65713e-17 0.92388 0.382683 37 -0.653281 0.653281 0.382683 38 -0.92388 1.13143e-16 0.382683 39 -0.653281 -0.653281 0.382683 40 -1.69714e-16 -0.92388 0.382683 41 0.653281 -0.653281 0.382683 42 0.707107 0 0.707107 43 0.5 0.5 0.707107 44 4.32978e-17 0.707107 0.707107 45 -0.5 0.5 0.707107 46 -0.707107 8.65956e-17 0.707107 47 -0.5 -0.5 0.707107 48 -1.29893e-16 -0.707107 0.707107 49 0.5 -0.5 0.707107 50 0.382683 0 0.92388 51 0.270598 0.270598 0.92388 52 2.34326e-17 0.382683 0.92388 53 -0.270598 0.270598 0.92388 54 -0.382683 4.68652e-17 0.92388 55 -0.270598 -0.270598 0.92388 56 -7.02978e-17 -0.382683 0.92388 57 0.270598 -0.270598 0.92388 58 0 0 1 TEST23 SPHERE_GRID_T6_ELEMENT sets up a grid of T6 triangles on a sphere. SPHERE_GRID_T6_ELEMENT_NUM returns the number of elements in the grid SPHERE_GRID_T6_NODE_NUM returns the number of nodes in the grid SPHERE_GRID_T6_NODE_XYZ returns the coordinates of nodes in the grid. Expected number of nodes = 44 Expected number of elements = 18 The elements and their nodes: 1 10 8 1 9 3 4 2 12 10 1 11 5 6 3 8 12 1 13 7 2 4 8 10 20 9 15 14 5 22 20 10 21 15 16 6 10 12 22 11 17 16 7 24 22 12 23 17 18 8 12 8 24 13 19 18 9 20 24 8 25 19 14 10 20 22 32 21 27 26 11 34 32 22 33 27 28 12 22 24 34 23 29 28 13 36 34 24 35 29 30 14 24 20 36 25 31 30 15 32 36 20 37 31 26 16 32 34 44 33 39 38 17 34 36 44 35 41 40 18 36 32 44 37 43 42 The node coordinates: 1 0 0 -1 2 0.382683 0 -0.92388 3 0.191342 0.331414 -0.92388 4 -0.191342 0.331414 -0.92388 5 -0.382683 4.68652e-17 -0.92388 6 -0.191342 -0.331414 -0.92388 7 0.191342 -0.331414 -0.92388 8 0.707107 0 -0.707107 9 0.353553 0.612372 -0.707107 10 -0.353553 0.612372 -0.707107 11 -0.707107 8.65956e-17 -0.707107 12 -0.353553 -0.612372 -0.707107 13 0.353553 -0.612372 -0.707107 14 0.92388 0 -0.382683 15 0.46194 0.800103 -0.382683 16 -0.46194 0.800103 -0.382683 17 -0.92388 1.13143e-16 -0.382683 18 -0.46194 -0.800103 -0.382683 19 0.46194 -0.800103 -0.382683 20 1 0 6.12323e-17 21 0.5 0.866025 6.12323e-17 22 -0.5 0.866025 6.12323e-17 23 -1 1.22465e-16 6.12323e-17 24 -0.5 -0.866025 6.12323e-17 25 0.5 -0.866025 6.12323e-17 26 0.92388 0 0.382683 27 0.46194 0.800103 0.382683 28 -0.46194 0.800103 0.382683 29 -0.92388 1.13143e-16 0.382683 30 -0.46194 -0.800103 0.382683 31 0.46194 -0.800103 0.382683 32 0.707107 0 0.707107 33 0.353553 0.612372 0.707107 34 -0.353553 0.612372 0.707107 35 -0.707107 8.65956e-17 0.707107 36 -0.353553 -0.612372 0.707107 37 0.353553 -0.612372 0.707107 38 0.382683 0 0.92388 39 0.191342 0.331414 0.92388 40 -0.191342 0.331414 0.92388 41 -0.382683 4.68652e-17 0.92388 42 -0.191342 -0.331414 0.92388 43 0.191342 -0.331414 0.92388 44 0 0 1 TEST24 TRIANGLE_UNIT_SET sets up a quadrature in the unit triangle, A = 0 B = 0 Rule QUAD ERROR 1 1 0 2 1 0 3 1 0 4 1 0 5 1 0 6 1 1.11022e-16 7 1 0 8 1 9.99201e-16 9 1 2.22045e-16 10 1 0 11 1 6.66134e-16 12 1 1.9984e-15 13 1 1.11022e-15 14 1 0 15 1 2e-10 16 1 0 17 1 6.66134e-16 18 1 5.55112e-16 19 1 2.22045e-16 20 1 2.22045e-16 A = 0 B = 1 Rule QUAD ERROR 1 1 0 2 1 0 3 1 1.11022e-16 4 1 0 5 1 0 6 1 0 7 1 0 8 1 6.66134e-16 9 1 4.49929e-12 10 1 1.11022e-16 11 1 8.88178e-16 12 1 1.55431e-15 13 1 6.66134e-16 14 1 0 15 1 1.29636e-10 16 1 4.36073e-12 17 1 2.22045e-16 18 1 2.59701e-08 19 1 2.22045e-16 20 1 0 A = 0 B = 2 Rule QUAD ERROR 1 0.666667 0.333333 2 2 1 3 1 0 4 1 0 5 1 2.22045e-16 6 1 4.44089e-16 7 1 0 8 1 0 9 1 4.5508e-13 10 1 2.22045e-16 11 1 1.66533e-15 12 1 3.10862e-15 13 1 1.11022e-15 14 1 0 15 1 9.83384e-11 16 1 5.13056e-12 17 1 2.22045e-16 18 1 2.30576e-08 19 1 3.33067e-16 20 1 1.11022e-16 A = 0 B = 3 Rule QUAD ERROR 1 0.37037 0.62963 2 3.33333 2.33333 3 1.01852 0.0185185 4 0.833333 0.166667 5 1 2.22045e-16 6 1 7.77156e-16 7 1 1.11022e-16 8 1 1.11022e-15 9 1 2.70894e-14 10 1 2.22045e-16 11 1 2.33147e-15 12 1 4.44089e-15 13 1 1.77636e-15 14 1 1.11022e-16 15 1 7.88721e-11 16 1 3.68705e-12 17 1 4.44089e-16 18 1 1.27948e-08 19 1 4.44089e-16 20 1 2.22045e-16 A = 0 B = 4 Rule QUAD ERROR 1 0.185185 0.814815 2 5 4 3 0.99537 0.00462963 4 0.625 0.375 5 0.933333 0.0666667 6 0.958333 0.0416667 7 0.958333 0.0416667 8 1 1.9984e-15 9 1 6.66134e-16 10 1 4.44089e-16 11 1 3.10862e-15 12 1 5.32907e-15 13 1 2.22045e-15 14 1.08333 0.0833333 15 1 6.48875e-11 16 1 9.68114e-13 17 1 4.44089e-16 18 1 5.67996e-09 19 1 2.22045e-16 20 1 1.11022e-16 A = 0 B = 5 Rule QUAD ERROR 1 0.0864198 0.91358 2 7 6 3 0.923611 0.0763889 4 0.4375 0.5625 5 0.808889 0.191111 6 0.875 0.125 7 0.875 0.125 8 1.00527 0.00527393 9 0.985185 0.0148148 10 1 4.44089e-16 11 1 3.55271e-15 12 1 5.9952e-15 13 1 2.88658e-15 14 1.26389 0.263889 15 1 5.4827e-11 16 1 2.26719e-12 17 1 6.66134e-16 18 1 2.2063e-09 19 1 1.11022e-16 20 1 1.11022e-16 A = 0 B = 6 Rule QUAD ERROR 1 0.0384088 0.961591 2 9.33333 8.33333 3 0.819787 0.180213 4 0.291667 0.708333 5 0.660662 0.339338 6 0.766111 0.233889 7 0.766975 0.233025 8 1.01295 0.0129495 9 0.948082 0.0519177 10 0.995414 0.00458554 11 0.995927 0.00407261 12 1 6.88338e-15 13 1 2.88658e-15 14 1.53395 0.533951 15 1 4.75356e-11 16 1 5.61706e-12 17 1 1.33227e-15 18 1 7.83546e-10 19 1 2.22045e-16 20 1 1.11022e-16 A = 0 B = 7 Rule QUAD ERROR 1 0.0164609 0.983539 2 12 11 3 0.702418 0.297582 4 0.1875 0.8125 5 0.516101 0.483899 6 0.648333 0.351667 7 0.650926 0.349074 8 1.0169 0.0169016 9 0.889877 0.110123 10 0.981708 0.0182918 11 0.983203 0.0167971 12 1.0002 0.000198039 13 1 3.10862e-15 14 1.88241 0.882407 15 1 4.18393e-11 16 1 8.9202e-12 17 1 1.55431e-15 18 1 2.60878e-10 19 1 0 20 1 1.11022e-16 A = 0 B = 8 Rule QUAD ERROR 1 0.00685871 0.993141 2 15 14 3 0.585294 0.414706 4 0.117188 0.882812 5 0.389922 0.610078 6 0.533854 0.466146 7 0.538484 0.461516 8 1.01216 0.0121614 9 0.816165 0.183835 10 0.956602 0.0433976 11 0.959183 0.0408174 12 1.00044 0.000439377 13 0.999572 0.000427958 14 2.29996 1.29996 15 0.999733 0.000267057 16 1 1.21405e-11 17 1 1.77636e-15 18 1 8.27227e-11 19 1 3.33067e-16 20 1 2.22045e-16 A = 0 B = 9 Rule QUAD ERROR 1 0.00279429 0.997206 2 18.3333 17.3333 3 0.476896 0.523104 4 0.0716146 0.928385 5 0.287142 0.712858 6 0.429942 0.570058 7 0.436368 0.563632 8 0.996255 0.0037446 9 0.733569 0.266431 10 0.920084 0.079916 11 0.92354 0.0764596 12 1.00009 8.58161e-05 13 0.997886 0.00211358 14 2.7799 1.7799 15 0.998637 0.00136299 16 1 1.52949e-11 17 1.00009 9.14859e-05 18 1 2.52481e-11 19 1 6.66134e-16 20 1 1.11022e-16 A = 0 B = 10 Rule QUAD ERROR 1 0.00111772 0.998882 2 22 21 3 0.381514 0.618486 4 0.0429688 0.957031 5 0.207231 0.792769 6 0.339994 0.660006 7 0.34766 0.65234 8 0.968918 0.031082 9 0.647992 0.352008 10 0.873671 0.126329 11 0.877655 0.122345 12 0.998219 0.00178063 13 0.993975 0.00602508 14 3.31769 2.31769 15 0.995992 0.0040076 16 1 1.84178e-11 17 1.00049 0.000488606 18 0.99997 3.0264e-05 19 1 9.99201e-16 20 1 6.66134e-16 A = 1 B = 0 Rule QUAD ERROR 1 1 0 2 1 0 3 1 1.11022e-16 4 1 0 5 1 0 6 1 1.11022e-16 7 1 0 8 1 6.66134e-16 9 1 4.49929e-12 10 1 1.11022e-16 11 1 8.88178e-16 12 1 1.55431e-15 13 1 6.66134e-16 14 1 0 15 1 3.40727e-10 16 1 8.72202e-12 17 1 2.22045e-16 18 1 2.59701e-08 19 1 0 20 1 0 A = 1 B = 1 Rule QUAD ERROR 1 1.33333 0.333333 2 0 1 3 1 0 4 1 0 5 1 2.22045e-16 6 1 2.22045e-16 7 1 1.11022e-16 8 1 8.88178e-16 9 1 1.7542e-11 10 1 2.22045e-16 11 1 4.44089e-16 12 1 3.33067e-16 13 1 1.33227e-15 14 1 0 15 1 2.2353e-10 16 1 2.05103e-12 17 1 2.22045e-16 18 1 8.08229e-08 19 1 1.11022e-16 20 1 2.22045e-16 A = 1 B = 2 Rule QUAD ERROR 1 1.11111 0.111111 2 0 1 3 0.972222 0.0277778 4 1.25 0.25 5 1 4.44089e-16 6 1 2.22045e-16 7 1 1.11022e-16 8 1 4.44089e-16 9 1 1.4595e-11 10 1 2.22045e-16 11 1 4.44089e-16 12 1 8.88178e-16 13 1 8.88178e-16 14 1 0 15 1 1.76203e-10 16 1 1.09059e-11 17 1 0 18 1 1.16362e-07 19 1 0 20 1 2.22045e-16 A = 1 B = 3 Rule QUAD ERROR 1 0.740741 0.259259 2 0 1 3 1.06481 0.0648148 4 1.25 0.25 5 1.13333 0.133333 6 1.08333 0.0833333 7 1.08333 0.0833333 8 1 8.88178e-16 9 1 1.86805e-11 10 1 4.44089e-16 11 1 7.77156e-16 12 1 1.9984e-15 13 1 0 14 0.833333 0.166667 15 1 1.48797e-10 16 1 1.72808e-11 17 1 3.33067e-16 18 1 1.35627e-07 19 1 1.11022e-16 20 1 2.22045e-16 A = 1 B = 4 Rule QUAD ERROR 1 0.432099 0.567901 2 0 1 3 1.17477 0.174769 4 1.09375 0.09375 5 1.24444 0.244444 6 1.16667 0.166667 7 1.16667 0.166667 8 0.986815 0.0131848 9 1.03704 0.037037 10 1 4.44089e-16 11 1 1.88738e-15 12 1 2.88658e-15 13 1 8.88178e-16 14 0.631944 0.368056 15 1 1.25252e-10 16 1 2.03799e-11 17 1 2.22045e-16 18 1 1.52635e-07 19 1 2.22045e-16 20 1 0 A = 1 B = 5 Rule QUAD ERROR 1 0.230453 0.769547 2 0 1 3 1.23508 0.235082 4 0.875 0.125 5 1.25357 0.25357 6 1.20167 0.201667 7 1.19907 0.199074 8 0.982247 0.0177527 9 1.09649 0.0964938 10 1.01376 0.0137566 11 1.01222 0.0122178 12 1 3.77476e-15 13 1 1.55431e-15 14 0.453704 0.546296 15 1 1.05866e-10 16 1 2.11806e-11 17 1 2.22045e-16 18 1 1.69535e-07 19 1 1.11022e-16 20 1 2.22045e-16 A = 1 B = 6 Rule QUAD ERROR 1 0.115226 0.884774 2 0 1 3 1.23058 0.230581 4 0.65625 0.34375 5 1.16663 0.166625 6 1.17833 0.178333 7 1.17315 0.173148 8 0.999117 0.000882959 9 1.1518 0.151802 10 1.04339 0.0433862 11 1.04046 0.040463 12 0.999307 0.000693137 13 1 1.9984e-15 14 0.314352 0.685648 15 1 9.31069e-11 16 1 2.08102e-11 17 1 4.44089e-16 18 1 1.83598e-07 19 1 3.33067e-16 20 1 0 A = 1 B = 7 Rule QUAD ERROR 1 0.0548697 0.94513 2 0 1 3 1.17091 0.17091 4 0.46875 0.53125 5 1.02082 0.0208158 6 1.10625 0.10625 7 1.10069 0.100694 8 1.03586 0.0358626 9 1.18472 0.184724 10 1.08213 0.0821314 11 1.07928 0.0792842 12 0.999233 0.000767312 13 1.00171 0.00171183 14 0.212191 0.787809 15 1.00057 0.000566893 16 1 2.00624e-11 17 1 5.55112e-16 18 1 1.92514e-07 19 1 2.22045e-16 20 1 2.22045e-16 A = 1 B = 8 Rule QUAD ERROR 1 0.0251486 0.974851 2 0 1 3 1.07309 0.073089 4 0.322266 0.677734 5 0.852434 0.147566 6 1.00146 0.00145833 7 0.998007 0.00199331 8 1.08374 0.0837383 9 1.18784 0.187843 10 1.12094 0.120935 11 1.11957 0.119573 12 1.00203 0.0020304 13 1.00716 0.00715736 14 0.140223 0.859777 15 1.00253 0.00253185 16 1 1.94018e-11 17 0.999588 0.000411687 18 1 1.9552e-07 19 1 8.88178e-16 20 1 0 A = 1 B = 9 Rule QUAD ERROR 1 0.0111772 0.988823 2 0 1 3 0.953802 0.0461977 4 0.214844 0.785156 5 0.686695 0.313305 6 0.879682 0.120318 7 0.879906 0.120094 8 1.13294 0.132942 9 1.16146 0.161458 10 1.15215 0.152148 11 1.15297 0.152968 12 1.00942 0.00941807 13 1.01744 0.0174439 14 0.0909672 0.909033 15 1.00658 0.00658406 16 1 1.90594e-11 17 0.998106 0.00189411 18 1.00015 0.000151127 19 1 1.33227e-15 20 1 1.11022e-16 A = 2 B = 0 Rule QUAD ERROR 1 0.666667 0.333333 2 2 1 3 1 0 4 1 0 5 1 2.22045e-16 6 1 4.44089e-16 7 1 0 8 1 0 9 1 4.54969e-13 10 1 2.22045e-16 11 1 1.66533e-15 12 1 3.10862e-15 13 1 1.11022e-15 14 1 0 15 1 4.57924e-10 16 1 1.94944e-11 17 1 4.44089e-16 18 1 2.30576e-08 19 1 3.33067e-16 20 1 0 A = 2 B = 1 Rule QUAD ERROR 1 1.11111 0.111111 2 0 1 3 0.972222 0.0277778 4 1.25 0.25 5 1 2.22045e-16 6 1 2.22045e-16 7 1 2.22045e-16 8 1 4.44089e-16 9 1 1.45951e-11 10 1 2.22045e-16 11 1 5.55112e-16 12 1 8.88178e-16 13 1 8.88178e-16 14 1 0 15 1 2.94521e-10 16 1 1.12316e-11 17 1 0 18 1 1.16362e-07 19 1 1.11022e-16 20 1 2.22045e-16 A = 2 B = 2 Rule QUAD ERROR 1 1.11111 0.111111 2 0 1 3 0.763889 0.236111 4 1.875 0.875 5 0.8 0.2 6 0.875 0.125 7 0.875 0.125 8 1 4.44089e-16 9 1 4.10005e-12 10 1 2.22045e-16 11 1 0 12 1 2.66454e-15 13 1 8.88178e-16 14 1.25 0.25 15 1 2.31016e-10 16 1 1.84386e-12 17 1 2.22045e-16 18 1 1.90479e-07 19 1 2.22045e-16 20 1 0 A = 2 B = 3 Rule QUAD ERROR 1 0.864198 0.135802 2 0 1 3 0.729167 0.270833 4 2.1875 1.1875 5 0.808889 0.191111 6 0.875 0.125 7 0.875 0.125 8 1.00527 0.00527393 9 0.985185 0.0148148 10 1 5.55112e-16 11 1 2.22045e-16 12 1 3.9968e-15 13 1 1.11022e-16 14 1.26389 0.263889 15 1 2.0766e-10 16 1 9.53349e-12 17 1 0 18 1 2.38316e-07 19 1 3.33067e-16 20 1 0 A = 2 B = 4 Rule QUAD ERROR 1 0.576132 0.423868 2 0 1 3 0.819187 0.180813 4 2.1875 1.1875 5 0.949926 0.0500741 6 0.9625 0.0375 7 0.956019 0.0439815 8 0.958434 0.0415657 9 1.03012 0.0301235 10 0.962963 0.037037 11 0.983741 0.0162589 12 1 5.55112e-15 13 1 8.88178e-16 14 1.13426 0.134259 15 1 1.83407e-10 16 1 1.79778e-11 17 1 1.11022e-16 18 1 2.87364e-07 19 1 1.11022e-16 20 1 1.11022e-16 A = 2 B = 5 Rule QUAD ERROR 1 0.345679 0.654321 2 0 1 3 0.937114 0.0628858 4 1.96875 0.96875 5 1.08068 0.0806756 6 1.0675 0.0675 7 1.04028 0.0402778 8 0.881502 0.118498 9 1.14074 0.140741 10 0.909751 0.0902494 11 0.96727 0.0327297 12 1.00089 0.000891176 13 1 1.55431e-15 14 0.943056 0.0569444 15 1 1.50525e-10 16 1 2.24776e-11 17 1 2.22045e-16 18 1 3.4683e-07 19 1 4.44089e-16 20 1 3.33067e-16 A = 2 B = 6 Rule QUAD ERROR 1 0.192044 0.807956 2 0 1 3 1.02849 0.0284851 4 1.64062 0.640625 5 1.1335 0.133495 6 1.14479 0.144792 7 1.0897 0.0896991 8 0.811392 0.188608 9 1.28639 0.286387 10 0.866501 0.133499 11 0.970273 0.0297268 12 0.996224 0.00377556 13 0.994117 0.00588301 14 0.74267 0.25733 15 0.997732 0.00226757 16 1 2.38005e-11 17 1 4.44089e-16 18 1 4.11141e-07 19 1 8.88178e-16 20 1 0 A = 2 B = 7 Rule QUAD ERROR 1 0.100594 0.899406 2 0 1 3 1.07412 0.0741205 4 1.28906 0.289062 5 1.10343 0.103432 6 1.17562 0.175625 7 1.09561 0.0956147 8 0.769683 0.230317 9 1.43345 0.433446 10 0.84609 0.15391 11 1.00117 0.00116637 12 0.981678 0.0183219 13 0.979223 0.0207775 14 0.560892 0.439108 15 0.990804 0.00919627 16 1 2.30351e-11 17 1.00073 0.000729094 18 1 4.71514e-07 19 1 1.55431e-15 20 1 0 A = 2 B = 8 Rule QUAD ERROR 1 0.0502972 0.949703 2 0 1 3 1.07329 0.0732855 4 0.966797 0.0332031 5 1.01255 0.0125542 6 1.15872 0.158724 7 1.06264 0.0626366 8 0.761318 0.238682 9 1.55761 0.557612 10 0.849523 0.150477 11 1.057 0.0570024 12 0.957945 0.0420547 13 0.956754 0.043246 14 0.409352 0.590648 15 0.978158 0.0218419 16 1 2.11142e-11 17 1.00232 0.00232291 18 0.999302 0.000697842 19 1 1.33227e-15 20 1 3.33067e-16 A = 3 B = 0 Rule QUAD ERROR 1 0.37037 0.62963 2 3.33333 2.33333 3 1.01852 0.0185185 4 0.833333 0.166667 5 1 2.22045e-16 6 1 6.66134e-16 7 1 1.11022e-16 8 1 1.11022e-15 9 1 2.70894e-14 10 1 2.22045e-16 11 1 2.33147e-15 12 1 4.44089e-15 13 1 1.55431e-15 14 1 1.11022e-16 15 1 5.6686e-10 16 1 2.50031e-11 17 1 2.22045e-16 18 1 1.27948e-08 19 1 4.44089e-16 20 1 2.22045e-16 A = 3 B = 1 Rule QUAD ERROR 1 0.740741 0.259259 2 0 1 3 1.06481 0.0648148 4 1.25 0.25 5 1.13333 0.133333 6 1.08333 0.0833333 7 1.08333 0.0833333 8 1 8.88178e-16 9 1 1.86805e-11 10 1 4.44089e-16 11 1 8.88178e-16 12 1 1.9984e-15 13 1 2.22045e-16 14 0.833333 0.166667 15 1 3.58026e-10 16 1 2.06191e-11 17 1 3.33067e-16 18 1 1.35627e-07 19 1 2.22045e-16 20 1 1.11022e-16 A = 3 B = 2 Rule QUAD ERROR 1 0.864198 0.135802 2 0 1 3 0.729167 0.270833 4 2.1875 1.1875 5 0.808889 0.191111 6 0.875 0.125 7 0.875 0.125 8 1.00527 0.00527393 9 0.985185 0.0148148 10 1 5.55112e-16 11 1 0 12 1 4.21885e-15 13 1 1.11022e-16 14 1.26389 0.263889 15 1 2.62157e-10 16 1 1.70136e-11 17 1 0 18 1 2.38316e-07 19 1 2.22045e-16 20 1 0 A = 3 B = 3 Rule QUAD ERROR 1 0.768176 0.231824 2 0 1 3 0.516118 0.483882 4 2.91667 1.91667 5 0.594568 0.405432 6 0.738889 0.261111 7 0.756173 0.243827 8 1.0635 0.0634958 9 0.917202 0.0827984 10 1.05115 0.0511464 11 1.01215 0.0121547 12 1 6.43929e-15 13 1 8.88178e-16 14 1.51235 0.512346 15 1 2.48082e-10 16 1 4.53981e-12 17 1 2.22045e-16 18 1 2.91796e-07 19 1 3.33067e-16 20 1 2.22045e-16 A = 3 B = 4 Rule QUAD ERROR 1 0.576132 0.423868 2 0 1 3 0.48161 0.51839 4 3.28125 2.28125 5 0.591526 0.408474 6 0.729167 0.270833 7 0.761574 0.238426 8 1.05737 0.0573704 9 0.905679 0.094321 10 1.05367 0.0536659 11 0.999136 0.000864444 12 0.999505 0.000495098 13 1 1.77636e-15 14 1.57176 0.571759 15 1 2.4917e-10 16 1 8.9786e-12 17 1 0 18 1 3.40438e-07 19 1 1.11022e-16 20 1 0 A = 3 B = 5 Rule QUAD ERROR 1 0.384088 0.615912 2 0 1 3 0.544624 0.455376 4 3.28125 2.28125 5 0.694511 0.305489 6 0.802083 0.197917 7 0.818287 0.181713 8 0.972131 0.0278694 9 0.992329 0.00767078 10 0.992477 0.00752259 11 0.961034 0.0389662 12 1.01154 0.0115368 13 1.01166 0.0116576 14 1.48534 0.48534 15 1.00529 0.00529101 16 1 1.939e-11 17 1 2.22045e-16 18 1 4.14554e-07 19 1 1.11022e-16 20 1 0 A = 3 B = 6 Rule QUAD ERROR 1 0.23472 0.76528 2 0 1 3 0.635854 0.364146 4 3.00781 2.00781 5 0.800594 0.199406 6 0.903681 0.0963194 7 0.87305 0.12695 8 0.846318 0.153682 9 1.16057 0.16057 10 0.896388 0.103612 11 0.920984 0.0790163 12 1.03039 0.0303923 13 1.02893 0.0289321 14 1.30875 0.308749 15 1.01621 0.016209 16 1 2.58418e-11 17 0.999369 0.000630624 18 0.999999 5.16391e-07 19 1 4.44089e-16 20 1 3.33067e-16 A = 3 B = 7 Rule QUAD ERROR 1 0.134126 0.865874 2 0 1 3 0.718176 0.281824 4 2.57812 1.57812 5 0.86017 0.13983 6 0.995347 0.00465278 7 0.904171 0.095829 8 0.720497 0.279503 9 1.37501 0.375005 10 0.79708 0.20292 11 0.902218 0.0977818 12 1.04271 0.0427062 13 1.03934 0.0393434 14 1.09161 0.0916067 15 1.02814 0.0281417 16 1 2.87974e-11 17 0.999861 0.000139383 18 1.00188 0.00188073 19 1 6.66134e-16 20 1 2.22045e-16 A = 4 B = 0 Rule QUAD ERROR 1 0.185185 0.814815 2 5 4 3 0.99537 0.00462963 4 0.625 0.375 5 0.933333 0.0666667 6 0.958333 0.0416667 7 0.958333 0.0416667 8 1 1.9984e-15 9 1 6.66134e-16 10 1 4.44089e-16 11 1 3.10862e-15 12 1 5.32907e-15 13 1 2.22045e-15 14 1.08333 0.0833333 15 1 6.71276e-10 16 1 2.71948e-11 17 1 4.44089e-16 18 1 5.67996e-09 19 1 3.33067e-16 20 1 2.22045e-16 A = 4 B = 1 Rule QUAD ERROR 1 0.432099 0.567901 2 0 1 3 1.17477 0.174769 4 1.09375 0.09375 5 1.24444 0.244444 6 1.16667 0.166667 7 1.16667 0.166667 8 0.986815 0.0131848 9 1.03704 0.037037 10 1 4.44089e-16 11 1 1.88738e-15 12 1 2.88658e-15 13 1 6.66134e-16 14 0.631944 0.368056 15 1 4.29929e-10 16 1 2.33239e-11 17 1 2.22045e-16 18 1 1.52635e-07 19 1 2.22045e-16 20 1 1.11022e-16 A = 4 B = 2 Rule QUAD ERROR 1 0.576132 0.423868 2 0 1 3 0.819187 0.180813 4 2.1875 1.1875 5 0.949926 0.0500741 6 0.9625 0.0375 7 0.956019 0.0439815 8 0.958434 0.0415657 9 1.03012 0.0301235 10 0.962963 0.037037 11 0.983741 0.0162589 12 1 5.32907e-15 13 1 8.88178e-16 14 1.13426 0.134259 15 1 2.76231e-10 16 1 2.94873e-11 17 1 2.22045e-16 18 1 2.87364e-07 19 1 2.22045e-16 20 1 1.11022e-16 A = 4 B = 3 Rule QUAD ERROR 1 0.576132 0.423868 2 0 1 3 0.48161 0.51839 4 3.28125 2.28125 5 0.591526 0.408474 6 0.729167 0.270833 7 0.761574 0.238426 8 1.05737 0.0573704 9 0.905679 0.094321 10 1.05367 0.0536659 11 0.999136 0.000864444 12 0.999505 0.000495098 13 1 1.77636e-15 14 1.57176 0.571759 15 1 2.46721e-10 16 1 2.14381e-11 17 1 2.22045e-16 18 1 3.40438e-07 19 1 1.11022e-16 20 1 0 A = 4 B = 4 Rule QUAD ERROR 1 0.48011 0.51989 2 0 1 3 0.320698 0.679302 4 4.10156 3.10156 5 0.414538 0.585462 6 0.601563 0.398437 7 0.698785 0.301215 8 1.12991 0.129907 9 0.799671 0.200329 10 1.13775 0.137746 11 1.00434 0.00433838 12 0.983592 0.0164076 13 0.985455 0.0145449 14 1.85667 0.856674 15 0.992063 0.00793651 16 1 6.63858e-12 17 1 0 18 1 3.60003e-07 19 1 0 20 1 2.22045e-16 A = 4 B = 5 Rule QUAD ERROR 1 0.35208 0.64792 2 0 1 3 0.29363 0.70637 4 4.51172 3.51172 5 0.402555 0.597445 6 0.585521 0.414479 7 0.715439 0.284561 8 1.1061 0.106105 9 0.796322 0.203678 10 1.14249 0.14249 11 0.9744 0.0255996 12 0.982741 0.0172591 13 0.985505 0.0144954 14 1.96312 0.963124 15 0.986185 0.0138154 16 1 8.09952e-12 17 1.00018 0.000175988 18 1 4.03388e-07 19 1 4.44089e-16 20 1 0 A = 4 B = 6 Rule QUAD ERROR 1 0.23472 0.76528 2 0 1 3 0.332597 0.667403 4 4.51172 3.51172 5 0.467461 0.532539 6 0.643003 0.356997 7 0.750509 0.249491 8 0.998828 0.0011716 9 0.905044 0.094956 10 1.06917 0.0691732 11 0.916307 0.0836933 12 1.00786 0.00785815 13 1.01024 0.0102431 14 1.91031 0.910312 15 0.992343 0.00765656 16 1 1.99307e-11 17 0.997408 0.00259195 18 0.996691 0.00330926 19 1 2.22045e-16 20 1 0 A = 5 B = 0 Rule QUAD ERROR 1 0.0864198 0.91358 2 7 6 3 0.923611 0.0763889 4 0.4375 0.5625 5 0.808889 0.191111 6 0.875 0.125 7 0.875 0.125 8 1.00527 0.00527393 9 0.985185 0.0148148 10 1 4.44089e-16 11 1 3.55271e-15 12 1 5.9952e-15 13 1 2.44249e-15 14 1.26389 0.263889 15 1 7.67816e-10 16 1 2.87432e-11 17 1 6.66134e-16 18 1 2.2063e-09 19 1 1.11022e-16 20 1 1.11022e-16 A = 5 B = 1 Rule QUAD ERROR 1 0.230453 0.769547 2 0 1 3 1.23508 0.235082 4 0.875 0.125 5 1.25357 0.25357 6 1.20167 0.201667 7 1.19907 0.199074 8 0.982247 0.0177527 9 1.09649 0.0964938 10 1.01376 0.0137566 11 1.01222 0.0122178 12 1 3.55271e-15 13 1 1.77636e-15 14 0.453704 0.546296 15 1 5.22148e-10 16 1 1.96254e-11 17 1 3.33067e-16 18 1 1.69535e-07 19 1 2.22045e-16 20 1 1.11022e-16 A = 5 B = 2 Rule QUAD ERROR 1 0.345679 0.654321 2 0 1 3 0.937114 0.0628858 4 1.96875 0.96875 5 1.08068 0.0806756 6 1.0675 0.0675 7 1.04028 0.0402778 8 0.881502 0.118498 9 1.14074 0.140741 10 0.909751 0.0902494 11 0.96727 0.0327297 12 1.00089 0.000891176 13 1 1.33227e-15 14 0.943056 0.0569444 15 1 2.99839e-10 16 1 3.5927e-11 17 1 2.22045e-16 18 1 3.4683e-07 19 1 5.55112e-16 20 1 2.22045e-16 A = 5 B = 3 Rule QUAD ERROR 1 0.384088 0.615912 2 0 1 3 0.544624 0.455376 4 3.28125 2.28125 5 0.694511 0.305489 6 0.802083 0.197917 7 0.818287 0.181713 8 0.972131 0.0278694 9 0.992329 0.00767078 10 0.992477 0.00752259 11 0.961034 0.0389662 12 1.01154 0.0115368 13 1.01166 0.0116576 14 1.48534 0.48534 15 1.00794 0.00793651 16 1 3.62368e-11 17 1 2.22045e-16 18 1 4.14554e-07 19 1 2.22045e-16 20 1 1.11022e-16 A = 5 B = 4 Rule QUAD ERROR 1 0.35208 0.64792 2 0 1 3 0.29363 0.70637 4 4.51172 3.51172 5 0.402555 0.597445 6 0.585521 0.414479 7 0.715439 0.284561 8 1.1061 0.106105 9 0.796322 0.203678 10 1.14249 0.14249 11 0.9744 0.0255996 12 0.982741 0.0172591 13 0.985505 0.0144954 14 1.96312 0.963124 15 0.999118 0.000881834 16 1 2.4324e-11 17 1.00018 0.000175988 18 1 4.03388e-07 19 1 2.22045e-16 20 1 0 A = 5 B = 5 Rule QUAD ERROR 1 0.281664 0.718336 2 0 1 3 0.187868 0.812132 4 5.41406 4.41406 5 0.273552 0.726448 6 0.473229 0.526771 7 0.710487 0.289513 8 1.16618 0.166185 9 0.679036 0.320964 10 1.23958 0.239579 11 0.968337 0.0316626 12 0.950326 0.0496743 13 0.955574 0.044426 14 2.29237 1.29237 15 0.977562 0.0224378 16 1 8.46434e-12 17 1.00302 0.00301864 18 1.00398 0.00397668 19 1 4.44089e-16 20 1 0 A = 6 B = 0 Rule QUAD ERROR 1 0.0384088 0.961591 2 9.33333 8.33333 3 0.819787 0.180213 4 0.291667 0.708333 5 0.660662 0.339338 6 0.766111 0.233889 7 0.766975 0.233025 8 1.01295 0.0129495 9 0.948082 0.0519177 10 0.995414 0.00458554 11 0.995927 0.00407261 12 1 6.66134e-15 13 1 3.10862e-15 14 1.53395 0.533951 15 1 8.49704e-10 16 1 3.17828e-11 17 1 1.55431e-15 18 1 7.83546e-10 19 1 2.22045e-16 20 1 1.11022e-16 A = 6 B = 1 Rule QUAD ERROR 1 0.115226 0.884774 2 0 1 3 1.23058 0.230581 4 0.65625 0.34375 5 1.16663 0.166625 6 1.17833 0.178333 7 1.17315 0.173148 8 0.999117 0.000882959 9 1.1518 0.151802 10 1.04339 0.0433862 11 1.04046 0.040463 12 0.999307 0.000693137 13 1 2.22045e-15 14 0.314352 0.685648 15 1 6.33303e-10 16 1 1.1474e-11 17 1 4.44089e-16 18 1 1.83598e-07 19 1 3.33067e-16 20 1 0 A = 6 B = 2 Rule QUAD ERROR 1 0.192044 0.807956 2 0 1 3 1.02849 0.0284851 4 1.64062 0.640625 5 1.1335 0.133495 6 1.14479 0.144792 7 1.0897 0.0896991 8 0.811392 0.188608 9 1.28639 0.286387 10 0.866501 0.133499 11 0.970273 0.0297268 12 0.996224 0.00377556 13 0.994117 0.00588301 14 0.74267 0.25733 15 0.994709 0.00529101 16 1 3.57204e-11 17 1 3.33067e-16 18 1 4.11141e-07 19 1 7.77156e-16 20 1 0 A = 6 B = 3 Rule QUAD ERROR 1 0.23472 0.76528 2 0 1 3 0.635854 0.364146 4 3.00781 2.00781 5 0.800594 0.199406 6 0.903681 0.0963194 7 0.87305 0.12695 8 0.846318 0.153682 9 1.16057 0.16057 10 0.896388 0.103612 11 0.920984 0.0790163 12 1.03039 0.0303923 13 1.02893 0.0289321 14 1.30875 0.308749 15 1.01529 0.0152851 16 1 4.61646e-11 17 0.999369 0.000630624 18 0.999999 5.16391e-07 19 1 4.44089e-16 20 1 1.11022e-16 A = 6 B = 4 Rule QUAD ERROR 1 0.23472 0.76528 2 0 1 3 0.332597 0.667403 4 4.51172 3.51172 5 0.467461 0.532539 6 0.643003 0.356997 7 0.750509 0.249491 8 0.998828 0.0011716 9 0.905044 0.094956 10 1.06917 0.0691732 11 0.916307 0.0836933 12 1.00786 0.00785815 13 1.01024 0.0102431 14 1.91031 0.910312 15 1.02067 0.0206741 16 1 4.01833e-11 17 0.997408 0.00259195 18 0.996691 0.00330926 19 1 2.22045e-16 20 1 0 A = 7 B = 0 Rule QUAD ERROR 1 0.0164609 0.983539 2 12 11 3 0.702418 0.297582 4 0.1875 0.8125 5 0.516101 0.483899 6 0.648333 0.351667 7 0.650926 0.349074 8 1.0169 0.0169016 9 0.889877 0.110123 10 0.981708 0.0182918 11 0.983203 0.0167971 12 1.0002 0.000198039 13 1 3.10862e-15 14 1.88241 0.882407 15 1 9.11534e-10 16 1 3.75857e-11 17 1 1.33227e-15 18 1 2.60878e-10 19 1 0 20 1 0 A = 7 B = 1 Rule QUAD ERROR 1 0.0548697 0.94513 2 0 1 3 1.17091 0.17091 4 0.46875 0.53125 5 1.02082 0.0208158 6 1.10625 0.10625 7 1.10069 0.100694 8 1.03586 0.0358626 9 1.18472 0.184724 10 1.08213 0.0821314 11 1.07928 0.0792842 12 0.999233 0.000767312 13 1.00171 0.00171183 14 0.212191 0.787809 15 1.00227 0.00226757 16 1 1.08302e-12 17 1 6.66134e-16 18 1 1.92514e-07 19 1 2.22045e-16 20 1 0 A = 7 B = 2 Rule QUAD ERROR 1 0.100594 0.899406 2 0 1 3 1.07412 0.0741205 4 1.28906 0.289062 5 1.10343 0.103432 6 1.17562 0.175625 7 1.09561 0.0956147 8 0.769683 0.230317 9 1.43345 0.433446 10 0.84609 0.15391 11 1.00117 0.00116637 12 0.981678 0.0183219 13 0.979223 0.0207775 14 0.560892 0.439108 15 0.982951 0.0170488 16 1 2.97523e-11 17 1.00073 0.000729094 18 1 4.71514e-07 19 1 1.44329e-15 20 1 0 A = 7 B = 3 Rule QUAD ERROR 1 0.134126 0.865874 2 0 1 3 0.718176 0.281824 4 2.57812 1.57812 5 0.86017 0.13983 6 0.995347 0.00465278 7 0.904171 0.095829 8 0.720497 0.279503 9 1.37501 0.375005 10 0.79708 0.20292 11 0.902218 0.0977818 12 1.04271 0.0427062 13 1.03934 0.0393434 14 1.09161 0.0916067 15 1.01144 0.0114358 16 1 5.0437e-11 17 0.999861 0.000139383 18 1.00188 0.00188073 19 1 6.66134e-16 20 1 1.11022e-16 A = 8 B = 0 Rule QUAD ERROR 1 0.00685871 0.993141 2 15 14 3 0.585294 0.414706 4 0.117188 0.882812 5 0.389922 0.610078 6 0.533854 0.466146 7 0.538484 0.461516 8 1.01216 0.0121614 9 0.816165 0.183835 10 0.956602 0.0433976 11 0.959183 0.0408174 12 1.00044 0.000439377 13 0.999572 0.000427958 14 2.29996 1.29996 15 0.999433 0.000566894 16 1 4.67112e-11 17 1 1.55431e-15 18 1 8.27227e-11 19 1 3.33067e-16 20 1 3.33067e-16 A = 8 B = 1 Rule QUAD ERROR 1 0.0251486 0.974851 2 0 1 3 1.07309 0.073089 4 0.322266 0.677734 5 0.852434 0.147566 6 1.00146 0.00145833 7 0.998007 0.00199331 8 1.08374 0.0837383 9 1.18784 0.187843 10 1.12094 0.120935 11 1.11957 0.119573 12 1.00203 0.0020304 13 1.00716 0.00715736 14 0.140223 0.859777 15 1.00951 0.00951121 16 1 9.66871e-12 17 0.999588 0.000411687 18 1 1.9552e-07 19 1 8.88178e-16 20 1 2.22045e-16 A = 8 B = 2 Rule QUAD ERROR 1 0.0502972 0.949703 2 0 1 3 1.07329 0.0732855 4 0.966797 0.0332031 5 1.01255 0.0125542 6 1.15872 0.158724 7 1.06264 0.0626366 8 0.761318 0.238682 9 1.55761 0.557612 10 0.849523 0.150477 11 1.057 0.0570024 12 0.957945 0.0420547 13 0.956754 0.043246 14 0.409352 0.590648 15 0.968709 0.0312911 16 1 1.94103e-11 17 1.00232 0.00232291 18 0.999302 0.000697842 19 1 1.33227e-15 20 1 2.22045e-16 A = 9 B = 0 Rule QUAD ERROR 1 0.00279429 0.997206 2 18.3333 17.3333 3 0.476896 0.523104 4 0.0716146 0.928385 5 0.287142 0.712858 6 0.429942 0.570058 7 0.436368 0.563632 8 0.996255 0.0037446 9 0.733569 0.266431 10 0.920084 0.079916 11 0.92354 0.0764596 12 1.00009 8.58161e-05 13 0.997886 0.00211358 14 2.7799 1.7799 15 0.997194 0.00280647 16 1 5.92401e-11 17 1.00009 9.14859e-05 18 1 2.5248e-11 19 1 6.66134e-16 20 1 2.22045e-16 A = 9 B = 1 Rule QUAD ERROR 1 0.0111772 0.988823 2 0 1 3 0.953802 0.0461977 4 0.214844 0.785156 5 0.686695 0.313305 6 0.879682 0.120318 7 0.879906 0.120094 8 1.13294 0.132942 9 1.16146 0.161458 10 1.15215 0.152148 11 1.15297 0.152968 12 1.00942 0.00941807 13 1.01744 0.0174439 14 0.0909672 0.909033 15 1.02311 0.023112 16 1 1.93612e-11 17 0.998106 0.00189411 18 1.00015 0.000151127 19 1 1.55431e-15 20 1 1.11022e-16 A = 10 B = 0 Rule QUAD ERROR 1 0.00111772 0.998882 2 22 21 3 0.381514 0.618486 4 0.0429688 0.957031 5 0.207231 0.792769 6 0.339994 0.660006 7 0.34766 0.65234 8 0.968918 0.031082 9 0.647992 0.352008 10 0.873671 0.126329 11 0.877655 0.122345 12 0.998219 0.00178063 13 0.993975 0.00602508 14 3.31769 2.31769 15 0.99201 0.00799017 16 1 7.49603e-11 17 1.00049 0.000488606 18 0.99997 3.0264e-05 19 1 9.99201e-16 20 1 4.44089e-16 FEM2D_PACK_PRB: Normal end of execution. 05 July 2013 10:13:07 AM