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

version 1.14, 2002/09/02 05:16:07 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.13 2001/10/09 01:36:20 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
   
   #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 85  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 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 99  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 116  int estimate_length(VL vl,pointer p)
         return total_length;          return total_length;
 }  }
   
 #if PARI  #if defined(PARI)
 void PRINTBF(BF a)  void PRINTBF(BF a)
 {  {
         char *str;          char *str;
         char *GENtostr();  
         char *GENtostr0();          char *GENtostr0();
         void myoutbrute();          void myoutbrute();
   
Line 145  void PRINTBF(BF a)
Line 133  void PRINTBF(BF a)
 }  }
 #endif  #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)
 {  {
         if ( !q ) {          if ( !q ) {
Line 228  void PRINTNUM(Num q)
Line 157  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 PARI  #if defined(PARI)
                 case N_B:                  case N_B:
                         PRINTBF((BF)q); break;                          PRINTBF((BF)q); break;
 #endif  #endif
Line 259  void PRINTNUM(Num q)
Line 188  void PRINTNUM(Num q)
         }          }
 }  }
   
 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 244  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 268  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 291  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 314  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 385  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 748  void PRINTSF(unsigned int i)
Line 469  void PRINTSF(unsigned int i)
                 PUTS("0");                  PUTS("0");
         } else          } else
                 total_length += 15;                  total_length += 15;
   
   #ifndef CPRINT
   #define CPRINT
   #endif
 }  }
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.23
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>