login  home  contents  what's new  discussion  bug reports help  links  subscribe  changes  refresh  edit

 Submitted by : (unknown) at: 2007-11-17T22:29:19-08:00 (12 years ago) Name : Axiom Version : default friCAS-20090114 Axiom-20050901 OpenAxiom-20091012 OpenAxiom-20110220 OpenAxiom-Release-141 Category : Axiom Aldor Interface Axiom Compiler Axiom Library Axiom Interpreter Axiom Documentation Axiom User Interface building Axiom from source lisp system MathAction Doyen CD Reduce Axiom on Windows Axiom on Linux Severity : critical serious normal minor wishlist Status : open closed rejected not reproducible fix proposed fixed somewhere duplicate need more info Optional subject :   Optional comment :

Axiom provides functionality to evaluate operators using a given evaluation function. Unfortunately, it is broken.

axiom
g := operator 'g;
Type: BasicOperator
axiom
eval(g(1783), g, i+->i)
 (1)
Type: Expression(Integer)
axiom
eval(g(n), g, i+->i)
 (2)
Type: Expression(Integer)
axiom
eval(a*g(n), g, i+->i)
 (3)
Type: Expression(Integer)
axiom
eval(a*g(1783), g, i+->i)
 (4)
Type: Expression(Integer)

I guess that the problem is in 'smprep$FS':  smprep(lop, lexp, lfunc, p) == ..(v := mainVariable p) case "failed" => p::% ..symbolIfCan(k := v::K) case SY => p::% ..g := (op := operator k) .....(arg := [eval(a,lop,lexp,lfunc) for a in argument k]$List(%))
..q := map(eval(#1::%, lop, lexp, lfunc),
...........univariate(p, k))\$SparseUnivariatePolynomialFunctions2(MP, %)
..(n := position(name op, lop)) < minIndex lop => q g
..a:%  := 0
..f    := eval((lfunc.n) arg, lop, lexp, lfunc)
..e    := lexp.n
..while q ^= 0 repeat
....m  := degree q
....qr := divide(m, e)
....t1 := f ** (qr.quotient)::N
....t2 := g ** (qr.remainder)::N
....a  := a + leadingCoefficient(q) * t1 * t2
....q  := reductum q
..a


It seems that Axiom picks the wrong mainVariable in the broken case, namely a instead of g(1783). I badly need this fixed...

Martin

Status: open => fix proposed

Status: fix proposed => fixed somewhere

no patch available

From wh-sandbox --alfredo, Thu, 28 Aug 2008 18:49:01 -0700 reply