18

我想使用指数增长的值来测试一些代码的性能。因此,当将一个额外的数字添加到 numbers_size 时,增量乘以 10。这就是我到目前为止的做法,但它看起来有点 hacky。在不引入非标准库的情况下提出改进建议?

numbers_size = 100
increment = 100
numbers_range = 1000000000
while numbers_size < numbers_range:
    t = time.time()
    test( numbers_size )
    taken_t = time.time() - t
    print numbers_size, test, taken_t

    increment = 10 ** (len(str(numbers_size))-1)
    numbers_size += increment
4

8 回答 8

23

如果您将 numpy 视为标准之一;),您可以使用numpy.logspace因为这是它应该做的......(注意:100=10^2、1000000000=10^9)

for n in numpy.logspace(2,9,num=9-2, endpoint=False):
    test(n)

例子2(注:100=10^2, 1000000000=10^9,想走10x步,就是9-2+1分...):

In[14]: np.logspace(2,9,num=9-2+1,base=10,dtype='int')
Out[14]: 
array([       100,       1000,      10000,     100000,    1000000,
         10000000,  100000000, 1000000000])

示例 3:

In[10]: np.logspace(2,9,dtype='int')
Out[10]: 
array([       100,        138,        193,        268,        372,
              517,        719,       1000,       1389,       1930,
             2682,       3727,       5179,       7196,      10000,
            13894,      19306,      26826,      37275,      51794,
            71968,     100000,     138949,     193069,     268269,
           372759,     517947,     719685,    1000000,    1389495,
          1930697,    2682695,    3727593,    5179474,    7196856,
         10000000,   13894954,   19306977,   26826957,   37275937,
         51794746,   71968567,  100000000,  138949549,  193069772,
        268269579,  372759372,  517947467,  719685673, 1000000000])

在您的情况下,我们使用,endpoint=False因为您不想包含端点...(例如np.logspace(2,9,num=9-2, endpoint=False)

于 2015-08-10T12:45:39.853 回答
15

为什么不

for exponent in range(2, 10):
    test(10 ** exponent)

如果我正确地阅读了您的意图。

于 2012-07-12T01:12:33.800 回答
14

要生成与您的代码相同的数字:

numbers_sizes = (i*10**exp for exp in range(2, 9) for i in range(1, 10))
for n in numbers_sizes:
    test(n)
于 2012-07-12T01:35:34.730 回答
5

最简单的做法是使用指数的线性序列:

for e in range(1, 90):
    i = int(10**(e/10.0))
    test(i)

您可以将序列抽象为它自己的生成器:

def exponent_range(max, nsteps):
    max_e = math.log10(max)
    for e in xrange(1, nsteps+1):
        yield int(10**(e*max_e/nsteps))

for i in exponent_range(10**9, nsteps=100):
    test(i)
于 2012-07-12T01:11:43.623 回答
5

我喜欢 Ned Batcheldor 的回答,但我会让它更笼统一些:

def exp_range(start, end, mul):
    while start < end:
        yield start
        start *= mul

然后你的代码变成

for sz in exp_range(100, 1000000000, 10):
    t = time.time()
    test(sz)
    print sz, test(sz), time.time()-t
于 2012-07-12T01:25:28.567 回答
1

OP 写道“在不引入非标准库的情况下提出改进建议?”

为了完整起见,这里有一个生成指数范围的方法——每个元素都是一个比前一个更大的固定因子:

from math import exp
from math import log

def frange(start, stop, numelements):
    """range function for floats"""
    incr = (stop - start) / numelements
    return (start + x * incr for x in range(numelements))

def exprange(start, stop, numelements):
    """exponential range - each element is a fixed factor bigger than the previous"""
    return (exp(x) for x in frange(log(start), log(stop), numelements))

测试:

print(", ".join("%.3f" % x for x in exprange(3,81,6)))

输出:

3.000, 5.196, 9.000, 15.588, 27.000, 46.765
于 2018-06-18T06:13:38.890 回答
0

使用生成器表达式:

max_exponent = 100
for i in (10**n for n in xrange(1, max_exponent)):
    test(i)
于 2012-07-12T01:16:04.480 回答
-2

“没有正确阅读问题”和“不怎么做”的例子

for i in xrange(100, 1000000000, 100):
    # timer
    test(i)
    # whatever

就像它得到的一样简单......xrange相应地调整

于 2012-07-12T01:09:36.247 回答