大数质数生成算法

大数质数生成，RES用

import numpy as np

def is_prime(number):
    test = (number%np.arange(2,int(np.sqrt(number))))==0
    ones = np.ones(test.shape[0])
    return np.dot(test,ones)==0

def getprime(number1,number2):
   assert number > 3
   while True:
       integer = np.random.randint(number1,number2)
       if is_prime(integer):
           return integer