#define EXTERN #include "l.h" #include #ifndef DEFAULT #define DEFAULT '9' #endif char *noname = ""; char symname[] = SYMDEF; char thechar = '8'; char *thestring = "386"; /* * -H0 -T0x40004C -D0x10000000 is garbage unix * -H1 -T0xd0 -R4 is unix coff * -H2 -T4128 -R4096 is plan9 format * -H3 -Tx -Rx is MS-DOS .COM * -H4 -Tx -Rx is fake MS-DOS .EXE */ void main(int argc, char *argv[]) { int i, c; char *a; Binit(&bso, 1, OWRITE); cout = -1; listinit(); memset(debug, 0, sizeof(debug)); nerrors = 0; outfile = "8.out"; HEADTYPE = -1; INITTEXT = -1; INITDAT = -1; INITRND = -1; INITENTRY = 0; ARGBEGIN { default: c = ARGC(); if(c >= 0 && c < sizeof(debug)) debug[c]++; break; case 'o': /* output to (next arg) */ outfile = ARGF(); break; case 'E': a = ARGF(); if(a) INITENTRY = a; break; case 'H': a = ARGF(); if(a) HEADTYPE = atolwhex(a); break; case 'T': a = ARGF(); if(a) INITTEXT = atolwhex(a); break; case 'D': a = ARGF(); if(a) INITDAT = atolwhex(a); break; case 'R': a = ARGF(); if(a) INITRND = atolwhex(a); break; } ARGEND USED(argc); if(*argv == 0) { diag("usage: 8l [-options] objects\n"); errorexit(); } if(!debug['9'] && !debug['U'] && !debug['B']) debug[DEFAULT] = 1; if(HEADTYPE == -1) { if(debug['U']) HEADTYPE = 1; if(debug['B']) HEADTYPE = 2; if(debug['9']) HEADTYPE = 2; } switch(HEADTYPE) { default: diag("unknown -H option"); errorexit(); case 0: /* this is garbage */ HEADR = 20L+56L; if(INITTEXT == -1) INITTEXT = 0x40004CL; if(INITDAT == -1) INITDAT = 0x10000000L; if(INITRND == -1) INITRND = 0; break; case 1: /* is unix coff */ HEADR = 0xd0L; if(INITTEXT == -1) INITTEXT = 0xd0; if(INITDAT == -1) INITDAT = 0x400000; if(INITRND == -1) INITRND = 0; break; case 2: /* plan 9 */ HEADR = 32L; if(INITTEXT == -1) INITTEXT = 4096+32; if(INITDAT == -1) INITDAT = 0; if(INITRND == -1) INITRND = 4096; break; case 3: /* MS-DOS .COM */ HEADR = 0; if(INITTEXT == -1) INITTEXT = 0x0100; if(INITDAT == -1) INITDAT = 0; if(INITRND == -1) INITRND = 4; break; case 4: /* fake MS-DOS .EXE */ HEADR = 0x200; if(INITTEXT == -1) INITTEXT = 0x0100; if(INITDAT == -1) INITDAT = 0; if(INITRND == -1) INITRND = 4; HEADR += (INITTEXT & 0xFFFF); if(debug['v']) Bprint(&bso, "HEADR = 0x%ld\n", HEADR); break; } if(INITDAT != 0 && INITRND != 0) print("warning: -D0x%lux is ignored because of -R0x%lux\n", INITDAT, INITRND); if(debug['v']) Bprint(&bso, "HEADER = -H0x%ld -T0x%lux -D0x%lux -R0x%lux\n", HEADTYPE, INITTEXT, INITDAT, INITRND); Bflush(&bso); for(i=1; optab[i].as; i++) if(i != optab[i].as) { diag("phase error in optab: %d\n", i); errorexit(); } maxop = i; for(i=0; i= D_AL && i <= D_BH) reg[i] = (i-D_AL) & 7; if(i >= D_AX && i <= D_DI) reg[i] = (i-D_AX) & 7; if(i >= D_F0 && i <= D_F0+7) reg[i] = (i-D_F0) & 7; } zprg.link = P; zprg.pcond = P; zprg.back = 2; zprg.as = AGOK; zprg.from.type = D_NONE; zprg.from.index = D_NONE; zprg.from.scale = 1; zprg.to = zprg.from; pcstr = "%.6lux "; nuxiinit(); histgen = 0; textp = P; datap = P; edatap = P; pc = 0; dtype = 4; cout = create(outfile, 1, 0775); if(cout < 0) { diag("cannot create %s\n", outfile); errorexit(); } version = 0; cbp = buf.cbuf; cbc = sizeof(buf.cbuf); firstp = prg(); lastp = firstp; if(INITENTRY == 0) { INITENTRY = "_main"; if(debug['p']) INITENTRY = "_mainp"; if(!debug['l']) lookup(INITENTRY, 0)->type = SXREF; } else lookup(INITENTRY, 0)->type = SXREF; while(*argv) objfile(*argv++); if(!debug['l']) loadlib(); firstp = firstp->link; if(firstp == P) errorexit(); patch(); follow(); dodata(); dostkoff(); if(debug['p']) if(debug['1']) doprof1(); else doprof2(); span(); doinit(); asmb(); undef(); if(debug['v']) { Bprint(&bso, "%5.2f cpu time\n", cputime()); Bprint(&bso, "%ld symbols\n", nsymbol); Bprint(&bso, "%ld memory used\n", thunk); Bprint(&bso, "%d sizeof adr\n", sizeof(Adr)); Bprint(&bso, "%d sizeof prog\n", sizeof(Prog)); } Bflush(&bso); errorexit(); } void loadlib(void) { int i; long h; Sym *s; loop: xrefresolv = 0; for(i=0; ilink) if(s->type == SXREF) goto loop; } void errorexit(void) { if(nerrors) { if(cout >= 0) remove(outfile); exits("error"); } exits(0); } void objfile(char *file) { long off, esym, cnt, l; int f, work; Sym *s; char magbuf[SARMAG]; char name[100], pname[150]; struct ar_hdr arhdr; char *e, *start, *stop; if(file[0] == '-' && file[1] == 'l') { if(debug['9']) sprint(name, "/%s/lib/lib", thestring); else sprint(name, "/usr/%clib/lib", thechar); strcat(name, file+2); strcat(name, ".a"); file = name; } if(debug['v']) Bprint(&bso, "%5.2f ldobj: %s\n", cputime(), file); Bflush(&bso); f = open(file, 0); if(f < 0) { diag("cannot open file: %s\n", file); errorexit(); } l = read(f, magbuf, SARMAG); if(l != SARMAG || strncmp(magbuf, ARMAG, SARMAG)){ /* load it as a regular file */ l = seek(f, 0L, 2); seek(f, 0L, 0); ldobj(f, l, file); close(f); return; } l = read(f, &arhdr, SAR_HDR); if(l != SAR_HDR) { diag("%s: short read on archive file symbol header\n", file); goto out; } if(strncmp(arhdr.name, symname, strlen(symname))) { diag("%s: first entry not symbol header\n", file); goto out; } esym = SARMAG + SAR_HDR + atolwhex(arhdr.size); off = SARMAG + SAR_HDR; /* * just bang the whole symbol file into memory */ seek(f, off, 0); cnt = esym - off; start = malloc(cnt + 10); cnt = read(f, start, cnt); if(cnt <= 0){ close(f); return; } stop = &start[cnt]; memset(stop, 0, 10); work = 1; while(work) { if(debug['v']) Bprint(&bso, "%5.2f library pass: %s\n", cputime(), file); Bflush(&bso); work = 0; for(e = start; e < stop; e = strchr(e+5, 0) + 1) { s = lookup(e+5, 0); if(s->type != SXREF) continue; sprint(pname, "%s(%s)", file, s->name); if(debug['v']) Bprint(&bso, "%5.2f library: %s\n", cputime(), pname); Bflush(&bso); l = e[1] & 0xff; l |= (e[2] & 0xff) << 8; l |= (e[3] & 0xff) << 16; l |= (e[4] & 0xff) << 24; seek(f, l, 0); l = read(f, &arhdr, SAR_HDR); if(l != SAR_HDR) goto bad; if(strncmp(arhdr.fmag, ARFMAG, sizeof(arhdr.fmag))) goto bad; l = atolwhex(arhdr.size); ldobj(f, l, pname); if(s->type == SXREF) { diag("%s: failed to load: %s\n", file, s->name); errorexit(); } work = 1; xrefresolv = 1; } } return; bad: diag("%s: bad or out of date archive\n", file); out: close(f); } int zaddr(uchar *p, Adr *a, Sym *h[]) { int c, t, i; long l; Sym *s; Auto *u; t = p[0]; c = 1; if(t & T_INDEX) { a->index = p[c]; a->scale = p[c+1]; c += 2; } else { a->index = D_NONE; a->scale = 0; } a->offset = 0; if(t & T_OFFSET) { a->offset = p[c] | (p[c+1]<<8) | (p[c+2]<<16) | (p[c+3]<<24); c += 4; } a->sym = S; if(t & T_SYM) { a->sym = h[p[c]]; c++; } a->type = D_NONE; if(t & T_FCONST) { a->ieee.l = p[c] | (p[c+1]<<8) | (p[c+2]<<16) | (p[c+3]<<24); a->ieee.h = p[c+4] | (p[c+5]<<8) | (p[c+6]<<16) | (p[c+7]<<24); c += 8; a->type = D_FCONST; } else if(t & T_SCONST) { for(i=0; iscon[i] = p[c+i]; c += NSNAME; a->type = D_SCONST; } if(t & T_TYPE) { a->type = p[c]; c++; } s = a->sym; if(s == S) return c; t = a->type; if(t != D_AUTO && t != D_PARAM) return c; l = a->offset; for(u=curauto; u; u=u->link) { if(u->asym == s) if(u->type == t) { if(u->aoffset > l) u->aoffset = l; return c; } } while(nhunk < sizeof(Auto)) gethunk(); u = (Auto*)hunk; nhunk -= sizeof(Auto); hunk += sizeof(Auto); u->link = curauto; curauto = u; u->asym = s; u->aoffset = l; u->type = t; return c; } void addlib(long line) { char name[MAXHIST*NAMELEN], comp[4*NAMELEN], *p; int i; USED(line); if(histfrogp <= 0) return; if(histfrog[0]->name[1] == '/') { sprint(name, ""); i = 1; } else if(histfrog[0]->name[1] == '.') { sprint(name, "."); i = 0; } else { if(debug['9']) sprint(name, "/%s/lib", thestring); else sprint(name, "/usr/%clib", thechar); i = 0; } for(; iname+1); for(;;) { p = strstr(comp, "$O"); if(p == 0) break; memmove(p+1, p+2, strlen(p+2)+1); p[0] = thechar; } for(;;) { p = strstr(comp, "$M"); if(p == 0) break; memmove(p+strlen(thestring), p+2, strlen(p+2)+1); memmove(p, thestring, strlen(thestring)); if(strlen(comp) > NAMELEN) { diag("library component too long"); return; } } if(strlen(comp) > NAMELEN || strlen(name) + strlen(comp) + 3 >= sizeof(name)) { diag("library component too long"); return; } strcat(name, "/"); strcat(name, comp); } for(i=0; iname = malloc(2*(histfrogp+1) + 1); u->asym = s; u->type = type; u->aoffset = line; u->link = curhist; curhist = u; j = 1; for(i=0; ivalue; s->name[j+0] = k>>8; s->name[j+1] = k; j += 2; } } void histtoauto(void) { Auto *l; while(l = curhist) { curhist = l->link; l->link = curauto; curauto = l; } } void collapsefrog(Sym *s) { int i; /* * bad encoding of path components only allows * MAXHIST components. if there is an overflow, * first try to collapse xxx/.. */ for(i=1; iname+1, "..") == 0) { memmove(histfrog+i-1, histfrog+i+1, (histfrogp-i-1)*sizeof(histfrog[0])); histfrogp--; goto out; } /* * next try to collapse . */ for(i=0; iname+1, ".") == 0) { memmove(histfrog+i, histfrog+i+1, (histfrogp-i-1)*sizeof(histfrog[0])); goto out; } /* * last chance, just truncate from front */ memmove(histfrog+0, histfrog+1, (histfrogp-1)*sizeof(histfrog[0])); out: histfrog[histfrogp-1] = s; } void nopout(Prog *p) { p->as = ANOP; p->from.type = D_NONE; p->to.type = D_NONE; } uchar* readsome(int f, uchar *buf, uchar *good, uchar *stop, int max) { int n; n = stop - good; memmove(buf, good, stop - good); stop = buf + n; n = MAXIO - n; if(n > max) n = max; n = read(f, stop, n); if(n <= 0) return 0; return stop + n; } void ldobj(int f, long c, char *pn) { long ipc; Prog *p, *t; uchar *bloc, *bsize, *stop; int v, o, r, skip; Sym *h[NSYM], *s, *di; bsize = buf.xbuf; bloc = buf.xbuf; di = S; newloop: memset(h, 0, sizeof(h)); version++; histfrogp = 0; ipc = pc; skip = 0; loop: if(c <= 0) goto eof; r = bsize - bloc; if(r < 100 && r < c) { /* enough for largest prog */ bsize = readsome(f, buf.xbuf, bloc, bsize, c); if(bsize == 0) goto eof; bloc = buf.xbuf; goto loop; } o = bloc[0] | (bloc[1] << 8); if(o <= AXXX || o >= ALAST) { if(o < 0) goto eof; diag("%s: opcode out of range %d\n", pn, o); print(" probably not a .8 file\n"); errorexit(); } if(o == ANAME) { stop = memchr(&bloc[4], 0, bsize-&bloc[4]); if(stop == 0){ bsize = readsome(f, buf.xbuf, bloc, bsize, c); if(bsize == 0) goto eof; bloc = buf.xbuf; stop = memchr(&bloc[4], 0, bsize-&bloc[4]); if(stop == 0){ fprint(2, "%s: name too long\n", pn); errorexit(); } } v = bloc[2]; /* type */ o = bloc[3]; /* sym */ bloc += 4; c -= 4; r = 0; if(v == D_STATIC) r = version; s = lookup((char*)bloc, r); c -= &stop[1] - bloc; bloc = stop + 1; if(debug['W']) print(" ANAME %s\n", s->name); h[o] = s; if((v == D_EXTERN || v == D_STATIC) && s->type == 0) s->type = SXREF; if(v == D_FILE) { if(s->type != SFILE) { histgen++; s->type = SFILE; s->value = histgen; } if(histfrogp < MAXHIST) { histfrog[histfrogp] = s; histfrogp++; } else collapsefrog(s); } goto loop; } while(nhunk < sizeof(Prog)) gethunk(); p = (Prog*)hunk; nhunk -= sizeof(Prog); hunk += sizeof(Prog); p->as = o; p->line = bloc[2] | (bloc[3] << 8) | (bloc[4] << 16) | (bloc[5] << 24); p->back = 2; r = zaddr(bloc+6, &p->from, h) + 6; r += zaddr(bloc+r, &p->to, h); bloc += r; c -= r; if(debug['W']) print("%P\n", p); switch(p->as) { case AHISTORY: if(p->to.offset == -1) { addlib(p->line); histfrogp = 0; goto loop; } addhist(p->line, D_FILE); /* 'z' */ if(p->to.offset) addhist(p->to.offset, D_FILE1); /* 'Z' */ histfrogp = 0; goto loop; case AEND: histtoauto(); if(curtext != P) curtext->to.autom = curauto; curauto = 0; curtext = P; if(c) goto newloop; return; case AGLOBL: s = p->from.sym; if(s->type == 0 || s->type == SXREF) { s->type = SBSS; s->value = 0; } if(s->type != SBSS) { diag("%s: redefinition: %s in %s\n", pn, s->name, TNAME); s->type = SBSS; s->value = 0; } if(p->to.offset > s->value) s->value = p->to.offset; goto loop; case ADYNT: if(p->to.sym == S) { diag("DYNT without a sym\n%P\n", p); break; } di = p->to.sym; p->from.scale = 4; if(di->type == SXREF) { if(debug['z']) Bprint(&bso, "%P set to %d\n", p, dtype); di->type = SCONST; di->value = dtype; dtype += 4; } if(p->from.sym == S) break; p->from.offset = di->value; p->from.sym->type = SDATA; if(curtext == P) { diag("DYNT not in text: %P\n", p); break; } p->to.sym = curtext->from.sym; p->to.type = D_ADDR; p->to.index = D_EXTERN; goto data; case AINIT: if(p->from.sym == S) { diag("INIT without a sym\n%P\n", p); break; } if(di == S) { diag("INIT without previous DYNT\n%P\n", p); break; } p->from.offset = di->value; p->from.sym->type = SDATA; goto data; case ADATA: data: if(edatap == P) datap = p; else edatap->link = p; edatap = p; p->link = P; goto loop; case AGOK: diag("%s: GOK opcode in %s\n", pn, TNAME); pc++; goto loop; case ATEXT: if(curtext != P) { histtoauto(); curtext->to.autom = curauto; curauto = 0; } skip = 0; curtext = p; s = p->from.sym; if(s == S) { diag("%s: no TEXT symbol: %P\n", pn, p); errorexit(); } if(s->type != 0 && s->type != SXREF) { if(p->from.scale & DUPOK) { skip = 1; goto casdef; } diag("redefinition: %s\n%P\n", s->name, p); } s->type = STEXT; s->value = pc; lastp->link = p; lastp = p; p->pc = pc; pc++; if(textp == P) { textp = p; etextp = p; goto loop; } etextp->pcond = p; etextp = p; goto loop; case AFMOVF: case AFADDF: case AFSUBF: case AFSUBRF: case AFMULF: case AFDIVF: case AFDIVRF: case AFCOMF: case AFCOMFP: if(skip) goto casdef; if(p->from.type == D_FCONST) { /* size sb 9 max */ sprint(literal, "$%lux", ieeedtof(&p->from.ieee)); s = lookup(literal, 0); if(s->type == 0) { s->type = SBSS; s->value = 4; t = prg(); t->as = ADATA; t->line = p->line; t->from.type = D_EXTERN; t->from.sym = s; t->from.scale = 4; t->to = p->from; if(edatap == P) datap = t; else edatap->link = t; edatap = t; t->link = P; } p->from.type = D_EXTERN; p->from.sym = s; p->from.offset = 0; } goto casdef; case AFMOVD: case AFADDD: case AFSUBD: case AFSUBRD: case AFMULD: case AFDIVD: case AFDIVRD: case AFCOMD: case AFCOMDP: if(skip) goto casdef; if(p->from.type == D_FCONST) { /* size sb 18 max */ sprint(literal, "$%lux.%lux", p->from.ieee.l, p->from.ieee.h); s = lookup(literal, 0); if(s->type == 0) { s->type = SBSS; s->value = 8; t = prg(); t->as = ADATA; t->line = p->line; t->from.type = D_EXTERN; t->from.sym = s; t->from.scale = 8; t->to = p->from; if(edatap == P) datap = t; else edatap->link = t; edatap = t; t->link = P; } p->from.type = D_EXTERN; p->from.sym = s; p->from.offset = 0; } goto casdef; case AADDL: if(p->from.type == D_CONST) if(p->from.offset < 0) { p->as = ASUBL; p->from.offset = -p->from.offset; } goto casdef; case ASUBL: if(p->from.type == D_CONST) if(p->from.offset < 0) { p->as = AADDL; p->from.offset = -p->from.offset; } goto casdef; casdef: default: if(skip) nopout(p); if(p->to.type == D_BRANCH) p->to.offset += ipc; lastp->link = p; lastp = p; p->pc = pc; pc++; goto loop; } goto loop; eof: diag("truncated object file: %s\n", pn); } Sym* lookup(char *symb, int v) { Sym *s; char *p; long h; int l, c; h = v; for(p=symb; c = *p; p++) h = h+h+h + c; l = (p - symb) + 1; if(h < 0) h = ~h; h %= NHASH; for(s = hash[h]; s != S; s = s->link) if(s->version == v) if(memcmp(s->name, symb, l) == 0) return s; while(nhunk < sizeof(Sym)) gethunk(); s = (Sym*)hunk; nhunk -= sizeof(Sym); hunk += sizeof(Sym); s->name = malloc(l + 1); memmove(s->name, symb, l); s->link = hash[h]; s->type = 0; s->version = v; s->value = 0; hash[h] = s; nsymbol++; return s; } Prog* prg(void) { Prog *p; while(nhunk < sizeof(Prog)) gethunk(); p = (Prog*)hunk; nhunk -= sizeof(Prog); hunk += sizeof(Prog); *p = zprg; return p; } Prog* copyp(Prog *q) { Prog *p; p = prg(); *p = *q; return p; } Prog* appendp(Prog *q) { Prog *p; p = prg(); p->link = q->link; q->link = p; p->line = q->line; return p; } void gethunk(void) { char *h; long nh; nh = NHUNK; if(thunk >= 5L*NHUNK) { nh = 5L*NHUNK; if(thunk >= 25L*NHUNK) nh = 25L*NHUNK; } h = mysbrk(nh); if(h == (char*)-1) { diag("out of memory\n"); errorexit(); } hunk = h; nhunk = nh; thunk += nh; } void doprof1(void) { Sym *s; long n; Prog *p, *q; if(debug['v']) Bprint(&bso, "%5.2f profile 1\n", cputime()); Bflush(&bso); s = lookup("__mcount", 0); n = 1; for(p = firstp->link; p != P; p = p->link) { if(p->as == ATEXT) { q = prg(); q->line = p->line; q->link = datap; datap = q; q->as = ADATA; q->from.type = D_EXTERN; q->from.offset = n*4; q->from.sym = s; q->from.scale = 4; q->to = p->from; q->to.type = D_CONST; q = prg(); q->line = p->line; q->pc = p->pc; q->link = p->link; p->link = q; p = q; p->as = AADDL; p->from.type = D_CONST; p->from.offset = 1; p->to.type = D_EXTERN; p->to.sym = s; p->to.offset = n*4 + 4; n += 2; continue; } } q = prg(); q->line = 0; q->link = datap; datap = q; q->as = ADATA; q->from.type = D_EXTERN; q->from.sym = s; q->from.scale = 4; q->to.type = D_CONST; q->to.offset = n; s->type = SBSS; s->value = n*4; } void doprof2(void) { Sym *s2, *s4; Prog *p, *q, *ps2, *ps4; if(debug['v']) Bprint(&bso, "%5.2f profile 2\n", cputime()); Bflush(&bso); s2 = lookup("_profin", 0); s4 = lookup("_profout", 0); ps2 = P; ps4 = P; for(p = firstp; p != P; p = p->link) { if(p->as == ATEXT) { if(p->from.sym == s2) { p->from.scale = 1; ps2 = p; } if(p->from.sym == s4) { p->from.scale = 1; ps4 = p; } } } for(p = firstp; p != P; p = p->link) { if(p->as == ATEXT) { curtext = p; if(p->from.scale & NOPROF) { /* dont profile */ for(;;) { q = p->link; if(q == P) break; if(q->as == ATEXT) break; p = q; } continue; } /* * JMPL profin */ q = prg(); q->line = p->line; q->pc = p->pc; q->link = p->link; p->link = q; p = q; p->as = ACALL; p->to.type = D_BRANCH; p->pcond = ps2; p->to.sym = s2; continue; } if(p->as == ARET) { /* * RET */ q = prg(); q->as = ARET; q->from = p->from; q->to = p->to; q->link = p->link; p->link = q; /* * JAL profout */ p->as = ACALL; p->from = zprg.from; p->to = zprg.to; p->to.type = D_BRANCH; p->pcond = ps4; p->to.sym = s4; p = q; continue; } } } void nuxiinit(void) { int i, c; for(i=0; i<4; i++) { c = find1(0x04030201L, i+1); if(i < 2) inuxi2[i] = c; if(i < 1) inuxi1[i] = c; inuxi4[i] = c; fnuxi4[i] = c; fnuxi8[i] = c; fnuxi8[i+4] = c+4; } if(debug['v']) { Bprint(&bso, "inuxi = "); for(i=0; i<1; i++) Bprint(&bso, "%d", inuxi1[i]); Bprint(&bso, " "); for(i=0; i<2; i++) Bprint(&bso, "%d", inuxi2[i]); Bprint(&bso, " "); for(i=0; i<4; i++) Bprint(&bso, "%d", inuxi4[i]); Bprint(&bso, "\nfnuxi = "); for(i=0; i<4; i++) Bprint(&bso, "%d", fnuxi4[i]); Bprint(&bso, " "); for(i=0; i<8; i++) Bprint(&bso, "%d", fnuxi8[i]); Bprint(&bso, "\n"); } Bflush(&bso); } int find1(long l, int c) { char *p; int i; p = (char*)&l; for(i=0; i<4; i++) if(*p++ == c) return i; return 0; } int find2(long l, int c) { short *p; int i; p = (short*)&l; for(i=0; i<4; i+=2) { if(((*p >> 8) & 0xff) == c) return i; if((*p++ & 0xff) == c) return i+1; } return 0; } long ieeedtof(Ieee *e) { int exp; long v; if(e->h == 0) return 0; exp = (e->h>>20) & ((1L<<11)-1L); exp -= (1L<<10) - 2L; v = (e->h & 0xfffffL) << 3; v |= (e->l >> 29) & 0x7L; if((e->l >> 28) & 1) { v++; if(v & 0x800000L) { v = (v & 0x7fffffL) >> 1; exp++; } } if(exp <= -126 || exp >= 130) diag("double fp to single fp overflow\n"); v |= ((exp + 126) & 0xffL) << 23; v |= e->h & 0x80000000L; return v; } double ieeedtod(Ieee *ieeep) { Ieee e; double fr; int exp; if(ieeep->h & (1L<<31)) { e.h = ieeep->h & ~(1L<<31); e.l = ieeep->l; return -ieeedtod(&e); } if(ieeep->l == 0 && ieeep->h == 0) return 0; fr = ieeep->l & ((1L<<16)-1L); fr /= 1L<<16; fr += (ieeep->l>>16) & ((1L<<16)-1L); fr /= 1L<<16; fr += (ieeep->h & (1L<<20)-1L) | (1L<<20); fr /= 1L<<21; exp = (ieeep->h>>20) & ((1L<<11)-1L); exp -= (1L<<10) - 2L; return ldexp(fr, exp); }