===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/io/pexpr.c,v
retrieving revision 1.2
retrieving revision 1.34
diff -u -p -r1.2 -r1.34
--- OpenXM_contrib2/asir2000/io/pexpr.c	2000/02/07 03:14:39	1.2
+++ OpenXM_contrib2/asir2000/io/pexpr.c	2004/03/17 02:23:45	1.34
@@ -1,15 +1,67 @@
-/* $OpenXM: OpenXM_contrib2/asir2000/io/pexpr.c,v 1.1.1.1 1999/12/03 07:39:11 noro Exp $ */
+/*
+ * Copyright (c) 1994-2000 FUJITSU LABORATORIES LIMITED 
+ * All rights reserved.
+ * 
+ * FUJITSU LABORATORIES LIMITED ("FLL") hereby grants you a limited,
+ * non-exclusive and royalty-free license to use, copy, modify and
+ * redistribute, solely for non-commercial and non-profit purposes, the
+ * computer program, "Risa/Asir" ("SOFTWARE"), subject to the terms and
+ * conditions of this Agreement. For the avoidance of doubt, you acquire
+ * only a limited right to use the SOFTWARE hereunder, and FLL or any
+ * third party developer retains all rights, including but not limited to
+ * copyrights, in and to the SOFTWARE.
+ * 
+ * (1) FLL does not grant you a license in any way for commercial
+ * purposes. You may use the SOFTWARE only for non-commercial and
+ * non-profit purposes only, such as academic, research and internal
+ * business use.
+ * (2) The SOFTWARE is protected by the Copyright Law of Japan and
+ * international copyright treaties. If you make copies of the SOFTWARE,
+ * with or without modification, as permitted hereunder, you shall affix
+ * to all such copies of the SOFTWARE the above copyright notice.
+ * (3) An explicit reference to this SOFTWARE and its copyright owner
+ * shall be made on your publication or presentation in any form of the
+ * results obtained by use of the SOFTWARE.
+ * (4) In the event that you modify the SOFTWARE, you shall notify FLL by
+ * e-mail at risa-admin@sec.flab.fujitsu.co.jp of the detailed specification
+ * for such modification or the source code of the modified part of the
+ * SOFTWARE.
+ * 
+ * THE SOFTWARE IS PROVIDED AS IS WITHOUT ANY WARRANTY OF ANY KIND. FLL
+ * MAKES ABSOLUTELY NO WARRANTIES, EXPRESSED, IMPLIED OR STATUTORY, AND
+ * EXPRESSLY DISCLAIMS ANY IMPLIED WARRANTY OF MERCHANTABILITY, FITNESS
+ * FOR A PARTICULAR PURPOSE OR NONINFRINGEMENT OF THIRD PARTIES'
+ * RIGHTS. NO FLL DEALER, AGENT, EMPLOYEES IS AUTHORIZED TO MAKE ANY
+ * MODIFICATIONS, EXTENSIONS, OR ADDITIONS TO THIS WARRANTY.
+ * UNDER NO CIRCUMSTANCES AND UNDER NO LEGAL THEORY, TORT, CONTRACT,
+ * OR OTHERWISE, SHALL FLL BE LIABLE TO YOU OR ANY OTHER PERSON FOR ANY
+ * DIRECT, INDIRECT, SPECIAL, INCIDENTAL, PUNITIVE OR CONSEQUENTIAL
+ * DAMAGES OF ANY CHARACTER, INCLUDING, WITHOUT LIMITATION, DAMAGES
+ * ARISING OUT OF OR RELATING TO THE SOFTWARE OR THIS AGREEMENT, DAMAGES
+ * FOR LOSS OF GOODWILL, WORK STOPPAGE, OR LOSS OF DATA, OR FOR ANY
+ * DAMAGES, EVEN IF FLL SHALL HAVE BEEN INFORMED OF THE POSSIBILITY OF
+ * SUCH DAMAGES, OR FOR ANY CLAIM BY ANY OTHER PARTY. EVEN IF A PART
+ * 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.33 2004/03/17 02:10:31 noro Exp $
+*/
 #include "ca.h"
 #include "al.h"
 #include "parse.h"
 #include "comp.h"
 #include "base.h"
 
+#if defined(PARI)
+#include "genpari.h"
+#endif
+
 #ifndef FPRINT
 #define FPRINT
 #endif
 
 #define PRINTHAT (fortran_output?PUTS("**"):PUTS("^"))
+extern int outputstyle;
 
 #ifdef FPRINT
 FILE *asir_out;
@@ -17,6 +69,11 @@ FILE *asir_out;
 char DFORMAT[BUFSIZ];
 int hex_output;
 int fortran_output;
+int double_output;
+int real_digit;
+int real_binary;
+int print_quote;
+extern int asir_texmacs;
 
 #define TAIL
 #define PUTS(s) fputs(s,OUT)
@@ -29,7 +86,7 @@ int fortran_output;
 #define PRINTV printv
 #define PRINTEXPR printexpr
 #define PRINTNUM printnum
-#define PRINTP printp
+#define PRINTP asir_printp
 #define PRINTR printr
 #define PRINTLIST printlist
 #define PRINTVECT printvect
@@ -40,6 +97,8 @@ int fortran_output;
 #define PRINTUI printui
 #define PRINTGF2MAT printgf2mat
 #define PRINTGFMMAT printgfmmat
+#define PRINTBYTEARRAY printbytearray
+#define PRINTQUOTE printquote
 #define PRINTERR printerr
 #define PRINTLF printlf
 #define PRINTLOP printlop
@@ -47,6 +106,14 @@ int fortran_output;
 #define PRINTEOP printeop
 #define PRINTQOP printqop
 #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
@@ -55,7 +122,12 @@ static char *buf;
 extern char DFORMAT[BUFSIZ];
 extern int hex_output;
 extern int fortran_output;
+extern int double_output;
+extern int real_digit;
+extern int real_binary;
+extern int print_quote;
 
+
 #define TAIL while ( *OUT ) OUT++;
 #define PUTS(s) strcat(OUT,s)
 #define PRINTF sprintf
@@ -67,7 +139,7 @@ extern int fortran_output;
 #define PRINTV sprintv
 #define PRINTEXPR sprintexpr
 #define PRINTNUM sprintnum
-#define PRINTP sprintp
+#define PRINTP asir_sprintp
 #define PRINTR sprintr
 #define PRINTLIST sprintlist
 #define PRINTVECT sprintvect
@@ -78,6 +150,8 @@ extern int fortran_output;
 #define PRINTUI sprintui
 #define PRINTGF2MAT sprintgf2mat
 #define PRINTGFMMAT sprintgfmmat
+#define PRINTBYTEARRAY sprintbytearray
+#define PRINTQUOTE sprintquote
 #define PRINTERR sprinterr
 #define PRINTLF sprintlf
 #define PRINTLOP sprintlop
@@ -85,6 +159,14 @@ extern int fortran_output;
 #define PRINTEOP sprinteop
 #define PRINTQOP sprintqop
 #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();
@@ -102,11 +184,24 @@ void PRINTDP();
 void PRINTUI();
 void PRINTGF2MAT();
 void PRINTGFMMAT();
+void PRINTBYTEARRAY();
+void PRINTQUOTE();
 void PRINTERR();
 void PRINTCPLX();
 void PRINTLM();
 void PRINTLF();
 void PRINTUP2();
+void PRINTUP();
+void PRINTUM();
+void PRINTFOP();
+void PRINTEOP();
+void PRINTLOP();
+void PRINTQOP();
+void PRINTSF();
+void PRINTSYMBOL();
+void PRINTRANGE();
+void PRINTTB();
+void PRINTFNODE();
 
 #ifdef FPRINT
 void output_init() {
@@ -118,7 +213,12 @@ int mmono(p)
 P p;
 {
 	if ( NUM(p) )
+#if defined(INTERVAL)
+		if ( NID(p) != N_IP && NID(p) != N_IntervalDouble && NID(p) != N_IntervalQuad && NID(p) != N_IntervalBigFloat
+			&& compnum(CO,(Num)p,0) < 0 ) 
+#else
 		if ( compnum(CO,(Num)p,0) < 0 ) 
+#endif
 			return ( 1 );
 		else
 			return ( 0 );
@@ -128,11 +228,11 @@ P p;
 		return (mmono(COEF(DC(p))));
 }
 
-#if PARI
+#if defined(PARI)
 void printbf(a)
 BF a;
 {
-	sor(a->body,'g',-1,0);
+	sor(a->body,double_output ? 'f' : 'g',-1,0);
 }
 #endif
 #endif
@@ -144,63 +244,92 @@ char *s;
 	s[0] = 0; buf = s;
 }
 
-#if PARI
+#if defined(PARI)
+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();
 
-	str = GENtostr(a->body);
+	if ( double_output ) {
+		str = GENtostr0(a->body,myoutbrute);
+	} else {
+		str = GENtostr(a->body);
+	}
 	TAIL PRINTF(OUT,"%s",str);
 	free(str);
 }
 #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;
-		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)
@@ -234,9 +363,45 @@ N n;
 	}
 }
 
+#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)
 Num q;
 {
+	char real_format[20];
+
 	if ( !q ) {
 		PUTS("0");
 		return;
@@ -251,15 +416,191 @@ Num q;
 			}
 			break;
 		case N_R:
-			TAIL PRINTF(OUT,"%g",BDY((Real)q));
+			switch (printmode) {
+				case PRINTF_E:
+#if defined(INTERVAL)
+				case MID_PRINTF_E:
+#endif
+					TAIL PRINTF(OUT,"%.16e",BDY((Real)q));
+					break;
+				case PRINTF_G:
+#if defined(INTERVAL)
+				case MID_PRINTF_G:
+#endif
+				default:
+					if ( real_binary ) {
+						unsigned int *m;
+						unsigned int u,l,mask;
+						int i,expo;
+
+						m = (unsigned int *)&BDY((Real)q);
+#if defined(__i386__) || defined(MIPSEL) || defined(VISUAL) || defined(__alpha) || defined(__FreeBSD__) || defined(__NetBSD__) || defined(__x86_64)
+						u = m[1]; l = m[0];
+#else
+						u = m[0]; l = m[1];
+#endif
+						if (u&0x80000000) {
+							TAIL PRINTF(OUT,"-");
+						}
+						u &= 0x7fffffff;
+						expo = ((int)(u>>20));
+						u &= 0xfffff;
+						if ( expo == 2047 ) {
+							if ( u || l ) {
+								TAIL PRINTF(OUT,"NaN");
+							} else {
+								TAIL PRINTF(OUT,"Inf");
+							}
+						} else if ( expo == 0 ) {
+							if ( u || l ) {
+								TAIL PRINTF(OUT,"0b0.");
+								for ( i = 0, mask = 0x80000; i < 20; 
+									i++, mask >>= 1) {
+									TAIL
+									if ( u&mask )
+										PRINTF(OUT,"1");
+									else
+										PRINTF(OUT,"0");
+								}
+								for ( i = 0, mask = 0x80000000; i < 32;
+									i++, mask >>= 1) {
+									TAIL
+									if ( l&mask )
+										PRINTF(OUT,"1");
+									else
+										PRINTF(OUT,"0");
+								}
+								TAIL PRINTF(OUT,"*2^%d",-1022);
+							} else {
+								TAIL PRINTF(OUT,"0");
+							}
+						} else {
+							expo -= 1023;
+							TAIL PRINTF(OUT,"0b1.");
+							for ( i = 0, mask = 0x80000; i < 20; 
+								i++, mask >>= 1) {
+								TAIL
+								if ( u&mask )
+									PRINTF(OUT,"1");
+								else
+									PRINTF(OUT,"0");
+							}
+							for ( i = 0, mask = 0x80000000; i < 32;
+								i++, mask >>= 1) {
+								TAIL
+								if ( l&mask )
+									PRINTF(OUT,"1");
+								else
+									PRINTF(OUT,"0");
+							}
+							if ( expo ) {
+								TAIL PRINTF(OUT,"*2^%d",expo);
+							}
+						}
+					} else if ( real_digit ) {
+						sprintf(real_format,
+							double_output?"%%.%df":"%%.%dg",real_digit);
+						TAIL PRINTF(OUT,real_format,BDY((Real)q));
+					} else {
+						TAIL PRINTF(OUT,double_output?"%f":"%g",BDY((Real)q));
+					}
+					break;
+			}
 			break;
 		case N_A:
 			PUTS("("); PRINTR(ALG,(R)BDY((Alg)q)); PUTS(")");
 			break;
-#if PARI
+#if defined(PARI)
 		case N_B:
 			PRINTBF((BF)q); break;
 #endif
+#if defined(INTERVAL)
+		case N_IP:
+		case N_IntervalBigFloat:
+			switch ( outputstyle ) {
+				case 1:
+					PUTS("intval(");
+					PRINTNUM(INF((Itv)q));
+					PUTS(",");
+					PRINTNUM(SUP((Itv)q));
+					PUTS(")");
+					break;
+				case 0:
+				default:
+					PUTS("[");
+					PRINTNUM(INF((Itv)q));
+					PUTS(",");
+					PRINTNUM(SUP((Itv)q));
+					PUTS("]");
+					break;
+			}
+			break;
+		case N_IntervalDouble:
+			switch (printmode) {
+				case PRINTF_E:
+					switch ( outputstyle ) {
+						case 1:
+							TAIL PRINTF(OUT, "intval(%.16e,%.16e)",
+								INF((IntervalDouble)q),SUP((IntervalDouble)q));
+							break;
+						case 0:
+						default:
+							TAIL PRINTF(OUT, "[%.16e,%.16e]",
+								INF((IntervalDouble)q),SUP((IntervalDouble)q));
+							break;
+					}
+#if defined(ITVDEBUG)
+					printbin(INF((IntervalDouble)q));
+					printbin(SUP((IntervalDouble)q));
+#endif
+					break;
+				case MID_PRINTF_G:
+					switch ( outputstyle ) {
+						case 1:
+							TAIL PRINTF(OUT, "intvalm(%g,%g)",
+								(SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,
+								(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5);
+							break;
+						case 0:
+						default:
+							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:
+					switch ( outputstyle ) {
+						case 1:
+							TAIL PRINTF(OUT, "intvalm(%.16e,%.16e)",
+							(SUP((IntervalDouble)q)+INF((IntervalDouble)q))*0.5,
+							(SUP((IntervalDouble)q)-INF((IntervalDouble)q))*0.5);
+							break;
+						case 0:
+						default:
+							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:
+				default:
+					switch ( outputstyle ) {
+						case 1:
+							TAIL PRINTF(OUT, "intval(%g,%g)",
+								INF((IntervalDouble)q),SUP((IntervalDouble)q));
+							break;
+						case 0:
+						default:
+							TAIL PRINTF(OUT, "[%g,%g]",
+								INF((IntervalDouble)q),SUP((IntervalDouble)q));
+							break;
+					}
+				break;
+			}
+			break;
+#endif
 		case N_C:
 			PRINTCPLX((C)q); break;
 		case N_M:
@@ -275,6 +616,12 @@ Num q;
 		case N_GFPN:
 			PRINTUP((UP)(((GFPN)q)->body));
 			break;
+		case N_GFS:
+			TAIL PRINTF(OUT,"@_%d",CONT((GFS)q));
+			break;
+		case N_GFSN:
+			PRINTUM(BDY((GFSN)q));
+			break;
 	}
 }
 
@@ -285,7 +632,13 @@ C a;
 	if ( a->r )
 		PRINTNUM(a->r);
 	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) )
+#endif
 			PUTS("+");
 		PRINTNUM(a->i); PUTS("*@i");
 	}
@@ -362,8 +715,8 @@ V v;
 	else if ( (vid)v->attr == V_PF ) {
 		pf = ((PFINS)v->priv)->pf; ad = ((PFINS)v->priv)->ad;
 		if ( !strcmp(NAME(pf),"pow") ) {
-			PUTS("("); PRINTR(vl,(R)ad[0].arg); PUTS(")"); PRINTHAT; PUTS("(");
-			PRINTR(vl,(R)ad[1].arg); PUTS(")");
+			PUTS("(("); PRINTR(vl,(R)ad[0].arg); PUTS(")"); PRINTHAT; PUTS("(");
+			PRINTR(vl,(R)ad[1].arg); PUTS("))");
 		} else if ( !pf->argc ) {
 			TAIL PRINTF(OUT,"%s",NAME(pf));
 		} else {
@@ -426,11 +779,24 @@ VECT vect;
 	int i;
 	pointer *ptr;
 
-	PUTS("[ ");
-	for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) {
-		PRINTEXPR(vl,ptr[i]); PUTS(" ");
+	switch ( outputstyle ) {
+		case 1:
+			PUTS("vect(");
+			for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) {
+				if ( i != 0 ) PUTS(",");
+				PRINTEXPR(vl,ptr[i]);
+			}
+			PUTS(")");
+			break;
+		case 0:
+		default:
+			PUTS("[ ");
+			for ( i = 0, ptr = BDY(vect); i < vect->len; i++ ) {
+				PRINTEXPR(vl,ptr[i]); PUTS(" ");
+			}
+			PUTS("]");
+			break;
 	}
-	PUTS("]");
 }
 
 void PRINTMAT(vl,mat)
@@ -440,14 +806,32 @@ 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");
+	switch ( outputstyle ) {
+		case 1:
+			PUTS("mat(\n");
+			for ( i = 0, r = mat->row, c = mat->col; i < r; i++ ) {
+				if ( i != 0 ) PUTS(",\n");
+				PUTS("[ ");
+				for ( j = 0, ptr = BDY(mat)[i]; j < c; j++ ) {
+					if ( j != 0 ) PUTS(",");
+					PRINTEXPR(vl,ptr[j]);
+				}
+				PUTS(" ]");
+			}
+			PUTS(")");
+			break;
+		case 0:
+		default:
+			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");
+			}
+			break;
 	}
 }
 
@@ -552,7 +936,6 @@ VL vl;
 GFMMAT mat;
 {
 	int row,col,i,j;
-	unsigned int t;
 	unsigned int **b;
 
 	row = mat->row;
@@ -567,6 +950,39 @@ GFMMAT mat;
 	}
 }
 
+void PRINTBYTEARRAY(vl,array)
+VL vl;
+BYTEARRAY array;
+{
+	int len,i;
+	unsigned char *b;
+
+	len = array->len;
+	b = array->body;
+	PUTS("|");
+	for ( i = 0; i < len-1; i++ ) {
+		TAIL PRINTF(OUT,"%02x ",(unsigned int)b[i]);
+	}
+	TAIL PRINTF(OUT,"%02x",(unsigned int)b[i]);
+	PUTS("|");
+}
+
+void PRINTQUOTE(vl,quote)
+VL vl;
+QUOTE quote;
+{
+	LIST list;
+
+	if ( print_quote == 1 ) {
+		fnodetotree(BDY(quote),&list);
+		PRINTEXPR(vl,(Obj)list);
+	} else if ( print_quote == 2 ) {
+		PRINTFNODE(BDY(quote),0);
+	} else {
+		PUTS("<...quoted...>");
+	}
+}
+
 void PRINTERR(vl,e)
 VL vl;
 ERR e;
@@ -633,7 +1049,7 @@ F f;
 	}
 }
 
-PRINTFOP(vl,f)
+void PRINTFOP(vl,f)
 VL vl;
 F f;
 {
@@ -648,7 +1064,7 @@ F f;
 	}
 }
 
-PRINTEOP(vl,f)
+void PRINTEOP(vl,f)
 VL vl;
 F f;
 {
@@ -675,7 +1091,7 @@ F f;
 	PUTS(")");
 }
 
-PRINTLOP(vl,f)
+void PRINTLOP(vl,f)
 VL vl;
 F f;
 {
@@ -701,7 +1117,7 @@ F f;
 	PRINTEXPR(vl,(Obj)FPL(f)); PUTS(op); PUTS("0");
 }
 
-PRINTQOP(vl,f)
+void PRINTQOP(vl,f)
 VL vl;
 F f;
 {
@@ -712,7 +1128,7 @@ F f;
 	PRINTEXPR(vl,(Obj)FQMAT(f)); PUTS(")");
 }
 
-PRINTUP(n)
+void PRINTUP(n)
 UP n;
 {
 	int i,d;
@@ -744,4 +1160,200 @@ UP n;
 		}
 		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));
+	}
+}
+
+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: PUTS("!"); 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(")");
 }