[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.4 and 1.23

version 1.4, 2000/08/22 05:04:17 version 1.23, 2004/07/13 07:59:54
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.3 2000/08/21 08:31:38 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/io/cpexpr.c,v 1.22 2004/05/14 06:02:54 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"
   #if PARI
   #include "genpari.h"
   #endif
   
 extern int hex_output,fortran_output;  #ifndef CPRINT
   #define CPRINT
   #endif
   
 #define PRINTHAT (fortran_output?PUTS("**"):PUTS("^"))  extern int hex_output,fortran_output,double_output,real_digit;
   extern int hideargs,outputstyle;
   static int total_length;
   
 #define TAIL  #define TAIL
 #define PUTS(s) (total_length+=strlen(s))  #define PUTS(s) (total_length+=strlen(s))
Line 76  extern int hex_output,fortran_output;
Line 83  extern int hex_output,fortran_output;
 #define PRINTUI length_ui  #define PRINTUI length_ui
 #define PRINTGF2MAT length_gf2mat  #define PRINTGF2MAT length_gf2mat
 #define PRINTGFMMAT length_gfmmat  #define PRINTGFMMAT length_gfmmat
   #define PRINTBYTEARRAY length_bytearray
   #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 PRINTERR length_err  #define PRINTERR length_err
 #define PRINTLF length_lf  #define PRINTLF length_lf
 #define PRINTLOP length_lop  #define PRINTLOP length_lop
Line 83  extern int hex_output,fortran_output;
Line 97  extern int hex_output,fortran_output;
 #define PRINTEOP length_eop  #define PRINTEOP length_eop
 #define PRINTQOP length_qop  #define PRINTQOP length_qop
 #define PRINTUP length_up  #define PRINTUP length_up
   #define PRINTUM length_um
   #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 PRINTERR();  
 void PRINTCPLX();  
 void PRINTLM();  
 void PRINTLF();  
 void PRINTUP2();  
   
   /* special functions for estimating length */
   
 static int total_length;  static int total_length;
   
 int estimate_length(vl,p)  int estimate_length(VL vl,pointer p)
 VL vl;  
 pointer p;  
 {  {
         total_length = 0;          total_length = 0;
         PRINTEXPR(vl,p);          PRINTEXPR(vl,p);
         return total_length;          return total_length;
 }  }
   
 #if PARI  #if defined(PARI)
 void PRINTBF(a)  void PRINTBF(BF a)
 BF a;  
 {  {
         char *str;          char *str;
         char *GENtostr();          char *GENtostr0();
           void myoutbrute();
   
         str = GENtostr(a->body);          if ( double_output ) {
                   str = GENtostr0(a->body,myoutbrute);
           } else {
                   str = GENtostr(a->body);
           }
         total_length += strlen(str);          total_length += strlen(str);
         free(str);          free(str);
 }  }
 #endif  #endif
   
 void PRINTEXPR(vl,p)  void PRINTNUM(Num q)
 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;  
                 default:  
                         break;  
         }  
 }  
   
 void PRINTN(n)  
 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(q)  
 Num q;  
 {  
         if ( !q ) {          if ( !q ) {
                 PUTS("0");                  PUTS("0");
                 return;                  return;
Line 204  Num q;
Line 149  Num q;
                         }                          }
                         break;                          break;
                 case N_R:                  case N_R:
                         total_length += 20; /* XXX */                          if ( double_output )
                                   total_length += 400+real_digit; /* XXX */
                           else
                                   total_length += 20+real_digit; /* XXX */
                         break;                          break;
                 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 PARI  #if defined(PARI)
                 case N_B:                  case N_B:
                         PRINTBF((BF)q); break;                          PRINTBF((BF)q); break;
 #endif  #endif
Line 229  Num q;
Line 177  Num q;
                 case N_GFPN:                  case N_GFPN:
                         PRINTUP((UP)(((GFPN)q)->body));                          PRINTUP((UP)(((GFPN)q)->body));
                         break;                          break;
                   case N_GFS:
                           total_length += 13; /* XXX */
                           break;
                   case N_GFSN:
                           PRINTUM(BDY((GFSN)q));
                           break;
                 default:                  default:
                         break;                          break;
         }          }
 }  }
   
 void PRINTCPLX(a)  void PRINTV(VL vl,V v)
 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,p)  
 VL vl;  
 P p;  
 {  
         V v;  
         DCP dc;  
         int t;  
   
         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,v)  
 VL vl;  
 V v;  
 {  
         PF pf;          PF pf;
         PFAD ad;          PFAD ad;
         int i;          int i;
Line 319  V v;
Line 199  V v;
         else if ( (vid)v->attr == V_PF ) {          else if ( (vid)v->attr == V_PF ) {
                 pf = ((PFINS)v->priv)->pf; ad = ((PFINS)v->priv)->ad;                  pf = ((PFINS)v->priv)->pf; ad = ((PFINS)v->priv)->ad;
                 if ( !strcmp(NAME(pf),"pow") ) {                  if ( !strcmp(NAME(pf),"pow") ) {
                         PUTS("("); PRINTR(vl,(R)ad[0].arg); PUTS(")"); PRINTHAT; PUTS("(");                          PUTS("(("); PRINTR(vl,(R)ad[0].arg); PUTS(")"); PRINTHAT; PUTS("(");
                         PRINTR(vl,(R)ad[1].arg); PUTS(")");                          PRINTR(vl,(R)ad[1].arg); PUTS("))");
                 } else if ( !pf->argc )                  } else if ( !pf->argc )
                         PUTS(NAME(pf));                          PUTS(NAME(pf));
                 else {                  else {
Line 364  V v;
Line 244  V v;
         }          }
 }  }
   
 void PRINTR(vl,a)  void PRINTN(N n)
 VL vl;  
 R a;  
 {  {
         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,vect)  void PRINTSTR(STRING str)
 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,mat)  
 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,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;          char *p;
   
         for ( p = BDY(str); *p; p++ )          for ( p = BDY(str); *p; p++ )
Line 440  STRING str;
Line 268  STRING str;
                 }                  }
 }  }
   
 void PRINTCOMP(vl,c)  void PRINTDP(VL vl,DP d)
 VL vl;  
 COMP c;  
 {  {
         int n,i;          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;          MP m;
         DL dl;          DL dl;
   
Line 476  DP d;
Line 286  DP d;
         }          }
 }  }
   
 void PRINTUI(vl,u)  void PRINTUI(VL vl,USINT u)
 VL vl;  
 USINT u;  
 {  {
         total_length += 10;          total_length += 10;
 }  }
   
 void PRINTGF2MAT(vl,mat)  void PRINTGFMMAT(VL vl,GFMMAT 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;          int row,col,i,j;
         unsigned int t;  
         unsigned int **b;          unsigned int **b;
   
         row = mat->row;          row = mat->row;
Line 528  GFMMAT mat;
Line 308  GFMMAT mat;
         }          }
 }  }
   
 void PRINTERR(vl,e)  void PRINTBYTEARRAY(VL vl,BYTEARRAY array)
 VL vl;  
 ERR e;  
 {  {
         PUTS("error("); PRINTEXPR(vl,e->body); PUTS(")");          /* |xx xx ... xx| */
           total_length += 1+3*array->len;
 }  }
   
 void PRINTUP2(p)  extern int print_quote;
 UP2 p;  
   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...> */
           total_length += 20;
   }
   
   void PRINTSYMBOL(SYMBOL sym)
   {
           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)
   {
         int d,i;          int d,i;
   
         if ( !p ) {          if ( !p ) {
Line 568  UP2 p;
Line 385  UP2 p;
         }          }
 }  }
   
 void PRINTLF(vl,f)  void PRINTQOP(VL vl,F 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;  
         }  
 }  
   
 PRINTFOP(vl,f)  
 VL vl;  
 F f;  
 {  
         char *op;          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(")");  
         }  
 }  
   
 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(")");  
 }  
   
 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");  
 }  
   
 PRINTQOP(vl,f)  
 VL vl;  
 F f;  
 {  
         char *op;  
   
         op = FOP(f)==AL_EX?"ex":"all";          op = FOP(f)==AL_EX?"ex":"all";
         PUTS(op); PUTS(NAME(FQVR(f)));          PUTS(op); PUTS(NAME(FQVR(f)));
         total_length += 2;          total_length += 2;
         PRINTEXPR(vl,(Obj)FQMAT(f)); PUTS(")");          PRINTEXPR(vl,(Obj)FQMAT(f)); PUTS(")");
 }  }
   
 PRINTUP(n)  void PRINTUP(UP n)
 UP n;  
 {  {
         int i,d;          int i,d;
   
Line 710  UP n;
Line 428  UP n;
                 }                  }
                 PUTS(")");                  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; total_length += 13;
                   }
                   for ( i = d-1; i >= 0; i-- ) {
                           if ( n->c[i] ) {
                                   PUTS("+("); PRINTSF(n->c[i]); PUTS(")");
                                   if ( i >= 2 ) {
                                           PUTS("*@s"); PRINTHAT; total_length += 13;
                                   } else if ( i == 1 )
                                           PUTS("*@s");
                           }
                   }
                   PUTS(")");
           }
   }
   
   void PRINTSF(unsigned int i)
   {
           if ( !i ) {
                   PUTS("0");
           } else
                   total_length += 15;
   
   #ifndef CPRINT
   #define CPRINT
   #endif
 }  }
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.23
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>