OpenXM_contrib2/asir2000/lib/sp - diff

Return to sp CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / lib

version 1.10, 2002/06/21 00:34:21

version 1.13, 2004/04/13 07:43:20

Line 45

* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,

* 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.12 2004/01/07 08:23:11 noro Exp $

sp : functions related to algebraic number fields

Line 63 extern ASCENT,GCDTIME,UFTIME,RESTIME,SQTIME,PRINT$

extern Ord$

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)

{

RESTIME=UFTIME=GCDTIME=SQTIME=0;

Line 789 def norm_ch_lag(V,VM,P,P0) {

Line 803 def norm_ch_lag(V,VM,P,P0) {

def cr_gcda(P1,P2)

{

if ( !(V = var(P1)) || !var(P2) )

if ( !P1 )

return P2;

if ( !P2 )

return P1;

if ( !var(P1) || !var(P2) )

return 1;

V = var(P1);

EXT = union_sort(getalgtreep(P1),getalgtreep(P2));

if ( EXT == [] )

return gcd(P1,P2);