1----------------------------- GeomRefPyramid test ---------------------------------- 2geometric reference element : shape pyramid, 5 vertices, 5 sides, 8 sideOfSides 3 * vertex 1 -> (0 ,0 ,0) 4 * vertex 2 -> (1 ,0 ,0) 5 * vertex 3 -> (1 ,1 ,0) 6 * vertex 4 -> (0 ,1 ,0) 7 * vertex 5 -> (0 ,0 ,1) 8 * centroid : (0.4 ,0.4 ,0.2) 9 * measure : 0.333333333333 10 * side 1 -> shape quadrangle, vertices : 1 2 3 4 11 * side 2 -> shape triangle, vertices : 1 2 5 12 * side 3 -> shape triangle, vertices : 2 3 5 13 * side 4 -> shape triangle, vertices : 3 4 5 14 * side 5 -> shape triangle, vertices : 4 1 5 15 * side of side 1 -> [ 1 2 ] 16 * side of side 2 -> [ 1 4 ] 17 * side of side 3 -> [ 1 5 ] 18 * side of side 4 -> [ 2 3 ] 19 * side of side 5 -> [ 2 5 ] 20 * side of side 6 -> [ 3 4 ] 21 * side of side 7 -> [ 3 5 ] 22 * side of side 8 -> [ 4 5 ] 23 * numbers on side1 -> [ 1 4 6 -2 ] 24 * numbers on side2 -> [ 1 5 -3 ] 25 * numbers on side3 -> [ 4 7 -5 ] 26 * numbers on side4 -> [ 6 8 -7 ] 27 * numbers on side5 -> [ -2 3 -8 ] 28sideWithVertices(1,2,3) = 1 29sideWithVertices(1,2,5) = 2 30sideWithVertices(2,3,5) = 3 31sideWithVertices(3,4,5) = 4 32sideWithVertices(4,1,5) = 5 33sideWithVertices(1,2) = 1 34sideWithVertices(1,4) = 2 35sideWithVertices(1,5) = 3 36sideWithVertices(2,3) = 4 37sideWithVertices(2,5) = 5 38sideWithVertices(3,4) = 6 39sideWithVertices(3,5) = 7 40sideWithVertices(4,5) = 8 41 42----------------------------- interpolation test ---------------------------------- 43reference element pyramid_Lagrange_1 44- finite element interpolation : Lagrange of order 1 conforming in H1 45- geometric reference element : shape pyramid, 5 vertices, 5 sides, 8 sideOfSides 46- reference DoFs : 5 47 * 1 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 1, index 1, coords = [ 0 0 0 ] 48 * 2 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 2, index 1, coords = [ 1 0 0 ] 49 * 3 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 3, index 1, coords = [ 1 1 0 ] 50 * 4 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 4, index 1, coords = [ 0 1 0 ] 51 * 5 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 5, index 1, coords = [ 0 0 1 ] 52- Dof numbers on side of sides (RefElement pyramid_Lagrange_1) 53 . side of side 1 -> 1 2 54 . side of side 2 -> 1 4 55 . side of side 3 -> 1 5 56 . side of side 4 -> 2 3 57 . side of side 5 -> 2 5 58 . side of side 6 -> 3 4 59 . side of side 7 -> 3 5 60 . side of side 8 -> 4 5 61- Dof numbers on sides (RefElement pyramid_Lagrange_1) 62 . side 1 -> 1 2 3 4 63 . side 2 -> 1 2 5 64 . side 3 -> 2 3 5 65 . side 4 -> 3 4 5 66 . side 5 -> 4 1 5 67- Side of side shapes : 8 68segment_Lagrange_1, segment_Lagrange_1, segment_Lagrange_1, segment_Lagrange_1, segment_Lagrange_1, segment_Lagrange_1, segment_Lagrange_1, segment_Lagrange_1 69- Side shapes : 5 70quadrangle_Lagrange_1, triangle_Lagrange_1, triangle_Lagrange_1, triangle_Lagrange_1, triangle_Lagrange_1 71 72reference element pyramid_Lagrange_2 73- finite element interpolation : Lagrange of order 2 conforming in H1 74- geometric reference element : shape pyramid, 5 vertices, 5 sides, 8 sideOfSides 75- reference DoFs : 14 76 * 1 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 1, index 1, coords = [ 0 0 0 ] 77 * 2 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 2, index 1, coords = [ 1 0 0 ] 78 * 3 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 3, index 1, coords = [ 1 1 0 ] 79 * 4 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 4, index 1, coords = [ 0 1 0 ] 80 * 5 -> reference DoF : nodal value, sharable? yes, support : dim 0, on vertex 5, index 1, coords = [ 0 0 1 ] 81 * 6 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 1, index 1, coords = [ 0.5 0 0 ] 82 * 7 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 2, index 1, coords = [ 0 0.5 0 ] 83 * 8 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 3, index 1, coords = [ 0 0 0.5 ] 84 * 9 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 4, index 1, coords = [ 1 0.5 0 ] 85 * 10 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 5, index 1, coords = [ 0.5 0 0.5 ] 86 * 11 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 6, index 1, coords = [ 0.5 1 0 ] 87 * 12 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 7, index 1, coords = [ 0.5 0.5 0.5 ] 88 * 13 -> reference DoF : nodal value, sharable? yes, support : dim 0, on edge 8, index 1, coords = [ 0 0.5 0.5 ] 89 * 14 -> reference DoF : nodal value, sharable? yes, support : dim 0, on face 0, index 1, coords = [ 0.5 0.5 0 ] 90- Dof numbers on side of sides (RefElement pyramid_Lagrange_2) 91 . side of side 1 -> 1 2 6 92 . side of side 2 -> 1 4 7 93 . side of side 3 -> 1 5 8 94 . side of side 4 -> 2 3 9 95 . side of side 5 -> 2 5 10 96 . side of side 6 -> 3 4 11 97 . side of side 7 -> 3 5 12 98 . side of side 8 -> 4 5 13 99- Dof numbers on sides (RefElement pyramid_Lagrange_2) 100 . side 1 -> 1 2 3 4 6 9 11 7 14 101 . side 2 -> 1 2 5 6 10 8 102 . side 3 -> 2 3 5 9 12 10 103 . side 4 -> 3 4 5 11 13 12 104 . side 5 -> 4 1 5 7 8 13 105- Side of side shapes : 8 106segment_Lagrange_2, segment_Lagrange_2, segment_Lagrange_2, segment_Lagrange_2, segment_Lagrange_2, segment_Lagrange_2, segment_Lagrange_2, segment_Lagrange_2 107- Side shapes : 5 108quadrangle_Lagrange_2, triangle_Lagrange_2, triangle_Lagrange_2, triangle_Lagrange_2, triangle_Lagrange_2 109 110split Fast O1 to P1 : [ [ 2 4 5 1 ] [ 4 2 5 3 ] ] 111split Fast O2 to P1 : [ [ 6 7 8 1 ] [ 7 6 8 3 ] [ 2 14 10 6 ] [ 14 2 10 9 ] [ 9 11 12 14 ] [ 11 9 12 3 ] [ 14 4 13 7 ] [ 4 14 13 11 ] [ 10 13 5 8 ] [ 13 10 5 12 ] [ 13 10 14 12 ] [ 10 13 14 8 ] [ 11 12 14 13 ] [ 10 12 14 9 ] [ 8 10 14 6 ] [ 13 8 14 7 ] ] 112split Fast O1 to O1 : [ (9,[ 1 2 3 4 5 ]) ] 113split Fast O2 to O1 : [ (9,[ 1 6 14 7 8 ]) (9,[ 6 2 9 14 10 ]) (9,[ 14 9 3 11 12 ]) (9,[ 7 14 11 4 13 ]) (9,[ 8 10 12 13 5 ]) (9,[ 12 13 8 10 14 ]) (6,[ 11 12 14 13 ]) (6,[ 10 12 14 9 ]) (6,[ 8 10 14 6 ]) (6,[ 13 8 14 7 ]) ] 114