aboutsummaryrefslogtreecommitdiff
path: root/src/algebra/sf.spad.pamphlet
diff options
context:
space:
mode:
authordos-reis <gdr@axiomatics.org>2011-02-05 03:52:54 +0000
committerdos-reis <gdr@axiomatics.org>2011-02-05 03:52:54 +0000
commit5eebd42bedc639efb11826fcc56520afbb4637ca (patch)
treec692329216ca44718ddb41e58679d5c334187a4f /src/algebra/sf.spad.pamphlet
parentdd16b47e4904b1839b2918cf5b38c352979a6141 (diff)
downloadopen-axiom-5eebd42bedc639efb11826fcc56520afbb4637ca.tar.gz
* interp/g-opt.boot ($VMsideEffectFreeOperators): Remove CGREATERP
and GGREATERP. Add VM-level complex number operators. * interp/lisp-backend.boot: Expand them. * algebra/gaussian.spad.pamphlet (ComplexCategory): Use VM form instead of Lisp level forms. * algebra/pfr.spad.pamphlet (PartialFraction): Likewise. * algebra/si.spad.pamphlet (SingleInteger): Likewise. * algebra/symbol.spad.pamphlet (Symbol): Likewise.
Diffstat (limited to 'src/algebra/sf.spad.pamphlet')
-rw-r--r--src/algebra/sf.spad.pamphlet10
1 files changed, 6 insertions, 4 deletions
diff --git a/src/algebra/sf.spad.pamphlet b/src/algebra/sf.spad.pamphlet
index 745dc254..b78d9bdb 100644
--- a/src/algebra/sf.spad.pamphlet
+++ b/src/algebra/sf.spad.pamphlet
@@ -475,10 +475,12 @@ DoubleFloat(): Join(FloatingPointSystem, DifferentialRing, OpenMath,
manexp(x) ==
zero? x => [0,0]
- s := sign x; x := abs x
- if x > max()$% then return [s*mantissa(max())+1,exponent max()]
- me:Record(man:%,exp:Integer) := MANEXP(x)$Lisp
- two53:= base()**precision()
+ s := sign x
+ x := abs x
+ if x > max()$% then
+ return [s*mantissa(max())+1,exponent max()]
+ me: Record(man:%,exp:Integer) := MANEXP(x)$Lisp
+ two53 := base()**precision()
[s*wholePart(two53 * me.man ),me.exp-precision()]
-- rationalApproximation(y,d,b) ==