distance constraint energy calculation in wham