===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/builtin/gf.c,v
retrieving revision 1.3
retrieving revision 1.8
diff -u -p -r1.3 -r1.8
--- OpenXM_contrib2/asir2000/builtin/gf.c	2000/08/21 08:31:19	1.3
+++ OpenXM_contrib2/asir2000/builtin/gf.c	2001/06/25 10:01:27	1.8
@@ -23,7 +23,7 @@
  * shall be made on your publication or presentation in any form of the
  * results obtained by use of the SOFTWARE.
  * (4) In the event that you modify the SOFTWARE, you shall notify FLL by
- * e-mail at risa-admin@flab.fujitsu.co.jp of the detailed specification
+ * e-mail at risa-admin@sec.flab.fujitsu.co.jp of the detailed specification
  * for such modification or the source code of the modified part of the
  * SOFTWARE.
  * 
@@ -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/gf.c,v 1.2 2000/04/20 02:20:15 noro Exp $ 
+ * $OpenXM: OpenXM_contrib2/asir2000/builtin/gf.c,v 1.7 2001/06/25 04:11:41 noro Exp $ 
 */
 #include "ca.h"
 #include "parse.h"
@@ -66,6 +66,7 @@ void nullspace(UM **,UM,int,int,int *);
 void nullspace_lm(LM **,int,int *);
 void nullspace_gf2n(GF2N **,int,int *);
 void nullspace_gfpn(GFPN **,int,int *);
+void nullspace_gfs(GFS **,int,int *);
 void null_to_sol(int **,int *,int,int,UM *);
 
 void showgfmat(UM **,int);
@@ -73,6 +74,7 @@ void pwr_mod(P,P,V,P,int,N,P *);
 void rem_mod(P,P,V,P,int,P *);
 
 void Pnullspace(),Pgcda_mod(),Pftest(),Presfmain(),Ppwr_mod(),Puhensel();
+void Psfuhensel();
 
 void Pnullspace_ff();
 
@@ -82,6 +84,7 @@ void Plinear_form_to_vect(),Pvect_to_linear_form();
 void solve_linear_equation_gf2n(GF2N **,int,int,int *);
 void linear_form_to_array(P,VL,int,Num *);
 void array_to_linear_form(Num *,VL,int,P *);
+void sfuhensel(P,V,V,GFS *,NODE *);
 
 extern int current_ff;
 
@@ -94,6 +97,7 @@ struct ftab gf_tab[] = {
 	{"resfmain",Presfmain,4},
 	{"pwr_mod",Ppwr_mod,6},
 	{"uhensel",Puhensel,4},
+	{"sfuhensel",Psfuhensel,3},
 	{0,0,0},
 };
 
@@ -178,6 +182,48 @@ NODE *rp;
 	*rp = top;
 }
 
+void Psfuhensel(arg,rp)
+NODE arg;
+LIST *rp;
+{
+	P f;
+	V x,y;
+	NODE n,r;
+	GFS ev;
+
+	f = (P)ARG0(arg);
+	x = VR((P)ARG1(arg));
+	y = VR((P)ARG2(arg));
+	sfuhensel(f,x,y,&ev,&r);
+	MKNODE(n,ev,r);
+	MKLIST(*rp,n);
+}
+
+void sfuhensel(f,x,y,evp,rp)
+P f;
+V x,y;
+GFS *evp;
+NODE *rp;
+{
+	ML lift;
+	int i;
+	P s,u,sf;
+	NODE t,top;
+	struct oVL vl1,vl;
+
+	sfhensel(5,f,x,evp,&sf,&lift);
+
+	vl1.v = y; vl1.next = 0;
+	vl.v = x; vl.next = &vl1;
+
+	for ( i = lift->n-1, top = 0; i >= 0; i-- ) {
+		sfbmtop(lift->bound,lift->c[i],x,y,&s);
+		reorderp(CO,&vl,s,&u);
+		MKNODE(t,u,top); top = t;
+	}
+	*rp = top;
+}
+
 void Presfmain(arg,rp)
 NODE arg;
 LIST *rp;
@@ -809,6 +855,8 @@ LIST *rp;
 			nullspace_gf2n((GF2N **)w,n,ind); break;
 		case FF_GFPN:
 			nullspace_gfpn((GFPN **)w,n,ind); break;
+		case FF_GFS:
+			nullspace_gfs((GFS **)w,n,ind); break;
 		default:
 			error("nullspace_ff : current_ff is not set");
 	}
@@ -947,6 +995,49 @@ int *ind;
 		}
 	}
 }
+
+void nullspace_gfs(mat,n,ind)
+GFS **mat;
+int n;
+int *ind;
+{
+	int i,j,l,s;
+	GFS w,w1,h,inv;
+	GFS *t,*u;
+	GFS one;
+
+	bzero(ind,n*sizeof(int));
+	ind[0] = 0;
+	mqtogfs(ONEM,&one);
+
+	for ( i = j = 0; j < n; i++, j++ ) {
+		for ( ; j < n; j++ ) {
+			for ( l = i; l < n; l++ )
+				if ( mat[l][j] )
+					break;
+			if ( l < n ) {
+				t = mat[i]; mat[i] = mat[l]; mat[l] = t; break;
+			} else
+				ind[j] = 1;
+		}
+		if ( j == n )
+			break;
+		divgfs(one,mat[i][j],&inv);
+		for ( s = j, t = mat[i]; s < n; s++ ) {
+			mulgfs(t[s],inv,&w); t[s] = w;
+		}
+		for ( l = 0; l < n; l++ ) {
+			if ( l == i )
+				continue;
+			u = mat[l]; 
+			chsgngfs(u[j],&h);
+			for ( s = j; s < n; s++ ) {
+				mulgfs(h,t[s],&w); addgfs(w,u[s],&w1); u[s] = w1;
+			}
+		}
+	}
+}
+
 /* p = a(0)vl[0]+a(1)vl[1]+...+a(m-1)vl[m-1]+a(m) -> array = [a(0) a(1) ... a(m)] */
 
 void linear_form_to_array(p,vl,m,array)