version 1.10, 2002/06/21 00:34:21 |
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.9 2001/10/12 06:07:05 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 |
|
|
extern ASCENT,GCDTIME,UFTIME,RESTIME,SQTIME,PRINT$ |
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) |
{ |
{ |