| 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); |