| version 1.13, 2001/11/19 00:57:13 |
version 1.15, 2002/06/12 08:19:04 |
|
|
| * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
| * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
| * |
* |
| * $OpenXM: OpenXM_contrib2/asir2000/lib/gr,v 1.12 2001/11/01 10:00:19 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/lib/gr,v 1.14 2001/11/19 01:40:05 noro Exp $ |
| */ |
*/ |
| extern INIT_COUNT,ITOR_FAIL$ |
extern INIT_COUNT,ITOR_FAIL$ |
| extern REMOTE_MATRIX,REMOTE_NF,REMOTE_VARS$ |
extern REMOTE_MATRIX,REMOTE_NF,REMOTE_VARS$ |
| Line 352 def minipoly(G0,V,O,P,V0) |
|
| Line 352 def minipoly(G0,V,O,P,V0) |
|
| if ( !zero_dim(hmlist(G0,V,O),V,O) ) |
if ( !zero_dim(hmlist(G0,V,O),V,O) ) |
| error("tolex : ideal is not zero-dimensional!"); |
error("tolex : ideal is not zero-dimensional!"); |
| |
|
| |
Pin = P; |
| |
P = ptozp(P); |
| |
CP = sdiv(P,Pin); |
| G1 = cons(V0-P,G0); |
G1 = cons(V0-P,G0); |
| O1 = [[0,1],[O,length(V)]]; |
O1 = [[0,1],[O,length(V)]]; |
| V1 = cons(V0,V); |
V1 = cons(V0,V); |
| Line 372 def minipoly(G0,V,O,P,V0) |
|
| Line 375 def minipoly(G0,V,O,P,V0) |
|
| TL = cons(V0^J,TL); |
TL = cons(V0^J,TL); |
| NF = gennf(G1,TL,V1,O1,V0,1)[0]; |
NF = gennf(G1,TL,V1,O1,V0,1)[0]; |
| R = tolex_main(V1,O1,NF,[MP],M,MB); |
R = tolex_main(V1,O1,NF,[MP],M,MB); |
| return R[0]; |
return ptozp(subst(R[0],V0,CP*V0)); |
| } |
} |
| } |
} |
| |
|
| Line 1575 def gbcheck(B,V,O) |
|
| Line 1578 def gbcheck(B,V,O) |
|
| { |
{ |
| dp_ord(O); |
dp_ord(O); |
| D = map(dp_ptod,B,V); |
D = map(dp_ptod,B,V); |
| L = dp_gr_checklist(D); |
L = dp_gr_checklist(D,length(V)); |
| DP = L[0]; Plist = L[1]; |
DP = L[0]; Plist = L[1]; |
| for ( IL = [], I = size(DP)[0]-1; I >= 0; I-- ) |
for ( IL = [], I = size(DP)[0]-1; I >= 0; I-- ) |
| IL = cons(I,IL); |
IL = cons(I,IL); |