From 5eebd42bedc639efb11826fcc56520afbb4637ca Mon Sep 17 00:00:00 2001 From: dos-reis Date: Sat, 5 Feb 2011 03:52:54 +0000 Subject: * 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. --- src/algebra/sf.spad.pamphlet | 10 ++++++---- 1 file changed, 6 insertions(+), 4 deletions(-) (limited to 'src/algebra/sf.spad.pamphlet') 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) == -- cgit v1.2.3