big fix of barrier hessian calculation

Should remove nan and be more accurate, but now it seems that
convergence breaking doesn’t work.