wd.u.inner.Rd
Williamson-Downs bounds inner
wd.u.inner(y, delt, Y1t, ddid)
the value to obtain the DoTT for