Hi,

a complete Julia newbie here... I spent a couple of days learning the 
syntax and main aspects of Julia, and since I heard many good things about 
it, I decided to try a little program to see how it compares against the 
other ones I regularly use: Fortran and Python.

I wrote a minimal program to solve the 3D heat equation in a cube of 
100x100x100 points in the three languages and the time it takes to run in 
each one is:

Fortran: ~7s
Python: ~33s
Julia:    ~80s

The code runs for 1000 iterations, and I'm being nice to Julia, since the 
programs in Fortran and Python write 100 HDF5 files with the complete 100^3 
data (every 10 iterations).

I attach the code (and you can also get it at: http://pastebin.com/y5HnbWQ1)

Am I doing something obviously wrong? Any suggestions on how to improve its 
speed?

Thanks a lot,
Ángel de Vicente

Attachment: heat3D.jl
Description: Binary data

Reply via email to