d0p1/ack
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

improved conversion of extended to string

Browse source
This commit is contained in:
ceriel 1990-07-20 09:39:43 +00:00
parent b23125354b
commit cf28f45a59
2 changed files with 66 additions and 17 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

41
lang/cem/libcc.ansi/stdlib/ext_comp.c
View file

@ -548,6 +548,7 @@ _ext_str_cvt(struct EXTEND *e, int ndigit, int *decpt, int *sign, int ecvtflag)
static char buf[NDIGITS+1]; static char buf[NDIGITS+1];
register char *p = buf; register char *p = buf;
register char *pe; register char *pe;
int findex = 0;
if (ndigit < 0) ndigit = 0; if (ndigit < 0) ndigit = 0;
if (ndigit > NDIGITS) ndigit = NDIGITS; if (ndigit > NDIGITS) ndigit = NDIGITS;
@ -566,27 +567,31 @@ _ext_str_cvt(struct EXTEND *e, int ndigit, int *decpt, int *sign, int ecvtflag)
while(cmp_ext(e,pp) >= 0) pp++; while(cmp_ext(e,pp) >= 0) pp++;
pp--; pp--;
mul_ext(e,&r_big_ten_powers[pp-big_ten_powers],e); findex = pp - big_ten_powers;
*decpt += (pp - big_ten_powers) * TP; mul_ext(e,&r_big_ten_powers[findex],e);
*decpt += findex * TP;
pp = &ten_powers[1]; pp = &ten_powers[1];
while(pp < &ten_powers[TP] && cmp_ext(e, pp) >= 0) pp++; while(pp < &ten_powers[TP] && cmp_ext(e, pp) >= 0) pp++;
pp--; pp--;
mul_ext(e, &r_ten_powers[pp-ten_powers], e); findex = pp - ten_powers;
*decpt += pp - ten_powers; *decpt += findex;
if (cmp_ext(e, &ten_powers[0]) < 0) { if (cmp_ext(e, &ten_powers[0]) < 0) {
pp = &r_big_ten_powers[1]; pp = &r_big_ten_powers[1];
while(cmp_ext(e,pp) < 0) pp++; while(cmp_ext(e,pp) < 0) pp++;
pp--; pp--;
mul_ext(e, &big_ten_powers[pp - r_big_ten_powers], e); findex = pp - r_big_ten_powers;
*decpt -= (pp - r_big_ten_powers) * TP; mul_ext(e, &big_ten_powers[findex], e);
*decpt -= findex * TP;
/* here, value >= 10 ** -28 */ /* here, value >= 10 ** -28 */
mul_ext(e, &ten_powers[1], e); mul_ext(e, &ten_powers[1], e);
(*decpt)--; (*decpt)--;
pp = &r_ten_powers[0]; pp = &r_ten_powers[0];
while(cmp_ext(e, pp) < 0) pp++; while(cmp_ext(e, pp) < 0) pp++;
mul_ext(e, &ten_powers[pp - r_ten_powers], e); findex = pp - r_ten_powers;
*decpt -= pp - r_ten_powers; mul_ext(e, &ten_powers[findex], e);
*decpt -= findex;
findex = 0;
} }
(*decpt)++; /* because now value in [1.0, 10.0) */ (*decpt)++; /* because now value in [1.0, 10.0) */
} }
@ -596,6 +601,26 @@ _ext_str_cvt(struct EXTEND *e, int ndigit, int *decpt, int *sign, int ecvtflag)
if (pe > &buf[NDIGITS]) pe = &buf[NDIGITS]; if (pe > &buf[NDIGITS]) pe = &buf[NDIGITS];
} }
while (p <= pe) { while (p <= pe) {
if (findex) {
struct EXTEND tc, oldtc;
int count = 0;
oldtc.exp = 0;
oldtc.sign = 0;
oldtc.m1 = 0;
oldtc.m2 = 0;
tc = ten_powers[findex];
while (cmp_ext(e, &tc) >= 0) {
oldtc = tc;
add_ext(&tc, &ten_powers[findex], &tc);
count++;
}
*p++ = count + '0';
oldtc.sign = 1;
add_ext(e, &oldtc, e);
findex--;
continue;
}
if (e->exp >= 0 && e->m1 != 0) { if (e->exp >= 0 && e->m1 != 0) {
struct EXTEND x; struct EXTEND x;

42
lang/cem/libcc/gen/ext_comp.c
View file

@ -554,6 +554,7 @@ _ext_str_cvt(e, ndigit, decpt, sign, ecvtflag)
static char buf[NDIGITS+1]; static char buf[NDIGITS+1];
register char *p = buf; register char *p = buf;
register char *pe; register char *pe;
int findex = 0;
if (ndigit < 0) ndigit = 0; if (ndigit < 0) ndigit = 0;
if (ndigit > NDIGITS) ndigit = NDIGITS; if (ndigit > NDIGITS) ndigit = NDIGITS;
@ -572,30 +573,33 @@ _ext_str_cvt(e, ndigit, decpt, sign, ecvtflag)
while(cmp_ext(e,pp) >= 0) pp++; while(cmp_ext(e,pp) >= 0) pp++;
pp--; pp--;
mul_ext(e,&r_big_ten_powers[pp-big_ten_powers],e); findex = pp - big_ten_powers;
*decpt += (pp - big_ten_powers)* mul_ext(e,&r_big_ten_powers[findex],e);
(sizeof(ten_powers)/sizeof(ten_powers[0])); *decpt += findex * (sizeof(ten_powers)/sizeof(ten_powers[0]));
pp = &ten_powers[1]; pp = &ten_powers[1];
while(pp<&ten_powers[(sizeof(ten_powers)/sizeof(ten_powers[0]))] && while(pp<&ten_powers[(sizeof(ten_powers)/sizeof(ten_powers[0]))] &&
cmp_ext(e, pp) >= 0) pp++; cmp_ext(e, pp) >= 0) pp++;
pp--; pp--;
mul_ext(e, &r_ten_powers[pp-ten_powers], e); findex = pp - ten_powers;
*decpt += pp - ten_powers; *decpt += findex;
if (cmp_ext(e, &ten_powers[0]) < 0) { if (cmp_ext(e, &ten_powers[0]) < 0) {
pp = &r_big_ten_powers[1]; pp = &r_big_ten_powers[1];
while(cmp_ext(e,pp) < 0) pp++; while(cmp_ext(e,pp) < 0) pp++;
pp--; pp--;
mul_ext(e,&big_ten_powers[pp-r_big_ten_powers],e); findex = pp - r_big_ten_powers;
*decpt -= (pp - r_big_ten_powers)* mul_ext(e,&big_ten_powers[findex],e);
*decpt -= findex *
(sizeof(ten_powers)/sizeof(ten_powers[0])); (sizeof(ten_powers)/sizeof(ten_powers[0]));
/* here, value >= 10 ** -28 */ /* here, value >= 10 ** -28 */
mul_ext(e, &ten_powers[1], e); mul_ext(e, &ten_powers[1], e);
(*decpt)--; (*decpt)--;
pp = &r_ten_powers[0]; pp = &r_ten_powers[0];
while(cmp_ext(e, pp) < 0) pp++; while(cmp_ext(e, pp) < 0) pp++;
mul_ext(e, &ten_powers[pp-r_ten_powers], e); findex = -(pp - r_ten_powers);
*decpt -= pp - r_ten_powers; mul_ext(e, &ten_powers[-findex], e);
*decpt += findex;
findex = 0;
} }
(*decpt)++; /* because now value in [1.0, 10.0) */ (*decpt)++; /* because now value in [1.0, 10.0) */
} }
@ -605,6 +609,26 @@ _ext_str_cvt(e, ndigit, decpt, sign, ecvtflag)
if (pe > &buf[NDIGITS]) pe = &buf[NDIGITS]; if (pe > &buf[NDIGITS]) pe = &buf[NDIGITS];
} }
while (p <= pe) { while (p <= pe) {
if (findex) {
struct EXTEND tc, oldtc;
int count = 0;
oldtc.exp = 0;
oldtc.sign = 0;
oldtc.m1 = 0;
oldtc.m2 = 0;
tc = ten_powers[findex];
while (cmp_ext(e, &tc) >= 0) {
oldtc = tc;
add_ext(&tc, &ten_powers[findex], &tc);
count++;
}
*p++ = count + '0';
oldtc.sign = 1;
add_ext(e, &oldtc, e);
findex--;
continue;
}
if (e->exp >= 0 && e->m1 != 0) { if (e->exp >= 0 && e->m1 != 0) {
struct EXTEND x; struct EXTEND x;