Computing the Fourier transform of a distribution of electron charge (represented by atomic positions in space)
The most immediate optimization is to use coalesced reads by switching to array of structs form to structure of arrays.