=================================================================== RCS file: /home/cvs/OpenXM_contrib2/asir2000/lib/primdec_mod,v retrieving revision 1.9 retrieving revision 1.10 diff -u -p -r1.9 -r1.10 --- OpenXM_contrib2/asir2000/lib/primdec_mod 2003/04/24 07:54:15 1.9 +++ OpenXM_contrib2/asir2000/lib/primdec_mod 2003/05/07 06:26:51 1.10 @@ -1,4 +1,4 @@ -/* $OpenXM: OpenXM_contrib2/asir2000/lib/primdec_mod,v 1.8 2003/04/21 02:02:16 noro Exp $ */ +/* $OpenXM: OpenXM_contrib2/asir2000/lib/primdec_mod,v 1.9 2003/04/24 07:54:15 noro Exp $ */ extern Hom,GBTime$ extern DIVLIST,INTIDEAL,ORIGINAL,ORIGINALDIMENSION,STOP,Trials,REM$ @@ -2204,7 +2204,7 @@ def partial_decomp(B,V) map(ox_cmo_rpc,ParallelMinipoly,"setmod_ff",characteristic_ff(),extdeg_ff()); map(ox_pop_cmo,ParallelMinipoly); } - B = map(ptosfp,B); + B = map(simp_ff,B); B = dp_gr_f_main(B,V,0,0); R = partial_decomp0(B,V,length(V)-1); if ( PartialDecompByLex ) { @@ -2377,7 +2377,7 @@ def minipoly_sf_by_buchberger(G,V,O,F,V0,Server) if ( Server ) ox_sync(0); Vc = cons(V0,setminus(vars(G),V)); - Gf = cons(ptosfp(V0-F),G); + Gf = cons(simp_ff(V0-F),G); Vf = append(V,Vc); Gelim = dp_gr_f_main(Gf,Vf,1,[[0,length(V)],[0,length(Vc)]]); for ( Gc = [], T = Gelim; T != []; T = cdr(T) ) { @@ -2412,19 +2412,19 @@ def minipoly_sf_0dim(G,V,O,F,V0,Server) for ( I = Len - 1, GI = []; I >= 0; I-- ) GI = cons(I,GI); MB = dp_mbase(HL); DIM = length(MB); UT = newvect(DIM); - U = dp_ptod(ptosfp(F),V); + U = dp_ptod(simp_ff(F),V); U = dp_nf_f(GI,U,PS,1); for ( I = 0; I < DIM; I++ ) UT[I] = [MB[I],dp_nf_f(GI,U*MB[I],PS,1)]; - T = dp_ptod(ptosfp(1),[V0]); - TT = dp_ptod(ptosfp(1),V); + T = dp_ptod(simp_ff(1),[V0]); + TT = dp_ptod(simp_ff(1),V); G = H = [[TT,T]]; for ( I = 1; ; I++ ) { if ( dp_gr_print() ) print(".",2); - T = dp_ptod(ptosfp(V0^I),[V0]); + T = dp_ptod(simp_ff(V0^I),[V0]); TT = dp_nf_tab_f(H[0][0],UT); H = cons([TT,T],H); L = dp_lnf_f([TT,T],G); @@ -2443,7 +2443,7 @@ def minipoly_sf_rat(G,V,F,V0) Vc = setminus(vars(G),V); Gf = cons(V0-F,G); Vf = append(V,[V0]); - G3 = dp_gr_f_main(map(ptosfp,Gf),Vf,0,3); + G3 = dp_gr_f_main(map(simp_ff,Gf),Vf,0,3); for ( T = G3; T != []; T = cdr(T) ) { Vt = setminus(vars(car(T)),Vc); if ( Vt == [V0] ) @@ -2826,7 +2826,7 @@ def henleq_gsl_sfrat(L,B,Vc,Eval) X = map(subst,X,V0,V0-E0); if ( zerovector(RESTA*X+RESTB) ) { if ( dp_gr_print() ) print("end",0); - return [X,ptosfp(1)]; + return [X,simp_ff(1)]; } else return 0; } else if ( COUNT == CCC ) { @@ -2889,7 +2889,7 @@ def henleq_gsl_sfrat_higher(L,B,Vc,Eval) X = map(mshift,X,Vc,E,-1); if ( zerovector(RESTA*X+RESTB) ) { if ( dp_gr_print() ) print("end",0); - return [X,ptosfp(1)]; + return [X,simp_ff(1)]; } else return 0; } else if ( COUNT == CCC ) { @@ -2999,7 +2999,7 @@ def polyvtoratv_higher(Vect,Vc,K) def polytorat_gcd(F,V,K) { if ( deg(F,V) < K ) - return [F,ptosfp(1)]; + return [F,simp_ff(1)]; F1 = Mod^(K*2); F2 = F; B1 = 0; B2 = 1; while ( 1 ) { @@ -3031,7 +3031,7 @@ def polytorat_gcd(F,V,K) def polytorat(F,V,Mat,K) { if ( deg(F,V) < K ) - return [F,ptosfp(1)]; + return [F,simp_ff(1)]; for ( I = 0; I < K; I++ ) for ( J = 0; J < K; J++ ) Mat[I][J] = coef(F,I+K-J); @@ -3053,7 +3053,7 @@ def polytorat_higher(F,V,K) { if ( K < 2 ) return 0; if ( homogeneous_deg(F) < K ) - return [F,ptosfp(1)]; + return [F,simp_ff(1)]; D = create_icpoly(V,K); C = extract_coef(D*F,V,K,2*K); Vc = vars(C);