modified H calculation to use dual coordinate descent