EQ2EMu/EQ2/source/depends/boost_1_72_0/libs/histogram/benchmark/histogram_filling_numpy.py

from __future__ import print_function
import numpy as np
# pip install fast-histogram
from fast_histogram import histogram1d
import timeit

x = np.random.rand(1 << 20)
nrepeat = 10

print(timeit.timeit("np.histogram(x, bins=100, range=(0, 1))",
                    "from __main__ import x, np", number=nrepeat) / (nrepeat * len(x)) / 1e-9)
                    
print(timeit.timeit("histogram1d(x, bins=100, range=(0, 1))",
                    "from __main__ import x, histogram1d", number=nrepeat) / (nrepeat * len(x)) / 1e-9)