diff options
| author | dos-reis <gdr@axiomatics.org> | 2011-03-12 17:04:43 +0000 |
|---|---|---|
| committer | dos-reis <gdr@axiomatics.org> | 2011-03-12 17:04:43 +0000 |
| commit | 52fddea19454dc2b9bcb54c6edd5a4cd4f5765a9 (patch) | |
| tree | 4b13ccc6d57c7c1ee615c83615b246c98eae388a /src/algebra/numquad.spad.pamphlet | |
| parent | 42d38bee45a64edfc12641053e58581c20584363 (diff) | |
| download | open-axiom-52fddea19454dc2b9bcb54c6edd5a4cd4f5765a9.tar.gz | |
* src/algebra/: Systematically use negative? when comparing for
less than 0.
Diffstat (limited to 'src/algebra/numquad.spad.pamphlet')
| -rw-r--r-- | src/algebra/numquad.spad.pamphlet | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/src/algebra/numquad.spad.pamphlet b/src/algebra/numquad.spad.pamphlet index 82c82440..c63bf852 100644 --- a/src/algebra/numquad.spad.pamphlet +++ b/src/algebra/numquad.spad.pamphlet @@ -322,10 +322,10 @@ NumericalQuadrature(): Exports == Implementation where if (a = b) then output("romberg: integration limits are equal = ",a::E) return([0.0,0.0,1,true]) - if (epsrel < 0.0) then + if negative? epsrel then output("romberg: eps_r < 0.0 eps_r = ",epsrel::E) return([0.0,0.0,0,false]) - if (epsabs < 0.0) then + if negative? epsabs then output("romberg: eps_a < 0.0 eps_a = ",epsabs::E) return([0.0,0.0,0,false]) for n in 1..nmax repeat @@ -373,10 +373,10 @@ NumericalQuadrature(): Exports == Implementation where if (a = b) then output("simpson: integration limits are equal = ",a::E) return([0.0,0.0,1,true]) - if (epsrel < 0.0) then + if negative? epsrel then output("simpson: eps_r < 0.0 : eps_r = ",epsrel::E) return([0.0,0.0,0,false]) - if (epsabs < 0.0) then + if negative? epsabs then output("simpson: eps_a < 0.0 : eps_a = ",epsabs::E) return([0.0,0.0,0,false]) for n in 1..nmax repeat @@ -415,10 +415,10 @@ NumericalQuadrature(): Exports == Implementation where if (a = b) then output("trapezoidal: integration limits are equal = ",a::E) return([0.0,0.0,1,true]) - if (epsrel < 0.0) then + if negative? epsrel then output("trapezoidal: eps_r < 0.0 : eps_r = ",epsrel::E) return([0.0,0.0,0,false]) - if (epsabs < 0.0) then + if negative? epsabs then output("trapezoidal: eps_a < 0.0 : eps_a = ",epsabs::E) return([0.0,0.0,0,false]) for n in 1..nmax repeat |