implement the option 'density_tol' controlling density cutoff to keep numerical stability of LRT