性能追求第二部分：Perl 与 Python

import argparse
import time
import math
import numpy as np
import os
from numba import njit
from joblib import parallel, delayed

parser = argparse.argumentparser()
parser.add_argument("--workers", type=int, default=8)
parser.add_argument("--arraysize", type=int, default=100_000_000)
args = parser.parse_args()
# set the number of threads to 1 for different libraries
print("=" * 80)
print(
    f"\nstarting the benchmark for {args.arraysize} elements "
    f"using {args.workers} threads/workers\n"
)

# generate the data structures for the benchmark
array0 = [np.random.rand() for _ in range(args.arraysize)]
array1 = array0.copy()
array2 = array0.copy()
array_in_np = np.array(array1)
array_in_np_copy = array_in_np.copy()

  for i in range(len(array0)):
    array0[i] = math.cos(math.sin(math.sqrt(array0[i])))

np.sqrt(array_in_np, out=array_in_np)
np.sin(array_in_np, out=array_in_np)
np.cos(array_in_np, out=array_in_np)

def compute_inplace_with_joblib(chunk):
    return np.cos(np.sin(np.sqrt(chunk))) #parallel function for joblib

chunks = np.array_split(array1, args.workers)  # split the array into chunks
numresults = parallel(n_jobs=args.workers)(
        delayed(compute_inplace_with_joblib)(chunk) for chunk in chunks
    )# process each chunk in a separate thread
array1 = np.concatenate(numresults)  # concatenate the results

@njit
def compute_inplace_with_numba(array):
    np.sqrt(array,array)
    np.sin(array,array)
    np.cos(array,array)
    ## njit will compile this function to machine code
compute_inplace_with_numba(array_in_np_copy)

in place in (  base python): 11.42 seconds
in place in (python joblib): 4.59 seconds
in place in ( python numba): 2.62 seconds
in place in ( python numpy): 0.92 seconds

in place in (  base python): 11.89 seconds
in place in (python joblib): 4.42 seconds
in place in ( python numpy): 0.93 seconds
in place in ( python numba): 0.49 seconds

n<-50000000
x<-runif(n)
start_time <- sys.time()
result <- cos(sin(sqrt(x)))
end_time <- sys.time()

# calculate the time taken
time_taken <- end_time - start_time

# print the time taken
print(sprintf("time in base r: %.2f seconds", time_taken))

Time in base R: 1.30 seconds