[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.10 and 1.21

version 1.10, 2001/03/09 01:44:10 version 1.21, 2003/02/14 22:29:15
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/spexpr.c,v 1.9 2000/12/22 10:03:31 saito Exp $   * $OpenXM: OpenXM_contrib2/asir2000/io/spexpr.c,v 1.20 2002/09/11 07:23:27 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "al.h"  #include "al.h"
Line 64  FILE *asir_out;
Line 64  FILE *asir_out;
 char DFORMAT[BUFSIZ];  char DFORMAT[BUFSIZ];
 int hex_output;  int hex_output;
 int fortran_output;  int fortran_output;
   int double_output;
   int real_digit;
   int print_quote;
   
 #define TAIL  #define TAIL
 #define PUTS(s) fputs(s,OUT)  #define PUTS(s) fputs(s,OUT)
Line 76  int fortran_output;
Line 79  int fortran_output;
 #define PRINTV printv  #define PRINTV printv
 #define PRINTEXPR printexpr  #define PRINTEXPR printexpr
 #define PRINTNUM printnum  #define PRINTNUM printnum
 #define PRINTP printp  #define PRINTP asir_printp
 #define PRINTR printr  #define PRINTR printr
 #define PRINTLIST printlist  #define PRINTLIST printlist
 #define PRINTVECT printvect  #define PRINTVECT printvect
Line 88  int fortran_output;
Line 91  int fortran_output;
 #define PRINTGF2MAT printgf2mat  #define PRINTGF2MAT printgf2mat
 #define PRINTGFMMAT printgfmmat  #define PRINTGFMMAT printgfmmat
 #define PRINTBYTEARRAY printbytearray  #define PRINTBYTEARRAY printbytearray
   #define PRINTQUOTE printquote
 #define PRINTERR printerr  #define PRINTERR printerr
 #define PRINTLF printlf  #define PRINTLF printlf
 #define PRINTLOP printlop  #define PRINTLOP printlop
Line 95  int fortran_output;
Line 99  int fortran_output;
 #define PRINTEOP printeop  #define PRINTEOP printeop
 #define PRINTQOP printqop  #define PRINTQOP printqop
 #define PRINTUP printup  #define PRINTUP printup
   #define PRINTUM printum
   #define PRINTSF printsf
 #endif  #endif
   
 #ifdef SPRINT  #ifdef SPRINT
Line 105  extern int hex_output;
Line 111  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 print_quote;
   
 #define TAIL while ( *OUT ) OUT++;  #define TAIL while ( *OUT ) OUT++;
 #define PUTS(s) strcat(OUT,s)  #define PUTS(s) strcat(OUT,s)
Line 117  extern int real_digit;
Line 124  extern int real_digit;
 #define PRINTV sprintv  #define PRINTV sprintv
 #define PRINTEXPR sprintexpr  #define PRINTEXPR sprintexpr
 #define PRINTNUM sprintnum  #define PRINTNUM sprintnum
 #define PRINTP sprintp  #define PRINTP asir_sprintp
 #define PRINTR sprintr  #define PRINTR sprintr
 #define PRINTLIST sprintlist  #define PRINTLIST sprintlist
 #define PRINTVECT sprintvect  #define PRINTVECT sprintvect
Line 129  extern int real_digit;
Line 136  extern int real_digit;
 #define PRINTGF2MAT sprintgf2mat  #define PRINTGF2MAT sprintgf2mat
 #define PRINTGFMMAT sprintgfmmat  #define PRINTGFMMAT sprintgfmmat
 #define PRINTBYTEARRAY sprintbytearray  #define PRINTBYTEARRAY sprintbytearray
   #define PRINTQUOTE sprintquote
 #define PRINTERR sprinterr  #define PRINTERR sprinterr
 #define PRINTLF sprintlf  #define PRINTLF sprintlf
 #define PRINTLOP sprintlop  #define PRINTLOP sprintlop
Line 136  extern int real_digit;
Line 144  extern int real_digit;
 #define PRINTEOP sprinteop  #define PRINTEOP sprinteop
 #define PRINTQOP sprintqop  #define PRINTQOP sprintqop
 #define PRINTUP sprintup  #define PRINTUP sprintup
   #define PRINTUM sprintum
   #define PRINTSF sprintsf
 #endif  #endif
   
 void PRINTEXPR();  void PRINTEXPR();
Line 154  void PRINTUI();
Line 164  void PRINTUI();
 void PRINTGF2MAT();  void PRINTGF2MAT();
 void PRINTGFMMAT();  void PRINTGFMMAT();
 void PRINTBYTEARRAY();  void PRINTBYTEARRAY();
   void PRINTQUOTE();
 void PRINTERR();  void PRINTERR();
 void PRINTCPLX();  void PRINTCPLX();
 void PRINTLM();  void PRINTLM();
 void PRINTLF();  void PRINTLF();
 void PRINTUP2();  void PRINTUP2();
   void PRINTUP();
   void PRINTUM();
   void PRINTFOP();
   void PRINTEOP();
   void PRINTLOP();
   void PRINTQOP();
   void PRINTSF();
   
 #ifdef FPRINT  #ifdef FPRINT
 void output_init() {  void output_init() {
Line 171  P p;
Line 189  P p;
 {  {
         if ( NUM(p) )          if ( NUM(p) )
 #if defined(INTERVAL)  #if defined(INTERVAL)
                 if ( NID(p) != N_IP && compnum(CO,(Num)p,0) < 0 )                  if ( NID(p) != N_IP && NID(p) != N_IntervalDouble && NID(p) != N_IntervalQuad && NID(p) != N_IntervalBigFloat
                           && compnum(CO,(Num)p,0) < 0 )
 #else  #else
                 if ( compnum(CO,(Num)p,0) < 0 )                  if ( compnum(CO,(Num)p,0) < 0 )
 #endif  #endif
Line 184  P p;
Line 203  P p;
                 return (mmono(COEF(DC(p))));                  return (mmono(COEF(DC(p))));
 }  }
   
 #if PARI  #if defined(PARI)
 void printbf(a)  void printbf(a)
 BF a;  BF a;
 {  {
           void sor();
   
         sor(a->body,double_output ? 'f' : 'g',-1,0);          sor(a->body,double_output ? 'f' : 'g',-1,0);
 }  }
 #endif  #endif
Line 200  char *s;
Line 221  char *s;
         s[0] = 0; buf = s;          s[0] = 0; buf = s;
 }  }
   
 #if PARI  #if defined(PARI)
 #include "genpari.h"  #include "genpari.h"
   
 void myoutbrute(g)  void myoutbrute(g)
Line 269  Obj p;
Line 290  Obj p;
                         PRINTGFMMAT(vl,(GFMMAT)p); break;                          PRINTGFMMAT(vl,(GFMMAT)p); break;
                 case O_BYTEARRAY:                  case O_BYTEARRAY:
                         PRINTBYTEARRAY(vl,(BYTEARRAY)p); break;                          PRINTBYTEARRAY(vl,(BYTEARRAY)p); break;
                   case O_QUOTE:
                           PRINTQUOTE(vl,(QUOTE)p); break;
                 default:                  default:
                         break;                          break;
         }          }
Line 305  N n;
Line 328  N n;
         }          }
 }  }
   
 extern int      printmode;  #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)  void PRINTNUM(q)
 Num q;  Num q;
 {  {
Line 331  Num q;
Line 386  Num q;
 #if defined(INTERVAL)  #if defined(INTERVAL)
                                 case MID_PRINTF_E:                                  case MID_PRINTF_E:
 #endif  #endif
                                         TAIL PRINTF(OUT,"%.15e",BDY((Real)q));                                          TAIL PRINTF(OUT,"%.16e",BDY((Real)q));
                                         break;                                          break;
                                 case PRINTF_G:                                  case PRINTF_G:
 #if defined(INTERVAL)  #if defined(INTERVAL)
                                 case MID_PRINTF_G:                                  case MID_PRINTF_G:
 #endif  #endif
                                 default:                                  default:
                                 if ( real_digit ) {                                          if ( real_digit ) {
                                         sprintf(real_format,                                                  sprintf(real_format,
                                                 double_output?"%%.%df":"%%.%dg",real_digit);                                                          double_output?"%%.%df":"%%.%dg",real_digit);
                                         TAIL PRINTF(OUT,real_format,BDY((Real)q));                                                  TAIL PRINTF(OUT,real_format,BDY((Real)q));
                                 } else {                                          } else {
                                         TAIL PRINTF(OUT,double_output?"%f":"%g",BDY((Real)q));                                                  TAIL PRINTF(OUT,double_output?"%f":"%g",BDY((Real)q));
                                 }                                          }
                                 break;                                          break;
                         }                          }
                         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
 #if defined(INTERVAL)  #if defined(INTERVAL)
                 case N_IP:                  case N_IP:
                 case N_IF:                  case N_IntervalBigFloat:
                         PUTS("[");                          PUTS("[");
                         PRINTNUM(INF((Itv)q));                          PRINTNUM(INF((Itv)q));
                         PUTS(",");                          PUTS(",");
                         PRINTNUM(SUP((Itv)q));                          PRINTNUM(SUP((Itv)q));
                         PUTS("]");                          PUTS("]");
                         break;                          break;
                 case N_ID:                  case N_IntervalDouble:
                         switch (printmode) {                          switch (printmode) {
                                 case PRINTF_E:                                  case PRINTF_E:
                                         TAIL PRINTF(OUT, "[%.16e,%.16e]",INF((ItvD)q),SUP((ItvD)q));                                          TAIL PRINTF(OUT, "[%.16e,%.16e]",INF((IntervalDouble)q),SUP((IntervalDouble)q));
 #if defined(ITVDEBUG)  #if defined(ITVDEBUG)
                                         printbin(INF((ItvD)q));                                          printbin(INF((IntervalDouble)q));
                                         printbin(SUP((ItvD)q));                                          printbin(SUP((IntervalDouble)q));
 #endif  #endif
                                         break;                                          break;
                                 case MID_PRINTF_G:                                  case MID_PRINTF_G:
                                         TAIL PRINTF(OUT, "<%g,%g>", (SUP((ItvD)q)+INF((ItvD)q))*0.5,(SUP((ItvD)q)-INF((ItvD)q))*0.5);                                          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:                                  case MID_PRINTF_E:
                                         TAIL PRINTF(OUT, "<%.16e,%.16e>", (SUP((ItvD)q)+INF((ItvD)q))*0.5,(SUP((ItvD)q)-INF((ItvD)q))*0.5);                                          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:                                  case PRINTF_G:
                                 default:                                  default:
                                         TAIL PRINTF(OUT, "[%g,%g]",INF((ItvD)q),SUP((ItvD)q));                                          TAIL PRINTF(OUT, "[%g,%g]",INF((IntervalDouble)q),SUP((IntervalDouble)q));
                                 break;                                  break;
                         }                          }
                         break;                          break;
Line 401  Num q;
Line 456  Num q;
                 case N_GFPN:                  case N_GFPN:
                         PRINTUP((UP)(((GFPN)q)->body));                          PRINTUP((UP)(((GFPN)q)->body));
                         break;                          break;
                   case N_GFS:
                           TAIL PRINTF(OUT,"@_%d",CONT((GFS)q));
                           break;
                   case N_GFSN:
                           PRINTUM(BDY((GFSN)q));
                           break;
         }          }
 }  }
   
Line 411  C a;
Line 472  C a;
         if ( a->r )          if ( a->r )
                 PRINTNUM(a->r);                  PRINTNUM(a->r);
         if ( a->i ) {          if ( a->i ) {
   #if defined(INTERVAL)
                   if ( a->r && ((compnum(0,a->i,0) > 0)
                           || NID(a->i) == N_IP || NID(a->i) == N_IntervalDouble
                           || NID(a->i) == N_IntervalQuad || NID(a->i) == N_IntervalBigFloat) )
   #else
                 if ( a->r && (compnum(0,a->i,0) > 0) )                  if ( a->r && (compnum(0,a->i,0) > 0) )
   #endif
                         PUTS("+");                          PUTS("+");
                 PRINTNUM(a->i); PUTS("*@i");                  PRINTNUM(a->i); PUTS("*@i");
         }          }
Line 488  V v;
Line 555  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 ) {
                         TAIL PRINTF(OUT,"%s",NAME(pf));                          TAIL PRINTF(OUT,"%s",NAME(pf));
                 } else {                  } else {
Line 678  VL vl;
Line 745  VL vl;
 GFMMAT mat;  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 698  VL vl;
Line 764  VL vl;
 BYTEARRAY array;  BYTEARRAY array;
 {  {
         int len,i;          int len,i;
         unsigned int t;  
         unsigned char *b;          unsigned char *b;
   
         len = array->len;          len = array->len;
Line 711  BYTEARRAY array;
Line 776  BYTEARRAY array;
         PUTS("|");          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)  void PRINTERR(vl,e)
 VL vl;  VL vl;
 ERR e;  ERR e;
Line 777  F f;
Line 856  F f;
         }          }
 }  }
   
 PRINTFOP(vl,f)  void PRINTFOP(vl,f)
 VL vl;  VL vl;
 F f;  F f;
 {  {
Line 792  F f;
Line 871  F f;
         }          }
 }  }
   
 PRINTEOP(vl,f)  void PRINTEOP(vl,f)
 VL vl;  VL vl;
 F f;  F f;
 {  {
Line 819  F f;
Line 898  F f;
         PUTS(")");          PUTS(")");
 }  }
   
 PRINTLOP(vl,f)  void PRINTLOP(vl,f)
 VL vl;  VL vl;
 F f;  F f;
 {  {
Line 845  F f;
Line 924  F f;
         PRINTEXPR(vl,(Obj)FPL(f)); PUTS(op); PUTS("0");          PRINTEXPR(vl,(Obj)FPL(f)); PUTS(op); PUTS("0");
 }  }
   
 PRINTQOP(vl,f)  void PRINTQOP(vl,f)
 VL vl;  VL vl;
 F f;  F f;
 {  {
Line 856  F f;
Line 935  F f;
         PRINTEXPR(vl,(Obj)FQMAT(f)); PUTS(")");          PRINTEXPR(vl,(Obj)FQMAT(f)); PUTS(")");
 }  }
   
 PRINTUP(n)  void PRINTUP(n)
 UP n;  UP n;
 {  {
         int i,d;          int i,d;
Line 887  UP n;
Line 966  UP n;
                         }                          }
                 }                  }
                 PUTS(")");                  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.10  
changed lines
  Added in v.1.21

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>