version 1.3, 2000/08/22 05:04:21 |
version 1.5, 2001/05/09 01:41:42 |
|
|
* 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.4 2001/03/08 07:49:13 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); |
Line 467 def lnf_insert(P,L,V) |
|
Line 469 def lnf_insert(P,L,V) |
|
def c_z_ff(F,E) |
def c_z_ff(F,E) |
{ |
{ |
Type = field_type_ff(); |
Type = field_type_ff(); |
if ( Type == 1 || Type == 3 ) |
if ( Type == 1 || Type == 3 || Type == 4 ) |
return c_z_lm(F,E); |
return c_z_lm(F,E); |
else |
else |
return c_z_gf2n(F,E); |
return c_z_gf2n(F,E); |