[BACK]Return to spexpr.c CVS log [TXT][DIR] Up to [local] / OpenXM_contrib2 / asir2000 / io

Diff for /OpenXM_contrib2/asir2000/io/spexpr.c between version 1.23 and 1.43

version 1.23, 2003/06/07 16:40:25 version 1.43, 2019/11/12 10:52:05
Line 44 
Line 44 
  * OF THE SOFTWARE HAS BEEN DEVELOPED BY A THIRD PARTY, THE THIRD PARTY   * OF THE SOFTWARE HAS BEEN DEVELOPED BY A THIRD PARTY, THE THIRD PARTY
  * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,   * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
  * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.   * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
  * $OpenXM: OpenXM_contrib2/asir2000/io/pexpr.c,v 1.22 2003/05/29 16:45:01 saito Exp $   * $OpenXM: OpenXM_contrib2/asir2000/io/spexpr.c,v 1.42 2018/03/29 07:50:06 ohara Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "al.h"  #include "al.h"
Line 56 
Line 56 
 #define SPRINT  #define SPRINT
 #endif  #endif
   
 #define PRINTHAT (fortran_output?PUTS("**"):PUTS("^"))  
 extern int outputstyle;  extern int outputstyle;
   
 #ifdef FPRINT  
 FILE *asir_out;  
 #define OUT asir_out  
 char DFORMAT[BUFSIZ];  
 int hex_output;  
 int fortran_output;  
 int double_output;  
 int real_digit;  
 int print_quote;  
   
 #define TAIL  
 #define PUTS(s) fputs(s,OUT)  
 #define PRINTF fprintf  
 #define PRINTN printn  
 #define PRINTBF printbf  
 #define PRINTCPLX printcplx  
 #define PRINTLM printlm  
 #define PRINTUP2 printup2  
 #define PRINTV printv  
 #define PRINTEXPR printexpr  
 #define PRINTNUM printnum  
 #define PRINTP asir_printp  
 #define PRINTR printr  
 #define PRINTLIST printlist  
 #define PRINTVECT printvect  
 #define PRINTMAT printmat  
 #define PRINTSTR printstr  
 #define PRINTCOMP printcomp  
 #define PRINTDP printdp  
 #define PRINTUI printui  
 #define PRINTGF2MAT printgf2mat  
 #define PRINTGFMMAT printgfmmat  
 #define PRINTBYTEARRAY printbytearray  
 #define PRINTQUOTE printquote  
 #define PRINTERR printerr  
 #define PRINTLF printlf  
 #define PRINTLOP printlop  
 #define PRINTFOP printfop  
 #define PRINTEOP printeop  
 #define PRINTQOP printqop  
 #define PRINTUP printup  
 #define PRINTUM printum  
 #define PRINTSF printsf  
 #endif  
   
 #ifdef SPRINT  
 static char *buf;  static char *buf;
 #define OUT buf  #define OUT buf
 extern char DFORMAT[BUFSIZ];  extern char DFORMAT[BUFSIZ];
Line 112  extern int hex_output;
Line 65  extern int hex_output;
 extern int fortran_output;  extern int fortran_output;
 extern int double_output;  extern int double_output;
 extern int real_digit;  extern int real_digit;
   extern int real_binary;
 extern int print_quote;  extern int print_quote;
   extern int printmode;
   extern int hideargs;
   
   
 #define TAIL while ( *OUT ) OUT++;  #define TAIL while ( *OUT ) OUT++;
 #define PUTS(s) strcat(OUT,s)  #define PUTS(s) strcat(OUT,s)
 #define PRINTF sprintf  #define PRINTF sprintf
Line 131  extern int print_quote;
Line 86  extern int print_quote;
 #define PRINTLIST sprintlist  #define PRINTLIST sprintlist
 #define PRINTVECT sprintvect  #define PRINTVECT sprintvect
 #define PRINTMAT sprintmat  #define PRINTMAT sprintmat
   #define PRINTIMAT sprintimat
 #define PRINTSTR sprintstr  #define PRINTSTR sprintstr
 #define PRINTCOMP sprintcomp  #define PRINTCOMP sprintcomp
 #define PRINTDP sprintdp  #define PRINTDP sprintdp
   #define PRINTDPM sprintdpm
 #define PRINTUI sprintui  #define PRINTUI sprintui
 #define PRINTGF2MAT sprintgf2mat  #define PRINTGF2MAT sprintgf2mat
 #define PRINTGFMMAT sprintgfmmat  #define PRINTGFMMAT sprintgfmmat
 #define PRINTBYTEARRAY sprintbytearray  #define PRINTBYTEARRAY sprintbytearray
 #define PRINTQUOTE sprintquote  #define PRINTQUOTE sprintquote
   #define PRINTQUOTEARG sprintquotearg
   #define PRINTNBP sprintnbp
 #define PRINTERR sprinterr  #define PRINTERR sprinterr
 #define PRINTLF sprintlf  #define PRINTLF sprintlf
 #define PRINTLOP sprintlop  #define PRINTLOP sprintlop
Line 148  extern int print_quote;
Line 107  extern int print_quote;
 #define PRINTUP sprintup  #define PRINTUP sprintup
 #define PRINTUM sprintum  #define PRINTUM sprintum
 #define PRINTSF sprintsf  #define PRINTSF sprintsf
 #endif  #define PRINTSYMBOL sprintsymbol
   #define PRINTRANGE sprintrange
   #define PRINTTB sprinttb
   #define PRINTDPV sprintdpv
   #define PRINTFNODE sprintfnode
   #define PRINTFNODENODE sprintfnodenode
   #define PRINTFARGS sprintfargs
   
 void PRINTEXPR();  
 void PRINTNUM();  
 void PRINTN();  
 void PRINTV();  
 void PRINTP();  
 void PRINTR();  
 void PRINTLIST();  
 void PRINTVECT();  
 void PRINTMAT();  
 void PRINTSTR();  
 void PRINTCOMP();  
 void PRINTDP();  
 void PRINTUI();  
 void PRINTGF2MAT();  
 void PRINTGFMMAT();  
 void PRINTBYTEARRAY();  
 void PRINTQUOTE();  
 void PRINTERR();  
 void PRINTCPLX();  
 void PRINTLM();  
 void PRINTLF();  
 void PRINTUP2();  
 void PRINTUP();  
 void PRINTUM();  
 void PRINTFOP();  
 void PRINTEOP();  
 void PRINTLOP();  
 void PRINTQOP();  
 void PRINTSF();  
   
 #ifdef FPRINT  
 void output_init() {  
         OUT = stdout;  
         sprintf(DFORMAT,"%%0%dd",DLENGTH);  
 }  
   
 int mmono(p)  
 P p;  
 {  
         if ( NUM(p) )  
 #if defined(INTERVAL)  #if defined(INTERVAL)
                 if ( NID(p) != N_IP && NID(p) != N_IntervalDouble && NID(p) != N_IntervalQuad && NID(p) != N_IntervalBigFloat  #define PRINTBF4ITV sprintbf4itv
                         && compnum(CO,(Num)p,0) < 0 )  
 #else  
                 if ( compnum(CO,(Num)p,0) < 0 )  
 #endif  #endif
                         return ( 1 );  
                 else  
                         return ( 0 );  
         else if ( NEXT(DC(p)) )  
                 return ( 0 );  
         else  
                 return (mmono(COEF(DC(p))));  
 }  
   
 #if defined(PARI)  #include "pexpr_body.c"
 void printbf(a)  
 BF a;  
 {  
         void sor();  
   
         sor(a->body,double_output ? 'f' : 'g',-1,0);  /* special functions for string output */
 }  
 #endif  
 #endif  
   
 #ifdef SPRINT  
 void soutput_init(s)  void soutput_init(s)
 char *s;  char *s;
 {  {
         s[0] = 0; buf = s;    s[0] = 0; buf = s;
 }  }
   
 #if defined(PARI)  void sprintbf(BF a)
 #include "genpari.h"  
   
 void myoutbrute(g)  
 GEN g;  
 {  {
         bruteall(g,'f',-1,1);    int dprec;
     char fbuf[BUFSIZ];
     char *s;
     dprec = a->body->_mpfr_prec*0.30103;
     sprintf(fbuf,"%%.%dR%c",dprec,(double_output==1)?'f':(double_output==2)?'e':'g');
     mpfr_asprintf(&s,fbuf,a->body);
     TAIL PUTS(s);
     mpfr_free_str(s);
 }  }
   
 void sprintbf(a)  
 BF a;  
 {  
         char *str;  
         char *GENtostr();  
         char *GENtostr0();  
   
         if ( double_output ) {  
                 str = GENtostr0(a->body,myoutbrute);  
         } else {  
                 str = GENtostr(a->body);  
         }  
         TAIL PRINTF(OUT,"%s",str);  
         free(str);  
 }  
 #endif  
 #endif  
   
 void PRINTEXPR(vl,p)  
 VL vl;  
 Obj p;  
 {  
         if ( !p ) {  
                 PRINTR(vl,(R)p);  
                 return;  
         }  
   
         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;  
                 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;  
                 default:  
                         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]);  
                 }  
         }  
 }  
   
 #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  
   
 #if defined(FPRINT)  
 #if 0  
 int     printmode = PRINTF_E;  
 #else  
 int     printmode = PRINTF_G;  
 #endif  
 #else  
 extern int printmode;  
 #endif  
   
 void PRINTNUM(q)  
 Num q;  
 {  
         char real_format[20];  
   
         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)  #if defined(INTERVAL)
                                 case MID_PRINTF_E:  void sprintbf4itv(BF a)
 #endif  
                                         TAIL PRINTF(OUT,"%.16e",BDY((Real)q));  
                                         break;  
                                 case PRINTF_G:  
 #if defined(INTERVAL)  
                                 case MID_PRINTF_G:  
 #endif  
                                 default:  
                                         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 0:  
                                         PUTS("[");  
                                         PRINTNUM(INF((Itv)q));  
                                         PUTS(",");  
                                         PRINTNUM(SUP((Itv)q));  
                                         PUTS("]");  
                                         break;  
                                 case 1:  
                                         PUTS("intval(");  
                                         PRINTNUM(INF((Itv)q));  
                                         PUTS(",");  
                                         PRINTNUM(SUP((Itv)q));  
                                         PUTS(")");  
                                         break;  
                                 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 0:  
                                                         TAIL PRINTF(OUT, "[%.16e,%.16e]",  
                                                                 INF((IntervalDouble)q),SUP((IntervalDouble)q));  
                                                         break;  
                                                 case 1:  
                                                         TAIL PRINTF(OUT, "intval(%.16e,%.16e)",  
                                                                 INF((IntervalDouble)q),SUP((IntervalDouble)q));  
                                                         break;  
                                                 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 0:  
                                                         TAIL PRINTF(OUT, "<%g,%g>",  
                                                                 (SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,  
                                                                 (SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5);  
                                                         break;  
                                                 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;  
                                                 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 0:  
                                                         TAIL PRINTF(OUT, "<%.16e,%.16e>",  
                                                         (SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,  
                                                         (SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5);  
                                                         break;  
                                                 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;  
                                                 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 0:  
                                                         TAIL PRINTF(OUT, "[%g,%g]",  
                                                                 INF((IntervalDouble)q),SUP((IntervalDouble)q));  
                                                         break;  
                                                 case 1:  
                                                         TAIL PRINTF(OUT, "intval(%g,%g)",  
                                                                 INF((IntervalDouble)q),SUP((IntervalDouble)q));  
                                                         break;  
                                                 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;  
         }  
 }  
   
 void PRINTCPLX(a)  
 C a;  
 {  {
         PUTS("(");    int dprec;
         if ( a->r )    char fbuf[BUFSIZ];
                 PRINTNUM(a->r);    char *s;
         if ( a->i ) {    dprec = a->body->_mpfr_prec*0.30103;
 #if defined(INTERVAL)    dprec += 1;
                 if ( a->r && ((compnum(0,a->i,0) > 0)    sprintf(fbuf,"%%.%dR%c",dprec,(double_output==1)?'f':(double_output==2)?'e':'g');
                         || NID(a->i) == N_IP || NID(a->i) == N_IntervalDouble    mpfr_asprintf(&s,fbuf,a->body);
                         || NID(a->i) == N_IntervalQuad || NID(a->i) == N_IntervalBigFloat) )    TAIL PUTS(s);
 #else    mpfr_free_str(s);
                 if ( a->r && (compnum(0,a->i,0) > 0) )  
 #endif  
                         PUTS("+");  
                 PRINTNUM(a->i); PUTS("*@i");  
         }  
         PUTS(")");  
 }  }
   #endif
 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("+");  
                                 }  
                         }  
                 }  
 }  
   
 int hideargs;  
   
 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 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 0:  
                         PUTS("[ ");  
                         for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) {  
                                 PRINTEXPR(vl,ptr[i]); PUTS(" ");  
                         }  
                         PUTS("]");  
                         break;  
                 case 1:  
                         PUTS("vect(");  
                         for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) {  
                                 if ( i != 0 ) PUTS(", ");  
                                 PRINTEXPR(vl,ptr[i]);  
                         }  
                         PUTS(")");  
                         break;  
                 default:  
                         PUTS("[ ");  
                         for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) {  
                                 PRINTEXPR(vl,ptr[i]); PUTS(" ");  
                         }  
                         PUTS("]");  
                         break;  
         }  
 }  
   
 void PRINTMAT(vl,mat)  
 VL vl;  
 MAT mat;  
 {  
         int i,j,r,c;  
         pointer *ptr;  
   
         switch ( outputstyle ) {  
                 case 0:  
                         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;  
                 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;  
                 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;  
         }  
 }  
   
 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 PRINTSTR(str)  
 STRING str;  
 {  
         char *p;  
   
         for ( p = BDY(str); *p; p++ )  
                 if ( *p == '"' )  
                         PUTS("\"");  
                 else {  
                         TAIL PRINTF(OUT,"%c",*p);  
                 }  
 }  
   
 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 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 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 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 ) {  
                 fnodetotree(BDY(quote),&list);  
                 PRINTEXPR(vl,(Obj)list);  
         } else {  
                 PUTS("<...quoted...>");  
         }  
 }  
   
 void PRINTERR(vl,e)  
 VL vl;  
 ERR e;  
 {  
         PUTS("error("); PRINTEXPR(vl,e->body); PUTS(")");  
 }  
   
 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 PRINTLF(vl,f)  
 VL vl;  
 F f;  
 {  
         switch ( FOP(f) ) {  
                 case AL_TRUE:  
                         TAIL PRINTF(OUT,"@true");  
                         break;  
                 case AL_FALSE:  
                         TAIL PRINTF(OUT,"@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 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));  
         }  
 }  
Legend:
Removed from v.1.23  
changed lines
  Added in v.1.43
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>