07-Jan-2022 20:01:54 fem2d_pack_test(): MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2 Test fem2d_pack(). fem2d_pack_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 # 1 2 3 4 5 6 1 1 3 11 2 7 6 2 13 11 3 12 7 8 3 3 5 13 4 9 8 4 15 13 5 14 9 10 5 11 13 21 12 17 16 6 23 21 13 22 17 18 7 13 15 23 14 19 18 8 25 23 15 24 19 20 9 21 23 31 22 27 26 10 33 31 23 32 27 28 11 23 25 33 24 29 28 12 35 33 25 34 29 30 13 31 33 41 32 37 36 14 43 41 33 42 37 38 15 33 35 43 34 39 38 16 45 43 35 44 39 40 17 41 43 51 42 47 46 18 53 51 43 52 47 48 19 43 45 53 44 49 48 20 55 53 45 54 49 50 21 51 53 61 52 57 56 22 63 61 53 62 57 58 23 53 55 63 54 59 58 24 65 63 55 64 59 60 Lower bandwidth ML = 10 Upper bandwidth MU = 10 Total bandwidth M = 21 fem2d_pack_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 # 1 2 3 4 5 6 1 1 3 11 2 7 6 2 13 11 3 12 7 8 3 3 5 13 4 9 8 4 15 13 5 14 9 10 5 11 13 21 12 17 16 6 23 21 13 22 17 18 7 13 15 23 14 19 18 8 25 23 15 24 19 20 9 21 23 31 22 27 26 10 33 31 23 32 27 28 11 23 25 33 24 29 28 12 35 33 25 34 29 30 13 31 33 41 32 37 36 14 43 41 33 42 37 38 15 33 35 43 34 39 38 16 45 43 35 44 39 40 17 41 43 51 42 47 46 18 53 51 43 52 47 48 19 43 45 53 44 49 48 20 55 53 45 54 49 50 21 51 53 61 52 57 56 22 63 61 53 62 57 58 23 53 55 63 54 59 58 24 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 var = Columns 1 through 13 1 2 3 4 5 6 7 8 9 10 11 12 13 Columns 14 through 26 14 15 16 17 18 19 20 21 22 23 24 25 26 Columns 27 through 39 27 28 29 30 31 32 33 34 35 36 37 38 39 Columns 40 through 52 40 41 42 43 44 45 46 47 48 49 50 51 52 Columns 53 through 65 53 54 55 56 57 58 59 60 61 62 63 64 65 Columns 66 through 78 66 67 68 69 70 71 72 73 74 75 76 77 78 Columns 79 through 91 79 80 81 82 83 84 85 86 87 88 89 90 91 Columns 92 through 104 92 93 94 95 96 97 98 99 100 101 102 103 104 Columns 105 through 117 105 106 107 108 109 110 111 112 113 114 115 116 117 Columns 118 through 130 118 119 120 121 122 123 124 125 126 127 128 129 130 Columns 131 through 143 131 132 133 134 135 136 137 138 139 140 141 142 143 Columns 144 through 151 144 145 146 147 148 149 150 151 Node U_Var V_Var P_Var 1 1 2 3 2 4 5 3 6 7 8 4 9 10 5 11 12 13 6 14 15 7 16 17 8 18 19 9 20 21 10 22 23 11 24 25 26 12 27 28 13 29 30 31 14 32 33 15 34 35 36 16 37 38 17 39 40 18 41 42 19 43 44 20 45 46 21 47 48 49 22 50 51 23 52 53 54 24 55 56 25 57 58 59 26 60 61 27 62 63 28 64 65 29 66 67 30 68 69 31 70 71 72 32 73 74 33 75 76 77 34 78 79 35 80 81 82 36 83 84 37 85 86 38 87 88 39 89 90 40 91 92 41 93 94 95 42 96 97 43 98 99 100 44 101 102 45 103 104 105 46 106 107 47 108 109 48 110 111 49 112 113 50 114 115 51 116 117 118 52 119 120 53 121 122 123 54 124 125 55 126 127 128 56 129 130 57 131 132 58 133 134 59 135 136 60 137 138 61 139 140 141 62 142 143 63 144 145 146 64 147 148 65 149 150 151 Lower bandwidth ML = 25 Upper bandwidth MU = 25 Total bandwidth M = 51 fem2d_pack_test03 BASIS_11_T3_TEST tests BASIS_11_T3. BASIS_11_T4_TEST tests BASIS_11_T4. BASIS_11_T6_TEST tests BASIS_11_T6. BASIS_11_T3_TEST: Verify basis functions for element T3. Number of nodes = 3 Physical Nodes: 1 2.000000 0.000000 2 4.000000 3.000000 3 0.000000 4.000000 The basis function values at basis nodes should form the identity matrix. 1.000000 -0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 1.000000 The X and Y derivatives should sum to 0. dPhidX sum dPhidY sum 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 BASIS_11_T4_TEST: Verify basis functions for element T4. Number of nodes = 4 Physical Nodes: 1 2.000000 0.000000 2 4.000000 3.000000 3 0.000000 4.000000 4 2.000000 2.333333 The basis function values at basis nodes should form the identity matrix. 1.000000 0.000000 0.000000 -0.000000 0.000000 1.000000 0.000000 -0.000000 0.000000 0.000000 1.000000 -0.000000 0.000000 0.000000 0.000000 1.000000 The X and Y derivatives should sum to 0. dPhidX sum dPhidY sum 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 2.081668e-17 1.734723e-17 BASIS_11_T6_TEST: Verify basis functions for element T6. Number of nodes = 6 Physical Nodes: 1 2.000000 0.000000 2 4.000000 3.000000 3 0.000000 4.000000 4 3.000000 1.500000 5 2.000000 3.500000 6 1.000000 2.000000 The basis function values at basis nodes should form the identity matrix. 1.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 0.000000 1.000000 -0.000000 -0.000000 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 -0.000000 1.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 1.000000 -0.000000 -0.000000 -0.000000 -0.000000 -0.000000 0.000000 1.000000 The X and Y derivatives should sum to 0. dPhidX sum dPhidY sum 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 fem2d_pack_test04 Test the computation of basis functions at the nodes defining the element. BASIS_MN_Q4_TEST: the Q4 element BASIS_MN_T3_TEST: the T3 element BASIS_MN_T4_TEST: the T4 element BASIS_MN_T6_TEST: the T6 element BASIS_MN_Q4_TEST: Verify basis functions for element Q4. Number of nodes = 4 Physical Nodes: I X Y 1 3.000000 1.000000 2 5.000000 1.000000 3 5.000000 4.000000 4 3.000000 4.000000 The basis function values at basis nodes should form the identity matrix. 1.000000 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 1.000000 The X and Y derivatives should sum to 0. dPhidX sum dPhidY sum: 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 BASIS_MN_T3_TEST: Verify basis functions for element T3. Number of nodes = 3 Physical Nodes: I X Y 1 2.000000 0.000000 2 4.000000 3.000000 3 0.000000 4.000000 The basis function values at basis nodes should form the identity matrix. 1.000000 0.000000 0.000000 -0.000000 1.000000 0.000000 0.000000 0.000000 1.000000 The X and Y derivatives should sum to 0. dPhidX sum dPhidY sum 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 BASIS_MN_T4_TEST: Verify basis functions for element T4. Number of nodes = 4 Physical Nodes: 1 2.000000 0.000000 2 4.000000 2.000000 3 0.000000 4.000000 4 2.000000 2.000000 The basis function values at basis nodes should form the identity matrix. 1.000000 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 -0.000000 0.000000 1.000000 The X and Y derivatives should sum to 0. dPhidX sum dPhidY sum 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 BASIS_MN_T6_TEST Verify basis functions for element T6. Number of nodes = 6 Physical Nodes: I X Y 1 2.000000 0.000000 2 4.000000 3.000000 3 0.000000 4.000000 4 3.000000 1.500000 5 2.000000 3.500000 6 1.000000 2.000000 The basis function values at basis nodes should form the identity matrix. 1.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 1.000000 -0.000000 -0.000000 0.000000 -0.000000 -0.000000 -0.000000 1.000000 -0.000000 -0.000000 -0.000000 -0.000000 0.000000 -0.000000 1.000000 0.000000 -0.000000 -0.000000 -0.000000 0.000000 -0.000000 1.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 1.000000 The X and Y derivatives should sum to 0. dPhidX sum dPhidY sum 0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 fem2d_pack_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 1.396743 -1.010324 0.380122 -0.918725 -0.870980 -0.532032 0.846724 1.488428 -1.188873 0.617845 -0.859062 -0.613237 -0.330279 0.943883 1.902478 -0.994043 0.842763 -0.612216 -0.222061 -0.111964 0.993712 2.187751 -0.531485 0.981281 -0.222056 0.222061 0.111964 0.993712 2.261946 0.092185 0.981281 0.222056 0.613237 0.330279 0.943883 2.102002 0.691453 0.842763 0.612216 0.870980 0.532032 0.846724 1.758902 1.093205 0.617845 0.859062 0.974455 0.707107 0.707107 1.354679 1.264479 0.380122 0.918725 -1.655277 -0.989949 0.989949 1.348065 -1.358487 0.537757 -1.655192 -1.630388 -0.744845 1.185414 1.649421 -1.553495 1.033098 -1.607084 -1.225047 -0.462391 1.321437 2.322298 -1.197669 1.562824 -1.178957 -0.458247 -0.156750 1.391197 2.798939 -0.521349 1.911701 -0.434357 0.458247 0.156750 1.391197 2.870455 0.309438 1.911701 0.434357 1.225047 0.462391 1.321437 2.494395 1.057265 1.562824 1.178957 1.630388 0.744845 1.185414 1.835013 1.521308 1.033098 1.607084 1.655277 0.989949 0.989949 1.284467 1.808422 0.537757 1.655192 -2.504522 -1.272792 1.272792 1.234249 -2.070535 0.475653 -2.433386 -2.717614 -0.957658 1.524104 1.762507 -2.340530 1.336678 -2.492964 -2.176853 -0.594502 1.698990 2.980593 -1.855541 2.391311 -1.903197 -0.840610 -0.201536 1.788682 3.886956 -0.798321 3.134629 -0.716097 0.840610 0.201536 1.788682 4.006382 0.538901 3.134629 0.716097 2.176853 0.594502 1.698990 3.251105 1.705491 2.391311 1.903197 2.717614 0.957658 1.524104 2.016078 2.345812 1.336678 2.492964 2.504522 1.272792 1.272792 1.042843 2.629324 0.475653 2.433386 -3.419607 -1.555635 1.555635 0.767728 -2.758061 0.036689 -3.110912 -4.116574 -1.170471 1.862793 1.607342 -3.144764 1.352322 -3.431020 -3.529009 -0.726614 2.076543 3.610862 -2.600555 3.222934 -2.759072 -1.409222 -0.246322 2.186167 5.185593 -1.122357 4.635066 -1.066143 1.409222 0.246322 2.186167 5.372364 0.824541 4.635066 1.066143 3.529009 0.726614 2.076543 4.008042 2.460057 3.222934 2.759072 4.116574 1.170471 1.862793 1.924943 3.213730 1.352322 3.431020 3.419607 1.555635 1.555635 0.393346 3.401743 0.036689 3.110912 -4.224014 -1.838478 1.838478 -0.182544 -3.319021 -0.893997 -3.546007 -5.744043 -1.383283 2.201483 1.049884 -3.876050 0.903470 -4.325786 -5.316174 -0.858726 2.454097 4.106977 -3.393173 3.935177 -3.715554 -2.202902 -0.291108 2.583652 6.657109 -1.485672 6.394405 -1.483086 2.202902 0.291108 2.583652 6.933709 1.159886 6.394405 1.483086 5.316174 0.858726 2.454097 4.656530 3.282002 3.935177 3.715554 5.744043 1.383283 2.201483 1.413110 4.032305 0.903470 4.325786 4.224014 1.838478 1.838478 -0.781344 4.010428 -0.893997 3.546007 -4.687378 -2.121320 2.121320 -0.979517 -4.163355 -2.354103 -3.615793 -7.450092 -1.596096 2.540173 1.319853 -4.708767 -0.163230 -5.078719 -7.543624 -0.990837 2.831650 5.198180 -3.550941 4.393914 -4.737266 -3.258940 -0.335893 2.981137 8.079616 -0.704037 8.390526 -1.965242 3.258940 0.335893 2.981137 7.648660 2.298268 8.390526 1.965242 7.543624 0.990837 2.831650 4.059436 4.029915 4.393914 4.737266 7.450092 1.596096 2.540173 -0.001466 4.177629 -0.163230 5.078719 4.687378 2.121320 2.121320 -2.199541 3.837779 -2.354103 3.615793 fem2d_pack_test06 DIV_Q4 estimates divergence and vorticity using 4 node quadrilateral elements. Original U, V data forms 21 rows and 13 columns. I, J, X, Y, U(I,J), V(I,J) Some of X: Col: 5 6 7 8 Row 5 : 0.333333 0.416667 0.5 0.583333 6 : 0.333333 0.416667 0.5 0.583333 7 : 0.333333 0.416667 0.5 0.583333 8 : 0.333333 0.416667 0.5 0.583333 Some of Y: Col: 5 6 7 8 Row 5 : 0.2 0.2 0.2 0.2 6 : 0.25 0.25 0.25 0.25 7 : 0.3 0.3 0.3 0.3 8 : 0.35 0.35 0.35 0.35 Some of U: Col: 5 6 7 8 Row 5 : 0.0666667 0.0833333 0.1 0.116667 6 : 0.0833333 0.104167 0.125 0.145833 7 : 0.1 0.125 0.15 0.175 8 : 0.116667 0.145833 0.175 0.204167 Some of V: Col: 5 6 7 8 Row 5 : 0.150537 0.21199 0.285952 0.371178 6 : 0.17274 0.233923 0.307439 0.391975 7 : 0.199758 0.260569 0.333487 0.417123 8 : 0.231492 0.291803 0.363945 0.446435 Some of XM: Col: 5 6 7 8 Row 5 : 0.375 0.458333 0.541667 0.625 6 : 0.375 0.458333 0.541667 0.625 7 : 0.375 0.458333 0.541667 0.625 8 : 0.375 0.458333 0.541667 0.625 Some of YM: Col: 5 6 7 8 Row 5 : 0.225 0.225 0.225 0.225 6 : 0.275 0.275 0.275 0.275 7 : 0.325 0.325 0.325 0.325 8 : 0.375 0.375 0.375 0.375 I, J, DIV(I,J), Exact Divergence, Difference 5 5 0.375000 0.225000 0.666367 0.666795 -0.000428 5 6 0.458333 0.225000 0.659194 0.659795 -0.000601 5 7 0.541667 0.225000 0.647832 0.648632 -0.000800 5 8 0.625000 0.225000 0.630884 0.631898 -0.001015 6 5 0.375000 0.275000 0.811631 0.812190 -0.000559 6 6 0.458333 0.275000 0.801938 0.802706 -0.000769 6 7 0.541667 0.275000 0.786965 0.787973 -0.001008 6 8 0.625000 0.275000 0.765027 0.766293 -0.001266 7 5 0.375000 0.325000 0.954681 0.955392 -0.000711 7 6 0.458333 0.325000 0.941922 0.942878 -0.000956 7 7 0.541667 0.325000 0.922702 0.923935 -0.001234 7 8 0.625000 0.325000 0.895062 0.896593 -0.001532 8 5 0.375000 0.375000 1.094645 1.095532 -0.000887 8 6 0.458333 0.375000 1.078186 1.079351 -0.001165 8 7 0.541667 0.375000 1.053981 1.055459 -0.001479 8 8 0.625000 0.375000 1.019819 1.021631 -0.001812 I, J, VORT(I,J), Exact Vorticity, Difference 5 5 0.375000 0.225000 0.360821 0.361326 -0.000505 5 6 0.458333 0.225000 0.426529 0.427360 -0.000832 5 7 0.541667 0.225000 0.476911 0.478187 -0.001276 5 8 0.625000 0.225000 0.503428 0.505273 -0.001845 6 5 0.375000 0.275000 0.356961 0.357532 -0.000570 6 6 0.458333 0.275000 0.420268 0.421177 -0.000909 6 7 0.541667 0.275000 0.467371 0.468735 -0.001364 6 8 0.625000 0.275000 0.489636 0.491575 -0.001939 7 5 0.375000 0.325000 0.351727 0.352375 -0.000648 7 6 0.458333 0.325000 0.412032 0.413033 -0.001001 7 7 0.541667 0.325000 0.455092 0.456559 -0.001467 7 8 0.625000 0.325000 0.472168 0.474218 -0.002049 8 5 0.375000 0.375000 0.344795 0.345532 -0.000737 8 6 0.458333 0.375000 0.401434 0.402540 -0.001106 8 7 0.541667 0.375000 0.439635 0.441219 -0.001584 8 8 0.625000 0.375000 0.450544 0.452718 -0.002174 fem2d_pack_test07 ELEMENTS_EPS creates an Encapsulated PostScript file containing an image of a Q4 mesh. fem2d_pack_test08 ELEMENTS_EPS creates an Encapsulated PostScript file containing an image of a T3 mesh. fem2d_pack_test09 ELEMENTS_EPS creates an Encapsulated PostScript file containing an image of a T4 mesh. fem2d_pack_test10 ELEMENTS_EPS creates an Encapsulated PostScript file containing an image of a T6 mesh. fem2d_pack_test105 GRID_NODES_01 creates a regular grid in the unit square. NUM_X = 5 NUM_Y = 3 NODE_NUM = 15 1 0.000000 0.000000 2 0.250000 0.000000 3 0.500000 0.000000 4 0.750000 0.000000 5 1.000000 0.000000 6 0.000000 0.500000 7 0.250000 0.500000 8 0.500000 0.500000 9 0.750000 0.500000 10 1.000000 0.500000 11 0.000000 1.000000 12 0.250000 1.000000 13 0.500000 1.000000 14 0.750000 1.000000 15 1.000000 1.000000 fem2d_pack_test11 Test the grid routines. GRID_TEST(): Test the grid routine for element "Q4" Number of elements in X direction = 3 Number of elements in Y direction = 2 Element order = 4 Nodes in grid = 12 GRID_PRINT: Element -> Node table. Element order = 4 Number of elements = 6 # 1 2 3 4 1 1 2 6 5 2 2 3 7 6 3 3 4 8 7 4 5 6 10 9 5 6 7 11 10 6 7 8 12 11 The grid width is 5 GRID_TEST(): Test the grid routine for element "Q8" Number of elements in X direction = 3 Number of elements in Y direction = 2 Element order = 8 Nodes in grid = 29 GRID_PRINT: Element -> Node table. Element order = 8 Number of elements = 6 # 1 2 3 4 5 6 7 8 1 1 3 14 12 2 9 13 8 2 3 5 16 14 4 10 15 9 3 5 7 18 16 6 11 17 10 4 12 14 25 23 13 20 24 19 5 14 16 27 25 15 21 26 20 6 16 18 29 27 17 22 28 21 The grid width is 13 GRID_TEST(): Test the grid routine for element "Q9" Number of elements in X direction = 3 Number of elements in Y direction = 2 Element order = 9 Nodes in grid = 35 GRID_PRINT: Element -> Node table. Element order = 9 Number of elements = 6 # 1 2 3 4 5 6 7 8 9 1 1 3 17 15 2 10 16 8 9 2 3 5 19 17 4 12 18 10 11 3 5 7 21 19 6 14 20 12 13 4 15 17 31 29 16 24 30 22 23 5 17 19 33 31 18 26 32 24 25 6 19 21 35 33 20 28 34 26 27 The grid width is 16 GRID_TEST(): Test the grid routine for element "Q12" Number of elements in X direction = 3 Number of elements in Y direction = 2 Element order = 12 Nodes in grid = 46 GRID_PRINT: Element -> Node table. Element order = 12 Number of elements = 6 # 1 2 3 4 5 6 7 8 9 10 11 12 1 1 2 3 4 11 12 15 16 19 20 21 22 2 4 5 6 7 12 13 16 17 22 23 24 25 3 7 8 9 10 13 14 17 18 25 26 27 28 4 19 20 21 22 29 30 33 34 37 38 39 40 5 22 23 24 25 30 31 34 35 40 41 42 43 6 25 26 27 28 31 32 35 36 43 44 45 46 The grid width is 21 GRID_TEST(): Test the grid routine for element "Q16" Number of elements in X direction = 3 Number of elements in Y direction = 2 Element order = 16 Nodes in grid = 70 GRID_PRINT: Element -> Node table. Element order = 16 Number of elements = 6 # 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 1 1 2 3 4 11 12 13 14 21 22 23 24 31 32 33 34 2 4 5 6 7 14 15 16 17 24 25 26 27 34 35 36 37 3 7 8 9 10 17 18 19 20 27 28 29 30 37 38 39 40 4 31 32 33 34 41 42 43 44 51 52 53 54 61 62 63 64 5 34 35 36 37 44 45 46 47 54 55 56 57 64 65 66 67 6 37 38 39 40 47 48 49 50 57 58 59 60 67 68 69 70 The grid width is 33 GRID_TEST(): Test the grid routine for element "QL" Number of elements in X direction = 3 Number of elements in Y direction = 2 Element order = 6 Nodes in grid = 21 GRID_PRINT: Element -> Node table. Element order = 6 Number of elements = 6 # 1 2 3 4 5 6 1 1 2 3 8 9 10 2 3 4 5 10 11 12 3 5 6 7 12 13 14 4 8 9 10 15 16 17 5 10 11 12 17 18 19 6 12 13 14 19 20 21 The grid width is 9 GRID_TEST(): Test the grid routine for element "T3" Number of elements in X direction = 3 Number of elements in Y direction = 2 Element order = 3 Nodes in grid = 12 GRID_PRINT: Element -> Node table. Element order = 3 Number of elements = 12 # 1 2 3 1 1 2 5 2 6 5 2 3 2 3 6 4 7 6 3 5 3 4 7 6 8 7 4 7 5 6 9 8 10 9 6 9 6 7 10 10 11 10 7 11 7 8 11 12 12 11 8 The grid width is 4 GRID_TEST(): Test the grid routine for element "T4" Number of elements in X direction = 3 Number of elements in Y direction = 2 Element order = -1 Nodes in grid = 24 GRID_PRINT: Element -> Node table. Element order = -1 Number of elements = 12 # 1 2 3 4 5 6 7 8 9 10 11 12 The grid width is 0 GRID_TEST(): Test the grid routine for element "T6" Number of elements in X direction = 3 Number of elements in Y direction = 2 Element order = 6 Nodes in grid = 35 GRID_PRINT: Element -> Node table. Element order = 6 Number of elements = 12 # 1 2 3 4 5 6 1 1 3 15 2 9 8 2 17 15 3 16 9 10 3 3 5 17 4 11 10 4 19 17 5 18 11 12 5 5 7 19 6 13 12 6 21 19 7 20 13 14 7 15 17 29 16 23 22 8 31 29 17 30 23 24 9 17 19 31 18 25 24 10 33 31 19 32 25 26 11 19 21 33 20 27 26 12 35 33 21 34 27 28 The grid width is 14 GRID_TEST(): Test the grid routine for element "T10" Number of elements in X direction = 3 Number of elements in Y direction = 2 Element order = 10 Nodes in grid = 70 GRID_PRINT: Element -> Node table. Element order = 10 Number of elements = 12 # 1 2 3 4 5 6 7 8 9 10 1 1 2 3 4 13 22 31 21 10 11 2 34 33 32 31 22 13 4 14 24 23 3 4 5 6 7 16 25 34 24 13 14 4 37 36 35 34 25 16 7 17 27 26 5 7 8 9 10 19 28 37 27 16 17 6 40 39 38 37 28 19 10 20 30 29 7 31 32 33 34 43 52 61 51 40 41 8 64 63 62 61 52 43 34 44 54 53 9 34 35 36 37 46 55 64 54 43 44 10 67 66 65 64 55 46 37 47 57 56 11 37 38 39 40 49 58 67 57 46 47 12 70 69 68 67 58 49 40 50 60 59 The grid width is 30 fem2d_pack_test12 INTERP_TEST tests the interpolating power of the element. INTERP_TEST for element "Q4". Element order = 4 Interpolate R^0 * S^0 (R,S,U): 0.000000 0.000000 1.000000 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 1.000000 (R,S): 0.814724 0.905792 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.126987 0.913376 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.632359 0.097540 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.278498 0.546882 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.957507 0.964889 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 Interpolate R^0 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 1.000000 (R,S): 0.157613 0.970593 (U,U*,Error): 0.970593 0.970593 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.957167 0.485376 (U,U*,Error): 0.485376 0.485376 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.800280 0.141886 (U,U*,Error): 0.141886 0.141886 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.421761 0.915736 (U,U*,Error): 0.915736 0.915736 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.792207 0.959492 (U,U*,Error): 0.959492 0.959492 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 Interpolate R^1 * S^0 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S): 0.655741 0.035712 (U,U*,Error): 0.655741 0.655741 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.849129 0.933993 (U,U*,Error): 0.849129 0.849129 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.678735 0.757740 (U,U*,Error): 0.678735 0.678735 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.743132 0.392227 (U,U*,Error): 0.743132 0.743132 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.655478 0.171187 (U,U*,Error): 0.655478 0.655478 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 Interpolate R^1 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S): 0.706046 0.031833 (U,U*,Error): 0.022475 0.022475 0.000000 (Ur,Ur*,Error): 0.031833 0.031833 0.000000 (Us,Us*,Error): 0.706046 0.706046 0.000000 (R,S): 0.276923 0.046171 (U,U*,Error): 0.012786 0.012786 0.000000 (Ur,Ur*,Error): 0.046171 0.046171 0.000000 (Us,Us*,Error): 0.276923 0.276923 0.000000 (R,S): 0.097132 0.823458 (U,U*,Error): 0.079984 0.079984 0.000000 (Ur,Ur*,Error): 0.823458 0.823458 0.000000 (Us,Us*,Error): 0.097132 0.097132 0.000000 (R,S): 0.694829 0.317099 (U,U*,Error): 0.220330 0.220330 0.000000 (Ur,Ur*,Error): 0.317099 0.317099 0.000000 (Us,Us*,Error): 0.694829 0.694829 0.000000 (R,S): 0.950222 0.034446 (U,U*,Error): 0.032731 0.032731 0.000000 (Ur,Ur*,Error): 0.034446 0.034446 0.000000 (Us,Us*,Error): 0.950222 0.950222 0.000000 INTERP_TEST for element "Q8". Element order = 8 Interpolate R^0 * S^0 (R,S,U): 0.000000 0.000000 1.000000 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 1.000000 (R,S,U): 0.500000 0.000000 1.000000 (R,S,U): 1.000000 0.500000 1.000000 (R,S,U): 0.500000 1.000000 1.000000 (R,S,U): 0.000000 0.500000 1.000000 (R,S): 0.438744 0.381558 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.765517 0.795200 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.186873 0.489764 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.445586 0.646313 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.709365 0.754687 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 Interpolate R^0 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 1.000000 (R,S,U): 0.500000 0.000000 0.000000 (R,S,U): 1.000000 0.500000 0.500000 (R,S,U): 0.500000 1.000000 1.000000 (R,S,U): 0.000000 0.500000 0.500000 (R,S): 0.276025 0.679703 (U,U*,Error): 0.679703 0.679703 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.655098 0.162612 (U,U*,Error): 0.162612 0.162612 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.118998 0.498364 (U,U*,Error): 0.498364 0.498364 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.959744 0.340386 (U,U*,Error): 0.340386 0.340386 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.585268 0.223812 (U,U*,Error): 0.223812 0.223812 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 Interpolate R^1 * S^0 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.500000 0.000000 0.500000 (R,S,U): 1.000000 0.500000 1.000000 (R,S,U): 0.500000 1.000000 0.500000 (R,S,U): 0.000000 0.500000 0.000000 (R,S): 0.751267 0.255095 (U,U*,Error): 0.751267 0.751267 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.505957 0.699077 (U,U*,Error): 0.505957 0.505957 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.890903 0.959291 (U,U*,Error): 0.890903 0.890903 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.547216 0.138624 (U,U*,Error): 0.547216 0.547216 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.149294 0.257508 (U,U*,Error): 0.149294 0.149294 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 Interpolate R^0 * S^2 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 1.000000 (R,S,U): 0.500000 0.000000 0.000000 (R,S,U): 1.000000 0.500000 0.250000 (R,S,U): 0.500000 1.000000 1.000000 (R,S,U): 0.000000 0.500000 0.250000 (R,S): 0.840717 0.254282 (U,U*,Error): 0.064659 0.064659 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.508564 0.508564 0.000000 (R,S): 0.814285 0.243525 (U,U*,Error): 0.059304 0.059304 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.487050 0.487050 0.000000 (R,S): 0.929264 0.349984 (U,U*,Error): 0.122489 0.122489 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.699968 0.699968 0.000000 (R,S): 0.196595 0.251084 (U,U*,Error): 0.063043 0.063043 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.502168 0.502168 0.000000 (R,S): 0.616045 0.473289 (U,U*,Error): 0.224002 0.224002 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 0.946578 0.946578 0.000000 Interpolate R^1 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.500000 0.000000 0.000000 (R,S,U): 1.000000 0.500000 0.500000 (R,S,U): 0.500000 1.000000 0.500000 (R,S,U): 0.000000 0.500000 0.000000 (R,S): 0.351660 0.830829 (U,U*,Error): 0.292169 0.292169 0.000000 (Ur,Ur*,Error): 0.830829 0.830829 0.000000 (Us,Us*,Error): 0.351660 0.351660 0.000000 (R,S): 0.585264 0.549724 (U,U*,Error): 0.321733 0.321733 0.000000 (Ur,Ur*,Error): 0.549724 0.549724 0.000000 (Us,Us*,Error): 0.585264 0.585264 0.000000 (R,S): 0.917194 0.285839 (U,U*,Error): 0.262170 0.262170 0.000000 (Ur,Ur*,Error): 0.285839 0.285839 0.000000 (Us,Us*,Error): 0.917194 0.917194 0.000000 (R,S): 0.757200 0.753729 (U,U*,Error): 0.570724 0.570724 0.000000 (Ur,Ur*,Error): 0.753729 0.753729 0.000000 (Us,Us*,Error): 0.757200 0.757200 0.000000 (R,S): 0.380446 0.567822 (U,U*,Error): 0.216025 0.216025 0.000000 (Ur,Ur*,Error): 0.567822 0.567822 0.000000 (Us,Us*,Error): 0.380446 0.380446 0.000000 Interpolate R^2 * S^0 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.500000 0.000000 0.250000 (R,S,U): 1.000000 0.500000 1.000000 (R,S,U): 0.500000 1.000000 0.250000 (R,S,U): 0.000000 0.500000 0.000000 (R,S): 0.075854 0.053950 (U,U*,Error): 0.005754 0.005754 0.000000 (Ur,Ur*,Error): 0.151709 0.151709 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.530798 0.779167 (U,U*,Error): 0.281746 0.281746 0.000000 (Ur,Ur*,Error): 1.061595 1.061595 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.934011 0.129906 (U,U*,Error): 0.872376 0.872376 0.000000 (Ur,Ur*,Error): 1.868021 1.868021 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.568824 0.469391 (U,U*,Error): 0.323560 0.323560 0.000000 (Ur,Ur*,Error): 1.137647 1.137647 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.011902 0.337123 (U,U*,Error): 0.000142 0.000142 0.000000 (Ur,Ur*,Error): 0.023804 0.023804 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 Interpolate R^1 * S^2 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.500000 0.000000 0.000000 (R,S,U): 1.000000 0.500000 0.250000 (R,S,U): 0.500000 1.000000 0.500000 (R,S,U): 0.000000 0.500000 0.000000 (R,S): 0.162182 0.794285 (U,U*,Error): 0.102319 0.102319 0.000000 (Ur,Ur*,Error): 0.630888 0.630888 0.000000 (Us,Us*,Error): 0.257638 0.257638 0.000000 (R,S): 0.311215 0.528533 (U,U*,Error): 0.086937 0.086937 0.000000 (Ur,Ur*,Error): 0.279347 0.279347 0.000000 (Us,Us*,Error): 0.328975 0.328975 0.000000 (R,S): 0.165649 0.601982 (U,U*,Error): 0.060028 0.060028 0.000000 (Ur,Ur*,Error): 0.362382 0.362382 0.000000 (Us,Us*,Error): 0.199435 0.199435 0.000000 (R,S): 0.262971 0.654079 (U,U*,Error): 0.112504 0.112504 0.000000 (Ur,Ur*,Error): 0.427819 0.427819 0.000000 (Us,Us*,Error): 0.344008 0.344008 0.000000 (R,S): 0.689215 0.748152 (U,U*,Error): 0.385775 0.385775 0.000000 (Ur,Ur*,Error): 0.559731 0.559731 0.000000 (Us,Us*,Error): 1.031274 1.031274 0.000000 Interpolate R^2 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.500000 0.000000 0.000000 (R,S,U): 1.000000 0.500000 0.500000 (R,S,U): 0.500000 1.000000 0.250000 (R,S,U): 0.000000 0.500000 0.000000 (R,S): 0.450542 0.083821 (U,U*,Error): 0.017015 0.017015 0.000000 (Ur,Ur*,Error): 0.075530 0.075530 0.000000 (Us,Us*,Error): 0.202988 0.202988 0.000000 (R,S): 0.228977 0.913337 (U,U*,Error): 0.047887 0.047887 0.000000 (Ur,Ur*,Error): 0.418266 0.418266 0.000000 (Us,Us*,Error): 0.052430 0.052430 0.000000 (R,S): 0.152378 0.825817 (U,U*,Error): 0.019175 0.019175 0.000000 (Ur,Ur*,Error): 0.251673 0.251673 0.000000 (Us,Us*,Error): 0.023219 0.023219 0.000000 (R,S): 0.538342 0.996135 (U,U*,Error): 0.288692 0.288692 0.000000 (Ur,Ur*,Error): 1.072523 1.072523 0.000000 (Us,Us*,Error): 0.289813 0.289813 0.000000 (R,S): 0.078176 0.442678 (U,U*,Error): 0.002705 0.002705 0.000000 (Ur,Ur*,Error): 0.069213 0.069213 0.000000 (Us,Us*,Error): 0.006111 0.006111 0.000000 INTERP_TEST for element "Q9". Element order = 9 Interpolate R^0 * S^0 (R,S,U): 0.000000 0.000000 1.000000 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 1.000000 (R,S,U): 0.500000 0.000000 1.000000 (R,S,U): 1.000000 0.500000 1.000000 (R,S,U): 0.500000 1.000000 1.000000 (R,S,U): 0.000000 0.500000 1.000000 (R,S,U): 0.500000 0.500000 1.000000 (R,S): 0.106653 0.961898 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.004634 0.774910 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.817303 0.868695 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.084436 0.399783 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.259870 0.800068 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 Interpolate R^0 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 1.000000 (R,S,U): 0.500000 0.000000 0.000000 (R,S,U): 1.000000 0.500000 0.500000 (R,S,U): 0.500000 1.000000 1.000000 (R,S,U): 0.000000 0.500000 0.500000 (R,S,U): 0.500000 0.500000 0.500000 (R,S): 0.431414 0.910648 (U,U*,Error): 0.910648 0.910648 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.181847 0.263803 (U,U*,Error): 0.263803 0.263803 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.145539 0.136069 (U,U*,Error): 0.136069 0.136069 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.869292 0.579705 (U,U*,Error): 0.579705 0.579705 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.549860 0.144955 (U,U*,Error): 0.144955 0.144955 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 Interpolate R^1 * S^0 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.500000 0.000000 0.500000 (R,S,U): 1.000000 0.500000 1.000000 (R,S,U): 0.500000 1.000000 0.500000 (R,S,U): 0.000000 0.500000 0.000000 (R,S,U): 0.500000 0.500000 0.500000 (R,S): 0.853031 0.622055 (U,U*,Error): 0.853031 0.853031 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.350952 0.513250 (U,U*,Error): 0.350952 0.350952 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.401808 0.075967 (U,U*,Error): 0.401808 0.401808 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.239916 0.123319 (U,U*,Error): 0.239916 0.239916 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.183908 0.239953 (U,U*,Error): 0.183908 0.183908 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 Interpolate R^0 * S^2 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 1.000000 (R,S,U): 0.500000 0.000000 0.000000 (R,S,U): 1.000000 0.500000 0.250000 (R,S,U): 0.500000 1.000000 1.000000 (R,S,U): 0.000000 0.500000 0.250000 (R,S,U): 0.500000 0.500000 0.250000 (R,S): 0.417267 0.049654 (U,U*,Error): 0.002466 0.002466 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.099309 0.099309 0.000000 (R,S): 0.902716 0.944787 (U,U*,Error): 0.892623 0.892623 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 1.889574 1.889574 0.000000 (R,S): 0.490864 0.489253 (U,U*,Error): 0.239368 0.239368 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.978505 0.978505 0.000000 (R,S): 0.337719 0.900054 (U,U*,Error): 0.810097 0.810097 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.800108 1.800108 0.000000 (R,S): 0.369247 0.111203 (U,U*,Error): 0.012366 0.012366 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.222406 0.222406 0.000000 Interpolate R^1 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.500000 0.000000 0.000000 (R,S,U): 1.000000 0.500000 0.500000 (R,S,U): 0.500000 1.000000 0.500000 (R,S,U): 0.000000 0.500000 0.000000 (R,S,U): 0.500000 0.500000 0.250000 (R,S): 0.780252 0.389739 (U,U*,Error): 0.304095 0.304095 0.000000 (Ur,Ur*,Error): 0.389739 0.389739 0.000000 (Us,Us*,Error): 0.780252 0.780252 0.000000 (R,S): 0.241691 0.403912 (U,U*,Error): 0.097622 0.097622 0.000000 (Ur,Ur*,Error): 0.403912 0.403912 0.000000 (Us,Us*,Error): 0.241691 0.241691 0.000000 (R,S): 0.096455 0.131973 (U,U*,Error): 0.012729 0.012729 0.000000 (Ur,Ur*,Error): 0.131973 0.131973 0.000000 (Us,Us*,Error): 0.096455 0.096455 0.000000 (R,S): 0.942051 0.956135 (U,U*,Error): 0.900727 0.900727 0.000000 (Ur,Ur*,Error): 0.956135 0.956135 0.000000 (Us,Us*,Error): 0.942051 0.942051 0.000000 (R,S): 0.575209 0.059780 (U,U*,Error): 0.034386 0.034386 0.000000 (Ur,Ur*,Error): 0.059780 0.059780 0.000000 (Us,Us*,Error): 0.575209 0.575209 0.000000 Interpolate R^2 * S^0 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.500000 0.000000 0.250000 (R,S,U): 1.000000 0.500000 1.000000 (R,S,U): 0.500000 1.000000 0.250000 (R,S,U): 0.000000 0.500000 0.000000 (R,S,U): 0.500000 0.500000 0.250000 (R,S): 0.234780 0.353159 (U,U*,Error): 0.055122 0.055122 0.000000 (Ur,Ur*,Error): 0.469560 0.469560 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.821194 0.015403 (U,U*,Error): 0.674360 0.674360 0.000000 (Ur,Ur*,Error): 1.642388 1.642388 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.043024 0.168990 (U,U*,Error): 0.001851 0.001851 0.000000 (Ur,Ur*,Error): 0.086048 0.086048 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.649115 0.731722 (U,U*,Error): 0.421351 0.421351 0.000000 (Ur,Ur*,Error): 1.298231 1.298231 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.647746 0.450924 (U,U*,Error): 0.419575 0.419575 0.000000 (Ur,Ur*,Error): 1.295492 1.295492 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 Interpolate R^1 * S^2 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.500000 0.000000 0.000000 (R,S,U): 1.000000 0.500000 0.250000 (R,S,U): 0.500000 1.000000 0.500000 (R,S,U): 0.000000 0.500000 0.000000 (R,S,U): 0.500000 0.500000 0.125000 (R,S): 0.547009 0.296321 (U,U*,Error): 0.048031 0.048031 0.000000 (Ur,Ur*,Error): 0.087806 0.087806 0.000000 (Us,Us*,Error): 0.324180 0.324180 0.000000 (R,S): 0.744693 0.188955 (U,U*,Error): 0.026589 0.026589 0.000000 (Ur,Ur*,Error): 0.035704 0.035704 0.000000 (Us,Us*,Error): 0.281427 0.281427 0.000000 (R,S): 0.686775 0.183511 (U,U*,Error): 0.023128 0.023128 0.000000 (Ur,Ur*,Error): 0.033676 0.033676 0.000000 (Us,Us*,Error): 0.252062 0.252062 0.000000 (R,S): 0.368485 0.625619 (U,U*,Error): 0.144224 0.144224 0.000000 (Ur,Ur*,Error): 0.391399 0.391399 0.000000 (Us,Us*,Error): 0.461062 0.461062 0.000000 (R,S): 0.780227 0.081126 (U,U*,Error): 0.005135 0.005135 0.000000 (Ur,Ur*,Error): 0.006581 0.006581 0.000000 (Us,Us*,Error): 0.126593 0.126593 0.000000 Interpolate R^2 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.500000 0.000000 0.000000 (R,S,U): 1.000000 0.500000 0.500000 (R,S,U): 0.500000 1.000000 0.250000 (R,S,U): 0.000000 0.500000 0.000000 (R,S,U): 0.500000 0.500000 0.125000 (R,S): 0.929386 0.775713 (U,U*,Error): 0.670028 0.670028 0.000000 (Ur,Ur*,Error): 1.441873 1.441873 0.000000 (Us,Us*,Error): 0.863758 0.863758 0.000000 (R,S): 0.486792 0.435859 (U,U*,Error): 0.103284 0.103284 0.000000 (Ur,Ur*,Error): 0.424345 0.424345 0.000000 (Us,Us*,Error): 0.236966 0.236966 0.000000 (R,S): 0.446784 0.306349 (U,U*,Error): 0.061152 0.061152 0.000000 (Ur,Ur*,Error): 0.273744 0.273744 0.000000 (Us,Us*,Error): 0.199616 0.199616 0.000000 (R,S): 0.508509 0.510772 (U,U*,Error): 0.132076 0.132076 0.000000 (Ur,Ur*,Error): 0.519464 0.519464 0.000000 (Us,Us*,Error): 0.258581 0.258581 0.000000 (R,S): 0.817628 0.794831 (U,U*,Error): 0.531357 0.531357 0.000000 (Ur,Ur*,Error): 1.299752 1.299752 0.000000 (Us,Us*,Error): 0.668515 0.668515 0.000000 Interpolate R^2 * S^2 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.500000 0.000000 0.000000 (R,S,U): 1.000000 0.500000 0.250000 (R,S,U): 0.500000 1.000000 0.250000 (R,S,U): 0.000000 0.500000 0.000000 (R,S,U): 0.500000 0.500000 0.062500 (R,S): 0.644318 0.378609 (U,U*,Error): 0.059509 0.059509 0.000000 (Ur,Ur*,Error): 0.184720 0.184720 0.000000 (Us,Us*,Error): 0.314356 0.314356 0.000000 (R,S): 0.811580 0.532826 (U,U*,Error): 0.186996 0.186996 0.000000 (Ur,Ur*,Error): 0.460820 0.460820 0.000000 (Us,Us*,Error): 0.701905 0.701905 0.000000 (R,S): 0.350727 0.939002 (U,U*,Error): 0.108460 0.108460 0.000000 (Ur,Ur*,Error): 0.618489 0.618489 0.000000 (Us,Us*,Error): 0.231012 0.231012 0.000000 (R,S): 0.875943 0.550156 (U,U*,Error): 0.232233 0.232233 0.000000 (Ur,Ur*,Error): 0.530247 0.530247 0.000000 (Us,Us*,Error): 0.844243 0.844243 0.000000 (R,S): 0.622475 0.587045 (U,U*,Error): 0.133532 0.133532 0.000000 (Ur,Ur*,Error): 0.429037 0.429037 0.000000 (Us,Us*,Error): 0.454931 0.454931 0.000000 INTERP_TEST for element "Q12". Element order = 12 Interpolate R^0 * S^0 (R,S,U): 0.000000 0.000000 1.000000 (R,S,U): 0.333333 0.000000 1.000000 (R,S,U): 0.666667 0.000000 1.000000 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 0.000000 0.333333 1.000000 (R,S,U): 1.000000 0.333333 1.000000 (R,S,U): 0.000000 0.666667 1.000000 (R,S,U): 1.000000 0.666667 1.000000 (R,S,U): 0.000000 1.000000 1.000000 (R,S,U): 0.333333 1.000000 1.000000 (R,S,U): 0.666667 1.000000 1.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.207742 0.301246 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.470923 0.230488 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.844309 0.194764 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.225922 0.170708 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.227664 0.435699 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 Interpolate R^0 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 0.333333 0.333333 (R,S,U): 1.000000 0.333333 0.333333 (R,S,U): 0.000000 0.666667 0.666667 (R,S,U): 1.000000 0.666667 0.666667 (R,S,U): 0.000000 1.000000 1.000000 (R,S,U): 0.333333 1.000000 1.000000 (R,S,U): 0.666667 1.000000 1.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.311102 0.923380 (U,U*,Error): 0.923380 0.923380 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.430207 0.184816 (U,U*,Error): 0.184816 0.184816 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.904881 0.979748 (U,U*,Error): 0.979748 0.979748 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.438870 0.111119 (U,U*,Error): 0.111119 0.111119 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.258065 0.408720 (U,U*,Error): 0.408720 0.408720 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 Interpolate R^1 * S^0 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.333333 (R,S,U): 0.666667 0.000000 0.666667 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 1.000000 0.333333 1.000000 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 1.000000 0.666667 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.333333 1.000000 0.333333 (R,S,U): 0.666667 1.000000 0.666667 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.594896 0.262212 (U,U*,Error): 0.594896 0.594896 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.602843 0.711216 (U,U*,Error): 0.602843 0.602843 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.221747 0.117418 (U,U*,Error): 0.221747 0.221747 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.296676 0.318778 (U,U*,Error): 0.296676 0.296676 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.424167 0.507858 (U,U*,Error): 0.424167 0.424167 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 Interpolate R^0 * S^2 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 0.333333 0.111111 (R,S,U): 1.000000 0.333333 0.111111 (R,S,U): 0.000000 0.666667 0.444444 (R,S,U): 1.000000 0.666667 0.444444 (R,S,U): 0.000000 1.000000 1.000000 (R,S,U): 0.333333 1.000000 1.000000 (R,S,U): 0.666667 1.000000 1.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.085516 0.262482 (U,U*,Error): 0.068897 0.068897 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.524964 0.524964 0.000000 (R,S): 0.801015 0.029220 (U,U*,Error): 0.000854 0.000854 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.058441 0.058441 0.000000 (R,S): 0.928854 0.730331 (U,U*,Error): 0.533383 0.533383 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.460662 1.460662 0.000000 (R,S): 0.488609 0.578525 (U,U*,Error): 0.334691 0.334691 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.157050 1.157050 0.000000 (R,S): 0.237284 0.458849 (U,U*,Error): 0.210542 0.210542 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.917698 0.917698 0.000000 Interpolate R^1 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 1.000000 0.333333 0.333333 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 1.000000 0.666667 0.666667 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.333333 1.000000 0.333333 (R,S,U): 0.666667 1.000000 0.666667 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.963089 0.546806 (U,U*,Error): 0.526622 0.526622 0.000000 (Ur,Ur*,Error): 0.546806 0.546806 0.000000 (Us,Us*,Error): 0.963089 0.963089 0.000000 (R,S): 0.521136 0.231594 (U,U*,Error): 0.120692 0.120692 0.000000 (Ur,Ur*,Error): 0.231594 0.231594 0.000000 (Us,Us*,Error): 0.521136 0.521136 0.000000 (R,S): 0.488898 0.624060 (U,U*,Error): 0.305102 0.305102 0.000000 (Ur,Ur*,Error): 0.624060 0.624060 0.000000 (Us,Us*,Error): 0.488898 0.488898 0.000000 (R,S): 0.679136 0.395515 (U,U*,Error): 0.268608 0.268608 0.000000 (Ur,Ur*,Error): 0.395515 0.395515 0.000000 (Us,Us*,Error): 0.679136 0.679136 0.000000 (R,S): 0.367437 0.987982 (U,U*,Error): 0.363021 0.363021 0.000000 (Ur,Ur*,Error): 0.987982 0.987982 0.000000 (Us,Us*,Error): 0.367437 0.367437 0.000000 Interpolate R^2 * S^0 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.111111 (R,S,U): 0.666667 0.000000 0.444444 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 1.000000 0.333333 1.000000 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 1.000000 0.666667 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.333333 1.000000 0.111111 (R,S,U): 0.666667 1.000000 0.444444 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.037739 0.885168 (U,U*,Error): 0.001424 0.001424 0.000000 (Ur,Ur*,Error): 0.075478 0.075478 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.913287 0.796184 (U,U*,Error): 0.834093 0.834093 0.000000 (Ur,Ur*,Error): 1.826574 1.826574 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.098712 0.261871 (U,U*,Error): 0.009744 0.009744 0.000000 (Ur,Ur*,Error): 0.197425 0.197425 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.335357 0.679728 (U,U*,Error): 0.112464 0.112464 0.000000 (Ur,Ur*,Error): 0.670714 0.670714 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.136553 0.721227 (U,U*,Error): 0.018647 0.018647 0.000000 (Ur,Ur*,Error): 0.273106 0.273106 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 Interpolate R^0 * S^3 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 0.333333 0.037037 (R,S,U): 1.000000 0.333333 0.037037 (R,S,U): 0.000000 0.666667 0.296296 (R,S,U): 1.000000 0.666667 0.296296 (R,S,U): 0.000000 1.000000 1.000000 (R,S,U): 0.333333 1.000000 1.000000 (R,S,U): 0.666667 1.000000 1.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.106762 0.653757 (U,U*,Error): 0.279415 0.279415 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.282196 1.282196 0.000000 (R,S): 0.494174 0.779052 (U,U*,Error): 0.472823 0.472823 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 1.820765 1.820765 0.000000 (R,S): 0.715037 0.903721 (U,U*,Error): 0.738078 0.738078 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 2.450133 2.450133 0.000000 (R,S): 0.890923 0.334163 (U,U*,Error): 0.037314 0.037314 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 0.334995 0.334995 0.000000 (R,S): 0.698746 0.197810 (U,U*,Error): 0.007740 0.007740 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.117386 0.117386 0.000000 Interpolate R^1 * S^2 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 1.000000 0.333333 0.111111 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 1.000000 0.666667 0.444444 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.333333 1.000000 0.333333 (R,S,U): 0.666667 1.000000 0.666667 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.030541 0.744074 (U,U*,Error): 0.016909 0.016909 0.000000 (Ur,Ur*,Error): 0.553647 0.553647 0.000000 (Us,Us*,Error): 0.045449 0.045449 0.000000 (R,S): 0.500022 0.479922 (U,U*,Error): 0.115168 0.115168 0.000000 (Ur,Ur*,Error): 0.230325 0.230325 0.000000 (Us,Us*,Error): 0.479944 0.479944 0.000000 (R,S): 0.904722 0.609867 (U,U*,Error): 0.336500 0.336500 0.000000 (Ur,Ur*,Error): 0.371937 0.371937 0.000000 (Us,Us*,Error): 1.103520 1.103520 0.000000 (R,S): 0.617666 0.859442 (U,U*,Error): 0.456234 0.456234 0.000000 (Ur,Ur*,Error): 0.738641 0.738641 0.000000 (Us,Us*,Error): 1.061697 1.061697 0.000000 (R,S): 0.805489 0.576722 (U,U*,Error): 0.267912 0.267912 0.000000 (Ur,Ur*,Error): 0.332608 0.332608 0.000000 (Us,Us*,Error): 0.929086 0.929086 0.000000 Interpolate R^2 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 1.000000 0.333333 0.333333 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 1.000000 0.666667 0.666667 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.333333 1.000000 0.111111 (R,S,U): 0.666667 1.000000 0.444444 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.182922 0.239932 (U,U*,Error): 0.008028 0.008028 0.000000 (Ur,Ur*,Error): 0.087778 0.087778 0.000000 (Us,Us*,Error): 0.033461 0.033461 0.000000 (R,S): 0.886512 0.028674 (U,U*,Error): 0.022535 0.022535 0.000000 (Ur,Ur*,Error): 0.050840 0.050840 0.000000 (Us,Us*,Error): 0.785903 0.785903 0.000000 (R,S): 0.489901 0.167927 (U,U*,Error): 0.040303 0.040303 0.000000 (Ur,Ur*,Error): 0.164535 0.164535 0.000000 (Us,Us*,Error): 0.240003 0.240003 0.000000 (R,S): 0.978681 0.712694 (U,U*,Error): 0.682630 0.682630 0.000000 (Ur,Ur*,Error): 1.395001 1.395001 0.000000 (Us,Us*,Error): 0.957816 0.957816 0.000000 (R,S): 0.500472 0.471088 (U,U*,Error): 0.117994 0.117994 0.000000 (Ur,Ur*,Error): 0.471533 0.471533 0.000000 (Us,Us*,Error): 0.250472 0.250472 0.000000 Interpolate R^3 * S^0 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.037037 (R,S,U): 0.666667 0.000000 0.296296 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 1.000000 0.333333 1.000000 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 1.000000 0.666667 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.333333 1.000000 0.037037 (R,S,U): 0.666667 1.000000 0.296296 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.059619 0.681972 (U,U*,Error): 0.000212 0.000212 0.000000 (Ur,Ur*,Error): 0.010663 0.010663 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.042431 0.071445 (U,U*,Error): 0.000076 0.000076 0.000000 (Ur,Ur*,Error): 0.005401 0.005401 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.521650 0.096730 (U,U*,Error): 0.141951 0.141951 0.000000 (Ur,Ur*,Error): 0.816356 0.816356 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.818149 0.817547 (U,U*,Error): 0.547642 0.547642 0.000000 (Ur,Ur*,Error): 2.008101 2.008101 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.722440 0.149865 (U,U*,Error): 0.377055 0.377055 0.000000 (Ur,Ur*,Error): 1.565757 1.565757 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 Interpolate R^1 * S^3 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 1.000000 0.333333 0.037037 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 1.000000 0.666667 0.296296 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.333333 1.000000 0.333333 (R,S,U): 0.666667 1.000000 0.666667 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.659605 0.518595 (U,U*,Error): 0.091996 0.091996 0.000000 (Ur,Ur*,Error): 0.139471 0.139471 0.000000 (Us,Us*,Error): 0.532184 0.532184 0.000000 (R,S): 0.972975 0.648991 (U,U*,Error): 0.265961 0.265961 0.000000 (Ur,Ur*,Error): 0.273349 0.273349 0.000000 (Us,Us*,Error): 1.229421 1.229421 0.000000 (R,S): 0.800331 0.453798 (U,U*,Error): 0.074792 0.074792 0.000000 (Ur,Ur*,Error): 0.093452 0.093452 0.000000 (Us,Us*,Error): 0.494442 0.494442 0.000000 (R,S): 0.432392 0.825314 (U,U*,Error): 0.243072 0.243072 0.000000 (Ur,Ur*,Error): 0.562157 0.562157 0.000000 (Us,Us*,Error): 0.883561 0.883561 0.000000 (R,S): 0.083470 0.133171 (U,U*,Error): 0.000197 0.000197 0.000000 (Ur,Ur*,Error): 0.002362 0.002362 0.000000 (Us,Us*,Error): 0.004441 0.004441 0.000000 Interpolate R^3 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 1.000000 0.333333 0.333333 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 1.000000 0.666667 0.666667 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.333333 1.000000 0.037037 (R,S,U): 0.666667 1.000000 0.296296 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.173389 0.390938 (U,U*,Error): 0.002038 0.002038 0.000000 (Ur,Ur*,Error): 0.035259 0.035259 0.000000 (Us,Us*,Error): 0.005213 0.005213 0.000000 (R,S): 0.831380 0.803364 (U,U*,Error): 0.461648 0.461648 0.000000 (Ur,Ur*,Error): 1.665838 1.665838 0.000000 (Us,Us*,Error): 0.574643 0.574643 0.000000 (R,S): 0.060471 0.399258 (U,U*,Error): 0.000088 0.000088 0.000000 (Ur,Ur*,Error): 0.004380 0.004380 0.000000 (Us,Us*,Error): 0.000221 0.000221 0.000000 (R,S): 0.526876 0.416799 (U,U*,Error): 0.060961 0.060961 0.000000 (Ur,Ur*,Error): 0.347108 0.347108 0.000000 (Us,Us*,Error): 0.146260 0.146260 0.000000 (R,S): 0.656860 0.627973 (U,U*,Error): 0.177975 0.177975 0.000000 (Ur,Ur*,Error): 0.812845 0.812845 0.000000 (Us,Us*,Error): 0.283412 0.283412 0.000000 INTERP_TEST for element "Q16". Element order = 16 Interpolate R^0 * S^0 (R,S,U): 0.000000 0.000000 1.000000 (R,S,U): 0.333333 0.000000 1.000000 (R,S,U): 0.666667 0.000000 1.000000 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 0.000000 0.333333 1.000000 (R,S,U): 0.333333 0.333333 1.000000 (R,S,U): 0.666667 0.333333 1.000000 (R,S,U): 1.000000 0.333333 1.000000 (R,S,U): 0.000000 0.666667 1.000000 (R,S,U): 0.333333 0.666667 1.000000 (R,S,U): 0.666667 0.666667 1.000000 (R,S,U): 1.000000 0.666667 1.000000 (R,S,U): 0.000000 1.000000 1.000000 (R,S,U): 0.333333 1.000000 1.000000 (R,S,U): 0.666667 1.000000 1.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.291984 0.431651 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.015487 0.984064 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.167168 0.106216 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.372410 0.198118 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.489688 0.339493 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 Interpolate R^0 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 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.000000 0.333333 0.333333 (R,S,U): 0.000000 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.000000 0.666667 0.666667 (R,S,U): 0.000000 1.000000 1.000000 (R,S,U): 0.333333 1.000000 1.000000 (R,S,U): 0.666667 1.000000 1.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.951630 0.920332 (U,U*,Error): 0.920332 0.920332 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.052677 0.737858 (U,U*,Error): 0.737858 0.737858 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.269119 0.422836 (U,U*,Error): 0.422836 0.422836 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.547871 0.942737 (U,U*,Error): 0.942737 0.942737 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.417744 0.983052 (U,U*,Error): 0.983052 0.983052 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 Interpolate R^1 * S^0 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.333333 (R,S,U): 0.666667 0.000000 0.666667 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 0.333333 0.333333 0.333333 (R,S,U): 0.666667 0.333333 0.666667 (R,S,U): 1.000000 0.333333 1.000000 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 0.333333 0.666667 0.333333 (R,S,U): 0.666667 0.666667 0.666667 (R,S,U): 1.000000 0.666667 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.333333 1.000000 0.333333 (R,S,U): 0.666667 1.000000 0.666667 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.301455 0.701099 (U,U*,Error): 0.301455 0.301455 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.666339 0.539126 (U,U*,Error): 0.666339 0.666339 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.698106 0.666528 (U,U*,Error): 0.698106 0.698106 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.178132 0.128014 (U,U*,Error): 0.178132 0.178132 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.999080 0.171121 (U,U*,Error): 0.999080 0.999080 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 Interpolate R^0 * S^2 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 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.000000 0.333333 0.111111 (R,S,U): 0.000000 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.000000 0.666667 0.444444 (R,S,U): 0.000000 1.000000 1.000000 (R,S,U): 0.333333 1.000000 1.000000 (R,S,U): 0.666667 1.000000 1.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.032601 0.561200 (U,U*,Error): 0.314945 0.314945 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.122400 1.122400 0.000000 (R,S): 0.881867 0.669175 (U,U*,Error): 0.447796 0.447796 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 1.338351 1.338351 0.000000 (R,S): 0.190433 0.368917 (U,U*,Error): 0.136099 0.136099 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 0.737833 0.737833 0.000000 (R,S): 0.460726 0.981638 (U,U*,Error): 0.963613 0.963613 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.963276 1.963276 0.000000 (R,S): 0.156405 0.855523 (U,U*,Error): 0.731919 0.731919 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 1.711046 1.711046 0.000000 Interpolate R^1 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 0.333333 0.333333 0.111111 (R,S,U): 0.666667 0.333333 0.222222 (R,S,U): 1.000000 0.333333 0.333333 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 0.333333 0.666667 0.222222 (R,S,U): 0.666667 0.666667 0.444444 (R,S,U): 1.000000 0.666667 0.666667 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.333333 1.000000 0.333333 (R,S,U): 0.666667 1.000000 0.666667 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.644765 0.376272 (U,U*,Error): 0.242607 0.242607 0.000000 (Ur,Ur*,Error): 0.376272 0.376272 0.000000 (Us,Us*,Error): 0.644765 0.644765 0.000000 (R,S): 0.190924 0.428253 (U,U*,Error): 0.081764 0.081764 0.000000 (Ur,Ur*,Error): 0.428253 0.428253 0.000000 (Us,Us*,Error): 0.190924 0.190924 0.000000 (R,S): 0.482022 0.120612 (U,U*,Error): 0.058137 0.058137 0.000000 (Ur,Ur*,Error): 0.120612 0.120612 0.000000 (Us,Us*,Error): 0.482022 0.482022 0.000000 (R,S): 0.589507 0.226188 (U,U*,Error): 0.133339 0.133339 0.000000 (Ur,Ur*,Error): 0.226188 0.226188 0.000000 (Us,Us*,Error): 0.589507 0.589507 0.000000 (R,S): 0.384619 0.582986 (U,U*,Error): 0.224228 0.224228 0.000000 (Ur,Ur*,Error): 0.582986 0.582986 0.000000 (Us,Us*,Error): 0.384619 0.384619 0.000000 Interpolate R^2 * S^0 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.111111 (R,S,U): 0.666667 0.000000 0.444444 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 0.333333 0.333333 0.111111 (R,S,U): 0.666667 0.333333 0.444444 (R,S,U): 1.000000 0.333333 1.000000 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 0.333333 0.666667 0.111111 (R,S,U): 0.666667 0.666667 0.444444 (R,S,U): 1.000000 0.666667 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.333333 1.000000 0.111111 (R,S,U): 0.666667 1.000000 0.444444 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.251806 0.290441 (U,U*,Error): 0.063406 0.063406 0.000000 (Ur,Ur*,Error): 0.503612 0.503612 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.617091 0.265281 (U,U*,Error): 0.380801 0.380801 0.000000 (Ur,Ur*,Error): 1.234182 1.234182 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.824376 0.982663 (U,U*,Error): 0.679596 0.679596 0.000000 (Ur,Ur*,Error): 1.648753 1.648753 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.730249 0.343877 (U,U*,Error): 0.533263 0.533263 0.000000 (Ur,Ur*,Error): 1.460498 1.460498 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.584069 0.107769 (U,U*,Error): 0.341137 0.341137 0.000000 (Ur,Ur*,Error): 1.168139 1.168139 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 Interpolate R^0 * S^3 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 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.000000 0.333333 0.037037 (R,S,U): 0.000000 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.000000 0.666667 0.296296 (R,S,U): 0.000000 1.000000 1.000000 (R,S,U): 0.333333 1.000000 1.000000 (R,S,U): 0.666667 1.000000 1.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.906308 0.879654 (U,U*,Error): 0.680668 0.680668 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 2.321372 2.321372 0.000000 (R,S): 0.817761 0.260728 (U,U*,Error): 0.017724 0.017724 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 0.203937 0.203937 0.000000 (R,S): 0.594356 0.022513 (U,U*,Error): 0.000011 0.000011 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.001520 0.001520 0.000000 (R,S): 0.425259 0.312719 (U,U*,Error): 0.030582 0.030582 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 0.293379 0.293379 0.000000 (R,S): 0.161485 0.178766 (U,U*,Error): 0.005713 0.005713 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 0.095872 0.095872 0.000000 Interpolate R^1 * S^2 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 0.333333 0.333333 0.037037 (R,S,U): 0.666667 0.333333 0.074074 (R,S,U): 1.000000 0.333333 0.111111 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 0.333333 0.666667 0.148148 (R,S,U): 0.666667 0.666667 0.296296 (R,S,U): 1.000000 0.666667 0.444444 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.333333 1.000000 0.333333 (R,S,U): 0.666667 1.000000 0.666667 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.422886 0.094229 (U,U*,Error): 0.003755 0.003755 0.000000 (Ur,Ur*,Error): 0.008879 0.008879 0.000000 (Us,Us*,Error): 0.079696 0.079696 0.000000 (R,S): 0.598524 0.470924 (U,U*,Error): 0.132734 0.132734 0.000000 (Ur,Ur*,Error): 0.221770 0.221770 0.000000 (Us,Us*,Error): 0.563719 0.563719 0.000000 (R,S): 0.695949 0.699888 (U,U*,Error): 0.340906 0.340906 0.000000 (Ur,Ur*,Error): 0.489843 0.489843 0.000000 (Us,Us*,Error): 0.974173 0.974173 0.000000 (R,S): 0.638531 0.033604 (U,U*,Error): 0.000721 0.000721 0.000000 (Ur,Ur*,Error): 0.001129 0.001129 0.000000 (Us,Us*,Error): 0.042914 0.042914 0.000000 (R,S): 0.068806 0.319600 (U,U*,Error): 0.007028 0.007028 0.000000 (Ur,Ur*,Error): 0.102144 0.102144 0.000000 (Us,Us*,Error): 0.043981 0.043981 0.000000 Interpolate R^2 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 0.333333 0.333333 0.037037 (R,S,U): 0.666667 0.333333 0.148148 (R,S,U): 1.000000 0.333333 0.333333 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 0.333333 0.666667 0.074074 (R,S,U): 0.666667 0.666667 0.296296 (R,S,U): 1.000000 0.666667 0.666667 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.333333 1.000000 0.111111 (R,S,U): 0.666667 1.000000 0.444444 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.530864 0.654446 (U,U*,Error): 0.184434 0.184434 0.000000 (Ur,Ur*,Error): 0.694844 0.694844 0.000000 (Us,Us*,Error): 0.281817 0.281817 0.000000 (R,S): 0.407619 0.819981 (U,U*,Error): 0.136243 0.136243 0.000000 (Ur,Ur*,Error): 0.668480 0.668480 0.000000 (Us,Us*,Error): 0.166153 0.166153 0.000000 (R,S): 0.718359 0.968649 (U,U*,Error): 0.499861 0.499861 0.000000 (Ur,Ur*,Error): 1.391676 1.391676 0.000000 (Us,Us*,Error): 0.516040 0.516040 0.000000 (R,S): 0.531334 0.325146 (U,U*,Error): 0.091794 0.091794 0.000000 (Ur,Ur*,Error): 0.345522 0.345522 0.000000 (Us,Us*,Error): 0.282316 0.282316 0.000000 (R,S): 0.105629 0.610959 (U,U*,Error): 0.006817 0.006817 0.000000 (Ur,Ur*,Error): 0.129070 0.129070 0.000000 (Us,Us*,Error): 0.011158 0.011158 0.000000 Interpolate R^3 * S^0 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.037037 (R,S,U): 0.666667 0.000000 0.296296 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 0.333333 0.333333 0.037037 (R,S,U): 0.666667 0.333333 0.296296 (R,S,U): 1.000000 0.333333 1.000000 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 0.333333 0.666667 0.037037 (R,S,U): 0.666667 0.666667 0.296296 (R,S,U): 1.000000 0.666667 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.333333 1.000000 0.037037 (R,S,U): 0.666667 1.000000 0.296296 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.778802 0.423453 (U,U*,Error): 0.472369 0.472369 0.000000 (Ur,Ur*,Error): 1.819599 1.819599 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.090823 0.266471 (U,U*,Error): 0.000749 0.000749 0.000000 (Ur,Ur*,Error): 0.024747 0.024747 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.153657 0.281005 (U,U*,Error): 0.003628 0.003628 0.000000 (Ur,Ur*,Error): 0.070831 0.070831 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.440085 0.527143 (U,U*,Error): 0.085233 0.085233 0.000000 (Ur,Ur*,Error): 0.581025 0.581025 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.457424 0.875372 (U,U*,Error): 0.095710 0.095710 0.000000 (Ur,Ur*,Error): 0.627711 0.627711 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 Interpolate R^1 * S^3 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 0.333333 0.333333 0.012346 (R,S,U): 0.666667 0.333333 0.024691 (R,S,U): 1.000000 0.333333 0.037037 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 0.333333 0.666667 0.098765 (R,S,U): 0.666667 0.666667 0.197531 (R,S,U): 1.000000 0.666667 0.296296 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.333333 1.000000 0.333333 (R,S,U): 0.666667 1.000000 0.666667 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.518052 0.943623 (U,U*,Error): 0.435280 0.435280 0.000000 (Ur,Ur*,Error): 0.840224 0.840224 0.000000 (Us,Us*,Error): 1.383858 1.383858 0.000000 (R,S): 0.637709 0.957694 (U,U*,Error): 0.560148 0.560148 0.000000 (Ur,Ur*,Error): 0.878376 0.878376 0.000000 (Us,Us*,Error): 1.754678 1.754678 0.000000 (R,S): 0.240707 0.676122 (U,U*,Error): 0.074399 0.074399 0.000000 (Ur,Ur*,Error): 0.309083 0.309083 0.000000 (Us,Us*,Error): 0.330111 0.330111 0.000000 (R,S): 0.289065 0.671808 (U,U*,Error): 0.087646 0.087646 0.000000 (Ur,Ur*,Error): 0.303205 0.303205 0.000000 (Us,Us*,Error): 0.391387 0.391387 0.000000 (R,S): 0.695140 0.067993 (U,U*,Error): 0.000219 0.000219 0.000000 (Ur,Ur*,Error): 0.000314 0.000314 0.000000 (Us,Us*,Error): 0.009641 0.009641 0.000000 Interpolate R^2 * S^2 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 0.333333 0.333333 0.012346 (R,S,U): 0.666667 0.333333 0.049383 (R,S,U): 1.000000 0.333333 0.111111 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 0.333333 0.666667 0.049383 (R,S,U): 0.666667 0.666667 0.197531 (R,S,U): 1.000000 0.666667 0.444444 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.333333 1.000000 0.111111 (R,S,U): 0.666667 1.000000 0.444444 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.254790 0.224040 (U,U*,Error): 0.003258 0.003258 0.000000 (Ur,Ur*,Error): 0.025578 0.025578 0.000000 (Us,Us*,Error): 0.029088 0.029088 0.000000 (R,S): 0.667833 0.844392 (U,U*,Error): 0.317998 0.317998 0.000000 (Ur,Ur*,Error): 0.952327 0.952327 0.000000 (Us,Us*,Error): 0.753199 0.753199 0.000000 (R,S): 0.344462 0.780520 (U,U*,Error): 0.072286 0.072286 0.000000 (Ur,Ur*,Error): 0.419701 0.419701 0.000000 (Us,Us*,Error): 0.185224 0.185224 0.000000 (R,S): 0.675332 0.006715 (U,U*,Error): 0.000021 0.000021 0.000000 (Ur,Ur*,Error): 0.000061 0.000061 0.000000 (Us,Us*,Error): 0.006125 0.006125 0.000000 (R,S): 0.602170 0.386771 (U,U*,Error): 0.054243 0.054243 0.000000 (Ur,Ur*,Error): 0.180160 0.180160 0.000000 (Us,Us*,Error): 0.280494 0.280494 0.000000 Interpolate R^3 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 0.333333 0.333333 0.012346 (R,S,U): 0.666667 0.333333 0.098765 (R,S,U): 1.000000 0.333333 0.333333 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 0.333333 0.666667 0.024691 (R,S,U): 0.666667 0.666667 0.197531 (R,S,U): 1.000000 0.666667 0.666667 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.333333 1.000000 0.037037 (R,S,U): 0.666667 1.000000 0.296296 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.915991 0.001151 (U,U*,Error): 0.000885 0.000885 0.000000 (Ur,Ur*,Error): 0.002897 0.002897 0.000000 (Us,Us*,Error): 0.768553 0.768553 0.000000 (R,S): 0.462449 0.424349 (U,U*,Error): 0.041968 0.041968 0.000000 (Ur,Ur*,Error): 0.272253 0.272253 0.000000 (Us,Us*,Error): 0.098899 0.098899 0.000000 (R,S): 0.460916 0.770160 (U,U*,Error): 0.075413 0.075413 0.000000 (Ur,Ur*,Error): 0.490847 0.490847 0.000000 (Us,Us*,Error): 0.097919 0.097919 0.000000 (R,S): 0.322472 0.784739 (U,U*,Error): 0.026315 0.026315 0.000000 (Ur,Ur*,Error): 0.244811 0.244811 0.000000 (Us,Us*,Error): 0.033533 0.033533 0.000000 (R,S): 0.471357 0.035763 (U,U*,Error): 0.003745 0.003745 0.000000 (Ur,Ur*,Error): 0.023837 0.023837 0.000000 (Us,Us*,Error): 0.104725 0.104725 0.000000 Interpolate R^2 * S^3 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 0.333333 0.333333 0.004115 (R,S,U): 0.666667 0.333333 0.016461 (R,S,U): 1.000000 0.333333 0.037037 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 0.333333 0.666667 0.032922 (R,S,U): 0.666667 0.666667 0.131687 (R,S,U): 1.000000 0.666667 0.296296 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.333333 1.000000 0.111111 (R,S,U): 0.666667 1.000000 0.444444 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.175874 0.721758 (U,U*,Error): 0.011630 0.011630 0.000000 (Ur,Ur*,Error): 0.132254 0.132254 0.000000 (Us,Us*,Error): 0.048340 0.048340 0.000000 (R,S): 0.473486 0.152721 (U,U*,Error): 0.000799 0.000799 0.000000 (Ur,Ur*,Error): 0.003373 0.003373 0.000000 (Us,Us*,Error): 0.015687 0.015687 0.000000 (R,S): 0.341125 0.607389 (U,U*,Error): 0.026075 0.026075 0.000000 (Ur,Ur*,Error): 0.152878 0.152878 0.000000 (Us,Us*,Error): 0.128790 0.128790 0.000000 (R,S): 0.191745 0.738427 (U,U*,Error): 0.014804 0.014804 0.000000 (Ur,Ur*,Error): 0.154411 0.154411 0.000000 (Us,Us*,Error): 0.060143 0.060143 0.000000 (R,S): 0.242850 0.917424 (U,U*,Error): 0.045539 0.045539 0.000000 (Ur,Ur*,Error): 0.375040 0.375040 0.000000 (Us,Us*,Error): 0.148914 0.148914 0.000000 Interpolate R^3 * S^2 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 0.333333 0.333333 0.004115 (R,S,U): 0.666667 0.333333 0.032922 (R,S,U): 1.000000 0.333333 0.111111 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 0.333333 0.666667 0.016461 (R,S,U): 0.666667 0.666667 0.131687 (R,S,U): 1.000000 0.666667 0.444444 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.333333 1.000000 0.037037 (R,S,U): 0.666667 1.000000 0.296296 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.269062 0.765500 (U,U*,Error): 0.011414 0.011414 0.000000 (Ur,Ur*,Error): 0.127267 0.127267 0.000000 (Us,Us*,Error): 0.029822 0.029822 0.000000 (R,S): 0.188662 0.287498 (U,U*,Error): 0.000555 0.000555 0.000000 (Ur,Ur*,Error): 0.008826 0.008826 0.000000 (Us,Us*,Error): 0.003861 0.003861 0.000000 (R,S): 0.091113 0.576209 (U,U*,Error): 0.000251 0.000251 0.000000 (Ur,Ur*,Error): 0.008269 0.008269 0.000000 (Us,Us*,Error): 0.000872 0.000872 0.000000 (R,S): 0.683363 0.546593 (U,U*,Error): 0.095342 0.095342 0.000000 (Ur,Ur*,Error): 0.418555 0.418555 0.000000 (Us,Us*,Error): 0.348858 0.348858 0.000000 (R,S): 0.425729 0.644443 (U,U*,Error): 0.032046 0.032046 0.000000 (Ur,Ur*,Error): 0.225817 0.225817 0.000000 (Us,Us*,Error): 0.099452 0.099452 0.000000 Interpolate R^3 * S^3 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 0.333333 0.333333 0.001372 (R,S,U): 0.666667 0.333333 0.010974 (R,S,U): 1.000000 0.333333 0.037037 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 0.333333 0.666667 0.010974 (R,S,U): 0.666667 0.666667 0.087791 (R,S,U): 1.000000 0.666667 0.296296 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.333333 1.000000 0.037037 (R,S,U): 0.666667 1.000000 0.296296 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.647618 0.679017 (U,U*,Error): 0.085035 0.085035 0.000000 (Ur,Ur*,Error): 0.393913 0.393913 0.000000 (Us,Us*,Error): 0.375697 0.375697 0.000000 (R,S): 0.635787 0.945174 (U,U*,Error): 0.217005 0.217005 0.000000 (Ur,Ur*,Error): 1.023952 1.023952 0.000000 (Us,Us*,Error): 0.688778 0.688778 0.000000 (R,S): 0.208935 0.709282 (U,U*,Error): 0.003255 0.003255 0.000000 (Ur,Ur*,Error): 0.046730 0.046730 0.000000 (Us,Us*,Error): 0.013765 0.013765 0.000000 (R,S): 0.236231 0.119396 (U,U*,Error): 0.000022 0.000022 0.000000 (Ur,Ur*,Error): 0.000285 0.000285 0.000000 (Us,Us*,Error): 0.000564 0.000564 0.000000 (R,S): 0.607304 0.450138 (U,U*,Error): 0.020429 0.020429 0.000000 (Ur,Ur*,Error): 0.100918 0.100918 0.000000 (Us,Us*,Error): 0.136154 0.136154 0.000000 INTERP_TEST for element "QL". Element order = 6 Interpolate R^0 * S^0 (R,S,U): 0.000000 0.000000 1.000000 (R,S,U): 0.500000 0.000000 1.000000 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 0.000000 1.000000 1.000000 (R,S,U): 0.500000 1.000000 1.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.458725 0.661945 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.770286 0.350218 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.662010 0.416159 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.841929 0.832917 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.256441 0.613461 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 Interpolate R^0 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.500000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 1.000000 1.000000 (R,S,U): 0.500000 1.000000 1.000000 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.582249 0.540739 (U,U*,Error): 0.540739 0.540739 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.869941 0.264779 (U,U*,Error): 0.264779 0.264779 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.318074 0.119215 (U,U*,Error): 0.119215 0.119215 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.939829 0.645552 (U,U*,Error): 0.645552 0.645552 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.479463 0.639317 (U,U*,Error): 0.639317 0.639317 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 Interpolate R^1 * S^0 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.500000 0.000000 0.500000 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.500000 1.000000 0.500000 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.544716 0.647311 (U,U*,Error): 0.544716 0.544716 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.543886 0.721047 (U,U*,Error): 0.543886 0.543886 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.522495 0.993705 (U,U*,Error): 0.522495 0.522495 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.218677 0.105798 (U,U*,Error): 0.218677 0.218677 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.109697 0.063591 (U,U*,Error): 0.109697 0.109697 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 Interpolate R^1 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.500000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.500000 1.000000 0.500000 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.404580 0.448373 (U,U*,Error): 0.181403 0.181403 0.000000 (Ur,Ur*,Error): 0.448373 0.448373 0.000000 (Us,Us*,Error): 0.404580 0.404580 0.000000 (R,S): 0.365816 0.763505 (U,U*,Error): 0.279302 0.279302 0.000000 (Ur,Ur*,Error): 0.763505 0.763505 0.000000 (Us,Us*,Error): 0.365816 0.365816 0.000000 (R,S): 0.627896 0.771980 (U,U*,Error): 0.484724 0.484724 0.000000 (Ur,Ur*,Error): 0.771980 0.771980 0.000000 (Us,Us*,Error): 0.627896 0.627896 0.000000 (R,S): 0.932854 0.972741 (U,U*,Error): 0.907425 0.907425 0.000000 (Ur,Ur*,Error): 0.972741 0.972741 0.000000 (Us,Us*,Error): 0.932854 0.932854 0.000000 (R,S): 0.192028 0.138874 (U,U*,Error): 0.026668 0.026668 0.000000 (Ur,Ur*,Error): 0.138874 0.138874 0.000000 (Us,Us*,Error): 0.192028 0.192028 0.000000 Interpolate R^2 * S^0 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.500000 0.000000 0.250000 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.500000 1.000000 0.250000 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.696266 0.093820 (U,U*,Error): 0.484787 0.484787 0.000000 (Ur,Ur*,Error): 1.392533 1.392533 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.525404 0.530344 (U,U*,Error): 0.276050 0.276050 0.000000 (Ur,Ur*,Error): 1.050809 1.050809 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.861140 0.484853 (U,U*,Error): 0.741562 0.741562 0.000000 (Ur,Ur*,Error): 1.722280 1.722280 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.393456 0.671431 (U,U*,Error): 0.154808 0.154808 0.000000 (Ur,Ur*,Error): 0.786913 0.786913 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.741258 0.520052 (U,U*,Error): 0.549463 0.549463 0.000000 (Ur,Ur*,Error): 1.482516 1.482516 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 Interpolate R^2 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.500000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.500000 1.000000 0.250000 (R,S,U): 1.000000 1.000000 1.000000 (R,S): 0.347713 0.149997 (U,U*,Error): 0.018135 0.018135 0.000000 (Ur,Ur*,Error): 0.104312 0.104312 0.000000 (Us,Us*,Error): 0.120904 0.120904 0.000000 (R,S): 0.586092 0.262145 (U,U*,Error): 0.090048 0.090048 0.000000 (Ur,Ur*,Error): 0.307283 0.307283 0.000000 (Us,Us*,Error): 0.343504 0.343504 0.000000 (R,S): 0.044454 0.754933 (U,U*,Error): 0.001492 0.001492 0.000000 (Ur,Ur*,Error): 0.067120 0.067120 0.000000 (Us,Us*,Error): 0.001976 0.001976 0.000000 (R,S): 0.242785 0.442402 (U,U*,Error): 0.026077 0.026077 0.000000 (Ur,Ur*,Error): 0.214818 0.214818 0.000000 (Us,Us*,Error): 0.058945 0.058945 0.000000 (R,S): 0.687796 0.359228 (U,U*,Error): 0.169938 0.169938 0.000000 (Ur,Ur*,Error): 0.494152 0.494152 0.000000 (Us,Us*,Error): 0.473063 0.473063 0.000000 INTERP_TEST for element "T3". Element order = 3 Interpolate R^0 * S^0 (R,S,U): 0.000000 0.000000 1.000000 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 0.000000 1.000000 1.000000 (R,S): 0.263660 0.605293 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.316584 0.295953 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.442305 0.019578 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.330858 0.424309 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.270270 0.197054 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 Interpolate R^0 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 1.000000 1.000000 (R,S): 0.178279 0.570079 (U,U*,Error): 0.570079 0.570079 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.112229 0.608817 (U,U*,Error): 0.608817 0.608817 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.230886 0.603208 (U,U*,Error): 0.603208 0.603208 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.191486 0.244923 (U,U*,Error): 0.244923 0.244923 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.377396 0.216019 (U,U*,Error): 0.216019 0.216019 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 Interpolate R^1 * S^0 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S): 0.209593 0.050696 (U,U*,Error): 0.209593 0.209593 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.327565 0.671264 (U,U*,Error): 0.327565 0.327565 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.561355 0.166499 (U,U*,Error): 0.561355 0.561355 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.768854 0.167254 (U,U*,Error): 0.768854 0.768854 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.138020 0.010128 (U,U*,Error): 0.138020 0.138020 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 INTERP_TEST - Warning! Skipping test for element "T4". INTERP_TEST for element "T6". Element order = 6 Interpolate R^0 * S^0 (R,S,U): 0.000000 0.000000 1.000000 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 0.000000 1.000000 1.000000 (R,S,U): 0.500000 0.000000 1.000000 (R,S,U): 0.500000 0.500000 1.000000 (R,S,U): 0.000000 0.500000 1.000000 (R,S): 0.485577 0.115719 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.588026 0.154752 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.199863 0.406955 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.251294 0.174416 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.210037 0.681476 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 Interpolate R^0 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 1.000000 1.000000 (R,S,U): 0.500000 0.000000 0.000000 (R,S,U): 0.500000 0.500000 0.500000 (R,S,U): 0.000000 0.500000 0.500000 (R,S): 0.534064 0.089951 (U,U*,Error): 0.089951 0.089951 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.111706 0.136293 (U,U*,Error): 0.136293 0.136293 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.321348 0.504823 (U,U*,Error): 0.504823 0.504823 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.189710 0.495006 (U,U*,Error): 0.495006 0.495006 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.147608 0.054974 (U,U*,Error): 0.054974 0.054974 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 Interpolate R^1 * S^0 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.500000 0.000000 0.500000 (R,S,U): 0.500000 0.500000 0.500000 (R,S,U): 0.000000 0.500000 0.000000 (R,S): 0.149287 0.439440 (U,U*,Error): 0.149287 0.149287 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.070391 0.303333 (U,U*,Error): 0.070391 0.070391 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.417209 0.184603 (U,U*,Error): 0.417209 0.417209 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.120986 0.011088 (U,U*,Error): 0.120986 0.120986 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.000522 0.865439 (U,U*,Error): 0.000522 0.000522 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 Interpolate R^0 * S^2 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 1.000000 1.000000 (R,S,U): 0.500000 0.000000 0.000000 (R,S,U): 0.500000 0.500000 0.250000 (R,S,U): 0.000000 0.500000 0.250000 (R,S): 0.387434 0.010050 (U,U*,Error): 0.000101 0.000101 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.020100 0.020100 0.000000 (R,S): 0.472320 0.520477 (U,U*,Error): 0.270896 0.270896 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.040953 1.040953 0.000000 (R,S): 0.198652 0.772157 (U,U*,Error): 0.596227 0.596227 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.544314 1.544314 0.000000 (R,S): 0.501906 0.099148 (U,U*,Error): 0.009830 0.009830 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.198295 0.198295 0.000000 (R,S): 0.425339 0.154822 (U,U*,Error): 0.023970 0.023970 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.309644 0.309644 0.000000 Interpolate R^1 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.500000 0.000000 0.000000 (R,S,U): 0.500000 0.500000 0.250000 (R,S,U): 0.000000 0.500000 0.000000 (R,S): 0.261360 0.414013 (U,U*,Error): 0.108206 0.108206 0.000000 (Ur,Ur*,Error): 0.414013 0.414013 0.000000 (Us,Us*,Error): 0.261360 0.261360 0.000000 (R,S): 0.246735 0.666416 (U,U*,Error): 0.164428 0.164428 0.000000 (Ur,Ur*,Error): 0.666416 0.666416 0.000000 (Us,Us*,Error): 0.246735 0.246735 0.000000 (R,S): 0.083483 0.625960 (U,U*,Error): 0.052257 0.052257 0.000000 (Ur,Ur*,Error): 0.625960 0.625960 0.000000 (Us,Us*,Error): 0.083483 0.083483 0.000000 (R,S): 0.339055 0.270248 (U,U*,Error): 0.091629 0.091629 0.000000 (Ur,Ur*,Error): 0.270248 0.270248 0.000000 (Us,Us*,Error): 0.339055 0.339055 0.000000 (R,S): 0.109248 0.017697 (U,U*,Error): 0.001933 0.001933 0.000000 (Ur,Ur*,Error): 0.017697 0.017697 0.000000 (Us,Us*,Error): 0.109248 0.109248 0.000000 Interpolate R^2 * S^0 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 0.000000 1.000000 0.000000 (R,S,U): 0.500000 0.000000 0.250000 (R,S,U): 0.500000 0.500000 0.250000 (R,S,U): 0.000000 0.500000 0.000000 (R,S): 0.230971 0.418554 (U,U*,Error): 0.053348 0.053348 0.000000 (Ur,Ur*,Error): 0.461942 0.461942 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.071687 0.419910 (U,U*,Error): 0.005139 0.005139 0.000000 (Ur,Ur*,Error): 0.143374 0.143374 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.016983 0.120860 (U,U*,Error): 0.000288 0.000288 0.000000 (Ur,Ur*,Error): 0.033966 0.033966 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.137289 0.515703 (U,U*,Error): 0.018848 0.018848 0.000000 (Ur,Ur*,Error): 0.274579 0.274579 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.155144 0.790595 (U,U*,Error): 0.024070 0.024070 0.000000 (Ur,Ur*,Error): 0.310289 0.310289 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 INTERP_TEST for element "T10". Element order = 10 Interpolate R^0 * S^0 (R,S,U): 0.000000 0.000000 1.000000 (R,S,U): 0.333333 0.000000 1.000000 (R,S,U): 0.666667 0.000000 1.000000 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 0.000000 0.333333 1.000000 (R,S,U): 0.333333 0.333333 1.000000 (R,S,U): 0.666667 0.333333 1.000000 (R,S,U): 0.000000 0.666667 1.000000 (R,S,U): 0.333333 0.666667 1.000000 (R,S,U): 0.000000 1.000000 1.000000 (R,S): 0.447709 0.370117 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.031991 0.614713 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.362411 0.049533 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.489570 0.192510 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.123084 0.205494 (U,U*,Error): 1.000000 1.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 Interpolate R^0 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 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.000000 0.666667 0.666667 (R,S,U): 0.333333 0.666667 0.666667 (R,S,U): 0.000000 1.000000 1.000000 (R,S): 0.146515 0.189072 (U,U*,Error): 0.189072 0.189072 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.042652 0.635198 (U,U*,Error): 0.635198 0.635198 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.281867 0.538597 (U,U*,Error): 0.538597 0.538597 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.304837 0.500884 (U,U*,Error): 0.500884 0.500884 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 (R,S): 0.535801 0.445183 (U,U*,Error): 0.445183 0.445183 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.000000 1.000000 0.000000 Interpolate R^1 * S^0 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.333333 (R,S,U): 0.666667 0.000000 0.666667 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 0.333333 0.333333 0.333333 (R,S,U): 0.666667 0.333333 0.666667 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 0.333333 0.666667 0.333333 (R,S,U): 0.000000 1.000000 0.000000 (R,S): 0.123932 0.490357 (U,U*,Error): 0.123932 0.123932 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.147002 0.126073 (U,U*,Error): 0.147002 0.147002 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.270294 0.208461 (U,U*,Error): 0.270294 0.270294 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.435020 0.359688 (U,U*,Error): 0.435020 0.435020 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.417029 0.205976 (U,U*,Error): 0.417029 0.417029 0.000000 (Ur,Ur*,Error): 1.000000 1.000000 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 Interpolate R^0 * S^2 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 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.000000 0.666667 0.444444 (R,S,U): 0.333333 0.666667 0.444444 (R,S,U): 0.000000 1.000000 1.000000 (R,S): 0.052067 0.917929 (U,U*,Error): 0.842593 0.842593 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.835858 1.835858 0.000000 (R,S): 0.105709 0.142041 (U,U*,Error): 0.020176 0.020176 0.000000 (Ur,Ur*,Error): 0.000000 -0.000000 0.000000 (Us,Us*,Error): 0.284082 0.284082 0.000000 (R,S): 0.166460 0.620959 (U,U*,Error): 0.385590 0.385590 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 1.241917 1.241917 0.000000 (R,S): 0.573710 0.052078 (U,U*,Error): 0.002712 0.002712 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.104156 0.104156 0.000000 (R,S): 0.068799 0.271338 (U,U*,Error): 0.073624 0.073624 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.542677 0.542677 0.000000 Interpolate R^1 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 0.333333 0.333333 0.111111 (R,S,U): 0.666667 0.333333 0.222222 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 0.333333 0.666667 0.222222 (R,S,U): 0.000000 1.000000 0.000000 (R,S): 0.737842 0.063405 (U,U*,Error): 0.046782 0.046782 0.000000 (Ur,Ur*,Error): 0.063405 0.063405 0.000000 (Us,Us*,Error): 0.737842 0.737842 0.000000 (R,S): 0.139559 0.065595 (U,U*,Error): 0.009154 0.009154 0.000000 (Ur,Ur*,Error): 0.065595 0.065595 0.000000 (Us,Us*,Error): 0.139559 0.139559 0.000000 (R,S): 0.015602 0.141061 (U,U*,Error): 0.002201 0.002201 0.000000 (Ur,Ur*,Error): 0.141061 0.141061 0.000000 (Us,Us*,Error): 0.015602 0.015602 0.000000 (R,S): 0.214441 0.486623 (U,U*,Error): 0.104352 0.104352 0.000000 (Ur,Ur*,Error): 0.486623 0.486623 0.000000 (Us,Us*,Error): 0.214441 0.214441 0.000000 (R,S): 0.177602 0.398589 (U,U*,Error): 0.070790 0.070790 0.000000 (Ur,Ur*,Error): 0.398589 0.398589 0.000000 (Us,Us*,Error): 0.177602 0.177602 0.000000 Interpolate R^2 * S^0 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.111111 (R,S,U): 0.666667 0.000000 0.444444 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 0.333333 0.333333 0.111111 (R,S,U): 0.666667 0.333333 0.444444 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 0.333333 0.666667 0.111111 (R,S,U): 0.000000 1.000000 0.000000 (R,S): 0.133931 0.030890 (U,U*,Error): 0.017938 0.017938 0.000000 (Ur,Ur*,Error): 0.267863 0.267863 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.060858 0.698694 (U,U*,Error): 0.003704 0.003704 0.000000 (Ur,Ur*,Error): 0.121717 0.121717 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.295534 0.332936 (U,U*,Error): 0.087340 0.087340 0.000000 (Ur,Ur*,Error): 0.591068 0.591068 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.532932 0.351802 (U,U*,Error): 0.284016 0.284016 0.000000 (Ur,Ur*,Error): 1.065864 1.065864 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.025228 0.842207 (U,U*,Error): 0.000636 0.000636 0.000000 (Ur,Ur*,Error): 0.050456 0.050456 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 Interpolate R^0 * S^3 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 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.000000 0.666667 0.296296 (R,S,U): 0.333333 0.666667 0.296296 (R,S,U): 0.000000 1.000000 1.000000 (R,S): 0.440967 0.145900 (U,U*,Error): 0.003106 0.003106 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.063860 0.063860 0.000000 (R,S): 0.347879 0.446027 (U,U*,Error): 0.088732 0.088732 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.596819 0.596819 0.000000 (R,S): 0.054239 0.177108 (U,U*,Error): 0.005555 0.005555 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.094101 0.094101 0.000000 (R,S): 0.662808 0.330829 (U,U*,Error): 0.036209 0.036209 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.328343 0.328343 0.000000 (R,S): 0.101514 0.881845 (U,U*,Error): 0.685767 0.685767 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 2.332951 2.332951 0.000000 Interpolate R^1 * S^2 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 0.333333 0.333333 0.037037 (R,S,U): 0.666667 0.333333 0.074074 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 0.333333 0.666667 0.148148 (R,S,U): 0.000000 1.000000 0.000000 (R,S): 0.011582 0.460018 (U,U*,Error): 0.002451 0.002451 0.000000 (Ur,Ur*,Error): 0.211616 0.211616 0.000000 (Us,Us*,Error): 0.010656 0.010656 0.000000 (R,S): 0.293083 0.000508 (U,U*,Error): 0.000000 0.000000 0.000000 (Ur,Ur*,Error): 0.000000 0.000000 0.000000 (Us,Us*,Error): 0.000298 0.000298 0.000000 (R,S): 0.287849 0.414523 (U,U*,Error): 0.049461 0.049461 0.000000 (Ur,Ur*,Error): 0.171829 0.171829 0.000000 (Us,Us*,Error): 0.238640 0.238640 0.000000 (R,S): 0.535160 0.236043 (U,U*,Error): 0.029817 0.029817 0.000000 (Ur,Ur*,Error): 0.055716 0.055716 0.000000 (Us,Us*,Error): 0.252641 0.252641 0.000000 (R,S): 0.818204 0.100222 (U,U*,Error): 0.008218 0.008218 0.000000 (Ur,Ur*,Error): 0.010044 0.010044 0.000000 (Us,Us*,Error): 0.164003 0.164003 0.000000 Interpolate R^2 * S^1 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.000000 (R,S,U): 0.666667 0.000000 0.000000 (R,S,U): 1.000000 0.000000 0.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 0.333333 0.333333 0.037037 (R,S,U): 0.666667 0.333333 0.148148 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 0.333333 0.666667 0.074074 (R,S,U): 0.000000 1.000000 0.000000 (R,S): 0.178117 0.359635 (U,U*,Error): 0.011410 0.011410 0.000000 (Ur,Ur*,Error): 0.128114 0.128114 0.000000 (Us,Us*,Error): 0.031726 0.031726 0.000000 (R,S): 0.056705 0.521886 (U,U*,Error): 0.001678 0.001678 0.000000 (Ur,Ur*,Error): 0.059187 0.059187 0.000000 (Us,Us*,Error): 0.003215 0.003215 0.000000 (R,S): 0.335849 0.175669 (U,U*,Error): 0.019815 0.019815 0.000000 (Ur,Ur*,Error): 0.117997 0.117997 0.000000 (Us,Us*,Error): 0.112795 0.112795 0.000000 (R,S): 0.791053 0.094846 (U,U*,Error): 0.059352 0.059352 0.000000 (Ur,Ur*,Error): 0.150057 0.150057 0.000000 (Us,Us*,Error): 0.625765 0.625765 0.000000 (R,S): 0.324609 0.531532 (U,U*,Error): 0.056008 0.056008 0.000000 (Ur,Ur*,Error): 0.345080 0.345080 0.000000 (Us,Us*,Error): 0.105371 0.105371 0.000000 Interpolate R^3 * S^0 (R,S,U): 0.000000 0.000000 0.000000 (R,S,U): 0.333333 0.000000 0.037037 (R,S,U): 0.666667 0.000000 0.296296 (R,S,U): 1.000000 0.000000 1.000000 (R,S,U): 0.000000 0.333333 0.000000 (R,S,U): 0.333333 0.333333 0.037037 (R,S,U): 0.666667 0.333333 0.296296 (R,S,U): 0.000000 0.666667 0.000000 (R,S,U): 0.333333 0.666667 0.037037 (R,S,U): 0.000000 1.000000 0.000000 (R,S): 0.087868 0.895988 (U,U*,Error): 0.000678 0.000678 0.000000 (Ur,Ur*,Error): 0.023162 0.023162 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.254454 0.263733 (U,U*,Error): 0.016475 0.016475 0.000000 (Ur,Ur*,Error): 0.194240 0.194240 0.000000 (Us,Us*,Error): 0.000000 0.000000 0.000000 (R,S): 0.561861 0.184194 (U,U*,Error): 0.177373 0.177373 0.000000 (Ur,Ur*,Error): 0.947065 0.947065 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.597211 0.299937 (U,U*,Error): 0.213002 0.213002 0.000000 (Ur,Ur*,Error): 1.069984 1.069984 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 (R,S): 0.134123 0.212602 (U,U*,Error): 0.002413 0.002413 0.000000 (Ur,Ur*,Error): 0.053967 0.053967 0.000000 (Us,Us*,Error): 0.000000 -0.000000 0.000000 fem2d_pack_test13(): Test 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 -6.31439e-16 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 -5.66214e-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 7.99361e-15 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 : 6.21725e-15 -1.18031e-14 6.21725e-15 1.2837e-16 -5.90153e-15 4 : 0 0 0 18 0 5 : 81 -40.5 9 24.75 -40.5 6 :-7.10543e-15 2.69784e-14 -1.42109e-14 1.20737e-14 -6.7446e-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 :-1.24345e-14 1.51461e-28 4.44089e-15 -1.1352e-14 2.02338e-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.24345e-14 -6.66134e-15 5.55112e-16 4.44089e-16 -2.81025e-15 4 : 0 0 -4.5 0 0 5 : 20.25 -4.5 -5.5 9 -4.5 6 :-1.77636e-14 1.24345e-14 -7.74072e-15 1.06581e-14 3.9443e-29 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 :-7.10543e-15 -1.77636e-15 4.67527e-15 -8.88178e-15 4.4964e-15 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 : 1.81105e-15 4 : -0 5 : 1 6 :-2.91742e-15 7 : 0 8 : -4.5 9 : -4.5 10 :-2.89892e-16 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 0 -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 fem2d_pack_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 fem2d_pack_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 fem2d_pack_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 ==> XSI ETA 0.894942 0.071453 2.861336 1.285811 0.894942 0.071453 0.242487 0.053754 1.538728 1.215018 0.242487 0.053754 0.441722 0.013283 1.896727 1.053133 0.441722 0.013283 0.102809 0.803342 2.008959 4.213367 0.102809 0.803342 0.093371 0.307367 1.494108 2.229468 0.093371 0.307367 0.456058 0.101669 2.013785 1.406678 0.456058 0.101669 0.004610 0.667907 1.677128 3.671629 0.004610 0.667907 0.297347 0.062045 1.656739 1.248181 0.297347 0.062045 0.298244 0.046351 1.642839 1.185405 0.298244 0.046351 0.494572 0.238574 2.227718 1.954296 0.494572 0.238574 fem2d_pack_test16(): For an order 6 triangle, REFERENCE_TO_PHYSICAL_T6 maps a reference point to a physical point. XSI ETA ==> X Y 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 2.000000 0.000000 0.000000 1.000000 0.000000 4.000000 0.500000 0.000000 1.000000 0.000000 0.500000 0.500000 1.000000 1.000000 0.000000 0.500000 0.000000 2.000000 0.250000 0.750000 0.500000 2.250000 0.750000 0.250000 1.500000 0.250000 0.400000 0.100000 0.800000 0.240000 0.300000 0.200000 0.600000 0.560000 0.200000 0.300000 0.400000 0.960000 0.100000 0.400000 0.200000 1.440000 0.100000 0.100000 0.200000 0.360000 0.200000 0.200000 0.400000 0.640000 0.300000 0.300000 0.600000 0.840000 0.400000 0.400000 0.800000 0.960000 fem2d_pack_test17() S_L2NORM computes the L2 norm of a scalar function S(X) over a region (of any dimension), assuming: * that there is a set of finite element basis functions PSI; * that the region is decomposed into a number of elements whose areas are known; * that the integral is to be computed by a quadrature rule applied in the same way to each element; * that the value of the basis functions is given at every quadrature node in every element; Our example will have one spatial dimension. XMIN = 0.000000 XMAX = 1.000000 The number of intervals will be 40 The number of basis functions is 41 The number of quadrature points is 5 The computed L2 norm is 0.522156 The exact value is 0.522184 fem2d_pack_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.000000 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 1.000000 The R and S derivatives should sum to 0. dTdR sum dTdS sum 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 SHAPE_TEST(): Verify shape functions of type "Q8" Element order = 8 Basis function values at basis nodes should form the identity matrix. 1.000000 -0.000000 -0.000000 -0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 1.000000 -0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 0.000000 0.000000 -0.000000 1.000000 -0.000000 -0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 -0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 1.000000 -0.000000 -0.000000 0.000000 0.000000 0.000000 -0.000000 -0.000000 -0.000000 1.000000 -0.000000 0.000000 -0.000000 -0.000000 -0.000000 0.000000 0.000000 0.000000 1.000000 The R and S derivatives should sum to 0. dTdR sum dTdS sum 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 SHAPE_TEST(): Verify shape functions of type "Q9" Element order = 9 Basis function values at basis nodes should form the identity matrix. 1.000000 -0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 0.000000 0.000000 1.000000 -0.000000 -0.000000 -0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 1.000000 0.000000 -0.000000 -0.000000 -0.000000 -0.000000 0.000000 0.000000 -0.000000 -0.000000 1.000000 0.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 1.000000 0.000000 -0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 1.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 1.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 1.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 0.000000 -0.000000 1.000000 The R and S derivatives should sum to 0. dTdR sum dTdS sum 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 SHAPE_TEST(): Verify shape functions of type "Q12" Element order = 12 Basis function values at basis nodes should form the identity matrix. 1.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 1.000000 0.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 1.000000 -0.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 1.000000 -0.000000 0.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 1.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 -0.000000 1.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 -0.000000 1.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 -0.000000 -0.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 1.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 1.000000 The R and S derivatives should sum to 0. dTdR sum dTdS sum 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 SHAPE_TEST(): Verify shape functions of type "Q16" Element order = 16 Basis function values at basis nodes should form the identity matrix. 1.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 1.000000 0.000000 -0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 1.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 1.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 1.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 1.000000 0.000000 -0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 1.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 1.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 1.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 1.000000 0.000000 -0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 1.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 1.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 1.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 1.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 1.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 1.000000 The R and S derivatives should sum to 0. dTdR sum dTdS sum 0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 -0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 -0.000000 SHAPE_TEST(): Verify shape functions of type "QL" Element order = 6 Basis function values at basis nodes should form the identity matrix. 1.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 1.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 -0.000000 1.000000 0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 1.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 -0.000000 1.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 1.000000 The R and S derivatives should sum to 0. dTdR sum dTdS sum 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 SHAPE_TEST(): Verify shape functions of type "T3" Element order = 3 Basis function values at basis nodes should form the identity matrix. 1.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 1.000000 The R and S derivatives should sum to 0. dTdR sum dTdS sum 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 SHAPE_TEST(): Verify shape functions of type "T4" Element order = -1 Basis function values at basis nodes should form the identity matrix. The R and S derivatives should sum to 0. dTdR sum dTdS sum SHAPE_TEST(): Verify shape functions of type "T6" Element order = 6 Basis function values at basis nodes should form the identity matrix. 1.000000 -0.000000 -0.000000 0.000000 0.000000 0.000000 0.000000 1.000000 -0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 1.000000 The R and S derivatives should sum to 0. dTdR sum dTdS sum 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 SHAPE_TEST(): Verify shape functions of type "T10" Element order = 10 Basis function values at basis nodes should form the identity matrix. 1.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 -0.000000 -0.000000 -0.000000 0.000000 0.000000 1.000000 -0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 0.000000 0.000000 1.000000 -0.000000 0.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 0.000000 -0.000000 0.000000 1.000000 -0.000000 0.000000 0.000000 -0.000000 -0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 1.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 -0.000000 1.000000 0.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 1.000000 -0.000000 0.000000 -0.000000 -0.000000 0.000000 -0.000000 0.000000 -0.000000 0.000000 0.000000 1.000000 The R and S derivatives should sum to 0. dTdR sum dTdS sum 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.000000 fem2d_pack_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.000000 0.000000 -1.000000 2 0.382683 0.000000 -0.923880 3 0.270598 0.270598 -0.923880 4 0.000000 0.382683 -0.923880 5 -0.270598 0.270598 -0.923880 6 -0.382683 0.000000 -0.923880 7 -0.270598 -0.270598 -0.923880 8 -0.000000 -0.382683 -0.923880 9 0.270598 -0.270598 -0.923880 10 0.707107 0.000000 -0.707107 11 0.500000 0.500000 -0.707107 12 0.000000 0.707107 -0.707107 13 -0.500000 0.500000 -0.707107 14 -0.707107 0.000000 -0.707107 15 -0.500000 -0.500000 -0.707107 16 -0.000000 -0.707107 -0.707107 17 0.500000 -0.500000 -0.707107 18 0.923880 0.000000 -0.382683 19 0.653281 0.653281 -0.382683 20 0.000000 0.923880 -0.382683 21 -0.653281 0.653281 -0.382683 22 -0.923880 0.000000 -0.382683 23 -0.653281 -0.653281 -0.382683 24 -0.000000 -0.923880 -0.382683 25 0.653281 -0.653281 -0.382683 26 1.000000 0.000000 0.000000 27 0.707107 0.707107 0.000000 28 0.000000 1.000000 0.000000 29 -0.707107 0.707107 0.000000 30 -1.000000 0.000000 0.000000 31 -0.707107 -0.707107 0.000000 32 -0.000000 -1.000000 0.000000 33 0.707107 -0.707107 0.000000 34 0.923880 0.000000 0.382683 35 0.653281 0.653281 0.382683 36 0.000000 0.923880 0.382683 37 -0.653281 0.653281 0.382683 38 -0.923880 0.000000 0.382683 39 -0.653281 -0.653281 0.382683 40 -0.000000 -0.923880 0.382683 41 0.653281 -0.653281 0.382683 42 0.707107 0.000000 0.707107 43 0.500000 0.500000 0.707107 44 0.000000 0.707107 0.707107 45 -0.500000 0.500000 0.707107 46 -0.707107 0.000000 0.707107 47 -0.500000 -0.500000 0.707107 48 -0.000000 -0.707107 0.707107 49 0.500000 -0.500000 0.707107 50 0.382683 0.000000 0.923880 51 0.270598 0.270598 0.923880 52 0.000000 0.382683 0.923880 53 -0.270598 0.270598 0.923880 54 -0.382683 0.000000 0.923880 55 -0.270598 -0.270598 0.923880 56 -0.000000 -0.382683 0.923880 57 0.270598 -0.270598 0.923880 58 0.000000 0.000000 1.000000 fem2d_pack_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.000000 0.000000 -1.000000 2 0.382683 0.000000 -0.923880 3 0.191342 0.331414 -0.923880 4 -0.191342 0.331414 -0.923880 5 -0.382683 0.000000 -0.923880 6 -0.191342 -0.331414 -0.923880 7 0.191342 -0.331414 -0.923880 8 0.707107 0.000000 -0.707107 9 0.353553 0.612372 -0.707107 10 -0.353553 0.612372 -0.707107 11 -0.707107 0.000000 -0.707107 12 -0.353553 -0.612372 -0.707107 13 0.353553 -0.612372 -0.707107 14 0.923880 0.000000 -0.382683 15 0.461940 0.800103 -0.382683 16 -0.461940 0.800103 -0.382683 17 -0.923880 0.000000 -0.382683 18 -0.461940 -0.800103 -0.382683 19 0.461940 -0.800103 -0.382683 20 1.000000 0.000000 0.000000 21 0.500000 0.866025 0.000000 22 -0.500000 0.866025 0.000000 23 -1.000000 0.000000 0.000000 24 -0.500000 -0.866025 0.000000 25 0.500000 -0.866025 0.000000 26 0.923880 0.000000 0.382683 27 0.461940 0.800103 0.382683 28 -0.461940 0.800103 0.382683 29 -0.923880 0.000000 0.382683 30 -0.461940 -0.800103 0.382683 31 0.461940 -0.800103 0.382683 32 0.707107 0.000000 0.707107 33 0.353553 0.612372 0.707107 34 -0.353553 0.612372 0.707107 35 -0.707107 0.000000 0.707107 36 -0.353553 -0.612372 0.707107 37 0.353553 -0.612372 0.707107 38 0.382683 0.000000 0.923880 39 0.191342 0.331414 0.923880 40 -0.191342 0.331414 0.923880 41 -0.382683 0.000000 0.923880 42 -0.191342 -0.331414 0.923880 43 0.191342 -0.331414 0.923880 44 0.000000 0.000000 1.000000 fem2d_pack_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.000000 0.000000 -1.000000 2 0.500000 0.000000 -0.866025 3 0.250000 0.433013 -0.866025 4 -0.250000 0.433013 -0.866025 5 -0.500000 0.000000 -0.866025 6 -0.250000 -0.433013 -0.866025 7 0.250000 -0.433013 -0.866025 8 0.866025 0.000000 -0.500000 9 0.433013 0.750000 -0.500000 10 -0.433013 0.750000 -0.500000 11 -0.866025 0.000000 -0.500000 12 -0.433013 -0.750000 -0.500000 13 0.433013 -0.750000 -0.500000 14 1.000000 0.000000 0.000000 15 0.500000 0.866025 0.000000 16 -0.500000 0.866025 0.000000 17 -1.000000 0.000000 0.000000 18 -0.500000 -0.866025 0.000000 19 0.500000 -0.866025 0.000000 20 0.866025 0.000000 0.500000 21 0.433013 0.750000 0.500000 22 -0.433013 0.750000 0.500000 23 -0.866025 0.000000 0.500000 24 -0.433013 -0.750000 0.500000 25 0.433013 -0.750000 0.500000 26 0.500000 0.000000 0.866025 27 0.250000 0.433013 0.866025 28 -0.250000 0.433013 0.866025 29 -0.500000 0.000000 0.866025 30 -0.250000 -0.433013 0.866025 31 0.250000 -0.433013 0.866025 32 0.000000 0.000000 1.000000 fem2d_pack_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.000000 0.000000 -1.000000 2 0.382683 0.000000 -0.923880 3 0.270598 0.270598 -0.923880 4 0.000000 0.382683 -0.923880 5 -0.270598 0.270598 -0.923880 6 -0.382683 0.000000 -0.923880 7 -0.270598 -0.270598 -0.923880 8 -0.000000 -0.382683 -0.923880 9 0.270598 -0.270598 -0.923880 10 0.707107 0.000000 -0.707107 11 0.500000 0.500000 -0.707107 12 0.000000 0.707107 -0.707107 13 -0.500000 0.500000 -0.707107 14 -0.707107 0.000000 -0.707107 15 -0.500000 -0.500000 -0.707107 16 -0.000000 -0.707107 -0.707107 17 0.500000 -0.500000 -0.707107 18 0.923880 0.000000 -0.382683 19 0.653281 0.653281 -0.382683 20 0.000000 0.923880 -0.382683 21 -0.653281 0.653281 -0.382683 22 -0.923880 0.000000 -0.382683 23 -0.653281 -0.653281 -0.382683 24 -0.000000 -0.923880 -0.382683 25 0.653281 -0.653281 -0.382683 26 1.000000 0.000000 0.000000 27 0.707107 0.707107 0.000000 28 0.000000 1.000000 0.000000 29 -0.707107 0.707107 0.000000 30 -1.000000 0.000000 0.000000 31 -0.707107 -0.707107 0.000000 32 -0.000000 -1.000000 0.000000 33 0.707107 -0.707107 0.000000 34 0.923880 0.000000 0.382683 35 0.653281 0.653281 0.382683 36 0.000000 0.923880 0.382683 37 -0.653281 0.653281 0.382683 38 -0.923880 0.000000 0.382683 39 -0.653281 -0.653281 0.382683 40 -0.000000 -0.923880 0.382683 41 0.653281 -0.653281 0.382683 42 0.707107 0.000000 0.707107 43 0.500000 0.500000 0.707107 44 0.000000 0.707107 0.707107 45 -0.500000 0.500000 0.707107 46 -0.707107 0.000000 0.707107 47 -0.500000 -0.500000 0.707107 48 -0.000000 -0.707107 0.707107 49 0.500000 -0.500000 0.707107 50 0.382683 0.000000 0.923880 51 0.270598 0.270598 0.923880 52 0.000000 0.382683 0.923880 53 -0.270598 0.270598 0.923880 54 -0.382683 0.000000 0.923880 55 -0.270598 -0.270598 0.923880 56 -0.000000 -0.382683 0.923880 57 0.270598 -0.270598 0.923880 58 0.000000 0.000000 1.000000 fem2d_pack_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.000000 0.000000 -1.000000 2 0.382683 0.000000 -0.923880 3 0.191342 0.331414 -0.923880 4 -0.191342 0.331414 -0.923880 5 -0.382683 0.000000 -0.923880 6 -0.191342 -0.331414 -0.923880 7 0.191342 -0.331414 -0.923880 8 0.707107 0.000000 -0.707107 9 0.353553 0.612372 -0.707107 10 -0.353553 0.612372 -0.707107 11 -0.707107 0.000000 -0.707107 12 -0.353553 -0.612372 -0.707107 13 0.353553 -0.612372 -0.707107 14 0.923880 0.000000 -0.382683 15 0.461940 0.800103 -0.382683 16 -0.461940 0.800103 -0.382683 17 -0.923880 0.000000 -0.382683 18 -0.461940 -0.800103 -0.382683 19 0.461940 -0.800103 -0.382683 20 1.000000 0.000000 0.000000 21 0.500000 0.866025 0.000000 22 -0.500000 0.866025 0.000000 23 -1.000000 0.000000 0.000000 24 -0.500000 -0.866025 0.000000 25 0.500000 -0.866025 0.000000 26 0.923880 0.000000 0.382683 27 0.461940 0.800103 0.382683 28 -0.461940 0.800103 0.382683 29 -0.923880 0.000000 0.382683 30 -0.461940 -0.800103 0.382683 31 0.461940 -0.800103 0.382683 32 0.707107 0.000000 0.707107 33 0.353553 0.612372 0.707107 34 -0.353553 0.612372 0.707107 35 -0.707107 0.000000 0.707107 36 -0.353553 -0.612372 0.707107 37 0.353553 -0.612372 0.707107 38 0.382683 0.000000 0.923880 39 0.191342 0.331414 0.923880 40 -0.191342 0.331414 0.923880 41 -0.382683 0.000000 0.923880 42 -0.191342 -0.331414 0.923880 43 0.191342 -0.331414 0.923880 44 0.000000 0.000000 1.000000 fem2d_pack_test24 TRIANGLE_UNIT_SET sets up a quadrature in the unit triangle, A = 0.000000 B = 0.000000 Rule QUAD ERROR 1 1.000000 0.000000 2 1.000000 0.000000 3 1.000000 0.000000 4 1.000000 0.000000 5 1.000000 0.000000 6 1.000000 0.000000 7 1.000000 0.000000 8 1.000000 0.000000 9 1.000000 0.000000 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 1.000000 0.000000 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 0.000000 B = 1.000000 Rule QUAD ERROR 1 1.000000 0.000000 2 1.000000 0.000000 3 1.000000 0.000000 4 1.000000 0.000000 5 1.000000 0.000000 6 1.000000 0.000000 7 1.000000 0.000000 8 1.000000 0.000000 9 1.000000 0.000000 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 1.000000 0.000000 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 0.000000 B = 2.000000 Rule QUAD ERROR 1 0.666667 0.333333 2 2.000000 1.000000 3 1.000000 0.000000 4 1.000000 0.000000 5 1.000000 0.000000 6 1.000000 0.000000 7 1.000000 0.000000 8 1.000000 0.000000 9 1.000000 0.000000 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 1.000000 0.000000 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 0.000000 B = 3.000000 Rule QUAD ERROR 1 0.370370 0.629630 2 3.333333 2.333333 3 1.018519 0.018519 4 0.833333 0.166667 5 1.000000 0.000000 6 1.000000 0.000000 7 1.000000 0.000000 8 1.000000 0.000000 9 1.000000 0.000000 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 1.000000 0.000000 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 0.000000 B = 4.000000 Rule QUAD ERROR 1 0.185185 0.814815 2 5.000000 4.000000 3 0.995370 0.004630 4 0.625000 0.375000 5 0.933333 0.066667 6 0.958333 0.041667 7 0.958333 0.041667 8 1.000000 0.000000 9 1.000000 0.000000 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 1.083333 0.083333 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 0.000000 B = 5.000000 Rule QUAD ERROR 1 0.086420 0.913580 2 7.000000 6.000000 3 0.923611 0.076389 4 0.437500 0.562500 5 0.808889 0.191111 6 0.875000 0.125000 7 0.875000 0.125000 8 1.005274 0.005274 9 0.985185 0.014815 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 1.263889 0.263889 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 0.000000 B = 6.000000 Rule QUAD ERROR 1 0.038409 0.961591 2 9.333333 8.333333 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.012949 0.012949 9 0.948082 0.051918 10 0.995414 0.004586 11 0.995927 0.004073 12 1.000000 0.000000 13 1.000000 0.000000 14 1.533951 0.533951 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 0.000000 B = 7.000000 Rule QUAD ERROR 1 0.016461 0.983539 2 12.000000 11.000000 3 0.702418 0.297582 4 0.187500 0.812500 5 0.516101 0.483899 6 0.648333 0.351667 7 0.650926 0.349074 8 1.016902 0.016902 9 0.889877 0.110123 10 0.981708 0.018292 11 0.983203 0.016797 12 1.000198 0.000198 13 1.000000 0.000000 14 1.882407 0.882407 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 0.000000 B = 8.000000 Rule QUAD ERROR 1 0.006859 0.993141 2 15.000000 14.000000 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.012161 0.012161 9 0.816165 0.183835 10 0.956602 0.043398 11 0.959183 0.040817 12 1.000439 0.000439 13 0.999572 0.000428 14 2.299961 1.299961 15 0.999733 0.000267 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 0.000000 B = 9.000000 Rule QUAD ERROR 1 0.002794 0.997206 2 18.333333 17.333333 3 0.476896 0.523104 4 0.071615 0.928385 5 0.287142 0.712858 6 0.429942 0.570058 7 0.436368 0.563632 8 0.996255 0.003745 9 0.733569 0.266431 10 0.920084 0.079916 11 0.923540 0.076460 12 1.000086 0.000086 13 0.997886 0.002114 14 2.779903 1.779903 15 0.998637 0.001363 16 1.000000 0.000000 17 1.000091 0.000091 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 0.000000 B = 10.000000 Rule QUAD ERROR 1 0.001118 0.998882 2 22.000000 21.000000 3 0.381514 0.618486 4 0.042969 0.957031 5 0.207231 0.792769 6 0.339994 0.660006 7 0.347660 0.652340 8 0.968918 0.031082 9 0.647992 0.352008 10 0.873671 0.126329 11 0.877655 0.122345 12 0.998219 0.001781 13 0.993975 0.006025 14 3.317690 2.317690 15 0.995992 0.004008 16 1.000000 0.000000 17 1.000489 0.000489 18 0.999970 0.000030 19 1.000000 0.000000 20 1.000000 0.000000 A = 1.000000 B = 0.000000 Rule QUAD ERROR 1 1.000000 0.000000 2 1.000000 0.000000 3 1.000000 0.000000 4 1.000000 0.000000 5 1.000000 0.000000 6 1.000000 0.000000 7 1.000000 0.000000 8 1.000000 0.000000 9 1.000000 0.000000 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 1.000000 0.000000 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 1.000000 B = 1.000000 Rule QUAD ERROR 1 1.333333 0.333333 2 0.000000 1.000000 3 1.000000 0.000000 4 1.000000 0.000000 5 1.000000 0.000000 6 1.000000 0.000000 7 1.000000 0.000000 8 1.000000 0.000000 9 1.000000 0.000000 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 1.000000 0.000000 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 1.000000 B = 2.000000 Rule QUAD ERROR 1 1.111111 0.111111 2 0.000000 1.000000 3 0.972222 0.027778 4 1.250000 0.250000 5 1.000000 0.000000 6 1.000000 0.000000 7 1.000000 0.000000 8 1.000000 0.000000 9 1.000000 0.000000 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 1.000000 0.000000 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 1.000000 B = 3.000000 Rule QUAD ERROR 1 0.740741 0.259259 2 0.000000 1.000000 3 1.064815 0.064815 4 1.250000 0.250000 5 1.133333 0.133333 6 1.083333 0.083333 7 1.083333 0.083333 8 1.000000 0.000000 9 1.000000 0.000000 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 0.833333 0.166667 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 1.000000 B = 4.000000 Rule QUAD ERROR 1 0.432099 0.567901 2 0.000000 1.000000 3 1.174769 0.174769 4 1.093750 0.093750 5 1.244444 0.244444 6 1.166667 0.166667 7 1.166667 0.166667 8 0.986815 0.013185 9 1.037037 0.037037 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 0.631944 0.368056 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 1.000000 B = 5.000000 Rule QUAD ERROR 1 0.230453 0.769547 2 0.000000 1.000000 3 1.235082 0.235082 4 0.875000 0.125000 5 1.253570 0.253570 6 1.201667 0.201667 7 1.199074 0.199074 8 0.982247 0.017753 9 1.096494 0.096494 10 1.013757 0.013757 11 1.012218 0.012218 12 1.000000 0.000000 13 1.000000 0.000000 14 0.453704 0.546296 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 1.000000 B = 6.000000 Rule QUAD ERROR 1 0.115226 0.884774 2 0.000000 1.000000 3 1.230581 0.230581 4 0.656250 0.343750 5 1.166625 0.166625 6 1.178333 0.178333 7 1.173148 0.173148 8 0.999117 0.000883 9 1.151802 0.151802 10 1.043386 0.043386 11 1.040463 0.040463 12 0.999307 0.000693 13 1.000000 0.000000 14 0.314352 0.685648 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 1.000000 B = 7.000000 Rule QUAD ERROR 1 0.054870 0.945130 2 0.000000 1.000000 3 1.170910 0.170910 4 0.468750 0.531250 5 1.020816 0.020816 6 1.106250 0.106250 7 1.100694 0.100694 8 1.035863 0.035863 9 1.184724 0.184724 10 1.082131 0.082131 11 1.079284 0.079284 12 0.999233 0.000767 13 1.001712 0.001712 14 0.212191 0.787809 15 1.000567 0.000567 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 1.000000 B = 8.000000 Rule QUAD ERROR 1 0.025149 0.974851 2 0.000000 1.000000 3 1.073089 0.073089 4 0.322266 0.677734 5 0.852434 0.147566 6 1.001458 0.001458 7 0.998007 0.001993 8 1.083738 0.083738 9 1.187843 0.187843 10 1.120935 0.120935 11 1.119573 0.119573 12 1.002030 0.002030 13 1.007157 0.007157 14 0.140223 0.859777 15 1.002532 0.002532 16 1.000000 0.000000 17 0.999588 0.000412 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 1.000000 B = 9.000000 Rule QUAD ERROR 1 0.011177 0.988823 2 0.000000 1.000000 3 0.953802 0.046198 4 0.214844 0.785156 5 0.686695 0.313305 6 0.879682 0.120318 7 0.879906 0.120094 8 1.132942 0.132942 9 1.161458 0.161458 10 1.152148 0.152148 11 1.152968 0.152968 12 1.009418 0.009418 13 1.017444 0.017444 14 0.090967 0.909033 15 1.006584 0.006584 16 1.000000 0.000000 17 0.998106 0.001894 18 1.000151 0.000151 19 1.000000 0.000000 20 1.000000 0.000000 A = 2.000000 B = 0.000000 Rule QUAD ERROR 1 0.666667 0.333333 2 2.000000 1.000000 3 1.000000 0.000000 4 1.000000 0.000000 5 1.000000 0.000000 6 1.000000 0.000000 7 1.000000 0.000000 8 1.000000 0.000000 9 1.000000 0.000000 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 1.000000 0.000000 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 2.000000 B = 1.000000 Rule QUAD ERROR 1 1.111111 0.111111 2 0.000000 1.000000 3 0.972222 0.027778 4 1.250000 0.250000 5 1.000000 0.000000 6 1.000000 0.000000 7 1.000000 0.000000 8 1.000000 0.000000 9 1.000000 0.000000 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 1.000000 0.000000 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 2.000000 B = 2.000000 Rule QUAD ERROR 1 1.111111 0.111111 2 0.000000 1.000000 3 0.763889 0.236111 4 1.875000 0.875000 5 0.800000 0.200000 6 0.875000 0.125000 7 0.875000 0.125000 8 1.000000 0.000000 9 1.000000 0.000000 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 1.250000 0.250000 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 2.000000 B = 3.000000 Rule QUAD ERROR 1 0.864198 0.135802 2 0.000000 1.000000 3 0.729167 0.270833 4 2.187500 1.187500 5 0.808889 0.191111 6 0.875000 0.125000 7 0.875000 0.125000 8 1.005274 0.005274 9 0.985185 0.014815 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 1.263889 0.263889 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 2.000000 B = 4.000000 Rule QUAD ERROR 1 0.576132 0.423868 2 0.000000 1.000000 3 0.819187 0.180813 4 2.187500 1.187500 5 0.949926 0.050074 6 0.962500 0.037500 7 0.956019 0.043981 8 0.958434 0.041566 9 1.030123 0.030123 10 0.962963 0.037037 11 0.983741 0.016259 12 1.000000 0.000000 13 1.000000 0.000000 14 1.134259 0.134259 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 2.000000 B = 5.000000 Rule QUAD ERROR 1 0.345679 0.654321 2 0.000000 1.000000 3 0.937114 0.062886 4 1.968750 0.968750 5 1.080676 0.080676 6 1.067500 0.067500 7 1.040278 0.040278 8 0.881502 0.118498 9 1.140741 0.140741 10 0.909751 0.090249 11 0.967270 0.032730 12 1.000891 0.000891 13 1.000000 0.000000 14 0.943056 0.056944 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 2.000000 B = 6.000000 Rule QUAD ERROR 1 0.192044 0.807956 2 0.000000 1.000000 3 1.028485 0.028485 4 1.640625 0.640625 5 1.133495 0.133495 6 1.144792 0.144792 7 1.089699 0.089699 8 0.811392 0.188608 9 1.286387 0.286387 10 0.866501 0.133499 11 0.970273 0.029727 12 0.996224 0.003776 13 0.994117 0.005883 14 0.742670 0.257330 15 0.997732 0.002268 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 2.000000 B = 7.000000 Rule QUAD ERROR 1 0.100594 0.899406 2 0.000000 1.000000 3 1.074121 0.074121 4 1.289062 0.289062 5 1.103432 0.103432 6 1.175625 0.175625 7 1.095615 0.095615 8 0.769683 0.230317 9 1.433446 0.433446 10 0.846090 0.153910 11 1.001166 0.001166 12 0.981678 0.018322 13 0.979223 0.020777 14 0.560892 0.439108 15 0.990804 0.009196 16 1.000000 0.000000 17 1.000729 0.000729 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 2.000000 B = 8.000000 Rule QUAD ERROR 1 0.050297 0.949703 2 0.000000 1.000000 3 1.073286 0.073286 4 0.966797 0.033203 5 1.012554 0.012554 6 1.158724 0.158724 7 1.062637 0.062637 8 0.761318 0.238682 9 1.557612 0.557612 10 0.849523 0.150477 11 1.057002 0.057002 12 0.957945 0.042055 13 0.956754 0.043246 14 0.409352 0.590648 15 0.978158 0.021842 16 1.000000 0.000000 17 1.002323 0.002323 18 0.999302 0.000698 19 1.000000 0.000000 20 1.000000 0.000000 A = 3.000000 B = 0.000000 Rule QUAD ERROR 1 0.370370 0.629630 2 3.333333 2.333333 3 1.018519 0.018519 4 0.833333 0.166667 5 1.000000 0.000000 6 1.000000 0.000000 7 1.000000 0.000000 8 1.000000 0.000000 9 1.000000 0.000000 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 1.000000 0.000000 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 3.000000 B = 1.000000 Rule QUAD ERROR 1 0.740741 0.259259 2 0.000000 1.000000 3 1.064815 0.064815 4 1.250000 0.250000 5 1.133333 0.133333 6 1.083333 0.083333 7 1.083333 0.083333 8 1.000000 0.000000 9 1.000000 0.000000 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 0.833333 0.166667 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 3.000000 B = 2.000000 Rule QUAD ERROR 1 0.864198 0.135802 2 0.000000 1.000000 3 0.729167 0.270833 4 2.187500 1.187500 5 0.808889 0.191111 6 0.875000 0.125000 7 0.875000 0.125000 8 1.005274 0.005274 9 0.985185 0.014815 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 1.263889 0.263889 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 3.000000 B = 3.000000 Rule QUAD ERROR 1 0.768176 0.231824 2 0.000000 1.000000 3 0.516118 0.483882 4 2.916667 1.916667 5 0.594568 0.405432 6 0.738889 0.261111 7 0.756173 0.243827 8 1.063496 0.063496 9 0.917202 0.082798 10 1.051146 0.051146 11 1.012155 0.012155 12 1.000000 0.000000 13 1.000000 0.000000 14 1.512346 0.512346 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 3.000000 B = 4.000000 Rule QUAD ERROR 1 0.576132 0.423868 2 0.000000 1.000000 3 0.481610 0.518390 4 3.281250 2.281250 5 0.591526 0.408474 6 0.729167 0.270833 7 0.761574 0.238426 8 1.057370 0.057370 9 0.905679 0.094321 10 1.053666 0.053666 11 0.999136 0.000864 12 0.999505 0.000495 13 1.000000 0.000000 14 1.571759 0.571759 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 3.000000 B = 5.000000 Rule QUAD ERROR 1 0.384088 0.615912 2 0.000000 1.000000 3 0.544624 0.455376 4 3.281250 2.281250 5 0.694511 0.305489 6 0.802083 0.197917 7 0.818287 0.181713 8 0.972131 0.027869 9 0.992329 0.007671 10 0.992477 0.007523 11 0.961034 0.038966 12 1.011537 0.011537 13 1.011658 0.011658 14 1.485340 0.485340 15 1.005291 0.005291 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 3.000000 B = 6.000000 Rule QUAD ERROR 1 0.234720 0.765280 2 0.000000 1.000000 3 0.635854 0.364146 4 3.007812 2.007812 5 0.800594 0.199406 6 0.903681 0.096319 7 0.873050 0.126950 8 0.846318 0.153682 9 1.160570 0.160570 10 0.896388 0.103612 11 0.920984 0.079016 12 1.030392 0.030392 13 1.028932 0.028932 14 1.308749 0.308749 15 1.016209 0.016209 16 1.000000 0.000000 17 0.999369 0.000631 18 0.999999 0.000001 19 1.000000 0.000000 20 1.000000 0.000000 A = 3.000000 B = 7.000000 Rule QUAD ERROR 1 0.134126 0.865874 2 0.000000 1.000000 3 0.718176 0.281824 4 2.578125 1.578125 5 0.860170 0.139830 6 0.995347 0.004653 7 0.904171 0.095829 8 0.720497 0.279503 9 1.375005 0.375005 10 0.797080 0.202920 11 0.902218 0.097782 12 1.042706 0.042706 13 1.039343 0.039343 14 1.091607 0.091607 15 1.028142 0.028142 16 1.000000 0.000000 17 0.999861 0.000139 18 1.001881 0.001881 19 1.000000 0.000000 20 1.000000 0.000000 A = 4.000000 B = 0.000000 Rule QUAD ERROR 1 0.185185 0.814815 2 5.000000 4.000000 3 0.995370 0.004630 4 0.625000 0.375000 5 0.933333 0.066667 6 0.958333 0.041667 7 0.958333 0.041667 8 1.000000 0.000000 9 1.000000 0.000000 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 1.083333 0.083333 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 4.000000 B = 1.000000 Rule QUAD ERROR 1 0.432099 0.567901 2 0.000000 1.000000 3 1.174769 0.174769 4 1.093750 0.093750 5 1.244444 0.244444 6 1.166667 0.166667 7 1.166667 0.166667 8 0.986815 0.013185 9 1.037037 0.037037 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 0.631944 0.368056 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 4.000000 B = 2.000000 Rule QUAD ERROR 1 0.576132 0.423868 2 0.000000 1.000000 3 0.819187 0.180813 4 2.187500 1.187500 5 0.949926 0.050074 6 0.962500 0.037500 7 0.956019 0.043981 8 0.958434 0.041566 9 1.030123 0.030123 10 0.962963 0.037037 11 0.983741 0.016259 12 1.000000 0.000000 13 1.000000 0.000000 14 1.134259 0.134259 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 4.000000 B = 3.000000 Rule QUAD ERROR 1 0.576132 0.423868 2 0.000000 1.000000 3 0.481610 0.518390 4 3.281250 2.281250 5 0.591526 0.408474 6 0.729167 0.270833 7 0.761574 0.238426 8 1.057370 0.057370 9 0.905679 0.094321 10 1.053666 0.053666 11 0.999136 0.000864 12 0.999505 0.000495 13 1.000000 0.000000 14 1.571759 0.571759 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 4.000000 B = 4.000000 Rule QUAD ERROR 1 0.480110 0.519890 2 0.000000 1.000000 3 0.320698 0.679302 4 4.101562 3.101562 5 0.414538 0.585462 6 0.601563 0.398437 7 0.698785 0.301215 8 1.129907 0.129907 9 0.799671 0.200329 10 1.137746 0.137746 11 1.004338 0.004338 12 0.983592 0.016408 13 0.985455 0.014545 14 1.856674 0.856674 15 0.992063 0.007937 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 4.000000 B = 5.000000 Rule QUAD ERROR 1 0.352080 0.647920 2 0.000000 1.000000 3 0.293630 0.706370 4 4.511719 3.511719 5 0.402555 0.597445 6 0.585521 0.414479 7 0.715439 0.284561 8 1.106105 0.106105 9 0.796322 0.203678 10 1.142490 0.142490 11 0.974400 0.025600 12 0.982741 0.017259 13 0.985505 0.014495 14 1.963124 0.963124 15 0.986185 0.013815 16 1.000000 0.000000 17 1.000176 0.000176 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 4.000000 B = 6.000000 Rule QUAD ERROR 1 0.234720 0.765280 2 0.000000 1.000000 3 0.332597 0.667403 4 4.511719 3.511719 5 0.467461 0.532539 6 0.643003 0.356997 7 0.750509 0.249491 8 0.998828 0.001172 9 0.905044 0.094956 10 1.069173 0.069173 11 0.916307 0.083693 12 1.007858 0.007858 13 1.010243 0.010243 14 1.910312 0.910312 15 0.992343 0.007657 16 1.000000 0.000000 17 0.997408 0.002592 18 0.996691 0.003309 19 1.000000 0.000000 20 1.000000 0.000000 A = 5.000000 B = 0.000000 Rule QUAD ERROR 1 0.086420 0.913580 2 7.000000 6.000000 3 0.923611 0.076389 4 0.437500 0.562500 5 0.808889 0.191111 6 0.875000 0.125000 7 0.875000 0.125000 8 1.005274 0.005274 9 0.985185 0.014815 10 1.000000 0.000000 11 1.000000 0.000000 12 1.000000 0.000000 13 1.000000 0.000000 14 1.263889 0.263889 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 5.000000 B = 1.000000 Rule QUAD ERROR 1 0.230453 0.769547 2 0.000000 1.000000 3 1.235082 0.235082 4 0.875000 0.125000 5 1.253570 0.253570 6 1.201667 0.201667 7 1.199074 0.199074 8 0.982247 0.017753 9 1.096494 0.096494 10 1.013757 0.013757 11 1.012218 0.012218 12 1.000000 0.000000 13 1.000000 0.000000 14 0.453704 0.546296 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 5.000000 B = 2.000000 Rule QUAD ERROR 1 0.345679 0.654321 2 0.000000 1.000000 3 0.937114 0.062886 4 1.968750 0.968750 5 1.080676 0.080676 6 1.067500 0.067500 7 1.040278 0.040278 8 0.881502 0.118498 9 1.140741 0.140741 10 0.909751 0.090249 11 0.967270 0.032730 12 1.000891 0.000891 13 1.000000 0.000000 14 0.943056 0.056944 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 5.000000 B = 3.000000 Rule QUAD ERROR 1 0.384088 0.615912 2 0.000000 1.000000 3 0.544624 0.455376 4 3.281250 2.281250 5 0.694511 0.305489 6 0.802083 0.197917 7 0.818287 0.181713 8 0.972131 0.027869 9 0.992329 0.007671 10 0.992477 0.007523 11 0.961034 0.038966 12 1.011537 0.011537 13 1.011658 0.011658 14 1.485340 0.485340 15 1.007937 0.007937 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 5.000000 B = 4.000000 Rule QUAD ERROR 1 0.352080 0.647920 2 0.000000 1.000000 3 0.293630 0.706370 4 4.511719 3.511719 5 0.402555 0.597445 6 0.585521 0.414479 7 0.715439 0.284561 8 1.106105 0.106105 9 0.796322 0.203678 10 1.142490 0.142490 11 0.974400 0.025600 12 0.982741 0.017259 13 0.985505 0.014495 14 1.963124 0.963124 15 0.999118 0.000882 16 1.000000 0.000000 17 1.000176 0.000176 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 5.000000 B = 5.000000 Rule QUAD ERROR 1 0.281664 0.718336 2 0.000000 1.000000 3 0.187868 0.812132 4 5.414062 4.414062 5 0.273552 0.726448 6 0.473229 0.526771 7 0.710487 0.289513 8 1.166185 0.166185 9 0.679036 0.320964 10 1.239579 0.239579 11 0.968337 0.031663 12 0.950326 0.049674 13 0.955574 0.044426 14 2.292374 1.292374 15 0.977562 0.022438 16 1.000000 0.000000 17 1.003019 0.003019 18 1.003977 0.003977 19 1.000000 0.000000 20 1.000000 0.000000 A = 6.000000 B = 0.000000 Rule QUAD ERROR 1 0.038409 0.961591 2 9.333333 8.333333 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.012949 0.012949 9 0.948082 0.051918 10 0.995414 0.004586 11 0.995927 0.004073 12 1.000000 0.000000 13 1.000000 0.000000 14 1.533951 0.533951 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 6.000000 B = 1.000000 Rule QUAD ERROR 1 0.115226 0.884774 2 0.000000 1.000000 3 1.230581 0.230581 4 0.656250 0.343750 5 1.166625 0.166625 6 1.178333 0.178333 7 1.173148 0.173148 8 0.999117 0.000883 9 1.151802 0.151802 10 1.043386 0.043386 11 1.040463 0.040463 12 0.999307 0.000693 13 1.000000 0.000000 14 0.314352 0.685648 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 6.000000 B = 2.000000 Rule QUAD ERROR 1 0.192044 0.807956 2 0.000000 1.000000 3 1.028485 0.028485 4 1.640625 0.640625 5 1.133495 0.133495 6 1.144792 0.144792 7 1.089699 0.089699 8 0.811392 0.188608 9 1.286387 0.286387 10 0.866501 0.133499 11 0.970273 0.029727 12 0.996224 0.003776 13 0.994117 0.005883 14 0.742670 0.257330 15 0.994709 0.005291 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 6.000000 B = 3.000000 Rule QUAD ERROR 1 0.234720 0.765280 2 0.000000 1.000000 3 0.635854 0.364146 4 3.007812 2.007812 5 0.800594 0.199406 6 0.903681 0.096319 7 0.873050 0.126950 8 0.846318 0.153682 9 1.160570 0.160570 10 0.896388 0.103612 11 0.920984 0.079016 12 1.030392 0.030392 13 1.028932 0.028932 14 1.308749 0.308749 15 1.015285 0.015285 16 1.000000 0.000000 17 0.999369 0.000631 18 0.999999 0.000001 19 1.000000 0.000000 20 1.000000 0.000000 A = 6.000000 B = 4.000000 Rule QUAD ERROR 1 0.234720 0.765280 2 0.000000 1.000000 3 0.332597 0.667403 4 4.511719 3.511719 5 0.467461 0.532539 6 0.643003 0.356997 7 0.750509 0.249491 8 0.998828 0.001172 9 0.905044 0.094956 10 1.069173 0.069173 11 0.916307 0.083693 12 1.007858 0.007858 13 1.010243 0.010243 14 1.910312 0.910312 15 1.020674 0.020674 16 1.000000 0.000000 17 0.997408 0.002592 18 0.996691 0.003309 19 1.000000 0.000000 20 1.000000 0.000000 A = 7.000000 B = 0.000000 Rule QUAD ERROR 1 0.016461 0.983539 2 12.000000 11.000000 3 0.702418 0.297582 4 0.187500 0.812500 5 0.516101 0.483899 6 0.648333 0.351667 7 0.650926 0.349074 8 1.016902 0.016902 9 0.889877 0.110123 10 0.981708 0.018292 11 0.983203 0.016797 12 1.000198 0.000198 13 1.000000 0.000000 14 1.882407 0.882407 15 1.000000 0.000000 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 7.000000 B = 1.000000 Rule QUAD ERROR 1 0.054870 0.945130 2 0.000000 1.000000 3 1.170910 0.170910 4 0.468750 0.531250 5 1.020816 0.020816 6 1.106250 0.106250 7 1.100694 0.100694 8 1.035863 0.035863 9 1.184724 0.184724 10 1.082131 0.082131 11 1.079284 0.079284 12 0.999233 0.000767 13 1.001712 0.001712 14 0.212191 0.787809 15 1.002268 0.002268 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 7.000000 B = 2.000000 Rule QUAD ERROR 1 0.100594 0.899406 2 0.000000 1.000000 3 1.074121 0.074121 4 1.289062 0.289062 5 1.103432 0.103432 6 1.175625 0.175625 7 1.095615 0.095615 8 0.769683 0.230317 9 1.433446 0.433446 10 0.846090 0.153910 11 1.001166 0.001166 12 0.981678 0.018322 13 0.979223 0.020777 14 0.560892 0.439108 15 0.982951 0.017049 16 1.000000 0.000000 17 1.000729 0.000729 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 7.000000 B = 3.000000 Rule QUAD ERROR 1 0.134126 0.865874 2 0.000000 1.000000 3 0.718176 0.281824 4 2.578125 1.578125 5 0.860170 0.139830 6 0.995347 0.004653 7 0.904171 0.095829 8 0.720497 0.279503 9 1.375005 0.375005 10 0.797080 0.202920 11 0.902218 0.097782 12 1.042706 0.042706 13 1.039343 0.039343 14 1.091607 0.091607 15 1.011436 0.011436 16 1.000000 0.000000 17 0.999861 0.000139 18 1.001881 0.001881 19 1.000000 0.000000 20 1.000000 0.000000 A = 8.000000 B = 0.000000 Rule QUAD ERROR 1 0.006859 0.993141 2 15.000000 14.000000 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.012161 0.012161 9 0.816165 0.183835 10 0.956602 0.043398 11 0.959183 0.040817 12 1.000439 0.000439 13 0.999572 0.000428 14 2.299961 1.299961 15 0.999433 0.000567 16 1.000000 0.000000 17 1.000000 0.000000 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 8.000000 B = 1.000000 Rule QUAD ERROR 1 0.025149 0.974851 2 0.000000 1.000000 3 1.073089 0.073089 4 0.322266 0.677734 5 0.852434 0.147566 6 1.001458 0.001458 7 0.998007 0.001993 8 1.083738 0.083738 9 1.187843 0.187843 10 1.120935 0.120935 11 1.119573 0.119573 12 1.002030 0.002030 13 1.007157 0.007157 14 0.140223 0.859777 15 1.009511 0.009511 16 1.000000 0.000000 17 0.999588 0.000412 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 8.000000 B = 2.000000 Rule QUAD ERROR 1 0.050297 0.949703 2 0.000000 1.000000 3 1.073286 0.073286 4 0.966797 0.033203 5 1.012554 0.012554 6 1.158724 0.158724 7 1.062637 0.062637 8 0.761318 0.238682 9 1.557612 0.557612 10 0.849523 0.150477 11 1.057002 0.057002 12 0.957945 0.042055 13 0.956754 0.043246 14 0.409352 0.590648 15 0.968709 0.031291 16 1.000000 0.000000 17 1.002323 0.002323 18 0.999302 0.000698 19 1.000000 0.000000 20 1.000000 0.000000 A = 9.000000 B = 0.000000 Rule QUAD ERROR 1 0.002794 0.997206 2 18.333333 17.333333 3 0.476896 0.523104 4 0.071615 0.928385 5 0.287142 0.712858 6 0.429942 0.570058 7 0.436368 0.563632 8 0.996255 0.003745 9 0.733569 0.266431 10 0.920084 0.079916 11 0.923540 0.076460 12 1.000086 0.000086 13 0.997886 0.002114 14 2.779903 1.779903 15 0.997194 0.002806 16 1.000000 0.000000 17 1.000091 0.000091 18 1.000000 0.000000 19 1.000000 0.000000 20 1.000000 0.000000 A = 9.000000 B = 1.000000 Rule QUAD ERROR 1 0.011177 0.988823 2 0.000000 1.000000 3 0.953802 0.046198 4 0.214844 0.785156 5 0.686695 0.313305 6 0.879682 0.120318 7 0.879906 0.120094 8 1.132942 0.132942 9 1.161458 0.161458 10 1.152148 0.152148 11 1.152968 0.152968 12 1.009418 0.009418 13 1.017444 0.017444 14 0.090967 0.909033 15 1.023112 0.023112 16 1.000000 0.000000 17 0.998106 0.001894 18 1.000151 0.000151 19 1.000000 0.000000 20 1.000000 0.000000 A = 10.000000 B = 0.000000 Rule QUAD ERROR 1 0.001118 0.998882 2 22.000000 21.000000 3 0.381514 0.618486 4 0.042969 0.957031 5 0.207231 0.792769 6 0.339994 0.660006 7 0.347660 0.652340 8 0.968918 0.031082 9 0.647992 0.352008 10 0.873671 0.126329 11 0.877655 0.122345 12 0.998219 0.001781 13 0.993975 0.006025 14 3.317690 2.317690 15 0.992010 0.007990 16 1.000000 0.000000 17 1.000489 0.000489 18 0.999970 0.000030 19 1.000000 0.000000 20 1.000000 0.000000 grid_t3_neighbor_test(): MATLAB version Test grid_t3_neighbor() Number of elements in X direction = 3 Number of elements in Y direction = 2 Element order = 3 Nodes in grid = 12 GRID_PRINT: Element -> Node table. Element order = 3 Number of elements = 12 # 1 2 3 1 1 2 5 2 6 5 2 3 2 3 6 4 7 6 3 5 3 4 7 6 8 7 4 7 5 6 9 8 10 9 6 9 6 7 10 10 11 10 7 11 7 8 11 12 12 11 8 The grid width is 4 Element neighbor array: 1: 2 -1 -1 2: 1 3 7 3: 4 2 -1 4: 3 5 9 5: 6 4 -1 6: 5 -1 11 7: 8 -1 2 8: 7 9 -1 9: 10 8 4 10: 9 11 -1 11: 12 10 6 12: 11 -1 -1 fem2d_pack_test(): Normal end of execution. 07-Jan-2022 20:01:56