Hello, I am a new user of Julia and I would like to introduce Julia to researchers and students in the university.
As to give a try, I translated a fortran code to Julia. This code implements the propagation of a 2D acoustic wave using finite differences. The code produces good numerical results, and I am using gadfly to display the results (currently commented) The not so good news is that it is very slow, 100 times slower than fortran. It spends most the time in the loops at line 121. I do not know why it is so slow. I guess that I must be doing something really wrong. (maybe in the array allocations !) Can someone take a look at the code and give me some hints. Thanks Alain
array-acoustic.jl
Description: Binary data
