1----------------------------- GeomRefTetrahedron test ---------------------------------- 2geometric reference element : shape tetrahedron, 4 vertices, 4 sides, 6 sideOfSides 3 * vertex 1 -> (1 ,0 ,0) 4 * vertex 2 -> (0 ,1 ,0) 5 * vertex 3 -> (0 ,0 ,1) 6 * vertex 4 -> (0 ,0 ,0) 7 * centroid : (0.25 ,0.25 ,0.25) 8 * measure : 0.166666666667 9 * side 1 -> shape triangle, vertices : 4 3 2 10 * side 2 -> shape triangle, vertices : 4 1 3 11 * side 3 -> shape triangle, vertices : 4 2 1 12 * side 4 -> shape triangle, vertices : 1 2 3 13 * side of side 1 -> [ 1 4 ] 14 * side of side 2 -> [ 2 4 ] 15 * side of side 3 -> [ 3 4 ] 16 * side of side 4 -> [ 2 3 ] 17 * side of side 5 -> [ 1 3 ] 18 * side of side 6 -> [ 1 2 ] 19 * numbers on side1 -> [ -3 -4 2 ] 20 * numbers on side2 -> [ -1 5 3 ] 21 * numbers on side3 -> [ -2 -6 1 ] 22 * numbers on side4 -> [ 6 4 -5 ] 23sideWithVertices(4,3,2) = 1 24sideWithVertices(4,1,3) = 2 25sideWithVertices(4,2,1) = 3 26sideWithVertices(1,2,3) = 4 27sideWithVertices(1,4) = 1 28sideWithVertices(2,4) = 2 29sideWithVertices(3,4) = 3 30sideWithVertices(2,3) = 4 31sideWithVertices(1,3) = 5 32sideWithVertices(1,2) = 6 33 34----------------------------- interpolation test ---------------------------------- 35reference element tetrahedron_Lagrange_1 36- finite element interpolation : Lagrange of order 1 conforming in H1 37- geometric reference element : shape tetrahedron, 4 vertices, 4 sides, 6 sideOfSides 38- reference DoFs : 4 39 * 1 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 1, index 1, coords = [ 1 0 0 ] 40 * 2 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 2, index 1, coords = [ 0 1 0 ] 41 * 3 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 3, index 1, coords = [ 0 0 1 ] 42 * 4 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 4, index 1, coords = [ 0 0 0 ] 43- Dof numbers on side of sides (RefElement tetrahedron_Lagrange_1) 44 . side of side 1 -> 1 4 45 . side of side 2 -> 2 4 46 . side of side 3 -> 3 4 47 . side of side 4 -> 2 3 48 . side of side 5 -> 1 3 49 . side of side 6 -> 1 2 50- Dof numbers on sides (RefElement tetrahedron_Lagrange_1) 51 . side 1 -> 4 3 2 52 . side 2 -> 4 1 3 53 . side 3 -> 4 2 1 54 . side 4 -> 1 2 3 55- Side of side shapes : 6 56segment_Lagrange_1, segment_Lagrange_1, segment_Lagrange_1, segment_Lagrange_1, segment_Lagrange_1, segment_Lagrange_1 57- Side shapes : 4 58triangle_Lagrange_1, triangle_Lagrange_1, triangle_Lagrange_1, triangle_Lagrange_1 59 60reference element tetrahedron_Lagrange_2 61- finite element interpolation : Lagrange of order 2 conforming in H1 62- geometric reference element : shape tetrahedron, 4 vertices, 4 sides, 6 sideOfSides 63- reference DoFs : 10 64 * 1 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 1, index 1, coords = [ 1 0 0 ] 65 * 2 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 2, index 1, coords = [ 0 1 0 ] 66 * 3 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 3, index 1, coords = [ 0 0 1 ] 67 * 4 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 4, index 1, coords = [ 0 0 0 ] 68 * 5 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 1, index 1, coords = [ 0.5 0 0 ] 69 * 6 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 2, index 1, coords = [ 0 0.5 0 ] 70 * 7 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 3, index 1, coords = [ 0 0 0.5 ] 71 * 8 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 4, index 1, coords = [ 0 0.5 0.5 ] 72 * 9 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 5, index 1, coords = [ 0.5 0 0.5 ] 73 * 10 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 6, index 1, coords = [ 0.5 0.5 0 ] 74- Dof numbers on side of sides (RefElement tetrahedron_Lagrange_2) 75 . side of side 1 -> 1 4 5 76 . side of side 2 -> 2 4 6 77 . side of side 3 -> 3 4 7 78 . side of side 4 -> 2 3 8 79 . side of side 5 -> 1 3 9 80 . side of side 6 -> 1 2 10 81- Dof numbers on sides (RefElement tetrahedron_Lagrange_2) 82 . side 1 -> 4 3 2 7 8 6 83 . side 2 -> 4 1 3 5 9 7 84 . side 3 -> 4 2 1 6 10 5 85 . side 4 -> 1 2 3 10 8 9 86- Side of side shapes : 6 87segment_Lagrange_2, segment_Lagrange_2, segment_Lagrange_2, segment_Lagrange_2, segment_Lagrange_2, segment_Lagrange_2 88- Side shapes : 4 89triangle_Lagrange_2, triangle_Lagrange_2, triangle_Lagrange_2, triangle_Lagrange_2 90 91reference element tetrahedron_Lagrange_3 92- finite element interpolation : Lagrange of order 3 conforming in H1 93- geometric reference element : shape tetrahedron, 4 vertices, 4 sides, 6 sideOfSides 94- reference DoFs : 20 95 * 1 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 1, index 1, coords = [ 1 0 0 ] 96 * 2 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 2, index 1, coords = [ 0 1 0 ] 97 * 3 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 3, index 1, coords = [ 0 0 1 ] 98 * 4 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 4, index 1, coords = [ 0 0 0 ] 99 * 5 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 1, index 1, coords = [ 0.66667 0 0 ] 100 * 6 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 2, index 1, coords = [ 0 0.66667 0 ] 101 * 7 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 3, index 1, coords = [ 0 0 0.66667 ] 102 * 8 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 4, index 1, coords = [ 0 0.66667 0.33333 ] 103 * 9 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 5, index 1, coords = [ 0.66667 0 0.33333 ] 104 * 10 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 6, index 1, coords = [ 0.66667 0.33333 0 ] 105 * 11 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 1, index 2, coords = [ 0.33333 0 0 ] 106 * 12 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 2, index 2, coords = [ 0 0.33333 0 ] 107 * 13 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 3, index 2, coords = [ 0 0 0.33333 ] 108 * 14 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 4, index 2, coords = [ 0 0.33333 0.66667 ] 109 * 15 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 5, index 2, coords = [ 0.33333 0 0.66667 ] 110 * 16 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 6, index 2, coords = [ 0.33333 0.66667 0 ] 111 * 17 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 1, index 1, coords = [ 0 0.33333 0.33333 ] 112 * 18 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 2, index 1, coords = [ 0.33333 0 0.33333 ] 113 * 19 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 3, index 1, coords = [ 0.33333 0.33333 0 ] 114 * 20 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 4, index 1, coords = [ 0.33333 0.33333 0.33333 ] 115- Dof numbers on side of sides (RefElement tetrahedron_Lagrange_3) 116 . side of side 1 -> 1 4 5 11 117 . side of side 2 -> 2 4 6 12 118 . side of side 3 -> 3 4 7 13 119 . side of side 4 -> 2 3 8 14 120 . side of side 5 -> 1 3 9 15 121 . side of side 6 -> 1 2 10 16 122- Dof numbers on sides (RefElement tetrahedron_Lagrange_3) 123 . side 1 -> 4 3 2 13 14 6 7 8 12 17 124 . side 2 -> 4 1 3 11 9 7 5 15 13 18 125 . side 3 -> 4 2 1 12 16 5 6 10 11 19 126 . side 4 -> 1 2 3 10 8 15 16 14 9 20 127- Side of side shapes : 6 128segment_Lagrange_3, segment_Lagrange_3, segment_Lagrange_3, segment_Lagrange_3, segment_Lagrange_3, segment_Lagrange_3 129- Side shapes : 4 130triangle_Lagrange_3, triangle_Lagrange_3, triangle_Lagrange_3, triangle_Lagrange_3 131 132reference element tetrahedron_Lagrange_4 133- finite element interpolation : Lagrange of order 4 conforming in H1 134- geometric reference element : shape tetrahedron, 4 vertices, 4 sides, 6 sideOfSides 135- reference DoFs : 35 136 * 1 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 1, index 1, coords = [ 1 0 0 ] 137 * 2 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 2, index 1, coords = [ 0 1 0 ] 138 * 3 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 3, index 1, coords = [ 0 0 1 ] 139 * 4 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 4, index 1, coords = [ 0 0 0 ] 140 * 5 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 1, index 1, coords = [ 0.75 0 0 ] 141 * 6 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 2, index 1, coords = [ 0 0.75 0 ] 142 * 7 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 3, index 1, coords = [ 0 0 0.75 ] 143 * 8 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 4, index 1, coords = [ 0 0.75 0.25 ] 144 * 9 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 5, index 1, coords = [ 0.75 0 0.25 ] 145 * 10 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 6, index 1, coords = [ 0.75 0.25 0 ] 146 * 11 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 1, index 2, coords = [ 0.5 0 0 ] 147 * 12 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 2, index 2, coords = [ 0 0.5 0 ] 148 * 13 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 3, index 2, coords = [ 0 0 0.5 ] 149 * 14 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 4, index 2, coords = [ 0 0.5 0.5 ] 150 * 15 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 5, index 2, coords = [ 0.5 0 0.5 ] 151 * 16 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 6, index 2, coords = [ 0.5 0.5 0 ] 152 * 17 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 1, index 3, coords = [ 0.25 0 0 ] 153 * 18 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 2, index 3, coords = [ 0 0.25 0 ] 154 * 19 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 3, index 3, coords = [ 0 0 0.25 ] 155 * 20 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 4, index 3, coords = [ 0 0.25 0.75 ] 156 * 21 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 5, index 3, coords = [ 0.25 0 0.75 ] 157 * 22 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 6, index 3, coords = [ 0.25 0.75 0 ] 158 * 23 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 1, index 1, coords = [ 0 0.25 0.25 ] 159 * 24 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 1, index 2, coords = [ 0 0.25 0.5 ] 160 * 25 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 1, index 3, coords = [ 0 0.5 0.25 ] 161 * 26 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 2, index 1, coords = [ 0.25 0 0.25 ] 162 * 27 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 2, index 2, coords = [ 0.5 0 0.25 ] 163 * 28 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 2, index 3, coords = [ 0.25 0 0.5 ] 164 * 29 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 3, index 1, coords = [ 0.25 0.25 0 ] 165 * 30 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 3, index 2, coords = [ 0.25 0.5 0 ] 166 * 31 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 3, index 3, coords = [ 0.5 0.25 0 ] 167 * 32 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 4, index 1, coords = [ 0.5 0.25 0.25 ] 168 * 33 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 4, index 2, coords = [ 0.25 0.5 0.25 ] 169 * 34 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 4, index 3, coords = [ 0.25 0.25 0.5 ] 170 * 35 -> reference DoF : nodal value, sharable? no, support : dim 0, on element 0, index 1, coords = [ 0.25 0.25 0.25 ] 171- Dof numbers on side of sides (RefElement tetrahedron_Lagrange_4) 172 . side of side 1 -> 1 4 5 11 17 173 . side of side 2 -> 2 4 6 12 18 174 . side of side 3 -> 3 4 7 13 19 175 . side of side 4 -> 2 3 8 14 20 176 . side of side 5 -> 1 3 9 15 21 177 . side of side 6 -> 1 2 10 16 22 178- Dof numbers on sides (RefElement tetrahedron_Lagrange_4) 179 . side 1 -> 4 3 2 19 20 6 13 14 12 7 8 18 23 24 25 180 . side 2 -> 4 1 3 17 9 7 11 15 13 5 21 19 26 27 28 181 . side 3 -> 4 2 1 18 22 5 12 16 11 6 10 17 29 30 31 182 . side 4 -> 1 2 3 10 8 21 16 14 15 22 20 9 32 33 34 183- Side of side shapes : 6 184segment_Lagrange_4, segment_Lagrange_4, segment_Lagrange_4, segment_Lagrange_4, segment_Lagrange_4, segment_Lagrange_4 185- Side shapes : 4 186triangle_Lagrange_4, triangle_Lagrange_4, triangle_Lagrange_4, triangle_Lagrange_4 187 188reference element tetrahedron_Lagrange_5 189- finite element interpolation : Lagrange of order 5 conforming in H1 190- geometric reference element : shape tetrahedron, 4 vertices, 4 sides, 6 sideOfSides 191- reference DoFs : 56 192 * 1 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 1, index 1, coords = [ 1 0 0 ] 193 * 2 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 2, index 1, coords = [ 0 1 0 ] 194 * 3 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 3, index 1, coords = [ 0 0 1 ] 195 * 4 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 4, index 1, coords = [ 0 0 0 ] 196 * 5 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 1, index 1, coords = [ 0.8 0 0 ] 197 * 6 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 2, index 1, coords = [ 0 0.8 0 ] 198 * 7 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 3, index 1, coords = [ 0 0 0.8 ] 199 * 8 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 4, index 1, coords = [ 0 0.8 0.2 ] 200 * 9 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 5, index 1, coords = [ 0.8 0 0.2 ] 201 * 10 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 6, index 1, coords = [ 0.8 0.2 0 ] 202 * 11 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 1, index 2, coords = [ 0.6 0 0 ] 203 * 12 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 2, index 2, coords = [ 0 0.6 0 ] 204 * 13 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 3, index 2, coords = [ 0 0 0.6 ] 205 * 14 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 4, index 2, coords = [ 0 0.6 0.4 ] 206 * 15 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 5, index 2, coords = [ 0.6 0 0.4 ] 207 * 16 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 6, index 2, coords = [ 0.6 0.4 0 ] 208 * 17 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 1, index 3, coords = [ 0.4 0 0 ] 209 * 18 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 2, index 3, coords = [ 0 0.4 0 ] 210 * 19 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 3, index 3, coords = [ 0 0 0.4 ] 211 * 20 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 4, index 3, coords = [ 0 0.4 0.6 ] 212 * 21 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 5, index 3, coords = [ 0.4 0 0.6 ] 213 * 22 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 6, index 3, coords = [ 0.4 0.6 0 ] 214 * 23 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 1, index 4, coords = [ 0.2 0 0 ] 215 * 24 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 2, index 4, coords = [ 0 0.2 0 ] 216 * 25 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 3, index 4, coords = [ 0 0 0.2 ] 217 * 26 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 4, index 4, coords = [ 0 0.2 0.8 ] 218 * 27 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 5, index 4, coords = [ 0.2 0 0.8 ] 219 * 28 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 6, index 4, coords = [ 0.2 0.8 0 ] 220 * 29 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 1, index 1, coords = [ 0 0.2 0.2 ] 221 * 30 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 1, index 2, coords = [ 0 0.2 0.6 ] 222 * 31 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 1, index 3, coords = [ 0 0.6 0.2 ] 223 * 32 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 1, index 4, coords = [ 0 0.2 0.4 ] 224 * 33 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 1, index 5, coords = [ 0 0.4 0.4 ] 225 * 34 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 1, index 6, coords = [ 0 0.4 0.2 ] 226 * 35 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 2, index 1, coords = [ 0.2 0 0.2 ] 227 * 36 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 2, index 2, coords = [ 0.6 0 0.2 ] 228 * 37 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 2, index 3, coords = [ 0.2 0 0.6 ] 229 * 38 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 2, index 4, coords = [ 0.4 0 0.2 ] 230 * 39 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 2, index 5, coords = [ 0.4 0 0.4 ] 231 * 40 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 2, index 6, coords = [ 0.2 0 0.4 ] 232 * 41 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 3, index 1, coords = [ 0.2 0.2 0 ] 233 * 42 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 3, index 2, coords = [ 0.2 0.6 0 ] 234 * 43 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 3, index 3, coords = [ 0.6 0.2 0 ] 235 * 44 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 3, index 4, coords = [ 0.2 0.4 0 ] 236 * 45 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 3, index 5, coords = [ 0.4 0.4 0 ] 237 * 46 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 3, index 6, coords = [ 0.4 0.2 0 ] 238 * 47 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 4, index 1, coords = [ 0.6 0.2 0.2 ] 239 * 48 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 4, index 2, coords = [ 0.2 0.6 0.2 ] 240 * 49 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 4, index 3, coords = [ 0.2 0.2 0.6 ] 241 * 50 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 4, index 4, coords = [ 0.4 0.4 0.2 ] 242 * 51 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 4, index 5, coords = [ 0.2 0.4 0.4 ] 243 * 52 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 4, index 6, coords = [ 0.4 0.2 0.4 ] 244 * 53 -> reference DoF : nodal value, sharable? no, support : dim 0, on element 0, index 1, coords = [ 0.4 0.2 0.2 ] 245 * 54 -> reference DoF : nodal value, sharable? no, support : dim 0, on element 0, index 2, coords = [ 0.2 0.4 0.2 ] 246 * 55 -> reference DoF : nodal value, sharable? no, support : dim 0, on element 0, index 3, coords = [ 0.2 0.2 0.4 ] 247 * 56 -> reference DoF : nodal value, sharable? no, support : dim 0, on element 0, index 4, coords = [ 0.2 0.2 0.2 ] 248- Dof numbers on side of sides (RefElement tetrahedron_Lagrange_5) 249 . side of side 1 -> 1 4 5 11 17 23 250 . side of side 2 -> 2 4 6 12 18 24 251 . side of side 3 -> 3 4 7 13 19 25 252 . side of side 4 -> 2 3 8 14 20 26 253 . side of side 5 -> 1 3 9 15 21 27 254 . side of side 6 -> 1 2 10 16 22 28 255- Dof numbers on sides (RefElement tetrahedron_Lagrange_5) 256 . side 1 -> 4 3 2 25 26 6 19 20 12 13 14 18 7 8 24 29 30 31 32 33 34 257 . side 2 -> 4 1 3 23 9 7 17 15 13 11 21 19 5 27 25 35 36 37 38 39 40 258 . side 3 -> 4 2 1 24 28 5 18 22 11 12 16 17 6 10 23 41 42 43 44 45 46 259 . side 4 -> 1 2 3 10 8 27 16 14 21 22 20 15 28 26 9 47 48 49 50 51 52 260- Side of side shapes : 6 261segment_Lagrange_5, segment_Lagrange_5, segment_Lagrange_5, segment_Lagrange_5, segment_Lagrange_5, segment_Lagrange_5 262- Side shapes : 4 263triangle_Lagrange_5, triangle_Lagrange_5, triangle_Lagrange_5, triangle_Lagrange_5 264 265test of P3 barycentric side dofs 266 BarycentricSideDofs = [ (1, 1, 1) ] 267 BarycentricSideDofMap ={ ((1, 1, 1),1) } 268test of P4 barycentric side dofs 269 BarycentricSideDofs = [ (2, 1, 1) (1, 2, 1) (1, 1, 2) ] 270 BarycentricSideDofMap ={ ((1, 1, 2),3) ((1, 2, 1),2) ((2, 1, 1),1) } 271test of P5 barycentric side dofs 272 BarycentricSideDofs = [ (3, 1, 1) (1, 3, 1) (1, 1, 3) (2, 2, 1) (1, 2, 2) (2, 1, 2) ] 273 BarycentricSideDofMap ={ ((1, 1, 3),3) ((1, 2, 2),5) ((1, 3, 1),2) ((2, 1, 2),6) ((2, 2, 1),4) ((3, 1, 1),1) } 274test of P6 barycentric side dofs 275 BarycentricSideDofs = [ (4, 1, 1) (1, 4, 1) (1, 1, 4) (3, 2, 1) (1, 3, 2) (2, 1, 3) (2, 3, 1) (1, 2, 3) (3, 1, 2) (2, 2, 2) ] 276 BarycentricSideDofMap ={ ((1, 1, 4),3) ((1, 2, 3),8) ((1, 3, 2),5) ((1, 4, 1),2) ((2, 1, 3),6) ((2, 2, 2),10) ((2, 3, 1),7) ((3, 1, 2),9) ((3, 2, 1),4) ((4, 1, 1),1) } 277test of P7 barycentric side dofs 278 BarycentricSideDofs = [ (5, 1, 1) (1, 5, 1) (1, 1, 5) (4, 2, 1) (1, 4, 2) (2, 1, 4) (3, 3, 1) (1, 3, 3) (3, 1, 3) (2, 4, 1) (1, 2, 4) (4, 1, 2) (3, 2, 2) (2, 3, 2) (2, 2, 3) ] 279 BarycentricSideDofMap ={ ((1, 1, 5),3) ((1, 2, 4),11) ((1, 3, 3),8) ((1, 4, 2),5) ((1, 5, 1),2) ((2, 1, 4),6) ((2, 2, 3),15) ((2, 3, 2),14) ((2, 4, 1),10) ((3, 1, 3),9) ((3, 2, 2),13) ((3, 3, 1),7) ((4, 1, 2),12) ((4, 2, 1),4) ((5, 1, 1),1) } 280 lt->sideDofsMap(1,1,2,3) = 1 281 lt->sideDofsMap(1,2,1,3) = 2 282 lt->sideDofsMap(10,2,1,3) = 4 283split P1 to P1 : [ [ 1 2 3 4 ] ] 284split P2 to P1 : [ [ 1 5 9 10 ] [ 2 6 8 10 ] [ 3 7 8 9 ] [ 4 5 6 7 ] [ 5 6 7 8 ] [ 5 6 8 10 ] [ 5 7 8 9 ] [ 5 8 9 10 ] ] 285split P3 to P1 : [ [ 1 5 9 10 ] [ 2 6 8 16 ] [ 3 7 14 15 ] [ 4 11 12 13 ] [ 5 9 10 11 ] [ 6 8 12 16 ] [ 7 13 14 15 ] [ 8 12 16 17 ] [ 8 14 17 20 ] [ 8 16 17 20 ] [ 9 10 11 18 ] [ 9 10 18 20 ] [ 9 15 18 20 ] [ 10 11 18 19 ] [ 10 16 19 20 ] [ 10 18 19 20 ] [ 11 12 13 17 ] [ 11 12 17 19 ] [ 11 13 17 18 ] [ 11 17 18 19 ] [ 12 16 17 19 ] [ 13 14 15 17 ] [ 13 15 17 18 ] [ 14 15 17 20 ] [ 15 17 18 20 ] [ 16 17 19 20 ] [ 17 18 19 20 ] ] 286split P4 to P1 : [ [ 1 5 9 10 ] [ 2 6 8 22 ] [ 3 7 20 21 ] [ 4 17 18 19 ] [ 5 9 10 11 ] [ 6 8 12 22 ] [ 7 13 20 21 ] [ 8 12 22 25 ] [ 8 14 25 33 ] [ 8 22 25 33 ] [ 9 10 11 27 ] [ 9 10 27 32 ] [ 9 15 27 32 ] [ 10 11 27 31 ] [ 10 16 31 32 ] [ 10 27 31 32 ] [ 11 17 26 29 ] [ 11 26 27 29 ] [ 11 27 29 31 ] [ 12 18 23 29 ] [ 12 22 25 30 ] [ 12 23 25 29 ] [ 12 25 29 30 ] [ 13 19 23 26 ] [ 13 20 21 24 ] [ 13 21 24 28 ] [ 13 23 24 26 ] [ 13 24 26 28 ] [ 14 20 24 34 ] [ 14 23 24 33 ] [ 14 23 25 33 ] [ 14 24 33 34 ] [ 15 21 28 34 ] [ 15 26 27 32 ] [ 15 26 28 32 ] [ 15 28 32 34 ] [ 16 22 30 33 ] [ 16 29 30 32 ] [ 16 29 31 32 ] [ 16 30 32 33 ] [ 17 18 19 23 ] [ 17 18 23 29 ] [ 17 19 23 26 ] [ 17 23 26 29 ] [ 20 21 24 34 ] [ 21 24 28 34 ] [ 22 25 30 33 ] [ 23 24 26 34 ] [ 23 24 33 34 ] [ 23 25 29 33 ] [ 23 26 29 35 ] [ 23 26 34 35 ] [ 23 29 33 35 ] [ 23 33 34 35 ] [ 24 26 28 34 ] [ 25 29 30 33 ] [ 26 27 29 32 ] [ 26 28 32 34 ] [ 26 29 32 35 ] [ 26 32 34 35 ] [ 27 29 31 32 ] [ 29 30 32 33 ] [ 29 32 33 35 ] [ 32 33 34 35 ] ] 287split P5 to P1 : [ [ 1 5 9 10 ] [ 2 6 8 28 ] [ 3 7 26 27 ] [ 4 23 24 25 ] [ 5 9 10 11 ] [ 6 8 12 28 ] [ 7 13 26 27 ] [ 8 12 28 31 ] [ 8 14 31 48 ] [ 8 28 31 48 ] [ 9 10 11 36 ] [ 9 10 36 47 ] [ 9 15 36 47 ] [ 10 11 36 43 ] [ 10 16 43 47 ] [ 10 36 43 47 ] [ 11 17 36 43 ] [ 12 18 31 42 ] [ 12 28 31 42 ] [ 13 19 30 37 ] [ 13 26 27 30 ] [ 13 27 30 37 ] [ 14 20 33 51 ] [ 14 31 33 48 ] [ 14 33 48 51 ] [ 15 21 39 52 ] [ 15 36 38 47 ] [ 15 38 39 47 ] [ 15 39 47 52 ] [ 16 22 45 50 ] [ 16 43 45 47 ] [ 16 45 47 50 ] [ 17 23 35 41 ] [ 17 35 38 41 ] [ 17 36 38 43 ] [ 17 38 41 46 ] [ 17 38 43 46 ] [ 18 24 29 41 ] [ 18 29 34 41 ] [ 18 31 34 42 ] [ 18 34 41 44 ] [ 18 34 42 44 ] [ 19 25 29 35 ] [ 19 29 32 35 ] [ 19 30 32 37 ] [ 19 32 35 40 ] [ 19 32 37 40 ] [ 20 26 30 49 ] [ 20 30 32 49 ] [ 20 32 33 49 ] [ 20 33 49 51 ] [ 21 27 37 49 ] [ 21 37 39 49 ] [ 21 39 49 52 ] [ 22 28 42 48 ] [ 22 42 44 48 ] [ 22 44 45 48 ] [ 22 45 48 50 ] [ 23 24 25 29 ] [ 23 24 29 41 ] [ 23 25 29 35 ] [ 23 29 35 41 ] [ 26 27 30 49 ] [ 27 30 37 49 ] [ 28 31 42 48 ] [ 29 32 33 51 ] [ 29 32 35 55 ] [ 29 32 51 55 ] [ 29 33 34 51 ] [ 29 34 41 54 ] [ 29 34 51 54 ] [ 29 35 41 56 ] [ 29 35 55 56 ] [ 29 41 54 56 ] [ 29 51 54 55 ] [ 29 54 55 56 ] [ 30 32 37 49 ] [ 31 33 34 48 ] [ 31 34 42 48 ] [ 32 33 49 51 ] [ 32 35 40 55 ] [ 32 37 40 49 ] [ 32 40 49 55 ] [ 32 49 51 55 ] [ 33 34 48 51 ] [ 34 41 44 54 ] [ 34 42 44 48 ] [ 34 44 48 54 ] [ 34 48 51 54 ] [ 35 38 39 52 ] [ 35 38 41 53 ] [ 35 38 52 53 ] [ 35 39 40 52 ] [ 35 40 52 55 ] [ 35 41 53 56 ] [ 35 52 53 55 ] [ 35 53 55 56 ] [ 36 38 43 47 ] [ 37 39 40 49 ] [ 38 39 47 52 ] [ 38 41 46 53 ] [ 38 43 46 47 ] [ 38 46 47 53 ] [ 38 47 52 53 ] [ 39 40 49 52 ] [ 40 49 52 55 ] [ 41 44 45 50 ] [ 41 44 50 54 ] [ 41 45 46 50 ] [ 41 46 50 53 ] [ 41 50 53 54 ] [ 41 53 54 56 ] [ 43 45 46 47 ] [ 44 45 48 50 ] [ 44 48 50 54 ] [ 45 46 47 50 ] [ 46 47 50 53 ] [ 47 50 52 53 ] [ 48 50 51 54 ] [ 49 51 52 55 ] [ 50 51 52 53 ] [ 50 51 53 54 ] [ 51 52 53 55 ] [ 51 53 54 55 ] [ 53 54 55 56 ] ] 288