In [61]:

    

import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline


    

In [62]:

    

L = 2 * np.pi
dx = 0.0122718463031
x = np.arange(0, L, dx)
A = np.linspace(0, L, 100000, endpoint=False)
A.max()


    

    
Out[62]:





6.2831224753265138

In [63]:

    

%%timeit
res1 = (A / dx).astype(int)


    

    




1000 loops, best of 3: 644 µs per loop

In [64]:

    

%%timeit
res2 = (A // dx).astype(int)


    

    




100 loops, best of 3: 5.51 ms per loop

In [65]:

    

%%timeit
res3 = np.floor_divide(A, dx).astype(int)


    

    




100 loops, best of 3: 5.56 ms per loop

In [66]:

    

results = [res1, res2, res3]
for r1 in results:
    plt.plot(r1)
    for r2 in results:
        if r2 is not r1:
            assert np.allclose(r1, r2)