Added equilibration to main.cpp. Fixed silly bug where I was choosing hbonded bonds to swap, not the opposite. Made HB swaps symmettrical i.e. keep adj[i][j] = -adj[j][i]