******** PBC SCF flags ********
N kpts = 64
Exchange divergence treatment (exxdiv) = ewald
    madelung (= occupied orbital energy shift) = 0.1052367084460809
    Total energy shift due to Ewald probe charge = -1/2 * Nelec*madelung = -1.68378733514
DF object = <pyscf.pbc.df.fft.FFTDF object at 0x7b069b54d550>
XC functionals = pbe
small_rho_cutoff = 1e-07
Uniform grid, mesh = [65 65 65]
Traceback (most recent call last):
  File "/home/chillenb/src/VASPsol/work/cuslab/usersjk.py", line 25, in <module>
    mf.run()
    ~~~~~~^^
  File "/home/chillenb/miniforge3/envs/condapyscf/lib/python3.13/site-packages/pyscf/lib/misc.py", line 638, in run
    self.kernel(*args)
    ~~~~~~~~~~~^^^^^^^
  File "<string>", line 2, in kernel
  File "/home/chillenb/miniforge3/envs/condapyscf/lib/python3.13/site-packages/pyscf/scf/hf.py", line 2009, in scf
    self.build(self.mol)
    ~~~~~~~~~~^^^^^^^^^^
  File "/home/chillenb/miniforge3/envs/condapyscf/lib/python3.13/site-packages/pyscf/pbc/dft/rks.py", line 308, in build
    if not numpy.all(self.rsjk.kpts == self.kpt):
                                       ^^^^^^^^
  File "/home/chillenb/miniforge3/envs/condapyscf/lib/python3.13/site-packages/pyscf/pbc/scf/hf.py", line 539, in kpt
    return self.with_df.kpts.reshape(3)
           ~~~~~~~~~~~~~~~~~~~~~~~~~^^^
ValueError: cannot reshape array of size 192 into shape (3,)