DoubleDouble.jl

DoubleDouble.jl copied to clipboard

With current definition of ` y # false ```

GregPlowman

I thought `

dpo

Correctly-rounded operations

1

It would be great to have correctly-rounded basic operations, i.e. be able to fix the rounding mode to `RoundDown` or `RoundUp` and do `a + b` for `a::Double{T}` and `b::Double{T}`.

dpsanders

No sqrt(Double(0.0)) and sign(Double(0.0))

2

Simon, presently sqrt(Double(0.0)) and sign(Double(0.0)) throw an error. Possible remedy for sqrt is to change ``` # Dekker sqrt2 function sqrt{T}(x::Double{T}) if x.hi

ivanslapnicar

compare with Bailey package

5

The most popular implementation of this sort of thing seems to be Bailey's [QD package](http://crd-legacy.lbl.gov/~dhbailey/mpdist/) (which also includes double-double arithmetic). Would be good to do some comparison with that package...

stevengj