
24 Sep
2005
24 Sep
'05
6:50 p.m.
What should we do about the accuracy of double_ operations? The implementation of plus,minus,times and divide mimics the behaviour of the runtime equivalent, double. This means that the mantissa is trunkated from 61 to 52 bit for every fundamental operation. The result of this, is that complex functions such as sine and exponential will differ from their runtime counterpart, unless a specialization is made for double_. The problem would disappear if we allow calculations with double_ to be more accurate than calculations with double. Is this a problem? Regards Peder