Question

我正在制作一个素数计算器，它完美地工作，但我希望通过使用多线程来加快它。我想知道你们中是否有人能指出某些资源（python文档不是很有用）或者在我的代码中给我一个例子。

import time
#Set variables
check2 = 0
check3 = 0
check5 = 0
check7 = 0
#get number
primenum = int(input("What number do you want to find out whether it is prime or not?\nIt must be a natural number\n**Calculations may take some time depending of the power of the computer and the size of the number**\n"))
#set divisor
primediv = primenum - 2
#assume it's prime until proven innocent
prime = 1
#Create variable for GUI
working = "Calculating"
work = 10000000
#set the number of divides to 0
divnum = 0
#start the clock
starttime = time.perf_counter()
#until it is not prime or divided by 1...
while prime == 1 and primediv >7:
    #does simple checks to deal with large numbers quickly
    #Does this by dividing with small numbers first
    if primenum != 0 and check2 == 0:
        primemod = primenum % 2
        check2 = 1
        print(working + ".")
        working = working +"."
    elif primenum != 0 and check3 == 0:
        primemod = primenum % 3
        check3 = 1
        print(working + ".")
        working = working +"."
    elif primenum != 0 and check5 == 0:
        primemod = primenum % 5
        check5 = 1
        print(working + ".")
        working = working + "."
    elif primenum != 0 and check7 == 0:
        primemod = primenum % 7
        check7 = 1
        print(working + ".")
        working = working + "."
    #divde and get remainder
    else:
        primemod = primenum % primediv
        #Working visuals
        if divnum == work:
            print(working + ".")
            working = working +"."
            work = work + 10000000
    #if the can't be devided evenly
    if primemod == 0:
        #then it isn't a prime
        prime = 0
    else:
        #if it can, keep going
        primediv = primediv - 2
        divnum = divnum + 1
#print results
if prime == 1:
    print("That number is prime")
    print ("It took ", time.perf_counter()-starttime, " seconds to perform that calculation\n")
else:
    print("That number is not prime")
    print ("It took ", time.perf_counter()-starttime, " seconds to perform that calculation\n")

Answer 1

最流行的Python（CPython和PyPy）实现带有全局解释器锁（“GIL”）。这样做的目的基本上是简化内存管理。它的作用是一次只有一个线程可以执行Python字节码。因此，对于在Python中执行所有计算的程序，线程化不会使其更快。

通常最好的方法是让程序更快地找到更好的算法。看，例如在这：

def prime_list(num):
    """Returns a list of all prime numbers up to and including num.

    :num: highest number to test
    :returns: a list of primes up to num

    """
    if num < 3:
        raise ValueError('this function only accepts arguments > 2')
    candidates = range(3, num+1, 2)  # (1)
    L = [c for c in candidates if all(c % p != 0 for p in range(3, c, 2))]  #(2)
    return [2] + L  # (3)

（1）偶数＆gt; 2可以除以2，因此它们不能是素数，不需要考虑。

（2）：对于奇数c为素数，必须确保c模数所有先前的奇数（p）必须为非零。

（3）：2也是素数。

一点点改进是p只需要达到sqrt(c)。

def prime_list2(num):
    if num < 3:
        raise ValueError('this function only accepts arguments > 2')
    candidates = range(3, num+1, 2)
    L = [c for c in candidates if all(c % p != 0 for p in
         range(3, int(math.sqrt(c))+1, 2))]
    return [2] + L

请注意，我使用列表推导而不是for-loops，因为它们通常更快。

最后，亚当史密斯提到的一种筛子;

def prime_list3(num):
    num += 1
    candidates = range(3, num, 2)
    results = [2]
    while len(candidates):
        t = candidates[0]
        results.append(t)
        candidates = [i for i in candidates if not i in range(t, num, t)]
    return results

要查看哪一个更快，我们需要运行测试;

首先是num=100。

In [8]: %timeit prime_list(100)
1000 loops, best of 3: 185 µs per loop

In [9]:  %timeit prime_list2(100)
10000 loops, best of 3: 156 µs per loop

In [10]: %timeit prime_list3(100)
1000 loops, best of 3: 313 µs per loop

然后1000：

In [11]: %timeit prime_list(1000)
100 loops, best of 3: 8.67 ms per loop

In [12]: %timeit prime_list2(1000)
1000 loops, best of 3: 1.94 ms per loop

In [13]: %timeit prime_list3(1000)
10 loops, best of 3: 21.6 ms per loop

然后10000;

In [14]: %timeit prime_list(10000)
1 loops, best of 3: 680 ms per loop

In [15]: %timeit prime_list2(10000)
10 loops, best of 3: 25.7 ms per loop

In [16]: %timeit prime_list3(10000)
1 loops, best of 3: 1.99 s per loop

在这些算法中，prime_list2似乎保持最佳状态。

使程序多线程

1 个答案: