TFA is a decent first cut. It may be inefficient (the author already mentioned and discarded a N^2 approach), but it’s easily understandable. The author also mentioned that he was planning a NumPy version, so maybe some of the topics you addressed will be covered.