diff options
| author | dos-reis <gdr@axiomatics.org> | 2011-03-12 19:17:32 +0000 |
|---|---|---|
| committer | dos-reis <gdr@axiomatics.org> | 2011-03-12 19:17:32 +0000 |
| commit | c93e84048964194e5674e859d6bd8827010f09f6 (patch) | |
| tree | 407ca23d232bf3f974550fce0175c49495ef49c0 /src/algebra/rdeef.spad.pamphlet | |
| parent | 52fddea19454dc2b9bcb54c6edd5a4cd4f5765a9 (diff) | |
| download | open-axiom-c93e84048964194e5674e859d6bd8827010f09f6.tar.gz | |
* src/algebra/: Systematically use negative? when comparing for
greater than 0.
Diffstat (limited to 'src/algebra/rdeef.spad.pamphlet')
| -rw-r--r-- | src/algebra/rdeef.spad.pamphlet | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/src/algebra/rdeef.spad.pamphlet b/src/algebra/rdeef.spad.pamphlet index c7d5e368..b028efda 100644 --- a/src/algebra/rdeef.spad.pamphlet +++ b/src/algebra/rdeef.spad.pamphlet @@ -272,7 +272,7 @@ ElementaryRischDE(R, F): Exports == Implementation where else data := data1::Data f' := nfprime := n * differentiate(f, x) p:P := 1 - for v in data | (m := n * v.coeff) > 0 repeat + for v in data | positive?(m := n * v.coeff) repeat p := p * v.argument ** (m::N) f' := f' - m * differentiate(v.argument::F, x) / (v.argument::F) rec := baseRDE(univariate(f', k), univariate(p::F * g, k)) @@ -286,7 +286,7 @@ ElementaryRischDE(R, F): Exports == Implementation where tab.f := data := makeData(f, x, k) else data := data1::Data p:P := 1 - for v in data | (m := n * v.coeff) > 0 repeat + for v in data | positive?(m := n * v.coeff) repeat p := p * v.argument ** (m::N) f := f - v.coeff * log(v.argument::F) nfp := nfp - m * differentiate(v.argument::F, x) / (v.argument::F) @@ -318,9 +318,9 @@ ElementaryRischDE(R, F): Exports == Implementation where one?(degree(num := univariate(fnum, u))) then cf := (leadingCoefficient num) / fden if (n := retractIfCan(cf)@Union(Z, "failed")) case Z then - if degree(numer logand, k) > 0 then + if positive? degree(numer logand, k) then disasters := concat([n::Z, numer logand], disasters) - if degree(denom logand, k) > 0 then + if positive? degree(denom logand, k) then disasters := concat([-(n::Z), denom logand], disasters) disasters @@ -343,7 +343,7 @@ ElementaryRischDE(R, F): Exports == Implementation where (u := cc exquo aa) case "failed" => "failed" primintfldpoly(u::UP, extint(#1, t'), t') n := degree(cc)::Z - (db := degree(bb)::Z) - if ((da := degree(aa)::Z) = db) and (da > 0) then + if ((da := degree(aa)::Z) = db) and positive? da then lk0 := tower(f0 := - (leadingCoefficient bb) / (leadingCoefficient aa)) lk1 := logdiff(twr, lk0) |