I have a set of k-points, e.g. k1,k2,k3,k4. and they are equivalent by symmetry.
Now I have calculated the momentum matrix element <i|p|j> at k1 point ONLY,
and then calculate the optical properties which, for example, depend on <i|p|j><j|p|i>
I have to make a summation on four k-points...