A convergent numerical method for \(\alpha \) -dissipative solutions of the Hunter–Saxton equation is derived. The method is based on applying a tailor-made projection operator to the initial data, and then solving exactly using the generalized method of characteristics. The projection step is the only step that introduces any approximation error. It is therefore crucial that its design ensures not only a good approximation of the initial data, but also that errors due to the energy dissipation at later times remain small. Furthermore, it is shown that the main quantity of interest, the wave profile, converges in \(L^{\infty }\) for all \(t \ge 0\) , while a subsequence of the energy density converges weakly for almost every time.
