/* $OpenXM: OpenXM_contrib2/asir2018/io/pexpr_body.c,v 1.5 2021/03/11 03:41:13 noro Exp $ */

#define PRINTHAT (fortran_output?PUTS("**"):PUTS("^"))

void PRINTEXPR(VL vl,Obj p);
void printbin(double z);
void PRINTCPLX(C a);
void PRINTP(VL vl,P p);
void PRINTR(VL vl,R a);
void PRINTVECT(VL vl,VECT vect);
void PRINTDPV(VL vl,DPV d);
void PRINTMAT(VL vl,MAT mat);
void PRINTIMAT(VL vl,IMAT xmat);
void PRINTLIST(VL vl,LIST list);
void PRINTCOMP(VL vl,COMP c);
void PRINTGF2MAT(VL vl,GF2MAT mat);
void PRINTERR(VL vl,ERR e);
void PRINTLF(VL vl,F f);
void PRINTFOP(VL vl,F f);
void PRINTEOP(VL vl,F f);
void PRINTLOP(VL vl,F f);
void PRINTRANGE(VL vl,RANGE p);
void PRINTFNODENODE(NODE n);
void PRINTFARGS(FNODE f);
void PRINTFNODE(FNODE f,int paren);
void PRINTNUM(Num q);
void PRINTZ(Z n);
void PRINTMPZ(mpz_t z);
void PRINTV(VL vl,V v);
void PRINTSTR(STRING str);
void PRINTDP(VL vl,DP d);
void PRINTDPM(VL vl,DPM d);
void PRINTUI(VL vl,USINT u);
void PRINTGFMMAT(VL vl,GFMMAT mat);
void PRINTBYTEARRAY(VL vl,BYTEARRAY array);
void PRINTQUOTE(VL vl,QUOTE quote);
void PRINTQUOTEARG(VL vl,QUOTEARG quote);
void PRINTNBP(VL vl,NBP p);
void PRINTTB(VL vl,TB p);
void PRINTUP2(UP2 p);
void PRINTQOP(VL vl,F f);
void PRINTUP(UP n);
void PRINTUM(UM n);
void PRINTSF(unsigned int i);
void PRINTSYMBOL(SYMBOL sym);
void PRINTBF(BF a);
void PRINTMPZ(mpz_t z);
#if defined(INTERVAL)
void PRINTBF4ITV(BF);
#endif

void PRINTEXPR(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_DPM:
        PRINTDPM(vl,(DPM)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(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 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 ( cmpz(DEG(dc),ONE) ) {
          PRINTHAT;
          if ( INT(DEG(dc)) && sgnz(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 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 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 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 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 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 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 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 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 vl,ERR e)
{
  PUTS("error("); PRINTEXPR(vl,e->body); PUTS(")");
}

void PRINTLF(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 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 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 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;
        default: 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;
    case I_PFDERIV:
      PUTS(((FUNC)FA0(f))->name);
      PUTS("{"); PRINTFARGS(FA2(f)); PUTS("}");
      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(Num q)
{
  char real_format[20];
  DAlg d;
  DP nm;
  Z num,den;

  if ( !q ) {
    PUTS("0");
    return;
  }
  switch ( NID(q) ) {
    case N_Q:
      nmq((Q)q,&num); PRINTZ(num);
      if ( !INT((Q)q) ) {
        PUTS("/"); dnq((Q)q,&den); PRINTZ(den);
      }
      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(__MINGW32__) || defined(__alpha) || defined(__FreeBSD__) || defined(__NetBSD__) || defined(__x86_64) || defined(__ARM_ARCH) || defined(ANDROID)
            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,"%%.%d%c",
              real_digit,(double_output==1)?'f':(double_output==2)?'e':'g');
            TAIL PRINTF(OUT,real_format,BDY((Real)q));
          } else {
            TAIL PRINTF(OUT,(double_output==1)?"%f":(double_output==2)?"%e":"%g",BDY((Real)q));
          }
          break;
      }
      break;
    case N_A:
      PUTS("("); PRINTR(ALG,(R)BDY((Alg)q)); PUTS(")");
      break;
    case N_B:
      PRINTBF((BF)q); break;
#if defined(INTERVAL)
    case N_IP:
      PUTS("[");
      PRINTNUM(INF((Itv)q));
      PUTS(",");
      PRINTNUM(SUP((Itv)q));
      PUTS("]");
      break;
    case N_IntervalBigFloat:
      switch ( outputstyle ) {
        case 1:
          PUTS("intval(");
          PRINTBF4ITV(INF((IntervalBigFloat)q));
          PUTS(",");
          PRINTBF4ITV(SUP((IntervalBigFloat)q));
          PUTS(")");
          break;
        case 0:
        default:
          PUTS("[");
          PRINTBF4ITV(INF((IntervalBigFloat)q));
          PUTS(",");
          PRINTBF4ITV(SUP((IntervalBigFloat)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:
      PRINTMPZ(BDY((LM)q)); break;
    case N_GF2N:
      PRINTUP2(BDY((GF2N)q));
      break;
    case N_GFPN:
      PRINTUP((UP)BDY((GFPN)q));
      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;
      den = d->dn;
      if ( sgnq((Q)den) < 0 ) PUTS("-");
      PUTS("(");
      PRINTDP(CO,((DAlg)q)->nm);
      PUTS(")");
      if ( !UNIQ(den) ) {
        PUTS("/");
        PRINTZ(den);
      }
      break;
    default:
      break;
  }
}

void PRINTV(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(STRING str)
{
  char *p;

  for ( p = BDY(str); *p; p++ )
    if ( *p == '"' )
      PUTS("\"");
    else {
      TAIL PRINTF(OUT,"%c",*p);
    }
}

void PRINTDP(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 PRINTDPM(VL vl,DPM d)
{
  int n,i;
  DMM 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:%d",dl->d[i],m->pos);
    PUTS(">>");
    if ( NEXT(m) )
      PUTS("+");
  }
}

void PRINTUI(VL vl,USINT u)
{
  TAIL PRINTF(OUT,"%u",BDY(u));
}

void PRINTGFMMAT(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 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 vl,QUOTE quote)
{
  LIST list;

  if ( print_quote == 1 ) {
    fnodetotree(BDY(quote),quote->pvs,&list);
    PRINTEXPR(vl,(Obj)list);
  } else if ( print_quote == 2 ) {
    PRINTFNODE(BDY(quote),0);
  } else {
    PUTS("<...quoted...>");
  }
}

void PRINTQUOTEARG(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(UP2 p)
{
  int d,i;

  if ( !p ) {
    TAIL PRINTF(OUT,"0");
  } else if ( hex_output ) {
    d = p->w;
    TAIL PRINTF(OUT,"0x%x",p->b[d-1]);
    if ( hex_output < 0 )
      for ( i = d - 2; i >= 0; i-- ) {
        TAIL PRINTF(OUT,"|%08x",p->b[i]);
      }
    else
      for ( i = d - 2; i >= 0; i-- ) {
      TAIL PRINTF(OUT,"%08x",p->b[i]);
    }
  } 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 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(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(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(unsigned int i)
{
  if ( !i ) {
    PUTS("0");
  } else {
    TAIL PRINTF(OUT,"@_%d",IFTOF(i));
  }
}

void PRINTSYMBOL(SYMBOL sym)
{
  PUTS(sym->name);
}

#endif /* CPRINT */