version 1.28, 2002/11/22 07:32:10 |
version 1.29, 2002/11/22 08:44:57 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/engine/Hgfs.c,v 1.27 2002/11/01 05:43:35 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/engine/Hgfs.c,v 1.28 2002/11/22 07:32:10 noro Exp $ */ |
|
|
#include "ca.h" |
#include "ca.h" |
#include "inline.h" |
#include "inline.h" |
Line 79 void gensqfrsfum(UM p,struct oDUM *dc) |
|
Line 79 void gensqfrsfum(UM p,struct oDUM *dc) |
|
{ |
{ |
int n,i,j,d,mod; |
int n,i,j,d,mod; |
UM t,s,g,f,f1,b; |
UM t,s,g,f,f1,b; |
|
GFS u,v; |
|
|
if ( (n = DEG(p)) == 1 ) { |
if ( (n = DEG(p)) == 1 ) { |
dc[0].f = UMALLOC(DEG(p)); cpyum(p,dc[0].f); dc[0].n = 1; |
dc[0].f = UMALLOC(DEG(p)); cpyum(p,dc[0].f); dc[0].n = 1; |
Line 117 void gensqfrsfum(UM p,struct oDUM *dc) |
|
Line 118 void gensqfrsfum(UM p,struct oDUM *dc) |
|
break; |
break; |
else { |
else { |
DEG(s) = DEG(t)/mod; |
DEG(s) = DEG(t)/mod; |
for ( j = 0; j <= DEG(t); j++ ) |
for ( j = 0; j <= DEG(t); j++ ) { |
COEF(s)[j] = COEF(t)[j*mod]; |
iftogfs(COEF(t)[j*mod],&u); |
|
pthrootgfs(u,&v); |
|
COEF(s)[j] = v?FTOIF(CONT(v)):0; |
|
} |
cpyum(s,b); d *= mod; |
cpyum(s,b); d *= mod; |
} |
} |
} |
} |