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/odealg.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/odealg.spad.pamphlet')
| -rw-r--r-- | src/algebra/odealg.spad.pamphlet | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/src/algebra/odealg.spad.pamphlet b/src/algebra/odealg.spad.pamphlet index 653dece6..1a3d355d 100644 --- a/src/algebra/odealg.spad.pamphlet +++ b/src/algebra/odealg.spad.pamphlet @@ -152,7 +152,7 @@ SystemODESolver(F, LO): Exports == Implementation where r := r - 1 c := c - 1 firstnonzero(m, r) ~= c => error "backsolve: undetermined system" - degree(eq := m(r, c)) > 0 => error "backsolve: pivot of order > 0" + positive? degree(eq := m(r, c)) => error "backsolve: pivot of order > 0" a := leadingCoefficient(eq)::F if part? then part(r + offset) := (v(r + offset) - applyLodo(m, r, part, n)) / a |