Hi All,

I was comparing the performance of Julia and Matlab in solving systems of 
linear equations and I got the following

Matlab:
>> A = randn(1000); A = A+A'+eye(1000); x = randn(1000,1);
>> tic; A\x; toc
Elapsed time is 0.034763 seconds.

Julia:
julia> A = randn(1000,1000); A = A+A'+eye(1000); x = randn(1000);

julia> @time A\x;
elapsed time: 0.192572124 seconds (8012424 bytes allocated)

I was wondering whether there is anything I can do to improve the 
performances of Julia in solving this kind of problems?

many thanks!
andrea

Reply via email to