relax somewhat the check for non-pd theta/eta using a tolerance of .Machine$double.eps^(3/4)