QR decomposition in julia 2x slower than octave

Not so for Matlab:

>> A=rand(5000, 5000);
>> tic; qr(A); toc
Elapsed time is 4.189508 seconds.
>> version

ans =

    '9.4.0.813654 (R2018a)'

>> tic; [q,r]=qr(A); toc
Elapsed time is 5.785541 seconds.