Open StevenCHowell opened 8 years ago
Calculate it before the first move, store the result, then calculate it after each subsequent move. That way you avoid doing a N^2 loop twice.
certainly, probably should store it in mvars
This will also require adding an input for the ionic conditions of the buffer.
This stills need to be added and tested. The planned design pattern is the screened Debye energy will be calculated before and after performing MC moves. Then accept new moves based on the Boltzmann factor, based on the energy change.