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/lingrob.spad.pamphlet | |
| parent | 039d244dbef31e4ca6d34d2cf42b58b41a572884 (diff) | |
| download | open-axiom-f51aee83708673ef9941174951bec9aee80cb03c.tar.gz | |
Avoid modulus bias in 'random()' usage.
Diffstat (limited to 'src/algebra/lingrob.spad.pamphlet')
| -rw-r--r-- | src/algebra/lingrob.spad.pamphlet | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/algebra/lingrob.spad.pamphlet b/src/algebra/lingrob.spad.pamphlet index a3018052..25ab9976 100644 --- a/src/algebra/lingrob.spad.pamphlet +++ b/src/algebra/lingrob.spad.pamphlet @@ -265,7 +265,7 @@ LinGroebnerPackage(lv,F) : C == T rval:List Z :=[] for ii in 1..(#lvar-1) repeat c:Z:=0 - while c=0 repeat c:=random()$Z rem 11 + while c=0 repeat c:=random(11)$Z rval:=concat(c,rval) nval:DPoly := (last.lvar)::DPoly - (+/[r*(vv)::DPoly for r in rval for vv in lvar]) @@ -311,7 +311,7 @@ LinGroebnerPackage(lv,F) : C == T xn:=lvar.last val := xn::DPoly nvar1:NNI:=(#lvar-1):NNI - ll: List Z :=[random()$Z rem 11 for i in 1..nvar1] + ll: List Z :=[random(11)$Z for i in 1..nvar1] val:=val+ +/[ll.i*(lvar.i)::DPoly for i in 1..nvar1] LL:=[elt(univariate(f,xn),val) for f in L] LL:= groebner(LL) |