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

Diff for /OpenXM_contrib2/asir2000/io/cpexpr.c between version 1.15 and 1.27

version 1.15, 2003/02/14 22:29:15 version 1.27, 2015/08/04 06:20:45
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/cpexpr.c,v 1.14 2002/09/02 05:16:07 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/io/cpexpr.c,v 1.26 2006/08/27 22:17:28 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "parse.h"  #include "parse.h"
 #include "al.h"  #include "al.h"
 #include "base.h"  #include "base.h"
   
   #ifndef CPRINT
   #define CPRINT
   #endif
   
 extern int hex_output,fortran_output,double_output,real_digit;  extern int hex_output,fortran_output,double_output,real_digit;
   extern int hideargs,outputstyle;
   static int total_length;
   
 #define PRINTHAT (fortran_output?PUTS("**"):PUTS("^"))  
   
 #define TAIL  #define TAIL
 #define PUTS(s) (total_length+=strlen(s))  #define PUTS(s) (total_length+=strlen(s))
 #define PRINTN length_n  #define PRINTN length_n
Line 78  extern int hex_output,fortran_output,double_output,rea
Line 82  extern int hex_output,fortran_output,double_output,rea
 #define PRINTGFMMAT length_gfmmat  #define PRINTGFMMAT length_gfmmat
 #define PRINTBYTEARRAY length_bytearray  #define PRINTBYTEARRAY length_bytearray
 #define PRINTQUOTE length_QUOTE  #define PRINTQUOTE length_QUOTE
   #define PRINTQUOTEARG length_QUOTEARG
   #define PRINTSYMBOL length_SYMBOL
   #define PRINTRANGE length_RANGE
   #define PRINTTB length_TB
   #define PRINTDPV length_DPV
   #define PRINTNBP length_nbp
 #define PRINTERR length_err  #define PRINTERR length_err
 #define PRINTLF length_lf  #define PRINTLF length_lf
 #define PRINTLOP length_lop  #define PRINTLOP length_lop
Line 87  extern int hex_output,fortran_output,double_output,rea
Line 97  extern int hex_output,fortran_output,double_output,rea
 #define PRINTUP length_up  #define PRINTUP length_up
 #define PRINTUM length_um  #define PRINTUM length_um
 #define PRINTSF length_sf  #define PRINTSF length_sf
   #define PRINTFARGS length_fargs
   #define PRINTFNODENODE length_fnodenode
   #define PRINTFNODE length_fnode
   
 void PRINTEXPR();  #include "pexpr_body.c"
 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 PRINTUP();  
 void PRINTUM();  
 void PRINTUP2();  
 void PRINTFOP();  
 void PRINTEOP();  
 void PRINTLOP();  
 void PRINTQOP();  
 void PRINTSF();  
   
   /* special functions for estimating length */
   
 static int total_length;  static int total_length;
   
 int estimate_length(VL vl,pointer p)  int estimate_length(VL vl,pointer p)
Line 127  int estimate_length(VL vl,pointer p)
Line 114  int estimate_length(VL vl,pointer p)
         return total_length;          return total_length;
 }  }
   
 #if defined(PARI)  
 void PRINTBF(BF a)  void PRINTBF(BF a)
 {  {
         char *str;          char *str;
         char *GENtostr();          int dprec;
         char *GENtostr0();          char fbuf[BUFSIZ];
         void myoutbrute();  
   
         if ( double_output ) {          dprec = a->body->_mpfr_prec*0.30103;
                 str = GENtostr0(a->body,myoutbrute);          sprintf(fbuf,"%%.%dR%c",dprec,double_output?'f':'g');
         } else {          mpfr_asprintf(&str,fbuf,a->body);
                 str = GENtostr(a->body);  
         }  
         total_length += strlen(str);          total_length += strlen(str);
         free(str);          mpfr_free_str(str);
 }  }
 #endif  
   
 void PRINTEXPR(VL vl,pointer p)  
 {  
         if ( !p ) {  
                 total_length++;  
                 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)  
 {  
         double ceil();  
   
         if ( !n )  
                 PUTS("0");  
         else if ( hex_output )  
                 total_length += 2+(int)(PL(n)*8);  
         else  
                 total_length += (int)(ceil(0.31*((double)(BSH*PL(n))))+1);  
 }  
   
 void PRINTNUM(Num q)  void PRINTNUM(Num q)
 {  {
           DAlg d;
           DP nm;
           Q dn;
   
         if ( !q ) {          if ( !q ) {
                 PUTS("0");                  PUTS("0");
                 return;                  return;
Line 228  void PRINTNUM(Num q)
Line 155  void PRINTNUM(Num q)
                 case N_A:                  case N_A:
                         PUTS("("); PRINTR(ALG,(R)BDY((Alg)q)); PUTS(")");                          PUTS("("); PRINTR(ALG,(R)BDY((Alg)q)); PUTS(")");
                         break;                          break;
 #if defined(PARI)  
                 case N_B:                  case N_B:
                         PRINTBF((BF)q); break;                          PRINTBF((BF)q); break;
 #endif  
                 case N_C:                  case N_C:
                         PRINTCPLX((C)q); break;                          PRINTCPLX((C)q); break;
                 case N_M:                  case N_M:
Line 254  void PRINTNUM(Num q)
Line 179  void PRINTNUM(Num q)
                 case N_GFSN:                  case N_GFSN:
                         PRINTUM(BDY((GFSN)q));                          PRINTUM(BDY((GFSN)q));
                         break;                          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;
                 default:                  default:
                         break;                          break;
         }          }
 }  }
   
 void PRINTCPLX(C a)  
 {  
         PUTS("(");  
         if ( a->r )  
                 PRINTNUM(a->r);  
         if ( a->i ) {  
                 if ( a->r && (compnum(0,a->i,0) > 0) )  
                         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 ( 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("+");  
                                 }  
                         }  
                 }  
 }  
   
 extern int hideargs;  
   
 void PRINTV(VL vl,V v)  void PRINTV(VL vl,V v)
 {  {
         PF pf;          PF pf;
Line 383  void PRINTV(VL vl,V v)
Line 253  void PRINTV(VL vl,V v)
         }          }
 }  }
   
 void PRINTR(VL vl,R a)  void PRINTN(N n)
 {  {
         if ( !a )          double ceil();
   
           if ( !n )
                 PUTS("0");                  PUTS("0");
           else if ( hex_output )
                   total_length += 2+(int)(PL(n)*8);
         else          else
                 switch (OID(a)) {                  total_length += (int)(ceil(0.31*((double)(BSH*PL(n))))+1);
                         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;  
   
         PUTS("[ ");  
         for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) {  
                 PRINTEXPR(vl,ptr[i]); PUTS(" ");  
         }  
         PUTS("]");  
 }  
   
 void PRINTMAT(VL vl,MAT mat)  
 {  
         int i,j,r,c;  
         pointer *ptr;  
   
         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");  
         }  
 }  
   
 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 PRINTSTR(STRING str)  void PRINTSTR(STRING str)
 {  {
         char *p;          char *p;
Line 450  void PRINTSTR(STRING str)
Line 277  void PRINTSTR(STRING str)
                 }                  }
 }  }
   
 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 PRINTDP(VL vl,DP d)  void PRINTDP(VL vl,DP d)
 {  {
         int n,i;          int n,i;
Line 487  void PRINTUI(VL vl,USINT u)
Line 300  void PRINTUI(VL vl,USINT u)
         total_length += 10;          total_length += 10;
 }  }
   
 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 PRINTGFMMAT(VL vl,GFMMAT mat)  void PRINTGFMMAT(VL vl,GFMMAT mat)
 {  {
         int row,col,i,j;          int row,col,i,j;
Line 533  void PRINTBYTEARRAY(VL vl,BYTEARRAY array)
Line 323  void PRINTBYTEARRAY(VL vl,BYTEARRAY array)
         total_length += 1+3*array->len;          total_length += 1+3*array->len;
 }  }
   
   extern int print_quote;
   
 void PRINTQUOTE(VL vl,QUOTE quote)  void PRINTQUOTE(VL vl,QUOTE quote)
 {  {
           LIST list;
   
           if ( print_quote == 2 ) {
                   PRINTFNODE(BDY(quote),0);
           } else if ( print_quote == 1 ) {
                   /* XXX */
                   fnodetotree(BDY(quote),&list);
                   PRINTEXPR(vl,(Obj)list);
           } else {
                   /* <...quoted...> */
                   total_length += 20;
           }
   }
   
   void PRINTQUOTEARG(VL vl,QUOTEARG quote)
   {
           /* XXX */
         /* <...quoted...> */          /* <...quoted...> */
         total_length += 20;          total_length += 20;
 }  }
   
 void PRINTERR(VL vl,ERR e)  void PRINTSYMBOL(SYMBOL sym)
 {  {
         PUTS("error("); PRINTEXPR(vl,e->body); PUTS(")");          total_length += strlen(sym->name);
 }  }
   
   void PRINTTB(VL vl,TB p)
   {
           int i;
   
           for ( i = 0; i < p->next; i++ ) {
                   total_length += strlen(p->body[i]);
           }
   }
   
 void PRINTUP2(UP2 p)  void PRINTUP2(UP2 p)
 {  {
         int d,i;          int d,i;
Line 576  void PRINTUP2(UP2 p)
Line 394  void PRINTUP2(UP2 p)
         }          }
 }  }
   
 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 PRINTQOP(VL vl,F f)  void PRINTQOP(VL vl,F f)
 {  {
         char *op;          char *op;
Line 742  void PRINTUM(UM n)
Line 472  void PRINTUM(UM n)
         }          }
 }  }
   
   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);
                           PRINTEXPR(vl,(Obj)m->c);
                           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 PRINTSF(unsigned int i)  void PRINTSF(unsigned int i)
 {  {
         if ( !i ) {          if ( !i ) {
                 PUTS("0");                  PUTS("0");
         } else          } else
                 total_length += 15;                  total_length += 15;
   
   #ifndef CPRINT
   #define CPRINT
   #endif
 }  }
Legend:
Removed from v.1.15  
changed lines
  Added in v.1.27
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>