Modified bigvertexgibbs() to be able to handle large beta.
Modified -m15 (temperature-finding mode) to use rolling average.
Modified -m18 to show entire "communication matrix".