Space group symmetry is exploited and implemented in density functional calculations of extended systems with periodic boundary conditions. Our scheme for reducing the number of two-electron integrals employs the entire set of operations of the space group, including glide plains and screw axes. Speedups observed for the Fock matrix formation in simple 3D systems range from 2X to 9X for the near field Coulomb part and from 3X to 8X for the Hartree–Fock-type exchange, the slowest steps of the procedure, thus leading to a substantial reduction of the computational time. The relatively small speedup factors in special cases are attributed to the highly symmetric positions atoms occupy in crystals, including the ones tested here, as well as to the choice of the smallest possible unit cells. For quasi-1D systems with most atoms staying invariant only under identity, the speedup factors often exceed one order of magnitude reaching almost 70X (near-field Coulomb) and 57X (HFx) for the largest tested (16,7) single-walled nanotube with 278 symmetry operations.