Dlh upperbound optimization

Have optimized the upper bound on DLH by 3 to 4 times. The last two commits titled "Optimized x(1-x/c)" and "Optimized upper bound" are to be merged.
Maybe I'm misreading this, but it looks like the last two terms both have divisor `(wdf_lower + 0.5)`.