I used the following script, and the raw results are attached.


    import pyperf
    runner = pyperf.Runner()

    ms = [12, 16, 24, 32, 48, 64, 96, 128, 
          192, 256, 384, 512, 768, 1024, 1536]

    ns = [2000, 3000, 4000, 6000, 8000,
          12000, 16000, 24000, 32000, 48000,
          64000, 96000]

    for n, m in product(ns, ms):
        runner.timeit(f"needle={m}, haystack={n}",
                      setup="needle = zipf_string(m); haystack = zipf_string(n)",
                      stmt="needle in haystack",
                      globals=locals() | globals())