Hi: I am interested in calling a lapack routine, dgesvd, that makes a singular value decomposition of a matrix. The library LinearAlgebra.jl already does it, but it is a generic library that does not take full advantage of the concrete CPU you use. For Intel CPU’s, there is MKL, that can be used …

AOCL Julia packages are now live! o GitHub - amd/AOCL.jl o GitHub - amd/AOCL_jll.jl

[image] aasdelat: function my_dgesvd( JOBU, JOBVT, M, N, A, LDA, S, U, LDU, VT, LDVT, & WORK, LWORK, INFO ) ................. end function my_dgesvd I have to say that there is a mistake in the Fortran code: it should say “subroutine” where it says “function”. Once correc…

ENVIRONMENT First of all, I will tell more about my environment: $ uname -a Linux ryzen-casa 6.5.0-27-generic #28~22.04.1-Ubuntu SMP PREEMPT_DYNAMIC Fri Mar 15 10:51:06 UTC 2 x86_64 x86_64 x86_64 GNU/Linux I have installed the AOCL libraries from the binaries: aocl-linux-aocc-4.2.0_1_amd64.deb. I …

[image] aasdelat: However, from a shell from Julia: Try like this: julia> ENV["LD_LIBRARY_PATH"] "/usr/local/lib/x86_64-unknown-linux-gnu/:/usr/local/lib/:/usr/local/lib/x86_64-unknown-linux-gnu/:/usr/local/lib/" Mind trying LinearAlgebra after using MKL? That’s likely to do better at smal…

[image] Elrod: Try like this: … julia> ENV["LD_LIBRARY_PATH"] "/opt/AMD/aocl/aocl-linux-aocc-4.2.0/aocc/lib:/opt/AMD/aocc-compiler-4.2.0/ompd:/opt/AMD/aocc-compiler-4.2.0/lib:/opt/AMD/aocc-compiler-4.2.0/lib32:/usr/lib/x86_64-linux-gnu:/usr/lib64:/usr/lib32:/usr/lib:" This way shows it!

When I sue the MKL.jl library, it uses more cpu threads, but I set: N_THREADS=30 export JULIA_NUM_THREADS=1 export OMP_NUM_THREADS=$N_THREADS export OPENBLAS_NUM_THREADS=$N_THREADS export MKL_NUM_THREADS=$N_THREADS export VECLIB_MAXIMUM_THREADS=$N_THREADS export NUMEXPR_NUM_THREADS=$N_THREADS The …

About the strange behavior of the dgesvd function, in order to test if it is a Julia or a aocl problem, I have made a fortran program that uses that function, and the result is that it behaves the same: Starts using all threads, then uses only one thread for a long time, and, finally, all threads ag…

There has been a misleading, because I have been using the dgesvd while I have should be using dgesdd. As stated in the LAPACK documentation for the singular value decomposition : There are two types of driver routines for the SVD. Originally LAPACK had just the simple driver described below, and t…

Pls contact toolchainsupport@amd.com to request for AMD optimized libraries.

new here in julia. i have tested AOCL, Openblas, intelMKL in my ryzen 3900x home computer. this the benchmark result by fortran and juliat version. as you can see AOCL actualy quite good. but still mkl wins eventhough in AMD. julia blas mkl is faster in dgetrf/dgetri than fortran but same performanc…

I contacted toolchainsupport@amd.com one week ago, but I still haven’t received any answer.

AOCL (not MKL) acceleration on AMD Ryzen CPU's

General Usage Performance

agodbole August 10, 2026, 10:54am 30

AOCL.jl is now registered (under the JuliaLinearAlgebra organization).

refer: AOCL.jl v5.3.1

Topic		Replies	Views
OpenBLAS is faster than Intel MKL on AMD Hardware (Ryzen) Performance blas , lapack	40	37588	June 19, 2020
Acceleration of Intel MKL on AMD Ryzen CPU's Performance performance , mkl , linearalgebra	34	7692	May 9, 2024
LU factorization performance issue New to Julia linearalgebra	30	1044	June 6, 2022
How to call blas getrf, getri properly ? i want to create a benchmark inverse matrix using gauss, crout, native julia inv(A) and BLAS direct New to Julia question , blas , mkl , benchmark , openblas	5	270	May 7, 2025
SVD 2x slower than in Matlab and how to get best performance on Windows10 Performance	20	3791	October 9, 2024

AOCL (not MKL) acceleration on AMD Ryzen CPU's

Related topics