version 1.9, 2001/10/12 06:07:05 |
version 1.11, 2003/02/07 09:25:58 |
|
|
* 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/sp,v 1.8 2000/08/22 05:04:23 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/lib/sp,v 1.10 2002/06/21 00:34:21 noro Exp $ |
*/ |
*/ |
/* |
/* |
sp : functions related to algebraic number fields |
sp : functions related to algebraic number fields |
Line 63 extern ASCENT,GCDTIME,UFTIME,RESTIME,SQTIME,PRINT$ |
|
Line 63 extern ASCENT,GCDTIME,UFTIME,RESTIME,SQTIME,PRINT$ |
|
extern Ord$ |
extern Ord$ |
extern USE_PARI_FACTOR$ |
extern USE_PARI_FACTOR$ |
|
|
|
/* gen_sp can handle non-monic poly */ |
|
|
|
def gen_sp(P) |
|
{ |
|
P = ptozp(P); |
|
V = var(P); |
|
D = deg(P,V); |
|
LC = coef(P,D,V); |
|
F = LC^(D-1)*subst(P,V,V/LC); |
|
/* F must be monic */ |
|
L = sp(F); |
|
return cons(map(subst,car(L),V,LC*V),cdr(L)); |
|
} |
|
|
def sp(P) |
def sp(P) |
{ |
{ |
RESTIME=UFTIME=GCDTIME=SQTIME=0; |
RESTIME=UFTIME=GCDTIME=SQTIME=0; |
Line 1327 def resfctr(F,L,V,N) |
|
Line 1341 def resfctr(F,L,V,N) |
|
N = ptozp(N); |
N = ptozp(N); |
V0 = var(N); |
V0 = var(N); |
DN = diff(N,V0); |
DN = diff(N,V0); |
|
LC = coef(N,deg(N,V0),V0); |
|
LCD = coef(DN,deg(DN,V0),V0); |
for ( I = 0, J = 2, Len = deg(N,V0)+1; I < 5; J++ ) { |
for ( I = 0, J = 2, Len = deg(N,V0)+1; I < 5; J++ ) { |
M = prime(J); |
M = prime(J); |
|
if ( !(LC%M) || !(LCD%M)) |
|
continue; |
G = gcd(N,DN,M); |
G = gcd(N,DN,M); |
if ( !deg(G,V0) ) { |
if ( !deg(G,V0) ) { |
I++; |
I++; |