===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/io/pexpr_body.c,v
retrieving revision 1.6
retrieving revision 1.16
diff -u -p -r1.6 -r1.16
--- OpenXM_contrib2/asir2000/io/pexpr_body.c	2004/12/17 03:09:08	1.6
+++ OpenXM_contrib2/asir2000/io/pexpr_body.c	2015/08/06 10:01:52	1.16
@@ -1,4 +1,4 @@
-/* $OpenXM: OpenXM_contrib2/asir2000/io/pexpr_body.c,v 1.5 2004/12/02 13:48:43 noro Exp $ */
+/* $OpenXM: OpenXM_contrib2/asir2000/io/pexpr_body.c,v 1.15 2015/08/04 06:20:45 noro Exp $ */
 
 #define PRINTHAT (fortran_output?PUTS("**"):PUTS("^"))
 
@@ -11,6 +11,7 @@ void PRINTR();
 void PRINTLIST();
 void PRINTVECT();
 void PRINTMAT();
+void PRINTIMAT(); /* IMAT */
 void PRINTSTR();
 void PRINTCOMP();
 void PRINTDP();
@@ -20,6 +21,7 @@ void PRINTGFMMAT();
 void PRINTBYTEARRAY();
 void PRINTQUOTE();
 void PRINTQUOTEARG();
+void PRINTNBP();
 void PRINTERR();
 void PRINTCPLX();
 void PRINTLM();
@@ -60,6 +62,10 @@ Obj p;
 				PRINTVECT(vl,(VECT)p); break;
 			case O_MAT:
 				PRINTMAT(vl,(MAT)p); break;
+/* IMAT */
+			case O_IMAT:
+				PRINTIMAT(vl,(IMAT)p); break;
+/* IMAT */
 			case O_STR:
 				PRINTSTR((STRING)p); break;
 			case O_COMP:
@@ -84,6 +90,8 @@ Obj p;
 				PRINTQUOTE(vl,(QUOTE)p); break;
 			case O_QUOTEARG:
 				PRINTQUOTEARG(vl,(QUOTEARG)p); break;
+			case O_NBP:
+				PRINTNBP(vl,(NBP)p); break;
 			case O_SYMBOL:
 				PRINTSYMBOL((SYMBOL)p); break;
 			case O_RANGE:
@@ -290,6 +298,49 @@ MAT mat;
 	}
 }
 
+/* IMAT */
+#if !defined(CPRINT)
+void PRINTIMAT(vl,xmat)
+VL vl;
+IMAT xmat;
+{
+	int i,j,c;
+	int row, col, cr;
+	IMATC Im;
+	IENT ent;
+
+	row = xmat->row;
+	col = xmat->col;
+	Im = (pointer)xmat->root;
+	if ( xmat->clen == 0 ) {
+		for(j = 0; j< row; j++) {
+			PUTS("( ");
+			for(i = 0; i < col; i++) PUTS("0 ");
+			PUTS(")");
+			if ( j < row - 1 ) PUTS("\n");
+		}
+		return;
+	}
+	c = -1;
+	GetNextIent(&Im, &ent, &c);
+	for(j = 0; j < row; j++) {
+		PUTS("( ");
+		for(i = 0; i < col; i++) {
+			cr = j * col + i;
+			if( ent.cr == cr ) {
+				PRINTEXPR(vl, (pointer)ent.body); PUTS(" ");
+				GetNextIent(&Im, &ent, &c);
+			} else {
+				PUTS("0 ");
+			}
+		}
+		PUTS(")");
+		if ( j < row - 1 ) PUTS("\n");
+	}
+}
+#endif
+/* IMAT */
+
 void PRINTLIST(vl,list)
 VL vl;
 LIST list;
@@ -506,6 +557,14 @@ void PRINTFNODE(FNODE f,int paren)
 			PUTS(((ARF)FA0(f))->name);
 			PRINTFNODE((FNODE)FA2(f),1);
 			break;
+		case I_NARYOP:
+			t = (NODE)FA1(f);
+			PRINTFNODE((FNODE)BDY(t),1);
+			for ( t = NEXT(t); t; t = NEXT(t) )	{
+				PUTS(((ARF)FA0(f))->name);
+				PRINTFNODE((FNODE)BDY(t),1);
+			}
+			break;
 		case I_COP:
 			switch( (cid)FA0(f) ) {
 				case C_EQ: opname = ("=="); break;
@@ -557,7 +616,7 @@ void PRINTFNODE(FNODE f,int paren)
 		/* lists */
 		case I_LIST: PUTS("["); PRINTFNODENODE((NODE)FA0(f)); PUTS("]"); break;
 		/* function */
-		case I_FUNC:
+		case I_FUNC: case I_FUNC_QARG:
 			if ( !strcmp(((FUNC)FA0(f))->name,"@pi") ) PUTS("@pi");
 			else if ( !strcmp(((FUNC)FA0(f))->name,"@e") ) PUTS("@e");
 			else {
@@ -630,7 +689,7 @@ Num q;
 						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)
+#if defined(__i386__) || defined(MIPSEL) || defined(VISUAL) || defined(__MINGW32__) || defined(__MINGW64__) || defined(__alpha) || defined(__FreeBSD__) || defined(__NetBSD__) || defined(__x86_64)
 						u = m[1]; l = m[0];
 #else
 						u = m[0]; l = m[1];
@@ -706,10 +765,8 @@ Num q;
 		case N_A:
 			PUTS("("); PRINTR(ALG,(R)BDY((Alg)q)); PUTS(")");
 			break;
-#if defined(PARI)
 		case N_B:
 			PRINTBF((BF)q); break;
-#endif
 #if defined(INTERVAL)
 		case N_IP:
 		case N_IntervalBigFloat:
@@ -831,6 +888,12 @@ Num q;
 				PRINTN(NM(dn));
 			}
 			break;
+		case N_GZ:
+			mpz_out_str(OUT,10,BDY((GZ)q));
+			break;
+		case N_GQ:
+			mpq_out_str(OUT,10,BDY((GQ)q));
+			break;
 	}
 }
 
@@ -1016,6 +1079,32 @@ VL vl;
 QUOTEARG quote;
 {
 	PUTS("<...quoted...>");
+}
+
+void PRINTNBP(VL vl,NBP p)
+{
+	NODE t;
+	NBM m;
+	int d,i;
+	unsigned int *b;
+	if ( !p ) PUTS("0");
+	else {
+		for ( t = BDY(p); t; t = NEXT(t) ) {
+			m = (NBM)BDY(t);
+			PUTS("(");
+			PRINTEXPR(vl,(Obj)m->c);
+			PUTS(")*");
+			d = m->d;
+			b = m->b;
+			if ( d ) 
+				for ( i = 0; i < d; i++ ) {
+					if ( NBM_GET(b,i) ) PUTS("x");
+					else PUTS("y");
+				}
+			else PUTS("1");
+			if ( NEXT(t) ) PUTS("+");
+		}
+	}
 }
 
 void PRINTTB(VL vl,TB p)