We review the various numerical methods that have been developed for calculating the reflection and transmission of ocean waves at a land-fast ice boundary, including recent developments. While an integral form of the solution, found by the Wiener-Hopf technique, has been known for many years, direct numerical computation of this exact solution has been thought to be prohibitively difficult. Instead, several numerical "matching" procedures have been developed, including some that are only approximate, along with asymptotic solutions based on the integral form. Recently it has been discovered that direct calculation of the integral form is feasible, actually requiring less computation than the matching methods. We outline the actual computations required and contrast each method, and provide examples of computation from the integral form.