diff options
| author | Gabriel Dos Reis <gdr@axiomatics.org> | 2015-12-31 14:56:12 -0800 |
|---|---|---|
| committer | Gabriel Dos Reis <gdr@axiomatics.org> | 2015-12-31 14:56:12 -0800 |
| commit | f51aee83708673ef9941174951bec9aee80cb03c (patch) | |
| tree | 4d8cfccd653dfa09819234bc61a8775091e4f1e6 /src/algebra/idecomp.spad.pamphlet | |
| parent | 039d244dbef31e4ca6d34d2cf42b58b41a572884 (diff) | |
| download | open-axiom-f51aee83708673ef9941174951bec9aee80cb03c.tar.gz | |
Avoid modulus bias in 'random()' usage.
Diffstat (limited to 'src/algebra/idecomp.spad.pamphlet')
| -rw-r--r-- | src/algebra/idecomp.spad.pamphlet | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/algebra/idecomp.spad.pamphlet b/src/algebra/idecomp.spad.pamphlet index d5c5d4e4..51ff3a62 100644 --- a/src/algebra/idecomp.spad.pamphlet +++ b/src/algebra/idecomp.spad.pamphlet @@ -146,7 +146,7 @@ IdealDecompositionPackage(vl,nv) : C == T -- take away nv, now doesn't pv:DPoly:=0 pw:DPoly:=0 while not one? degree(lf,y) repeat - val:=random()$Z rem 23 + val:=random(23)$Z pv:=px+val*py pw:=px-val*py Id:=groebner([(univariate(h,x)).pv for h in Id]) @@ -316,7 +316,7 @@ IdealDecompositionPackage(vl,nv) : C == T -- take away nv, now doesn't ---- put the ideal in general position ---- genPosLastVar(J:FIdeal,truelist:List OV):GenPos == x := last truelist ;lv1:List OV :=remove(x,truelist) - ranvals:List(Z):=[(random()$Z rem 23) for vv in lv1] + ranvals:List(Z):=[random(23)$Z for vv in lv1] val:=+/[rv*(vv::DPoly) for vv in lv1 for rv in ranvals] val:=val+(x::DPoly) [ranvals,groebnerIdeal(groebner([(univariate(p,x)).val |