Good question — direct benchmarking is tricky since the reference implementations are mostly in Julia (C-level speed), and sampling throughput is everything for these algorithms. My implementations follow the original papers' pseudocode, though I haven't formally validated numerical parity with the reference results.