`cupy.linalg.eigh` may be a drop-in replacement for `numpy.linalg.eigh` and yield big boost: https://docs.cupy.dev/en/stable/reference/generated/cupy.linalg.eigh.html @braceal