The electrostatic energy is calculated in the usual way: real space + reciprocal space + self. The only difference is that the reciprocal lattice vectors are set up along the cell unit vectors. The code for this is in energy_routines.f90 Subroutine Ewald_Reciprocal_Lattice_Vector_Setup.
