The algorithm implementation exists in pure Python currently. We should port it to rust to be more performant, but maybe more importantly to leverage it directly from other rust functions in the future.