version 1.3, 2000/08/22 05:04:21 |
version 1.4, 2001/03/08 07:49:13 |
|
|
* 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/fff,v 1.2 2000/08/21 08:31:41 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/lib/fff,v 1.3 2000/08/22 05:04:21 noro Exp $ |
*/ |
*/ |
/* |
/* |
fff : Univariate factorizer over a finite field. |
fff : Univariate factorizer over a finite field. |
|
|
F2 = ugcd(F,F1); |
F2 = ugcd(F,F1); |
/* FLAT = H/gcd(H,H') : square free part of H */ |
/* FLAT = H/gcd(H,H') : square free part of H */ |
FLAT = sdiv(F,F2); |
FLAT = sdiv(F,F2); |
|
FLAT /= LCOEF(FLAT); |
I = 0; |
I = 0; |
/* square free factorization of H */ |
/* square free factorization of H */ |
while ( deg(FLAT,V) ) { |
while ( deg(FLAT,V) ) { |
|
|
FLAT1 = simp_ff(1); |
FLAT1 = simp_ff(1); |
else |
else |
FLAT1 = ugcd(F,FLAT); |
FLAT1 = ugcd(F,FLAT); |
|
FLAT1 /= LCOEF(FLAT1); |
G = sdiv(FLAT,FLAT1); |
G = sdiv(FLAT,FLAT1); |
FLAT = FLAT1; |
FLAT = FLAT1; |
L = cons([G,I],L); |
L = cons([G,I],L); |