There seems to be an error in the EQeq charges algorithm. The original code:
unitcell = obuc->GetCellMatrix(); fourier = (2 * PI * unitcell.inverse()).transpose();
should be changed to:
unitcell = obuc->GetCellMatrix().transpose(); fourier = (2 * PI * unitcell.inverse());