As of now, BLAS.set_num_threads work in 2 different ways for OpenBLAS and MKL, the former sets both BLAS and LAPACK routines, while the latter only sets BLAS ones. It seems that community expects MKL to also behave like OpenBLAS (JuliaLinearAlgebra/MKL.jl#174, https://discourse.julialang.org/t/correct-way-to-set-mkl-threads-at-runtime/128470) and this seems to be creating ambiguities.

Would it be possible to use mkl_ instead of mkl_domain as the default way, and then wrap the domain specific directly in MKL.jl? This would make things easier for end-users but still allow experts a more fine grained control