Edit:
In reference to Sam Nead's answer:
For this I get 0.479646 + 2.91505 I.
In checking this with SnapPy (which gives a cusp shape of -2.4902446675 + 2.9794470665 I), I noticed that, while I got the same tetrahedra shapes as Takahashi, SnapPy disagrees with those (my y is the difference between SnapPy's first two shapes).
Further, https://knotinfo.math.indiana.edu gives translations yielding a cusp shape of 2.49024 - 2.97945 I. I am assuming this is simply a choice of orientation (the real part of the meridian differs by a sign). When I followed Yokota, I got SnapPy's result.
I know that there can exist inequivalent cusp tilings, but it seems to me that the tetrahedra shapes should be invariant, and certainly the cusp shape is an invariant. I am now pretty much thoroughly confused.
Further edit: