===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/engine/gfs.c,v
retrieving revision 1.12
retrieving revision 1.15
diff -u -p -r1.12 -r1.15
--- OpenXM_contrib2/asir2000/engine/gfs.c	2002/09/27 08:40:49	1.12
+++ OpenXM_contrib2/asir2000/engine/gfs.c	2003/01/06 01:16:39	1.15
@@ -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/engine/gfs.c,v 1.11 2002/09/27 04:24:04 noro Exp $
+ * $OpenXM: OpenXM_contrib2/asir2000/engine/gfs.c,v 1.14 2002/12/18 06:15:40 noro Exp $
 */
 #include "ca.h"
 #include "inline.h"
@@ -661,15 +661,15 @@ int mulremum_enc(int p,int n,UM dp,int a,int b)
 
 void gfs_galois_action(GFS a,Q e,GFS *c)
 {
-	Q q;
+	Q p;
 	int i,k;
 	GFS t,s;
 
 	t = a;
 	k = QTOS(e);
-	STOQ(current_gfs_q,q);
+	STOQ(current_gfs_p,p);
 	for ( i = 0; i < k; i++ ) {
-		pwrgfs(t,q,&s); t = s;
+		pwrgfs(t,p,&s); t = s;
 	}
 	*c = t;
 }
@@ -944,6 +944,21 @@ int cmpgfs(GFS a,GFS b)
 		}
 }
 
+void pthrootgfs(GFS a,GFS *b)
+{
+	Q p;
+	int e,i;
+	GFS t,s;
+
+	STOQ(characteristic_sf(),p);
+	e = extdeg_sf()-1;
+	t = a;
+	for ( i = 0; i < e; i++ ) {
+		pwrgfs(t,p,&s); t = s;
+	}
+	*b = t;
+}
+
 void randomgfs(GFS *r)
 {
 	unsigned int t;
@@ -1117,7 +1132,12 @@ int _itosf(int n)
 {
 	int i;
 
+	/* XXX */
+#if 0
 	n %= current_gfs_p;
+#else
+	n %= current_gfs_q;
+#endif
 	if ( !n )
 		return 0;
 	i = !current_gfs_ntoi ? n : current_gfs_ntoi[n];