File: [local] / OpenXM_contrib2 / asir2000 / io / pexpr_body.c (download)
Revision 1.14, Fri Dec 20 02:02:24 2013 UTC (10 years, 9 months ago) by noro
Branch: MAIN
CVS Tags: RELEASE_1_3_1_13b Changes since 1.13: +7 -1
lines
Modified nd_f4 and nd_f4_trace so that it partially use GMP.
|
/* $OpenXM: OpenXM_contrib2/asir2000/io/pexpr_body.c,v 1.14 2013/12/20 02:02:24 noro Exp $ */
#define PRINTHAT (fortran_output?PUTS("**"):PUTS("^"))
void PRINTEXPR();
void PRINTNUM();
void PRINTN();
void PRINTV();
void PRINTP();
void PRINTR();
void PRINTLIST();
void PRINTVECT();
void PRINTMAT();
void PRINTIMAT(); /* IMAT */
void PRINTSTR();
void PRINTCOMP();
void PRINTDP();
void PRINTUI();
void PRINTGF2MAT();
void PRINTGFMMAT();
void PRINTBYTEARRAY();
void PRINTQUOTE();
void PRINTQUOTEARG();
void PRINTNBP();
void PRINTERR();
void PRINTCPLX();
void PRINTLM();
void PRINTLF();
void PRINTUP2();
void PRINTUP();
void PRINTUM();
void PRINTFOP();
void PRINTEOP();
void PRINTLOP();
void PRINTQOP();
void PRINTSF();
void PRINTSYMBOL();
void PRINTRANGE();
void PRINTTB();
void PRINTDPV();
void PRINTFNODE();
void PRINTBF();
void PRINTDAlg();
void PRINTEXPR(vl,p)
VL vl;
Obj p;
{
if ( !p ) {
PRINTR(vl,(R)p);
} else
switch ( OID(p) ) {
case O_N:
PRINTNUM((Num)p); break;
case O_P:
PRINTP(vl,(P)p); break;
case O_R:
PRINTR(vl,(R)p); break;
case O_LIST:
PRINTLIST(vl,(LIST)p); break;
case O_VECT:
PRINTVECT(vl,(VECT)p); break;
case O_MAT:
PRINTMAT(vl,(MAT)p); break;
/* IMAT */
case O_IMAT:
PRINTIMAT(vl,(IMAT)p); break;
/* IMAT */
case O_STR:
PRINTSTR((STRING)p); break;
case O_COMP:
PRINTCOMP(vl,(COMP)p); break;
case O_DP:
PRINTDP(vl,(DP)p); break;
case O_USINT:
PRINTUI(vl,(USINT)p); break;
case O_GF2MAT:
PRINTGF2MAT(vl,(GF2MAT)p); break;
case O_ERR:
PRINTERR(vl,(ERR)p); break;
case O_MATHCAP:
PRINTLIST(vl,((MATHCAP)p)->body); break;
case O_F:
PRINTLF(vl,(F)p); break;
case O_GFMMAT:
PRINTGFMMAT(vl,(GFMMAT)p); break;
case O_BYTEARRAY:
PRINTBYTEARRAY(vl,(BYTEARRAY)p); break;
case O_QUOTE:
PRINTQUOTE(vl,(QUOTE)p); break;
case O_QUOTEARG:
PRINTQUOTEARG(vl,(QUOTEARG)p); break;
case O_NBP:
PRINTNBP(vl,(NBP)p); break;
case O_SYMBOL:
PRINTSYMBOL((SYMBOL)p); break;
case O_RANGE:
PRINTRANGE(vl,(RANGE)p); break;
case O_TB:
PRINTTB(vl,(TB)p); break;
case O_DPV:
PRINTDPV(vl,(DPV)p); break;
default:
break;
}
}
#if defined(FPRINT)
#if defined(ITVDEBUG)
void printbin(double z)
{
int i, j, mask;
union {
double x;
char c[8];
} a;
a.x = z;
for(i=7;i>=0;i--) {
mask = 0x80;
for(j=0;j<8;j++) {
if (a.c[i] & mask) fprintf(stderr,"1");
else fprintf(stderr,"0");
mask >>= 1;
}
}
fprintf(stderr,"\n");
}
#endif
#endif
void PRINTCPLX(a)
C a;
{
PUTS("(");
if ( a->r )
PRINTNUM(a->r);
if ( a->i ) {
#if defined(INTERVAL)
if ( a->r && ((compnum(0,a->i,0) > 0)
|| NID(a->i) == N_IP || NID(a->i) == N_IntervalDouble
|| NID(a->i) == N_IntervalQuad || NID(a->i) == N_IntervalBigFloat) )
#else
if ( a->r && (compnum(0,a->i,0) > 0) )
#endif
PUTS("+");
PRINTNUM(a->i); PUTS("*@i");
}
PUTS(")");
}
void PRINTP(vl,p)
VL vl;
P p;
{
V v;
DCP dc;
if ( !p )
PUTS("0");
else if ( NUM(p) )
PRINTNUM((Num)p);
else
for ( dc = DC(p), v = VR(p); dc; dc = NEXT(dc) ) {
if ( !DEG(dc) )
PRINTP(vl,COEF(dc));
else {
if ( NUM(COEF(dc)) && UNIQ((Q)COEF(dc)) ) {
;
} else if ( NUM(COEF(dc)) && MUNIQ((Q)COEF(dc)) )
PUTS("-");
else if ( NUM(COEF(dc)) || !NEXT(DC(COEF(dc)))) {
PRINTP(vl,COEF(dc)); PUTS("*");
} else {
PUTS("("); PRINTP(vl,COEF(dc)); PUTS(")*");
}
PRINTV(vl,v);
if ( cmpq(DEG(dc),ONE) ) {
PRINTHAT;
if ( INT(DEG(dc)) && SGN(DEG(dc))>0 )
PRINTNUM((Num)DEG(dc));
else {
PUTS("("); PRINTNUM((Num)DEG(dc)); PUTS(")");
}
}
}
if ( NEXT(dc) ) {
P t;
t = COEF(NEXT(dc));
if (!DEG(NEXT(dc))) {
if ( NUM(t) ) {
if ( !mmono(t) )
PUTS("+");
} else {
if (!mmono(COEF(DC(t))))
PUTS("+");
}
} else {
if ( !mmono(t) )
PUTS("+");
}
}
}
}
void PRINTR(vl,a)
VL vl;
R a;
{
if ( !a )
PUTS("0");
else
switch (OID(a)) {
case O_N: case O_P:
PRINTP(vl,(P)a); break;
default:
PUTS("("); PRINTP(vl,NM((R)a)); PUTS(")/("); PRINTP(vl,DN((R)a)); PUTS(")");
break;
}
}
void PRINTVECT(vl,vect)
VL vl;
VECT vect;
{
int i;
pointer *ptr;
switch ( outputstyle ) {
case 1:
PUTS("vect(");
for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) {
if ( i != 0 ) PUTS(",");
PRINTEXPR(vl,ptr[i]);
}
PUTS(")");
break;
case 0:
default:
PUTS("[ ");
for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) {
PRINTEXPR(vl,ptr[i]); PUTS(" ");
}
PUTS("]");
break;
}
}
void PRINTDPV(vl,d)
VL vl;
DPV d;
{
int i;
DP *ptr;
PUTS("modvect(");
for ( i = 0, ptr = BDY(d); i < d->len; i++ ) {
if ( i != 0 ) PUTS(",");
PRINTEXPR(vl,(Obj)ptr[i]);
}
PUTS(")");
}
void PRINTMAT(vl,mat)
VL vl;
MAT mat;
{
int i,j,r,c;
pointer *ptr;
switch ( outputstyle ) {
case 1:
PUTS("mat(\n");
for ( i = 0, r = mat->row, c = mat->col; i < r; i++ ) {
if ( i != 0 ) PUTS(",\n");
PUTS("[ ");
for ( j = 0, ptr = BDY(mat)[i]; j < c; j++ ) {
if ( j != 0 ) PUTS(",");
PRINTEXPR(vl,ptr[j]);
}
PUTS(" ]");
}
PUTS(")");
break;
case 0:
default:
for ( i = 0, r = mat->row, c = mat->col; i < r; i++ ) {
PUTS("[ ");
for ( j = 0, ptr = BDY(mat)[i]; j < c; j++ ) {
PRINTEXPR(vl,ptr[j]); PUTS(" ");
}
PUTS("]");
if ( i < r - 1 )
PUTS("\n");
}
break;
}
}
/* IMAT */
#if !defined(CPRINT)
void PRINTIMAT(vl,xmat)
VL vl;
IMAT xmat;
{
int i,j,c;
int row, col, cr;
IMATC Im;
IENT ent;
row = xmat->row;
col = xmat->col;
Im = (pointer)xmat->root;
if ( xmat->clen == 0 ) {
for(j = 0; j< row; j++) {
PUTS("( ");
for(i = 0; i < col; i++) PUTS("0 ");
PUTS(")");
if ( j < row - 1 ) PUTS("\n");
}
return;
}
c = -1;
GetNextIent(&Im, &ent, &c);
for(j = 0; j < row; j++) {
PUTS("( ");
for(i = 0; i < col; i++) {
cr = j * col + i;
if( ent.cr == cr ) {
PRINTEXPR(vl, (pointer)ent.body); PUTS(" ");
GetNextIent(&Im, &ent, &c);
} else {
PUTS("0 ");
}
}
PUTS(")");
if ( j < row - 1 ) PUTS("\n");
}
}
#endif
/* IMAT */
void PRINTLIST(vl,list)
VL vl;
LIST list;
{
NODE tnode;
PUTS("[");
for ( tnode = (NODE)list->body; tnode; tnode = NEXT(tnode) ) {
PRINTEXPR(vl,tnode->body);
if ( NEXT(tnode) )
PUTS(",");
}
PUTS("]");
}
void PRINTCOMP(vl,c)
VL vl;
COMP c;
{
int n,i;
n = getcompsize((int)c->type);
PUTS("{");
for ( i = 0; i < n; i++ ) {
PRINTEXPR(vl,(pointer)c->member[i]);
if ( i < n-1 )
PUTS(",");
}
PUTS("}");
}
void PRINTGF2MAT(vl,mat)
VL vl;
GF2MAT mat;
{
int row,col,w,i,j,k,m;
unsigned int t;
unsigned int **b;
row = mat->row;
col = mat->col;
w = (col+BSH-1)/BSH;
b = mat->body;
for ( i = 0; i < row; i++ ) {
for ( j = 0, m = 0; j < w; j++ ) {
t = b[i][j];
for ( k = 0; m < col && k < BSH; k++, m++ )
if ( t & (1<<k) )
PUTS("1");
else
PUTS("0");
}
PUTS("\n");
}
}
void PRINTERR(vl,e)
VL vl;
ERR e;
{
PUTS("error("); PRINTEXPR(vl,e->body); PUTS(")");
}
void PRINTLF(vl,f)
VL vl;
F f;
{
switch ( FOP(f) ) {
case AL_TRUE:
PUTS("@true");
break;
case AL_FALSE:
PUTS("@false");
break;
case AL_OR: case AL_AND:
PRINTFOP(vl,f); break;
case AL_NOT: case AL_IMPL: case AL_REPL: case AL_EQUIV:
PRINTEOP(vl,f); break;
case AL_EQUAL: case AL_NEQ: case AL_LESSP:
case AL_GREATERP: case AL_LEQ: case AL_GEQ:
PRINTLOP(vl,f); break;
case AL_EX: case AL_ALL:
PRINTQOP(vl,f); break;
default:
break;
}
}
void PRINTFOP(vl,f)
VL vl;
F f;
{
char *op;
NODE n;
op = FOP(f)==AL_OR?" @|| ":" @&& ";
n = FJARG(f);
PUTS("("); PRINTEXPR(vl,BDY(n)); PUTS(")");
for ( n = NEXT(n); n; n = NEXT(n) ) {
PUTS(op); PUTS("("); PRINTEXPR(vl,BDY(n)); PUTS(")");
}
}
void PRINTEOP(vl,f)
VL vl;
F f;
{
oFOP op;
char *sop;
if ( (op = FOP(f)) == AL_NOT ) {
PUTS("(@! "); PRINTEXPR(vl,(Obj)FARG(f)); PUTS(")"); return;
}
switch ( op ) {
case AL_IMPL:
sop = " @impl "; break;
case AL_REPL:
sop = " @repl "; break;
case AL_EQUIV:
sop = " @equiv "; break;
default:
break;
}
PUTS("(");
PRINTEXPR(vl,(Obj)FLHS(f));
PUTS(sop);
PRINTEXPR(vl,(Obj)FRHS(f));
PUTS(")");
}
void PRINTLOP(vl,f)
VL vl;
F f;
{
char *op;
switch ( FOP(f) ) {
case AL_EQUAL:
op = " @== "; break;
case AL_NEQ:
op = " @!= "; break;
case AL_LESSP:
op = " @< "; break;
case AL_GREATERP:
op = " @> "; break;
case AL_LEQ:
op = " @<= "; break;
case AL_GEQ:
op = " @>= "; break;
default:
error("PRINTLOP : invalid operator");
break;
}
PRINTEXPR(vl,(Obj)FPL(f)); PUTS(op); PUTS("0");
}
void PRINTRANGE(VL vl,RANGE p)
{
PUTS("range(");
PRINTEXPR(vl,p->start);
PUTS(",");
PRINTEXPR(vl,p->end);
PUTS(")");
}
void PRINTFNODENODE(NODE n)
{
for ( ; n; n = NEXT(n) ) {
PRINTFNODE((FNODE)BDY(n),0);
if ( NEXT(n) ) PUTS(",");
}
}
void PRINTFARGS(FNODE f)
{
NODE n;
if ( f->id == I_LIST ) {
n = (NODE)FA0(f);
PRINTFNODENODE(n);
} else
PRINTFNODE(f,0);
}
void PRINTFNODE(FNODE f,int paren)
{
NODE n,t,t0;
char vname[BUFSIZ],prefix[BUFSIZ];
char *opname,*vname_conv,*prefix_conv;
Obj obj;
int i,len,allzero,elen,elen2;
C cplx;
char *r;
FNODE fi,f2;
if ( !f ) {
PUTS("(0)");
return;
}
if ( paren ) PUTS("(");
switch ( f->id ) {
/* unary operators */
case I_NOT: PUTS("!"); PRINTFNODE((FNODE)FA0(f),1); break;
case I_PAREN: PRINTFNODE((FNODE)FA0(f),0); break;
case I_MINUS: PUTS("-"); PRINTFNODE((FNODE)FA0(f),1); break;
/* binary operators */
/* arg list */
/* I_AND, I_OR => FA0(f), FA1(f) */
/* otherwise => FA1(f), FA2(f) */
case I_BOP:
PRINTFNODE((FNODE)FA1(f),1);
PUTS(((ARF)FA0(f))->name);
PRINTFNODE((FNODE)FA2(f),1);
break;
case I_NARYOP:
t = (NODE)FA1(f);
PRINTFNODE((FNODE)BDY(t),1);
for ( t = NEXT(t); t; t = NEXT(t) ) {
PUTS(((ARF)FA0(f))->name);
PRINTFNODE((FNODE)BDY(t),1);
}
break;
case I_COP:
switch( (cid)FA0(f) ) {
case C_EQ: opname = ("=="); break;
case C_NE: opname = ("!="); break;
case C_GT: opname = (">"); break;
case C_LT: opname = ("<"); break;
case C_GE: opname = (">="); break;
case C_LE: opname = ("<="); break;
}
PRINTFNODE((FNODE)FA1(f),1);
PUTS(opname);
PRINTFNODE((FNODE)FA2(f),1);
break;
case I_LOP:
switch( (lid)FA0(f) ) {
case L_EQ: opname = ("@=="); break;
case L_NE: opname = ("@!="); break;
case L_GT: opname = ("@>"); break;
case L_LT: opname = ("@<"); break;
case L_GE: opname = ("@>="); break;
case L_LE: opname = ("@<="); break;
case L_AND: opname = ("@&&"); break;
case L_OR: opname = ("@||"); break;
case L_NOT: opname = ("@!"); break;
}
if ( (lid)FA0(f)==L_NOT ) {
PUTS(opname); PRINTFNODE((FNODE)FA1(f),1);
} else {
PRINTFNODE((FNODE)FA1(f),1);
PUTS(opname);
PRINTFNODE((FNODE)FA2(f),1);
}
break;
case I_AND:
PRINTFNODE((FNODE)FA0(f),1);
PUTS("&&");
PRINTFNODE((FNODE)FA1(f),1);
break;
case I_OR:
PRINTFNODE((FNODE)FA0(f),1);
PUTS("!!");
PRINTFNODE((FNODE)FA1(f),1);
break;
/* ternary operators */
case I_CE:
PRINTFNODE((FNODE)FA0(f),1); PUTS("?"); PRINTFNODE((FNODE)FA1(f),1);
PUTS(":"); PRINTFNODE((FNODE)FA2(f),1);
break;
/* lists */
case I_LIST: PUTS("["); PRINTFNODENODE((NODE)FA0(f)); PUTS("]"); break;
/* function */
case I_FUNC: case I_FUNC_QARG:
if ( !strcmp(((FUNC)FA0(f))->name,"@pi") ) PUTS("@pi");
else if ( !strcmp(((FUNC)FA0(f))->name,"@e") ) PUTS("@e");
else {
PUTS(((FUNC)FA0(f))->name);
PUTS("("); PRINTFARGS(FA1(f)); PUTS(")");
}
break;
/* XXX */
case I_CAR: PUTS("car("); PRINTFNODE(FA0(f),0); PUTS(")"); break;
case I_CDR: PUTS("cdr("); PRINTFNODE(FA0(f),0); PUTS(")"); break;
/* exponent vector */
case I_EV: PUTS("<<"); PRINTFNODENODE((NODE)FA0(f)); PUTS(">>"); break;
/* string */
case I_STR: PUTS((char *)FA0(f)); break;
/* internal object */
case I_FORMULA: obj = (Obj)FA0(f); PRINTEXPR(CO,obj); break;
/* program variable */
case I_PVAR:
if ( FA1(f) )
error("printfnode : not implemented yet");
GETPVNAME(FA0(f),opname);
PUTS(opname);
break;
default: error("printfnode : not implemented yet");
}
if ( paren ) PUTS(")");
}
/* functions not used in cpexpr.c */
#if !defined(CPRINT)
void PRINTNUM(q)
Num q;
{
char real_format[20];
DAlg d;
DP nm;
Q dn;
if ( !q ) {
PUTS("0");
return;
}
switch ( NID(q) ) {
case N_Q:
if ( SGN((Q)q) == -1 )
PUTS("-");
PRINTN(NM((Q)q));
if ( !INT((Q)q) ) {
PUTS("/"); PRINTN(DN((Q)q));
}
break;
case N_R:
switch (printmode) {
case PRINTF_E:
#if defined(INTERVAL)
case MID_PRINTF_E:
#endif
TAIL PRINTF(OUT,"%.16e",BDY((Real)q));
break;
case PRINTF_G:
#if defined(INTERVAL)
case MID_PRINTF_G:
#endif
default:
if ( real_binary ) {
unsigned int *m;
unsigned int u,l,mask;
int i,expo;
m = (unsigned int *)&BDY((Real)q);
#if defined(__i386__) || defined(MIPSEL) || defined(VISUAL) || defined(__alpha) || defined(__FreeBSD__) || defined(__NetBSD__) || defined(__x86_64)
u = m[1]; l = m[0];
#else
u = m[0]; l = m[1];
#endif
if (u&0x80000000) {
TAIL PRINTF(OUT,"-");
}
u &= 0x7fffffff;
expo = ((int)(u>>20));
u &= 0xfffff;
if ( expo == 2047 ) {
if ( u || l ) {
TAIL PRINTF(OUT,"NaN");
} else {
TAIL PRINTF(OUT,"Inf");
}
} else if ( expo == 0 ) {
if ( u || l ) {
TAIL PRINTF(OUT,"0b0.");
for ( i = 0, mask = 0x80000; i < 20;
i++, mask >>= 1) {
TAIL
if ( u&mask )
PRINTF(OUT,"1");
else
PRINTF(OUT,"0");
}
for ( i = 0, mask = 0x80000000; i < 32;
i++, mask >>= 1) {
TAIL
if ( l&mask )
PRINTF(OUT,"1");
else
PRINTF(OUT,"0");
}
TAIL PRINTF(OUT,"*2^%d",-1022);
} else {
TAIL PRINTF(OUT,"0");
}
} else {
expo -= 1023;
TAIL PRINTF(OUT,"0b1.");
for ( i = 0, mask = 0x80000; i < 20;
i++, mask >>= 1) {
TAIL
if ( u&mask )
PRINTF(OUT,"1");
else
PRINTF(OUT,"0");
}
for ( i = 0, mask = 0x80000000; i < 32;
i++, mask >>= 1) {
TAIL
if ( l&mask )
PRINTF(OUT,"1");
else
PRINTF(OUT,"0");
}
if ( expo ) {
TAIL PRINTF(OUT,"*2^%d",expo);
}
}
} else if ( real_digit ) {
sprintf(real_format,
double_output?"%%.%df":"%%.%dg",real_digit);
TAIL PRINTF(OUT,real_format,BDY((Real)q));
} else {
TAIL PRINTF(OUT,double_output?"%f":"%g",BDY((Real)q));
}
break;
}
break;
case N_A:
PUTS("("); PRINTR(ALG,(R)BDY((Alg)q)); PUTS(")");
break;
#if defined(PARI)
case N_B:
PRINTBF((BF)q); break;
#endif
#if defined(INTERVAL)
case N_IP:
case N_IntervalBigFloat:
switch ( outputstyle ) {
case 1:
PUTS("intval(");
PRINTNUM(INF((Itv)q));
PUTS(",");
PRINTNUM(SUP((Itv)q));
PUTS(")");
break;
case 0:
default:
PUTS("[");
PRINTNUM(INF((Itv)q));
PUTS(",");
PRINTNUM(SUP((Itv)q));
PUTS("]");
break;
}
break;
case N_IntervalDouble:
switch (printmode) {
case PRINTF_E:
switch ( outputstyle ) {
case 1:
TAIL PRINTF(OUT, "intval(%.16e,%.16e)",
INF((IntervalDouble)q),SUP((IntervalDouble)q));
break;
case 0:
default:
TAIL PRINTF(OUT, "[%.16e,%.16e]",
INF((IntervalDouble)q),SUP((IntervalDouble)q));
break;
}
#if defined(ITVDEBUG)
printbin(INF((IntervalDouble)q));
printbin(SUP((IntervalDouble)q));
#endif
break;
case MID_PRINTF_G:
switch ( outputstyle ) {
case 1:
TAIL PRINTF(OUT, "intvalm(%g,%g)",
(SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,
(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5);
break;
case 0:
default:
TAIL PRINTF(OUT, "<%g,%g>",
(SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,
(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5);
break;
}
break;
case MID_PRINTF_E:
switch ( outputstyle ) {
case 1:
TAIL PRINTF(OUT, "intvalm(%.16e,%.16e)",
(SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,
(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5);
break;
case 0:
default:
TAIL PRINTF(OUT, "<%.16e,%.16e>",
(SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,
(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5);
break;
}
break;
case PRINTF_G:
default:
switch ( outputstyle ) {
case 1:
TAIL PRINTF(OUT, "intval(%g,%g)",
INF((IntervalDouble)q),SUP((IntervalDouble)q));
break;
case 0:
default:
TAIL PRINTF(OUT, "[%g,%g]",
INF((IntervalDouble)q),SUP((IntervalDouble)q));
break;
}
break;
}
break;
#endif
case N_C:
PRINTCPLX((C)q); break;
case N_M:
TAIL PRINTF(OUT,"%d",CONT((MQ)q)); break;
case N_LM:
PRINTN(((LM)q)->body); break;
case N_GF2N:
if ( hex_output )
PRINTN((N)(((GF2N)q)->body));
else
PRINTUP2(((GF2N)q)->body);
break;
case N_GFPN:
PRINTUP((UP)(((GFPN)q)->body));
break;
case N_GFS:
TAIL PRINTF(OUT,"@_%d",CONT((GFS)q));
break;
case N_GFSN:
PRINTUM(BDY((GFSN)q));
break;
case N_DA:
d = (DAlg)q;
nm = d->nm;
dn = d->dn;
if ( SGN((Q)dn) == -1 ) PUTS("-");
PUTS("(");
PRINTDP(CO,((DAlg)q)->nm);
PUTS(")");
if ( !UNIN(NM(dn)) ) {
PUTS("/");
PRINTN(NM(dn));
}
break;
case N_GZ:
mpz_out_str(OUT,10,BDY((GZ)q));
break;
case N_GQ:
mpq_out_str(OUT,10,BDY((GQ)q));
break;
}
}
void PRINTN(n)
N n;
{
register int i,*ptr;
N tn;
if ( !n ) {
PUTS("0");
return;
}
if ( hex_output ) {
ptr = BD(n);
TAIL PRINTF(OUT,"0x%x",ptr[PL(n)-1]);
if ( hex_output < 0 )
for ( i = PL(n) - 2; i >= 0; i-- ) {
TAIL PRINTF(OUT,"|%08x",ptr[i]);
}
else
for ( i = PL(n) - 2; i >= 0; i-- ) {
TAIL PRINTF(OUT,"%08x",ptr[i]);
}
} else {
ntobn(DBASE,n,&tn);
ptr = BD(tn);
TAIL PRINTF(OUT,"%d",ptr[PL(tn) - 1]);
for ( i = PL(tn) - 2; i >= 0; i-- ) {
TAIL PRINTF(OUT,DFORMAT,ptr[i]);
}
}
}
void PRINTV(vl,v)
VL vl;
V v;
{
PF pf;
PFAD ad;
int i;
if ( NAME(v) )
PUTS(NAME(v));
else if ( (vid)v->attr == V_PF ) {
pf = ((PFINS)v->priv)->pf; ad = ((PFINS)v->priv)->ad;
if ( !strcmp(NAME(pf),"pow") ) {
PUTS("(("); PRINTR(vl,(R)ad[0].arg); PUTS(")"); PRINTHAT; PUTS("(");
PRINTR(vl,(R)ad[1].arg); PUTS("))");
} else if ( !pf->argc ) {
TAIL PRINTF(OUT,"%s",NAME(pf));
} else {
if ( hideargs ) {
for ( i = 0; i < pf->argc; i++ )
if ( ad[i].d )
break;
if ( i < pf->argc ) {
TAIL PRINTF(OUT,"%s{%d",NAME(pf),ad[0].d);
for ( i = 1; i < pf->argc; i++ ) {
TAIL PRINTF(OUT,",%d",ad[i].d);
}
PUTS("}");
} else {
TAIL PRINTF(OUT,"%s",NAME(pf));
}
} else {
for ( i = 0; i < pf->argc; i++ )
if ( ad[i].d )
break;
if ( i < pf->argc ) {
TAIL PRINTF(OUT,"%s{%d",NAME(pf),ad[0].d);
for ( i = 1; i < pf->argc; i++ ) {
TAIL PRINTF(OUT,",%d",ad[i].d);
}
PUTS("}(");
} else {
TAIL PRINTF(OUT,"%s(",NAME(pf));
}
PRINTR(vl,(R)ad[0].arg);
for ( i = 1; i < pf->argc; i++ ) {
PUTS(","); PRINTR(vl,(R)ad[i].arg);
}
PUTS(")");
}
}
}
}
void PRINTSTR(str)
STRING str;
{
char *p;
for ( p = BDY(str); *p; p++ )
if ( *p == '"' )
PUTS("\"");
else {
TAIL PRINTF(OUT,"%c",*p);
}
}
void PRINTDP(vl,d)
VL vl;
DP d;
{
int n,i;
MP m;
DL dl;
for ( n = d->nv, m = BDY(d); m; m = NEXT(m) ) {
PUTS("("); PRINTEXPR(vl,(pointer)m->c); PUTS(")*<<");
for ( i = 0, dl = m->dl; i < n-1; i++ ) {
TAIL PRINTF(OUT,"%d,",dl->d[i]);
}
TAIL PRINTF(OUT,"%d",dl->d[i]);
PUTS(">>");
if ( NEXT(m) )
PUTS("+");
}
}
void PRINTUI(vl,u)
VL vl;
USINT u;
{
TAIL PRINTF(OUT,"%u",BDY(u));
}
void PRINTGFMMAT(vl,mat)
VL vl;
GFMMAT mat;
{
int row,col,i,j;
unsigned int **b;
row = mat->row;
col = mat->col;
b = mat->body;
for ( i = 0; i < row; i++ ) {
PUTS("[");
for ( j = 0; j < col; j++ ) {
TAIL PRINTF(OUT,"%8d",b[i][j]);
}
PUTS("]\n");
}
}
void PRINTBYTEARRAY(vl,array)
VL vl;
BYTEARRAY array;
{
int len,i;
unsigned char *b;
len = array->len;
b = array->body;
PUTS("|");
for ( i = 0; i < len-1; i++ ) {
TAIL PRINTF(OUT,"%02x ",(unsigned int)b[i]);
}
TAIL PRINTF(OUT,"%02x",(unsigned int)b[i]);
PUTS("|");
}
void PRINTQUOTE(vl,quote)
VL vl;
QUOTE quote;
{
LIST list;
if ( print_quote == 1 ) {
fnodetotree(BDY(quote),&list);
PRINTEXPR(vl,(Obj)list);
} else if ( print_quote == 2 ) {
PRINTFNODE(BDY(quote),0);
} else {
PUTS("<...quoted...>");
}
}
void PRINTQUOTEARG(vl,quote)
VL vl;
QUOTEARG quote;
{
PUTS("<...quoted...>");
}
void PRINTNBP(VL vl,NBP p)
{
NODE t;
NBM m;
int d,i;
unsigned int *b;
if ( !p ) PUTS("0");
else {
for ( t = BDY(p); t; t = NEXT(t) ) {
m = (NBM)BDY(t);
PUTS("(");
PRINTEXPR(vl,(Obj)m->c);
PUTS(")*");
d = m->d;
b = m->b;
if ( d )
for ( i = 0; i < d; i++ ) {
if ( NBM_GET(b,i) ) PUTS("x");
else PUTS("y");
}
else PUTS("1");
if ( NEXT(t) ) PUTS("+");
}
}
}
void PRINTTB(VL vl,TB p)
{
int i;
#if defined(SPRINT)
char *ptr,*e;
for ( i = 0, ptr = buf+strlen(OUT); i < p->next; i++ )
for ( e = p->body[i]; *e; *ptr++ = *e++ );
*ptr = 0;
#else
for ( i = 0; i < p->next; i++ ) {
PUTS(p->body[i]);
}
#endif
}
void PRINTUP2(p)
UP2 p;
{
int d,i;
if ( !p ) {
TAIL PRINTF(OUT,"0");
} else {
d = degup2(p);
TAIL PRINTF(OUT,"(");
if ( !d ) {
TAIL PRINTF(OUT,"1");
} else if ( d == 1 ) {
TAIL PRINTF(OUT,"@");
} else {
PUTS("@"); PRINTHAT; TAIL PRINTF(OUT,"%d",d);
}
for ( i = d-1; i >= 0; i-- ) {
if ( p->b[i/BSH] & (1<<(i%BSH)) )
if ( !i ) {
TAIL PRINTF(OUT,"+1");
} else if ( i == 1 ) {
TAIL PRINTF(OUT,"+@");
} else {
PUTS("+@"); PRINTHAT; TAIL PRINTF(OUT,"%d",i);
}
}
TAIL PRINTF(OUT,")");
}
}
void PRINTQOP(vl,f)
VL vl;
F f;
{
char *op;
op = FOP(f)==AL_EX?"ex":"all";
TAIL PRINTF(OUT,"%s(%s,",op,NAME(FQVR(f)));
PRINTEXPR(vl,(Obj)FQMAT(f)); PUTS(")");
}
void PRINTUP(n)
UP n;
{
int i,d;
if ( !n )
PUTS("0");
else if ( !n->d )
PRINTNUM(n->c[0]);
else {
d = n->d;
PUTS("(");
if ( !d ) {
PRINTNUM(n->c[d]);
} else if ( d == 1 ) {
PRINTNUM(n->c[d]);
PUTS("*@p");
} else {
PRINTNUM(n->c[d]);
PUTS("*@p"); PRINTHAT; TAIL PRINTF(OUT,"%d",d);
}
for ( i = d-1; i >= 0; i-- ) {
if ( n->c[i] ) {
PUTS("+("); PRINTNUM(n->c[i]); PUTS(")");
if ( i >= 2 ) {
PUTS("*@p"); PRINTHAT; TAIL PRINTF(OUT,"%d",i);
} else if ( i == 1 )
PUTS("*@p");
}
}
PUTS(")");
}
}
void PRINTUM(n)
UM n;
{
int i,d;
if ( !n )
PUTS("0");
else if ( !n->d )
PRINTSF(n->c[0]);
else {
d = n->d;
PUTS("(");
if ( !d ) {
PRINTSF(n->c[d]);
} else if ( d == 1 ) {
PRINTSF(n->c[d]);
PUTS("*@s");
} else {
PRINTSF(n->c[d]);
PUTS("*@s"); PRINTHAT; TAIL PRINTF(OUT,"%d",d);
}
for ( i = d-1; i >= 0; i-- ) {
if ( n->c[i] ) {
PUTS("+("); PRINTSF(n->c[i]); PUTS(")");
if ( i >= 2 ) {
PUTS("*@s"); PRINTHAT; TAIL PRINTF(OUT,"%d",i);
} else if ( i == 1 )
PUTS("*@s");
}
}
PUTS(")");
}
}
void PRINTSF(i)
unsigned int i;
{
if ( !i ) {
PUTS("0");
} else {
TAIL PRINTF(OUT,"@_%d",IFTOF(i));
}
}
void PRINTSYMBOL(SYMBOL sym)
{
PUTS(sym->name);
}
#endif /* CPRINT */