I’m analyzing the performance of julia linear solver (A,B) vs the Matlab solver.

Julia is fast only on the second call of the solver and only for small (A) matrix. Over size of A 500x500, julia solver is always slower than Matlab.

Anyone did a similar performance comparison ?

I looked for a chart of the julia’s solver (A,B) but i found nothing.

Which file of julia source contain the solver algorithm ? I can reconstruct a flow chart of the algorithm.

My testing platform is :

Julia Version 0.5.0

Commit 3c9d753 (2016-09-19 18:14 UTC)

Platform Info:

System: NT (x86_64-w64-mingw32)

CPU: Intel® Core™ i7-7500U CPU @ 2.70GHz

WORD_SIZE: 64

BLAS: libopenblas (USE64BITINT DYNAMIC_ARCH NO_AFFINITY Prescott)

LAPACK: libopenblas64_

LIBM: libopenlibm

LLVM: libLLVM-3.7.1 (ORCJIT, broadwell)