# Linear solver \\(A, B) performance vs Matlab A\\b

**URL:** <https://discourse.julialang.org/t/linear-solver-a-b-performance-vs-matlab-a-b/2082>\
**Category:** General Usage\
**Created:** [February 13, 2017, 2:23pm UTC](https://discourse.julialang.org/t/linear-solver-a-b-performance-vs-matlab-a-b/2082 "2017-02-13T14:23:03Z")\
**Posts on this page:** 1\
**Showing post:** 12

<div class="post-metadata">

**Author:** ![Tecnezio](https://avatars.discourse-cdn.com/v4/letter/t/ed655f/32.png) [@Tecnezio](https://discourse.julialang.org/u/Tecnezio)\
**Post date:** [February 14, 2017, 2:53pm UTC](https://discourse.julialang.org/t/linear-solver-a-b-performance-vs-matlab-a-b/2082/12 "2017-02-14T14:53:16Z")

</div>

I’ve done it !  
Here the post about my test:

> [@Call Intel MKL as external library](https://discourse.julialang.org/t/call-intel-mkl-as-external-library/2097/3):
>
> No, i haven’t set anything… But i loaded the wrong version of the script… and now it work! Here the working version: # Test MKL solve linear systems N = 10000; A = rand(N,N); AA = copy(A); b = rand(N,1); c = copy(b); function juliaSol(A,b) \(A,b); end @time sol = juliaSol(A,b); # Test MKL const global librt = Libdl.find\_library(["libmkl\_rt"], ["/opt/intel/mkl/lib"]) # Open librt Libdl.dlopen(librt) # LU FACORIZATION function luFactMKL(A::StridedMatrix{Float64}) m, n = size(A) lda = …

---

_[View the full topic](https://discourse.julialang.org/t/linear-solver-a-b-performance-vs-matlab-a-b/2082)._
