===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/builtin/poly.c,v
retrieving revision 1.19
retrieving revision 1.25
diff -u -p -r1.19 -r1.25
--- OpenXM_contrib2/asir2000/builtin/poly.c	2003/06/19 07:08:18	1.19
+++ OpenXM_contrib2/asir2000/builtin/poly.c	2016/03/31 05:30:32	1.25
@@ -45,7 +45,7 @@
  * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
  * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
  *
- * $OpenXM: OpenXM_contrib2/asir2000/builtin/poly.c,v 1.18 2003/01/16 00:33:27 noro Exp $ 
+ * $OpenXM: OpenXM_contrib2/asir2000/builtin/poly.c,v 1.24 2011/07/21 04:43:26 noro Exp $ 
 */
 #include "ca.h"
 #include "parse.h"
@@ -53,10 +53,11 @@
 
 void Pranp();
 
-void Pmul_trunc();
+void Pheadsgn();
+void Pmul_trunc(),Pquo_trunc();
 void Pumul(),Pumul_ff(),Pusquare(),Pusquare_ff(),Putmul(),Putmul_ff();
 void Pkmul(),Pksquare(),Pktmul();
-void Pord(), Pcoef0(), Pcoef(), Pdeg(), Pmindeg(), Psetmod();
+void Pord(), Premove_vars(), Pcoef0(), Pcoef(), Pdeg(), Pmindeg(), Psetmod();
 void Pcoef_gf2n();
 void getcoef(), getdeglist(), mergedeglist(), change_mvar(), restore_mvar();
 
@@ -112,6 +113,8 @@ void field_order_ff(N *);
 int current_ff;
 
 struct ftab poly_tab[] = {
+	{"headsgn",Pheadsgn,1},
+	{"quo_trunc",Pquo_trunc,2},
 	{"mul_trunc",Pmul_trunc,3},
 	{"homogeneous_deg",Phomogeneous_deg,-2},
 	{"homogeneous_part",Phomogeneous_part,-3},
@@ -121,6 +124,8 @@ struct ftab poly_tab[] = {
 	{"p_mag",Pp_mag,1},
 	{"maxblen",Pmaxblen,1},
 	{"ord",Pord,-1},
+	{"remove_vars",Premove_vars,1},
+	{"delete_vars",Premove_vars,1},
 	{"coef0",Pcoef0,-3},
 	{"coef",Pcoef,-3},
 	{"coef_gf2n",Pcoef_gf2n,2},
@@ -225,6 +230,14 @@ struct ftab poly_tab[] = {
 	{0,0,0},
 };
 
+void Pheadsgn(NODE arg,Q *rp)
+{
+	int s;
+
+	s = headsgn((P)ARG0(arg));
+	STOQ(s,*rp);
+}
+
 void Pmul_trunc(NODE arg,P *rp)
 {
 	P p1,p2,p,h;
@@ -252,6 +265,37 @@ void Pmul_trunc(NODE arg,P *rp)
 	mulp_trunc(vl,p1,p2,vn,rp);
 }
 
+void Pquo_trunc(NODE arg,P *rp)
+{
+	P p1,p2,p,h;
+	VL vl0,vl1,vl2,tvl,vl;
+	VN vn;
+	int i,n;
+
+	p1 = (P)ARG0(arg);
+	p2 = (P)ARG1(arg);
+	if ( !p1 )
+		*rp = 0;
+	else if ( NUM(p2) )
+		divsp(CO,p1,p2,rp);
+	else {
+		get_vars((Obj)p1,&vl1); get_vars((Obj)p2,&vl2); mergev(CO,vl1,vl2,&vl);
+		for ( tvl = vl, n = 0; tvl; tvl = NEXT(tvl), n++ );
+		vn = (VN) ALLOCA((n+1)*sizeof(struct oVN));
+		for ( i = 0, tvl = vl; i < n; tvl = NEXT(tvl), i++ ) {
+			vn[i].v = tvl->v;
+			vn[i].n = 0;
+		}
+		vn[i].v = 0;
+		vn[i].n = 0;
+		for ( h = p2, i = 0; OID(h) == O_P; h = COEF(DC(h)) ) {
+			for ( ; vn[i].v != VR(h); i++ );
+			vn[i].n = QTOS(DEG(DC(h)));
+		}
+		quop_trunc(vl,p1,p2,vn,rp);
+	}
+}
+
 void Phomogeneous_part(NODE arg,P *rp)
 {
 	if ( argc(arg) == 2 )
@@ -402,7 +446,10 @@ void Pp_mag(NODE arg,Q *rp)
 
 void Pord(NODE arg,LIST *listp)
 {
-	NODE n,tn;
+	NODE n,tn,p,opt;
+	char *key;
+	Obj value;
+	int overwrite=0;
 	LIST l;
 	VL vl,tvl,svl;
 	P t;
@@ -410,6 +457,18 @@ void Pord(NODE arg,LIST *listp)
 	V *va;
 	V v;
 
+	if ( current_option ) {
+		for ( opt = current_option; opt; opt = NEXT(opt) ) {
+			p = BDY((LIST)BDY(opt));
+			key = BDY((STRING)BDY(p));
+			value = (Obj)BDY(NEXT(p));
+			if ( !strcmp(key,"overwrite") && value ) {
+				overwrite = value ? 1 : 0;
+				break;
+			}
+		}
+	}
+
 	if ( argc(arg) ) {
 		asir_assert(ARG0(arg),O_LIST,"ord");
 		for ( vl = 0, i = 0, n = BDY((LIST)ARG0(arg));
@@ -423,22 +482,29 @@ void Pord(NODE arg,LIST *listp)
 				error("ord : invalid argument");
 			VR(tvl) = VR(t);
 		}
-		va = (V *)ALLOCA(i*sizeof(V));
-		for ( j = 0, svl = vl; j < i; j++, svl = NEXT(svl) )
-			va[j] = VR(svl);
-		for ( svl = CO; svl; svl = NEXT(svl) ) {
-			v = VR(svl);
-			for ( j = 0; j < i; j++ )
-				if ( v == va[j] )
-					break;
-			if ( j == i ) {
-				if ( !vl ) {
-					NEWVL(vl); tvl = vl;
-				} else {
-					NEWVL(NEXT(tvl)); tvl = NEXT(tvl);
+		if ( !overwrite ) {
+			va = (V *)ALLOCA(i*sizeof(V));
+			for ( j = 0, svl = vl; j < i; j++, svl = NEXT(svl) )
+				va[j] = VR(svl);
+			for ( svl = CO; svl; svl = NEXT(svl) ) {
+				v = VR(svl);
+				for ( j = 0; j < i; j++ )
+					if ( v == va[j] )
+						break;
+				if ( j == i ) {
+					if ( !vl ) {
+						NEWVL(vl); tvl = vl;
+					} else {
+						NEWVL(NEXT(tvl)); tvl = NEXT(tvl);
+					}
+					VR(tvl) = v;
 				}
-				VR(tvl) = v;
 			}
+		} else {
+			for ( svl = vl; svl; svl = NEXT(svl) ) {
+				if ( svl->v->attr == (pointer)V_PF )
+					((PFINS)svl->v->priv)->pf->ins = 0;
+			}
 		}
 		if ( vl )
 			NEXT(tvl) = 0;
@@ -448,6 +514,40 @@ void Pord(NODE arg,LIST *listp)
 		NEXTNODE(n,tn); MKV(VR(vl),t); BDY(tn) = (pointer)t;
 	}
 	NEXT(tn) = 0; MKLIST(l,n); *listp = l;
+}
+
+void Premove_vars(NODE arg,LIST *listp)
+{
+	NODE l,nd,tnd;
+	V *v,*va;
+	int n,na,i,j;
+	VL vl,vl1;
+	P t;
+	LIST list;
+
+	asir_assert(ARG0(arg),O_LIST,"remove_vars");
+	l = BDY((LIST)ARG0(arg)); n = length(l);
+	v = (V *)ALLOCA(n*sizeof(V));
+	for ( i = 0; i < n; i++, l = NEXT(l) ) 
+		if ( !(t = (P)BDY(l)) || (OID(t) != O_P) )
+			error("ord : invalid argument");
+		else v[i] = VR(t);
+
+	for ( na = 0, vl = CO; vl; vl = NEXT(vl), na++ );
+	va = (V *)ALLOCA(na*sizeof(V));
+	for ( i = 0, vl = CO; i < na; i++, vl = NEXT(vl) ) va[i] = VR(vl);
+	for ( i = 0; i < na; i++ )
+		for ( j = 0; j < n; j++ ) if ( va[i] == v[j] ) va[i] = 0;
+	for ( vl = 0, i = na-1; i >= 0; i-- )
+		if ( va[i] ) {
+			NEWVL(vl1); VR(vl1) = va[i]; NEXT(vl1) = vl; vl = vl1;
+		}
+	CO = vl;
+	for ( nd = 0, vl = CO; vl; vl = NEXT(vl) ) {
+		NEXTNODE(nd,tnd); MKV(VR(vl),t); BDY(tnd) = (pointer)t;
+	}
+	if ( nd ) NEXT(tnd) = 0;
+	MKLIST(list,nd); *listp = list;
 }
 
 void Pcoef0(NODE arg,Obj *rp)