• # Et par rapport a du pur Python ?

    Posté par . En réponse au journal Pythran 0.2 : Python peut-il être aussi rapide que du C ?. Évalué à 3.

    Et qu'en est-il du gain de vitesse de Pythran par rapport a du pur numpy (codé correctement). L'exemple donné reviens à:

    def run(xmin, xmax, ymin, ymax, step, range_, t):
 X,Y = meshgrid(linspace(xmin,xmax,step),linspace(ymin,ymax,step))
 pt = zeros((step,step,3))
 pt[:,:,0] = X*180/pi
 pt[:,:,1] = Y*180/pi
 for k in range(t.shape[0]):
 tmp = 6368.*arccos(cos(X)*cos(t[k,0])*cos(Y-t[k,1])+sin(X)*sin(t[k,1]))
 pt[:,:,2] += where(tmp < range_, t[k,2]/(1+tmp), 0)
 return pt

    Je n'ai pas Pythran pour comparer. Quel gain peut-on en attendre?

    A noter que cet exemple a été écrit vite-fait, il y a probablement moyen d'optimiser en faisant sauter la boucle sur k.