Lines Matching refs:adjexp
1819 mpd_ssize_t adjexp, etiny, shift;
1822 adjexp = mpd_adjexp(dec);
1823 if (adjexp > ctx->emax) {
1871 /* At this point adjexp=exp+digits-1 <= emax and exp > etop=emax-prec+1:
1880 if (!mpd_iszerocoeff(dec) && adjexp < ctx->emin) {
1886 else if (adjexp < ctx->emin) {
1901 /* At this point adjexp=exp+digits-1 < emin and exp < etiny=emin-prec+1:
1919 * (1) adjexp=exp+digits-1 < emin
2214 * At this point adjexp(a) == adjexp(b) and a->exp == b->exp,
3363 * adjexp(small) < exp if and only if:
4329 * Lower bound for log10(p / abs(r)): adjexp(p) - (adjexp(r) + 1)
4332 * 2) -p < adjexp(r) <= -1
4405 * MAX-EMAX+1 < log10(e^(0.1*10*t)) <= log10(e^(r*10^t)) < adjexp(e^(r*10^t))+1
4409 * adjexp(e^(r*10^t)) <= log10(e^(r*10^t)) <= log10(e^(-0.1*10^t)) < MIN-ETINY
4450 * 3) log10(9) - prec - t - 1 < log10(abs(r)) < adjexp(abs(r)) + 1
4451 * 4) - prec - t - 2 < adjexp(abs(r)) <= -1
4938 * p = ctx->prec+2-adjexp(lower), then the relative error of
4939 * the result is (using 10**adjexp(x) <= abs(x)):
5014 mpd_ssize_t adjexp, t;
5041 * 1) log10(x) < 0 iff adjexp(x) < 0
5042 * 2) 0 < x /\ x <= y ==> adjexp(x) <= adjexp(y)
5044 * 4) adjexp(x) <= log10(x) < adjexp(x) + 1
5046 * Case adjexp(x) >= 0:
5047 * 5) 2 * adjexp(x) < abs(log(x))
5048 * Case adjexp(x) > 0:
5049 * 6) adjexp(2 * adjexp(x)) <= adjexp(abs(log(x)))
5050 * Case adjexp(x) == 0:
5053 * Case adjexp(x) < 0:
5054 * 7) 2 * (-adjexp(x) - 1) < abs(log(x))
5055 * Case adjexp(x) < -1:
5056 * 8) adjexp(2 * (-adjexp(x) - 1)) <= adjexp(abs(log(x)))
5057 * Case adjexp(x) == -1:
5060 adjexp = mpd_adjexp(a);
5061 t = (adjexp < 0) ? -adjexp-1 : adjexp;
5065 mpd_setspecial(result, (adjexp<0), MPD_INF);
5158 mpd_ssize_t adjexp, t;
5184 adjexp = mpd_adjexp(a);
5185 if (adjexp < 0) {
5187 adjexp = -adjexp;
5189 _settriple(result, sign, adjexp, 0);
5195 * 1) log10(x) < 0 iff adjexp(x) < 0
5196 * 2) 0 < x /\ x <= y ==> adjexp(x) <= adjexp(y)
5197 * 3) adjexp(x) <= log10(x) < adjexp(x) + 1
5199 * Case adjexp(x) >= 0:
5200 * 4) adjexp(x) <= abs(log10(x))
5201 * Case adjexp(x) > 0:
5202 * 5) adjexp(adjexp(x)) <= adjexp(abs(log10(x)))
5203 * Case adjexp(x) == 0:
5206 * Case adjexp(x) < 0:
5207 * 6) -adjexp(x) - 1 < abs(log10(x))
5208 * Case adjexp(x) < -1:
5209 * 7) adjexp(-adjexp(x) - 1) <= adjexp(abs(log(x)))
5210 * Case adjexp(x) == -1:
5213 adjexp = mpd_adjexp(a);
5214 t = (adjexp < 0) ? -adjexp-1 : adjexp;
5217 mpd_setspecial(result, (adjexp<0), MPD_INF);
6685 * 3) adjexp(y)+1 > log10(abs(y)) >= extra
7513 /* s.digits < t.digits <= 2*k+5, |adjexp(s)-adjexp(t)| <= 1,