Written by Alex Ganose @utf and Janosh Riebesell @janosh. Published 2022-03-28. Last updated 2024-03-30.
-
Install Xcode command line tools:
xcode-select --install
-
Install dependencies using Homebrew:
brew install gcc openmpi scalapack fftw qd openblas
Optionally, add
hdf5
for HDF5 support in VASP. -
Compile VASP:
These instructions are for VASP 6.4.1 but should work with minor adjustments for other versions.
cd /path/to/vasp-6.x.y cp arch/makefile.include.gnu_omp makefile.include
Edit
makefile.include
in the VASPsrc
directory:-
Add to
CPP_OPTIONS
:-D_OPENMP \ -Dqd_emulate
-
Change all instances of
gcc
togcc-13
andg++
tog++-13
-
Add after
LLIBS = -lstdc++
to emulate quad precision:QD ?= /opt/homebrew/ LLIBS += -L$(QD)/lib -lqdmod -lqd INCS += -I$(QD)/include/qd
-
Set
SCALAPACK_ROOT ?= /opt/homebrew
-
Set
OPENBLAS_ROOT ?= /opt/homebrew/Cellar/openblas/0.3.20
(Double check this is the path on your system) -
Set
FFTW_ROOT ?= /opt/homebrew
-
(optional but recommended by VASP) For HDF5 support, add
CPP_OPTIONS+= -DVASP_HDF5 HDF5_ROOT ?= /opt/homebrew/ LLIBS += -L$(HDF5_ROOT)/lib -lhdf5_fortran INCS += -I$(HDF5_ROOT)/include
-
Append
getshmem.o
toOBJECTS_LIB
inmakefile.include
(VASP wiki)- OBJECTS_LIB = linpack_double.o + OBJECTS_LIB = linpack_double.o getshmem.o
-
In
src/parser/makefile
, change (as noted by @zhuligs):- ar vq libparser.a $(CPPOBJ_PARS) $(COBJ_PARS) locproj.tab.h + ar vq libparser.a $(CPPOBJ_PARS) $(COBJ_PARS)
Do not replace the tab at the beginning of the line with spaces!
-
In
src/lib/getshmem.c
, add the line#define SHM_NORESERVE 0
(VASP forum):/*output: shmem id */ #define SHM_NORESERVE 0 // this line was added void getshmem_C(size_t _size, int *_id)
-
In
makefile.include
, update the parser section (VASP forum):# For the parser library CXX_PARS = g++-13 - LLIBS = -lstdc++ + LIBS += parser + LLIBS = -Lparser -lparser -lstdc++ QD ?= /opt/homebrew LLIBS += -L$(QD)/lib -lqdmod -lqd INCS += -I$(QD)/include/qd
-
-
Finally, run:
make all
If a previous compilation failed, remember to run
make veryclean
to start from a clean slate. Fixesgfortran
errors likeFatal Error:
string.mod
not found
If successful, the VASP binaries will be in src/bin
. Test with make test
.
Confirmed working with VASP 6.4.1 on M1 Pro with Sonoma 14.2.1 and [email protected]
.
See makefile.include
.
Initial performance testing suggests optimal parameters for the M1 Pro chip with 8 performance, 2 efficiency cores and 16" MBP cooling are
export OMP_NUM_THREADS=1 # important
mpiexec -np 8 vasp_std
NCORE = 4 # in INCAR
n_proc | n_threads | n_core | elapsed (sec) | |
---|---|---|---|---|
0 | 1 | 1 | 2 | 93.3 |
1 | 1 | 1 | 4 | 92.8 |
2 | 1 | 2 | 2 | 82.8 |
3 | 1 | 2 | 4 | 82.7 |
4 | 2 | 1 | 2 | 42.8 |
5 | 2 | 1 | 4 | 42.9 |
6 | 2 | 2 | 2 | 52.9 |
7 | 2 | 2 | 4 | 52.7 |
8 | 4 | 1 | 2 | 32.9 |
9 | 4 | 1 | 4 | 32.9 |
10 | 4 | 2 | 2 | 52.9 |
11 | 4 | 2 | 4 | 53.0 |
12 | 8 | 1 | 2 | 32.8 |
13 | 8 | 1 | 4 | 22.8 |
14 | 8 | 2 | 2 | 62.8 |
15 | 8 | 2 | 4 | 62.9 |
Brings wall time for this Si2 relaxation down to 23 sec.
from time import perf_counter
from atomate2.vasp.jobs.core import RelaxMaker
from atomate2.vasp.powerups import update_user_incar_settings
from jobflow import run_locally
from pymatgen.core import Structure
start = perf_counter()
# FCC silicon structure
si_structure = Structure(
lattice=[[0, 2.73, 2.73], [2.73, 0, 2.73], [2.73, 2.73, 0]],
species=["Si", "Si"],
coords=[[0, 0, 0], [0.25, 0.25, 0.25]],
)
# relax job to optimize structure
relax_job = RelaxMaker().make(si_structure)
relax_job = update_user_incar_settings(relax_job, {"NCORE": 4})
# run job
run_locally(relax_job, create_folders=True, ensure_success=True)
print(f"Si relaxation took {perf_counter() - start:.3f} sec")
Having this same problem on my M2 Mac, appreciate any help!