===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/engine/Q.c,v
retrieving revision 1.1.1.1
retrieving revision 1.3
diff -u -p -r1.1.1.1 -r1.3
--- OpenXM_contrib2/asir2000/engine/Q.c	1999/12/03 07:39:08	1.1.1.1
+++ OpenXM_contrib2/asir2000/engine/Q.c	2000/05/29 08:54:46	1.3
@@ -1,4 +1,4 @@
-/* $OpenXM: OpenXM/src/asir99/engine/Q.c,v 1.1.1.1 1999/11/10 08:12:26 noro Exp $ */
+/* $OpenXM: OpenXM_contrib2/asir2000/engine/Q.c,v 1.2 2000/04/05 08:32:17 noro Exp $ */
 #include "ca.h"
 #include "base.h"
 #include "inline.h"
@@ -302,6 +302,8 @@ Q *nr;
 	}
 }
 
+/* t = [nC0 nC1 ... nCn] */
+
 void mkbc(n,t)
 int n;
 Q *t;
@@ -315,6 +317,50 @@ Q *t;
 	}
 	for ( ; i <= n; i++ )
 		t[i] = t[n-i];
+}
+
+/*
+ *  Dx^k*x^l = W(k,l,0)*x^l*Dx^k+W(k,l,1)*x^(l-1)*x^(k-1)*+...
+ *
+ *  t = [W(k,l,0) W(k,l,1) ... W(k,l,min(k,l)]
+ *  where W(k,l,i) = i! * kCi * lCi
+ */
+
+void mkwc(k,l,t)
+int k,l;
+Q *t;
+{
+	int i,n,up,low;
+	N nm,d,c;
+
+	n = MIN(k,l);
+	for ( t[0] = ONE, i = 1; i <= n; i++ ) {
+		DM(k-i+1,l-i+1,up,low);
+		if ( up ) {
+			nm = NALLOC(2); PL(nm) = 2; BD(nm)[0] = low; BD(nm)[1] = up;
+		} else {
+			nm = NALLOC(1); PL(nm) = 1; BD(nm)[0] = low;
+		}
+		kmuln(NM(t[i-1]),nm,&d); divin(d,i,&c); NTOQ(c,1,t[i]);
+	}
+}
+
+/* mod m table */
+/* XXX : should be optimized */
+
+void mkwcm(k,l,m,t)
+int k,l,m;
+int *t;
+{
+	int i,n;
+	Q *s;
+
+	n = MIN(k,l);
+	s = (Q *)ALLOCA((n+1)*sizeof(Q));
+	mkwc(k,l,s);
+	for ( i = 0; i <= n; i++ ) {
+		t[i] = rem(NM(s[i]),m);
+	}
 }
 
 void factorial(n,r)