So what are the results with libxsmm and current AMD BLAS, as that must be for small dimensions? The reason it's serial BLAS that mainly matters is that HPC co

So what are the results with libxsmm and current AMD BLAS, as that must be for small dimensions? The reason it's serial BLAS that mainly matters is that HPC codes are usually parallelized above the BLAS; why do you want the nesting? Swapping in threaded OpenBLAS or BLIS is something you might do with basically serial stuff like vanilla R, e.g. https://loveshack.fedorapeople.org/blas-subversion.html#_addendum_example https://loveshack.fedorapeople.org/blas-subversion.html#_add... OpenBLAS threading has been somewhat buggy, but the main problem with its OpenMP support currently seems to be that using OMP_PLACES kills it.