Cleaning up manual, .Machine$double.xmax new detault for non-relative
upper bound.
Changed bw.switch to 100 since we are now exclusively in scale factor
mode (now this is scaled by n but regardless 100 sd n^{-1/(2p+q)} is
`large' though I might want to revisit this)
Scaled bw.switch to truly represent a 100 sigma smoothing