21.12 tetra.input
)set expose add con DenavitHartenbergMatrix Bring DH matrices into the
environment
x1:DFLOAT := sqrt(2.0@DFLOAT/3.0@DFLOAT) Set up the coordinates of the
x2:DFLOAT := sqrt(3.0@DFLOAT)/6 corners of the tetrahedron.
p1 := point [-0.5@DFLOAT, -x2, 0.0@DFLOAT] Some needed points
p2 := point [0.5@DFLOAT, -x2, 0.0@DFLOAT]
p3 := point [0.0@DFLOAT, 2*x2, 0.0@DFLOAT]
p4 := point [0.0@DFLOAT, 0.0@DFLOAT, x1]
baseTriangle := [p2, p1, p3] The base of the tetrahedron
mt := [0.5@DFLOAT*(p2+p1), 0.5@DFLOAT*(p1+p3), 0.5@DFLOAT*(p3+p2)]
The middle triangle inscribed
in the base of the tetrahedron
bt1 := [mt.1, p1, mt.2] The bases of the triangles of
bt2 := [p2, mt.1, mt.3] the subdivided tetrahedron
bt3 := [mt.2, p3, mt.3]
bt4 := [0.5@DFLOAT*(p2+p4), 0.5@DFLOAT*(p1+p4), 0.5@DFLOAT*(p3+p4)]
tt1 := tri2tri(baseTriangle, bt1) Create the transformations
tt2 := tri2tri(baseTriangle, bt2) that bring the base of the
tt3 := tri2tri(baseTriangle, bt3) tetrahedron to the bases of
tt4 := tri2tri(baseTriangle, bt4) the subdivided tetrahedron
drawPyramid(n) == Draw a Sierpinsky tetrahedron
s := createThreeSpace() with levels of recursive
dh := rotatex(0.0@DFLOAT) subdivision
drawPyramidInner(s, n, dh)
makeViewport3D(s, "Sierpinsky Tetrahedron")
drawPyramidInner(s, n, dh) == Recursively draw a Sierpinsky
n = 0 => makeTetrahedron(s, dh, n) tetrahedron
drawPyramidInner(s, n-1, dh * tt1) Draw the 4 recursive pyramids
drawPyramidInner(s, n-1, dh * tt2)
drawPyramidInner(s, n-1, dh * tt3)
drawPyramidInner(s, n-1, dh * tt4)
makeTetrahedron(sp, dh, color) == Draw a tetrahedron into the
w1 := dh*p1 given space with the given
w2 := dh*p2 color, transforming it by
w3 := dh*p3 the given DH matrix
w4 := dh*p4
polygon(sp, [w1, w2, w4])
polygon(sp, [w1, w3, w4])
polygon(sp, [w2, w3, w4])
void()
Sierpinsky's Tetrahedron