Abstract
An efficient calculational algorithm is provided for the operator decomposition approach to non-adiabatic Berry phases for systems with periodic Hamiltonians.