Performances, la NumWorks détrône la HP Prime ! :o

[compsystems]

> The _python one doesn't use semicolons.
true

Code: Tout sélectionner

#cas
def fractal_python(W,H,N):
   local x,y,z,c,j,t,col
   for x in range(W):
      for y in range(H):
         z=0
         c=2.7*x/(W-1)-2.1-i*(1.87*y/(H-1)-.935)
         j=0
         while j<N and abs(z)<2:
            j=j+1
            z=z*z+c
         #end while
         t=255j/N
         col=RGB(IP(t),IP(.75*t),IP(.25*t))
         PIXON_P(x,y,col)
      #end for
   #end for
   FREEZE
   WAIT(0)
#end

[critor]

J'ai adapté vos dernières améliorations, Parisse.

Voici les codes et durées pour les paramètres {320,222,10} :

HP Prime - CAS+Python :

Code: Tout sélectionner

#cas
def fractal_python(W,H,N) :
local x,y,z,c,j,t,w,h,n;
w=2.7/(W-1)
h=1.87/(H-1)
n=255.0/N
for x in range(W):
  for y in range(H):
   z=0
   c=w*x-2.1-i*(h*y-.935)
   for j in range(N):
    if abs(z=z*z+c)>2:
     break
   t=n*j
   PIXON_P(x,y,RGB(IP(t),IP(.75*t),IP(.25*t)))
FREEZE
WAIT(0)
#end

3mins44.96s

HP Prime - CAS+HPPPL :

Code: Tout sélectionner

#cas
fractal_cas(W,H,N)
BEGIN
local x,y,z,c,j,t,w,h,n;
w:=2.7/(W-1);
h:=1.87/(H-1);
n:=255.0/N;
FOR x FROM 0 TO W-1 DO
  FOR y FROM 0 TO H-1 DO
   z:=0;
   c:=w*x-2.1-i*(h*y-.935);
   FOR j FROM 0 to N-1 DO
    IF abs(z:=z*z+c)>2 THEN
     break;
    END;
   END;
   t:=n*j;
   PIXON_P(x,y,RGB(IP(t),IP(.75*t),IP(.25*t)));
  END;
END;
FREEZE;
WAIT(0);
END;
#end

3mins46.72s

HP Prime - HPPPL :

Code: Tout sélectionner

EXPORT fractal(W,H,N)
BEGIN
LOCAL x,y,z,c,j,t,w,h;
LOCAL n;
w:=2.7/(W-1);
h:=1.87/(H-1);
n:=255.0/N;
FOR x FROM 0 TO W-1 DO
  FOR y FROM 0 TO H-1 DO
   z:=0;
   c:=w*x-2.1-i*(h*y-.935);
   FOR j FROM 0 to N-1 DO
    IF ABS(z:=z*z+c)>2 THEN
     break;
    END;
   END;
   t:=n*j;
   PIXON_P(x,y,RGB(IP(t),IP(.75*t),IP(.25*t)));
  END;
END;
FREEZE;
WAIT(0);
END;

1min43.18s

NumWorks :

Code: Tout sélectionner

import kandinsky
def mandelbrot(W,H,N) :
w=2.7/(W-1)
h=1.87/(H-1)
n=255.0/N
for x in range(W):
  for y in range(H):
   z=complex(0,0)
   c=complex(w*x-2.1,.935-h*y)
   for j in range(N):
    z=z*z+c
    if abs(z)>2:
     break
   t=n*j
   kandinsky.set_pixel(x,y,kandinsky.color(int(t),int(.75*t),int(.25*t)))

1min21.8s
Notons que la NumWorks est ici désavantagée, car elle refuse les affectations en tant qu'argument, optimisation donc non présente contrairement aux versions ci-dessus.

Dans tous les cas, cela ne change rien au rapport de performances.

[critor]

Edité le 1er temps du message précédent.
Même si le code était visible dans l'éditeur, je pense que la machine continuait à exécuter la version d'hier.

[parisse]

Bon, ca s'approche des temps que j'avais en memoire. On a bien gagne en remplacant le t qui etait un rationnel avant par un flottant : les couts multiplication rationnel par flottant sont evidemment bien plus importants que flottant par flottant, et c'etait un biais defavorable a HP CAS. J'ai trouve un autre biais defavorable a la Prime, il provient du calcul de la couleur car RGB et IP ne sont pas des instructions CAS, il y a donc des couts de conversion CAS->Home->CAS, dans mon propre code je faisais juste PIXON_P(x,y,255*20*j+256). Il faudrait comparer les 2 versions HP et Numworks avec ce code couleur pour eviter le biais. Je pense que cote Numworks on va gagner quelques pourcents, mais cote HP CAS on va gagner bien plus.

[parisse]

Voila, ca prend 152s sans RGB ni IP, un peu moins que les 170s que j'avais en memoire, le pre-calcul fait gagner environ 10%. Bien sur on reste plus lent que HP PPL et que la Numworks, mais dans un rapport 2, ce qui me semble un prix plus que raisonnable a payer pour beneficier du calcul exact et formel. Mes tests sur Mac indiquent un rapport 7 entre Python et Xcas pour les calculs en complexes, une partie de ce facteur est compensee par le processeur de la Prime plus puissant que celui de la Numworks.
Concernant la difference HP PPL vs Numworks, je pense que le processeur plus puissant de la HP est compense par les calculs en flottants qui sont faits en base 10 dans Home au lieu de la base 2 dans Python.

[critor]

Voici donc les performances avec les 4 programmes équivalents n'utilisant pas les fonctions de couleurs.

HP Prime - CAS+Python :

Code: Tout sélectionner

#cas
def fractal_python(W,H,N) :
local x,y,z,c,j,t,w,h,n;
w=2.7/(W-1)
h=1.87/(H-1)
n=255.0/N
for x in range(W):
  for y in range(H):
   z=0
   c=w*x-2.1-i*(h*y-.935)
   for j in range(N):
    z:=z*z+c;
    if abs(z)>2:
     break
   t=n*j
   PIXON_P(x,y,255*20*j+256)
FREEZE
WAIT(0)
#end

2mins42.90s

HP Prime - CAS+HPPPL :

Code: Tout sélectionner

#cas
fractal_cas(W,H,N)
BEGIN
local x,y,z,c,j,t,w,h,n;
w:=2.7/(W-1);
h:=1.87/(H-1);
n:=255.0/N;
FOR x FROM 0 TO W-1 DO
  FOR y FROM 0 TO H-1 DO
   z:=0;
   c:=w*x-2.1-i*(h*y-.935);
   FOR j FROM 0 to N-1 DO
    z:=z*z+c;
    IF abs(z)>2 THEN
     break;
    END;
   END;
   t:=n*j;
   PIXON_P(x,y,255*20*j+256);
  END;
END;
FREEZE;
WAIT(0);
END;
#end

2mins43.08s

HP Prime - HPPPL :

Code: Tout sélectionner

EXPORT fractal(W,H,N)
BEGIN
LOCAL x,y,z,c,j,t,w,h;
LOCAL n;
w:=2.7/(W-1);
h:=1.87/(H-1);
n:=255.0/N;
FOR x FROM 0 TO W-1 DO
  FOR y FROM 0 TO H-1 DO
   z:=0;
   c:=w*x-2.1-i*(h*y-.935);
   FOR j FROM 0 to N-1 DO
    z:=z*z+c;
    IF ABS(z)>2 THEN
     break;
    END;
   END;
   t:=n*j;
   PIXON_P(x,y,255*20*j+256);
  END;
END;
FREEZE;
WAIT(0);
END;

1min50.81s

NumWorks :

Code: Tout sélectionner

import kandinsky
def mandelbrot(W,H,N) :
w=2.7/(W-1)
h=1.87/(H-1)
n=255.0/N
for x in range(W):
  for y in range(H):
   z=complex(0,0)
   c=complex(w*x-2.1,.935-h*y)
   for j in range(N):
    z=z*z+c
    if abs(z)>2:
     break
   t=n*j
   kandinsky.set_pixel(x,y,255*20*j+256)

1min15.57s

Comme prévu, cela arrange les choses dans le contexte CAS.

[critor]

Petite astuce sur NumWorks, même si ça n'économise pas tout-à-fait une seconde :

Code: Tout sélectionner

from kandinsky import set_pixel
def mandelbrot(W,H,N) :
w=2.7/(W-1)
h=1.87/(H-1)
n=255.0/N
for x in range(W):
  for y in range(H):
   z=complex(0,0)
   c=complex(w*x-2.1,.935-h*y)
   for j in range(N):
    z=z*z+c
    if abs(z)>2:
     break
   t=n*j
   set_pixel(x,y,255*20*j+256)

[parisse]

remarque: le t=n*j est devenu inutile

[critor]

En effet, oubli de ma part. Merci.

[critor]

Allez, voici donc pour référence les performances des versions finales.

HP Prime - CAS+Python :

Code: Tout sélectionner

#cas
def fractal_python(W,H,N) :
local x,y,z,c,j,w,h,n;
w=2.7/(W-1)
h=1.87/(H-1)
n=255.0/N
for x in range(W):
  for y in range(H):
   z=0
   c=w*x-2.1-i*(h*y-.935)
   for j in range(N):
    z:=z*z+c;
    if abs(z)>2:
     break
   PIXON_P(x,y,255*20*j+256)
FREEZE
WAIT(0)
#end

2mins38.41s

HP Prime - CAS+HPPPL :

Code: Tout sélectionner

#cas
fractal_cas(W,H,N)
BEGIN
local x,y,z,c,j,w,h,n;
w:=2.7/(W-1);
h:=1.87/(H-1);
n:=255.0/N;
FOR x FROM 0 TO W-1 DO
  FOR y FROM 0 TO H-1 DO
   z:=0;
   c:=w*x-2.1-i*(h*y-.935);
   FOR j FROM 0 to N-1 DO
    z:=z*z+c;
    IF abs(z)>2 THEN
     break;
    END;
   END;
   PIXON_P(x,y,255*20*j+256);
  END;
END;
FREEZE;
WAIT(0);
END;
#end

2mins37.69s

HP Prime - HPPPL :

Code: Tout sélectionner

EXPORT fractal(W,H,N)
BEGIN
LOCAL x,y,z,c,j,w,h,n;
w:=2.7/(W-1);
h:=1.87/(H-1);
n:=255.0/N;
FOR x FROM 0 TO W-1 DO
  FOR y FROM 0 TO H-1 DO
   z:=0;
   c:=w*x-2.1-i*(h*y-.935);
   FOR j FROM 0 to N-1 DO
    z:=z*z+c;
    IF ABS(z)>2 THEN
     break;
    END;
   END;
   PIXON_P(x,y,255*20*j+256);
  END;
END;
FREEZE;
WAIT(0);
END;

1min37.10s

NumWorks :

Code: Tout sélectionner

from kandinsky import set_pixel
def mandelbrot(W,H,N) :
w=2.7/(W-1)
h=1.87/(H-1)
n=255.0/N
for x in range(W):
  for y in range(H):
   z=complex(0,0)
   c=complex(w*x-2.1,.935-h*y)
   for j in range(N):
    z=z*z+c
    if abs(z)>2:
     break
   set_pixel(x,y,255*20*j+256)

1min13.80s

Dommage que l'on n'arrive pas à chuter sous la minute avec la NumWorks.