=================================================================== RCS file: /home/cvs/OpenXM_contrib2/asir2000/io/pexpr.c,v retrieving revision 1.26 retrieving revision 1.33 diff -u -p -r1.26 -r1.33 --- OpenXM_contrib2/asir2000/io/pexpr.c 2003/12/24 08:00:38 1.26 +++ OpenXM_contrib2/asir2000/io/pexpr.c 2004/03/17 02:10:31 1.33 @@ -44,7 +44,7 @@ * OF THE SOFTWARE HAS BEEN DEVELOPED BY A THIRD PARTY, THE THIRD PARTY * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. - * $OpenXM: OpenXM_contrib2/asir2000/io/pexpr.c,v 1.25 2003/12/02 06:56:48 noro Exp $ + * $OpenXM: OpenXM_contrib2/asir2000/io/pexpr.c,v 1.32 2004/03/03 09:25:30 noro Exp $ */ #include "ca.h" #include "al.h" @@ -52,6 +52,10 @@ #include "comp.h" #include "base.h" +#if defined(PARI) +#include "genpari.h" +#endif + #ifndef FPRINT #define FPRINT #endif @@ -69,6 +73,7 @@ int double_output; int real_digit; int real_binary; int print_quote; +extern int asir_texmacs; #define TAIL #define PUTS(s) fputs(s,OUT) @@ -103,6 +108,12 @@ int print_quote; #define PRINTUP printup #define PRINTUM printum #define PRINTSF printsf +#define PRINTSYMBOL printsymbol +#define PRINTRANGE printrange +#define PRINTTB printtb +#define PRINTFNODE printfnode +#define PRINTFNODENODE printfnodenode +#define PRINTFARGS printfargs #endif #ifdef SPRINT @@ -113,6 +124,7 @@ extern int hex_output; extern int fortran_output; extern int double_output; extern int real_digit; +extern int real_binary; extern int print_quote; @@ -149,6 +161,12 @@ extern int print_quote; #define PRINTUP sprintup #define PRINTUM sprintum #define PRINTSF sprintsf +#define PRINTSYMBOL sprintsymbol +#define PRINTRANGE sprintrange +#define PRINTTB sprinttb +#define PRINTFNODE sprintfnode +#define PRINTFNODENODE sprintfnodenode +#define PRINTFARGS sprintfargs #endif void PRINTEXPR(); @@ -180,6 +198,10 @@ void PRINTEOP(); void PRINTLOP(); void PRINTQOP(); void PRINTSF(); +void PRINTSYMBOL(); +void PRINTRANGE(); +void PRINTTB(); +void PRINTFNODE(); #ifdef FPRINT void output_init() { @@ -210,8 +232,6 @@ P p; void printbf(a) BF a; { - void sor(); - sor(a->body,double_output ? 'f' : 'g',-1,0); } #endif @@ -225,19 +245,20 @@ char *s; } #if defined(PARI) -#include "genpari.h" - void myoutbrute(g) GEN g; { +# if PARI_VERSION_CODE > 131588 + brute(g, 'f',-1); +# else bruteall(g,'f',-1,1); +# endif } void sprintbf(a) BF a; { char *str; - char *GENtostr(); char *GENtostr0(); if ( double_output ) { @@ -251,53 +272,64 @@ BF a; #endif #endif +#define DATA_BEGIN 2 +#define DATA_END 5 + +extern FUNC user_print_function; + void PRINTEXPR(vl,p) VL vl; Obj p; { + if ( asir_texmacs && !user_print_function ) printf("\2verbatim:"); 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: + } else + 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; - } + 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; + case O_SYMBOL: + PRINTSYMBOL((SYMBOL)p); break; + case O_RANGE: + PRINTRANGE(vl,(RANGE)p); break; + case O_TB: + PRINTTB(vl,(TB)p); break; + default: + break; + } + if ( asir_texmacs && !user_print_function ) { putchar('\5'); fflush(stdout); } } void PRINTN(n) @@ -941,9 +973,11 @@ QUOTE quote; { LIST list; - if ( print_quote ) { + if ( print_quote == 1 ) { fnodetotree(BDY(quote),&list); PRINTEXPR(vl,(Obj)list); + } else if ( print_quote == 2 ) { + PRINTFNODE(BDY(quote),0); } else { PUTS("<...quoted...>"); } @@ -1170,4 +1204,156 @@ unsigned int i; } else { TAIL PRINTF(OUT,"@_%d",IFTOF(i)); } +} + +void PRINTSYMBOL(SYMBOL sym) +{ + PUTS(sym->name); +} + +void PRINTRANGE(VL vl,RANGE p) +{ + PUTS("range("); + PRINTEXPR(vl,p->start); + PUTS(","); + PRINTEXPR(vl,p->end); + PUTS(")"); +} + +void PRINTTB(VL vl,TB p) +{ + int i; + + for ( i = 0; i < p->next; i++ ) { + PUTS(p->body[i]); + } +} + +void PRINTFNODENODE(NODE n) +{ + for ( ; n; n = NEXT(n) ) { + PRINTFNODE((FNODE)BDY(n),0); + if ( NEXT(n) ) PUTS(","); + } +} + +void PRINTFARGS(FNODE f) +{ + NODE n; + + if ( f->id == I_LIST ) { + n = (NODE)FA0(f); + PRINTFNODENODE(n); + } else + PRINTFNODE(f,0); +} + +void PRINTFNODE(FNODE f,int paren) +{ + NODE n,t,t0; + char vname[BUFSIZ],prefix[BUFSIZ]; + char *opname,*vname_conv,*prefix_conv; + Obj obj; + int i,len,allzero,elen,elen2; + C cplx; + char *r; + FNODE fi,f2; + + if ( !f ) { + PUTS("(0)"); + return; + } + if ( paren ) PUTS("("); + switch ( f->id ) { + /* unary operators */ + case I_NOT: PRINTFNODE((FNODE)FA0(f),1); break; + case I_PAREN: PRINTFNODE((FNODE)FA0(f),0); break; + case I_MINUS: PUTS("-"); PRINTFNODE((FNODE)FA0(f),1); break; + /* binary operators */ + /* arg list */ + /* I_AND, I_OR => FA0(f), FA1(f) */ + /* otherwise => FA1(f), FA2(f) */ + case I_BOP: + PRINTFNODE((FNODE)FA1(f),1); + PUTS(((ARF)FA0(f))->name); + PRINTFNODE((FNODE)FA2(f),1); + break; + case I_COP: + switch( (cid)FA0(f) ) { + case C_EQ: opname = ("=="); break; + case C_NE: opname = ("!="); break; + case C_GT: opname = (">"); break; + case C_LT: opname = ("<"); break; + case C_GE: opname = (">="); break; + case C_LE: opname = ("<="); break; + } + PRINTFNODE((FNODE)FA1(f),1); + PUTS(opname); + PRINTFNODE((FNODE)FA2(f),1); + break; + case I_LOP: + switch( (lid)FA0(f) ) { + case L_EQ: opname = ("@=="); break; + case L_NE: opname = ("@!="); break; + case L_GT: opname = ("@>"); break; + case L_LT: opname = ("@<"); break; + case L_GE: opname = ("@>="); break; + case L_LE: opname = ("@<="); break; + case L_AND: opname = ("@&&"); break; + case L_OR: opname = ("@||"); break; + case L_NOT: opname = ("@!"); break; + } + if ( (lid)FA0(f)==L_NOT ) { + PUTS(opname); PRINTFNODE((FNODE)FA1(f),1); + } else { + PRINTFNODE((FNODE)FA1(f),1); + PUTS(opname); + PRINTFNODE((FNODE)FA2(f),1); + } + break; + case I_AND: + PRINTFNODE((FNODE)FA0(f),1); + PUTS("&&"); + PRINTFNODE((FNODE)FA1(f),1); + break; + case I_OR: + PRINTFNODE((FNODE)FA0(f),1); + PUTS("!!"); + PRINTFNODE((FNODE)FA1(f),1); + break; + /* ternary operators */ + case I_CE: + PRINTFNODE((FNODE)FA0(f),1); PUTS("?"); PRINTFNODE((FNODE)FA1(f),1); + PUTS(":"); PRINTFNODE((FNODE)FA2(f),1); + break; + /* lists */ + case I_LIST: PUTS("["); PRINTFNODENODE((NODE)FA0(f)); PUTS("]"); break; + /* function */ + case I_FUNC: + if ( !strcmp(((FUNC)FA0(f))->name,"@pi") ) PUTS("@pi"); + else if ( !strcmp(((FUNC)FA0(f))->name,"@e") ) PUTS("@e"); + else { + PUTS(((FUNC)FA0(f))->name); + PUTS("("); PRINTFARGS(FA1(f)); PUTS(")"); + } + break; + /* XXX */ + case I_CAR: PUTS("car("); PRINTFNODE(FA0(f),0); PUTS(")"); break; + case I_CDR: PUTS("cdr("); PRINTFNODE(FA0(f),0); PUTS(")"); break; + /* exponent vector */ + case I_EV: PUTS("<<"); PRINTFNODENODE((NODE)FA0(f)); PUTS(">>"); break; + /* string */ + case I_STR: PUTS((char *)FA0(f)); break; + /* internal object */ + case I_FORMULA: obj = (Obj)FA0(f); PRINTEXPR(CO,obj); break; + /* program variable */ + case I_PVAR: + if ( FA1(f) ) + error("printfnode : not implemented yet"); + GETPVNAME(FA0(f),opname); + PUTS(opname); + break; + default: error("printfnode : not implemented yet"); + } + if ( paren ) PUTS(")"); }