First I tried "Poly-> Put Models on Faces/Vertices", choosing the vertex option and applying to all, but you can't easily control orientation with this feature, and it ended up like this:
So I thought of another approach. First, set the edge length to 1 (Scale->Edge Length), and then measure the height of a tetrahedron using measurement mode (enter measurement mode, Shift-click on a vertex, then Shift-click on opposite face).
Go to "Scale->Measured Distance" and copy the value. You can exit measurement mode now.
Now create a triangular prism. Set the edge length to 1.
Go to "Scale->Stretch (Non-Uniform Scale)". Select "Set height to" and paste in the height of a tetrahedron from earlier, then add "* 2" at the end. Yep, it can do inline calculations.
Now we have a triangular prism twice the height of a tetrahedron. If we excavate a tetrahedron from each end, their vertices will touch in the middle.
First, put this prism in a memory slot, we'll need it later (type "m1").
Let's excavate a tetrahedron from just one end: select one of the triangular faces, "Poly->Augment Polyhedron" (or just hit "a"), select the excavate option and "selected face only".
Hit OK and hit Enter to confirm the preview. Select the triangle at the other end, the one that wasn't excavated, and put it in another memory slot (type "m2").
Now go back to a tetrahedron and use Poly->Augment Polyhedron again. Select the second memory slot from the drop-down, and "Excavate" and "All faces of same type".
Hit OK and Enter. You should get this:
Select one of the rectangular (yellow) faces and augment again, this time choosing the first memory slot. Hit OK and Enter again, and prepare for some magic:
