Mathemagix

[Bernard Mourrain <address@concealed>]

Hello,

Attached is a classical picture of a Mandelbrot set (400x400 pixel)
obtained with shape+axel. For each point of a grid, a color is computed 
according to the number of iterations of a polynomial map staying in a 
given disk starting from this point
See code shape/expl/mandelbrot.cpp or shape/mmx/mandelbrot.mmx.

Below is a comparison of the "same" code
(for an image of 200x200 pixels) in different environments:

time ./expl/mandelbrot (c++)

real    0m0.078s
user    0m0.017s
sys     0m0.014s

time python mmx/shape/mmx/mandelbrot.py

real    0m0.484s
user    0m0.467s
sys     0m0.012s


time maple mmx/shape/mmx/mandelbrot.mpl

real    0m7.884s
user    0m0.051s
sys     0m0.084s

time mmx-light mmx/shape/mmx/mandelbrot.mmx    (see file attached) 


real    0m21.581s
user    0m21.444s
sys     0m0.028s


Why is there such a huge difference for mmx-light ?
where is the weakness ? Any way to improve it ?

My best wishes of expansion and efficiency for mmx in 2010.

Bernard

Attachment: mandelbrot.mmx
Description: application/mmx