[Note/suggestions] Got Arpack, etc sort of working with MKL [tested on macos]


Update: With fflags from here and borrowing from the Arpack builder script, here’s a way to potentially get Arpack and packages that depend on it working with MKL( Pkg.test("Arpack") , Pkg.test("LightGraphs") or Pkg.test("Distributions") all pass) [basically, use the suggestions here to provide a arpack binary compiled with mkl to the julia package]:

  1. Clone arpack-ng, checkout "b095052372aa95d4281a645ee1e367c28255c947" according to the script used by Arpack.jl’s builder
  2. Build arpack library binary with mkl: mkdir build, cd to build folder and use the script (64 bit):
#!/usr/bin/env bash
source /opt/intel/mkl/bin/mklvars.sh intel64 ilp64


# Symbols that have float32, float64, complexf32, and complexf64 support
SDCZ_SYMBOLS="axpy copy gemv geqr2 lacpy lahqr lanhs larnv lartg \
              lascl laset scal trevc trmm trsen gbmv gbtrf gbtrs \
              gttrf gttrs pttrf pttrs"
# All symbols that have float32/float64 support (including the SDCZ_SYMBOLS above)
SD_SYMBOLS="${SDCZ_SYMBOLS} dot ger labad laev2 lamch lanst lanv2 \
            lapy2 larf larfg lasr nrm2 orm2r rot steqr swap"
# All symbols that have complexf32/complexf64 support (including the SDCZ_SYMBOLS above)
CZ_SYMBOLS="${SDCZ_SYMBOLS} dotc geru unm2r"
# Add in (s|d)*_64 symbol remappings:
for sym in ${SD_SYMBOLS}; do
    SYMBOL_DEFS="${SYMBOL_DEFS} -Ds${sym}=s${sym}_64 -Dd${sym}=d${sym}_64"
# Add in (c|z)*_64 symbol remappings:
for sym in ${CZ_SYMBOLS}; do
    SYMBOL_DEFS="${SYMBOL_DEFS} -Dc${sym}=c${sym}_64 -Dz${sym}=z${sym}_64"
# Add one-off symbol mappings; things that don't fit into any other bucket:
for sym in scnrm2 dznrm2 csscal zdscal dgetrf dgetrs; do
    SYMBOL_DEFS="${SYMBOL_DEFS} -D${sym}=${sym}_64"
# Set up not only lowercase symbol remappings, but uppercase as well:

FFLAGS="${FFLAGS} -fdefault-integer-8 ${SYMBOL_DEFS} -ff2c"


Followed by make clean all; The built library libarpack.2.0.0.dylib is in the build/lib folder. Then, copy and overwrite the libarpack.2.0.0.dylib in the arpack’s folder in .julia, e.g. current version lives in:


Build it in julia:

(v1.1) pkg> build Arpack

then, e.g.

using Pkg, Arpack, LinearAlgebra
A = Diagonal(1:4)
eigs(A, nev=3)

works, and all Pkg.test("Arpack") , Pkg.test("LightGraphs") , Pkg.test("Distributions") pass.

Now the tests are working with mkl, this whole process may be cleaned up a bit. May be working on other platforms? Do not have access to linux/windows to test. If the tests pass maybe it’s relatively safe to use? Out of caution, test your own code with official binary without mkl and compare. Note that to use Arpack in the official binary may need removing the self-build arpack binary in e.g. ~/.julia/packages/Arpack/UiiMc/deps/usr/lib/libarpack.2.0.0.dylib (official non-mil binary should download the missing binary automatically), and rebuild build Arpack in pkg.

MKL and libm compile notes: macOS