Issue with number of k points
Dear libvdwxc team,
running the following script (co_hbn_gr.py) on 16 cores with any k point grid up to 10x10x1 Monkhorst pack (I haven't tried with more points) except for 7x7x1 grid (16 points in the irreducible Brillouin zone) causes the segmentation fault on a first iteration of the self-consistent calculation. With 7x7x1 Monkhorst-Pack the code also runs normally on 32 cores.
The problem could be related to toolchain used to build GPAW (version 1.4.0) with your code (version 0.3.2). The above problems happen when GPAW is built with the Intel 2018 toolchain. With the foss 2018 toolchain, the error changes to numpy.linalg.linalg.LinAlgError about minor of matrix not being positive definite. When the functional is changed to PBE, the error doesn't occur.
Could you try and see if you reproduce this behaviour and, in case you don't, inform us which toolchain you used to build GPAW with?