Interference2DLattice becomes very slow for short decay lengths
Function Interference2DLattice::interferenceForXi computes a double sum in q space. The number of terms under the sum is inverse proportional to the x and y decay lengths. Therefore, for small decay lengths, the simulation becomes unacceptably slow.
Thanks to Sebastian Köhler from Lund University for reporting.