const { EPSILON: eps } = Number;

const u = eps/2;
const es = (eps**2)/2;
const ups = u + es;


/** 
 * Subtract one unit in the last place (ulp) from the given number
 * 
 * * subnormal numbers (and 0) are returned unaltered
 * @internal 
 */
function sub1Ulp(n: number) {
    return n > 0 ? n - n*ups : n + n*ups;
}


export { sub1Ulp }