Stick const in front of T and RHS. On Sat, Apr 25, 2015 at 11:32 AM, Tim Holy <[email protected]> wrote:
> Did you read through > http://docs.julialang.org/en/release-0.3/manual/performance-tips/? You > should > memorize :-) the sections up through the Tools section; the rest you can > consult as you discover you need them. > > --Tim > > On Saturday, April 25, 2015 01:03:38 AM Ángel de Vicente wrote: > > 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 > >
