| 82 | [(1,0,1), (1,1,0), (0,1,1), (1,0,0)]] |
| 83 | |
| 84 | def makeP2Tets( p1_tets ): |
| 85 | midpoint = lambda p0, p1: tuple((0.5*(p0[i]+p1[i]) for i in range(3))) |
| 86 | p2_tets = [] |
| 87 | for tet in p1_tets: |
| 88 | tet.append( midpoint(tet[0], tet[3]) ) |
| 89 | tet.append( midpoint(tet[1], tet[3]) ) |
| 90 | tet.append( midpoint(tet[2], tet[3]) ) |
| 91 | tet.append( midpoint(tet[0], tet[1]) ) |
| 92 | tet.append( midpoint(tet[0], tet[2]) ) |
| 93 | tet.append( midpoint(tet[1], tet[2]) ) |
| 94 | p2_tets.append(tet) |
| 95 | return p2_tets |
| 96 | for eltype in p1_tets: |
| 97 | points = p1_tets[eltype] |
| 98 | if order == 2: |