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/rderf.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/rderf.spad.pamphlet')
-rw-r--r-- | src/algebra/rderf.spad.pamphlet | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/algebra/rderf.spad.pamphlet b/src/algebra/rderf.spad.pamphlet index 1cd6563d..dbf1eddf 100644 --- a/src/algebra/rderf.spad.pamphlet +++ b/src/algebra/rderf.spad.pamphlet @@ -73,7 +73,7 @@ TranscendentalRischDE(F, UP): Exports == Implementation where (~b0?) and (baseCase? or degree(u.b) > max(0, n)) => answ := SPDEnocancel1(u.b, u.c, u.m, derivation) [[u.alpha * answ.ans + u.beta, answ.nosol]] - (n > 0) and (b0? or degree(u.b) < n) => + positive? n and (b0? or degree(u.b) < n) => uansw := SPDEnocancel2(u.b,u.c,u.m,n,leadingCoefficient dt,derivation) uansw case ans=> [[u.alpha * uansw.ans.ans + u.beta, uansw.ans.nosol]] [[uansw.eq.b, uansw.eq.c, uansw.eq.m, @@ -131,7 +131,7 @@ TranscendentalRischDE(F, UP): Exports == Implementation where q:UP := 0 while cc ~= 0 repeat negative? d or negative?(n := (degree cc)::Z - dtm1) or n > d => return [[q, true]] - if n > 0 then + if positive? n then r := monomial((leadingCoefficient cc) / (n * lt), n::N) cc := cc - bb * r - derivation r d := n - 1 |