This library lets you boostrap vector-valued statistics fast as it uses parallel processing. Ploting as well as computation of bias-corrected and accelerated confidence intervals are available. 
pip install speedboot
from speedboot import speedboot
sb_object = speedboot(data=n_sample, stats_fun=estimators)
sb_object.fit(R=999, par=True, seed=1)
sb_object.plot()
sb_object.emp_ci()See a quick demo in demo.ipynb. See toy simulations in simulations.ipynb.
Release history is available on PyPI.