===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/lib/bfct,v
retrieving revision 1.21
retrieving revision 1.23
diff -u -p -r1.21 -r1.23
--- OpenXM_contrib2/asir2000/lib/bfct	2002/01/30 02:12:58	1.21
+++ OpenXM_contrib2/asir2000/lib/bfct	2003/04/20 11:59:57	1.23
@@ -45,14 +45,47 @@
  * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
  * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
  *
- * $OpenXM: OpenXM_contrib2/asir2000/lib/bfct,v 1.20 2002/01/29 05:37:12 noro Exp $
+ * $OpenXM: OpenXM_contrib2/asir2000/lib/bfct,v 1.22 2003/04/20 08:54:28 noro Exp $
  */
 /* requires 'primdec' */
 
+#define TMP_S ssssssss
+#define TMP_T tttttttt
+#define TMP_Y1 yyyyyyyy1
+#define TMP_Y2 yyyyyyyy2
+
+extern LIBRARY_GR_LOADED$
+extern LIBRARY_PRIMDEC_LOADED$
+
+if(!LIBRARY_GR_LOADED) load("gr"); else ; LIBRARY_GR_LOADED = 1$
+if(!LIBRARY_PRIMDEC_LOADED) load("primdec"); else ; LIBRARY_PRIMDEC_LOADED = 1$
+
+/* toplevel */
+
+def bfunction(F)
+{
+	/* XXX */
+	F = replace_vars_f(F);
+
+	V = vars(F);
+	N = length(V);
+	D = newvect(N);
+
+	for ( I = 0; I < N; I++ )
+		D[I] = [deg(F,V[I]),V[I]];
+	qsort(D,compare_first);
+	for ( V = [], I = 0; I < N; I++ )
+		V = cons(D[I][1],V);
+	return bfct_via_gbfct_weight(F,V);
+}
+
 /* annihilating ideal of F^s */
 
 def ann(F)
 {
+	/* XXX */
+	F = replace_vars_f(F);
+
 	V = vars(F);
 	N = length(V);
 	D = newvect(N);
@@ -95,6 +128,9 @@ def ann(F)
 
 def ann0(F)
 {
+	/* XXX */
+	F = replace_vars_f(F);
+
 	V = vars(F);
 	N = length(V);
 	D = newvect(N);
@@ -367,6 +403,9 @@ def initial_part(F,V,MW,W)
 
 def bfct(F)
 {
+	/* XXX */
+	F = replace_vars_f(F);
+
 	V = vars(F);
 	N = length(V);
 	D = newvect(N);
@@ -390,6 +429,9 @@ def bfct(F)
 
 def bfct_via_gbfct(F)
 {
+	/* XXX */
+	F = replace_vars_f(F);
+
 	V = vars(F);
 	N = length(V);
 	D = newvect(N);
@@ -419,6 +461,10 @@ def bfct_via_gbfct(F)
 
 def bfct_via_gbfct_weight(F,V)
 {
+	/* XXX */
+	F = replace_vars_f(F);
+	V = replace_vars_v(V);
+
 	N = length(V);
 	D = newvect(N);
 	Wt = getopt(weight);
@@ -456,6 +502,10 @@ def bfct_via_gbfct_weight(F,V)
 
 def bfct_via_gbfct_weight_1(F,V)
 {
+	/* XXX */
+	F = replace_vars_f(F);
+	V = replace_vars_v(V);
+
 	N = length(V);
 	D = newvect(N);
 	Wt = getopt(weight);
@@ -491,6 +541,10 @@ def bfct_via_gbfct_weight_1(F,V)
 
 def bfct_via_gbfct_weight_2(F,V)
 {
+	/* XXX */
+	F = replace_vars_f(F);
+	V = replace_vars_v(V);
+
 	N = length(V);
 	D = newvect(N);
 	Wt = getopt(weight);
@@ -810,6 +864,33 @@ def w_tdeg(F,V,W)
 	for ( R = 0; T; T = cdr(T) ) {
 		D = dp_td(T);
 		if ( D > R ) R = D;
+	}
+	return R;
+}
+
+def replace_vars_f(F)
+{
+	return subst(F,s,TMP_S,t,TMP_T,y1,TMP_Y1,y2,TMP_Y2);
+}
+
+def replace_vars_v(V)
+{
+	V = replace_var(V,s,TMP_S);
+	V = replace_var(V,t,TMP_T);
+	V = replace_var(V,y1,TMP_Y1);
+	V = replace_var(V,y2,TMP_Y2);
+	return V;
+}
+
+def replace_var(V,X,Y)
+{
+	V = reverse(V);
+	for ( R = []; V != []; V = cdr(V) ) {
+		Z = car(V);
+		if ( Z == X )
+			R = cons(Y,R);
+		else
+			R = cons(Z,R);
 	}
 	return R;
 }