这个问题有一个关于pandas apply函数并行化的更常见版本-所以这是一个令人耳目一新的问题:)
首先,我想提一下swifter,因为您需要一个“打包”的解决方案,并且它在大多数有关熊猫并行化的问题上都出现了。
但是,我仍然想分享我的个人要点代码,因为在使用DataFrame几年后,我从未找到100%并行化解决方案(主要是apply函数),而且我总是不得不为我的“手册”代码。
多亏您的支持,我才得以通过其名称来支持任何(理论上)DataFrame方法(因此您不必保留isin,apply等的版本)的通用性。
我使用python 2.7和3.6在“ isin”,“ apply”和“ isna”函数上进行了测试。它不到20行,我遵循了“ subset”和“ njobs”之类的熊猫命名约定。
我还添加了一个与“ isin”的等效代码进行时间比较的功能,它的速度似乎比此要点要慢X2倍。
它包括2个功能:
df_multi_core-这是您调用的那个。它接受:
- 您的df对象
- 您要调用的函数名称
- 可以执行该功能的列的子集(有助于减少时间/内存)
- 并行运行的作业数(所有内核为-1或忽略)
- df函数接受的其他任何变形(例如“轴”)
_df_split-这是一个内部帮助器函数,必须全局定位到正在运行的模块(Pool.map是“与位置相关的”),否则我将在内部定位它。
这是我的要旨中的代码(我将在其中添加更多的pandas功能测试):
import pandas as pd
import numpy as np
import multiprocessing
from functools import partial
def _df_split(tup_arg, **kwargs):
split_ind, df_split, df_f_name = tup_arg
return (split_ind, getattr(df_split, df_f_name)(**kwargs))
def df_multi_core(df, df_f_name, subset=None, njobs=-1, **kwargs):
if njobs == -1:
njobs = multiprocessing.cpu_count()
pool = multiprocessing.Pool(processes=njobs)
try:
splits = np.array_split(df[subset], njobs)
except ValueError:
splits = np.array_split(df, njobs)
pool_data = [(split_ind, df_split, df_f_name) for split_ind, df_split in enumerate(splits)]
results = pool.map(partial(_df_split, **kwargs), pool_data)
pool.close()
pool.join()
results = sorted(results, key=lambda x:x[0])
results = pd.concat([split[1] for split in results])
return results
贝娄(Bellow)是并行isin的测试代码,比较了本机,多核要害和敏捷性能。在具有8个物理内核的I7机器上,我的速度提高了大约X4倍。我很想听听您从真实数据中得到什么!
from time import time
if __name__ == '__main__':
sep = '-' * 50
# isin test
N = 10000000
df = pd.DataFrame({'c1': np.random.randint(low=1, high=N, size=N), 'c2': np.arange(N)})
lookfor = np.random.randint(low=1, high=N, size=1000000)
print('{}\ntesting pandas isin on {}\n{}'.format(sep, df.shape, sep))
t1 = time()
print('result\n{}'.format(df.isin(lookfor).sum()))
t2 = time()
print('time for native implementation {}\n{}'.format(round(t2 - t1, 2), sep))
t3 = time()
res = df_multi_core(df=df, df_f_name='isin', subset=['c1'], njobs=-1, values=lookfor)
print('result\n{}'.format(res.sum()))
t4 = time()
print('time for multi core implementation {}\n{}'.format(round(t4 - t3, 2), sep))
t5 = time()
ddata = dd.from_pandas(df, npartitions=njobs)
res = ddata.map_partitions(lambda df: df.apply(apply_f, axis=1)).compute(scheduler='processes')
t6 = time()
print('result random sample\n{}'.format(res.sample(n=3, random_state=0)))
print('time for dask implementation {}\n{}'.format(round(t6 - t5, 2), sep))
--------------------------------------------------
testing pandas isin on (10000000, 2)
--------------------------------------------------
result
c1 953213
c2 951942
dtype: int64
time for native implementation 3.87
--------------------------------------------------
result
c1 953213
dtype: int64
time for multi core implementation 1.16
--------------------------------------------------
result
c1 953213
c2 951942
dtype: int64
time for dask implementation 2.88