GHC: Ticket #3676: realToFrac doesn't sanely convert between floating types
http://ghc.haskell.org/trac/ghc/ticket/3676
<p>
As far as I can tell, the only way to convert between floating types in Haskell is to use realToFrac. Unfortunately, this function does some insane mangling of values when converting. Some examples:
</p>
<pre class="wiki">realToFrac (0/0 :: Double) :: Double
--> -Infinity
realToFrac (-1/0 :: Float) :: Double
--> -3.402823669209385e38
realToFrac (-0 :: Double) :: CDouble
--> 0.0
</pre><p>
The last item illustrates an important point: it is impossible to convert a value from Double to CDouble without potentially changing it. This makes it difficult or impossible to use the FFI to call any functions with floating point parameters/return values.
</p>
<p>
Using a combination of unsafeCoerce (to shoehorn Double values in/out of CDoubles) and some functions written in C (to perform float<=>double), I was able to work around these problems, but that hardly seems like a nice solution.
</p>
en-usGHChttp://ghc.haskell.org/trac/ghc/chrome/site/ghc_logo.png
http://ghc.haskell.org/trac/ghc/ticket/3676
Trac 1.0.1iglooFri, 20 Nov 2009 22:01:42 GMTdescription changed
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:1
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:1
<ul>
<li><strong>description</strong>
modified (<a href="/trac/ghc/ticket/3676?action=diff&version=1">diff</a>)
</li>
</ul>
TicketiglooFri, 20 Nov 2009 22:09:57 GMTmilestone set
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:2
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:2
<ul>
<li><strong>milestone</strong>
set to <em>6.12.2</em>
</li>
</ul>
<p>
Thanks for the report.
</p>
TicketiglooSun, 22 Nov 2009 00:00:24 GMT
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:3
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:3
<p>
See also <a class="new ticket" href="http://ghc.haskell.org/trac/ghc/ticket/3070" title="bug: floor(0/0) should not be defined (new)">#3070</a>
</p>
TicketdraconxTue, 09 Feb 2010 17:24:54 GMTversion changed
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:4
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:4
<ul>
<li><strong>version</strong>
changed from <em>6.10.4</em> to <em>6.12.1</em>
</li>
</ul>
<p>
Still present in 6.12.1.
</p>
<p>
I just discovered that the issue is worse than originally described: the behaviour of realToFrac *changes* if it gets inlined. Consider the following:
</p>
<pre class="wiki">module Main where
{-# NOINLINE holycow #-}
holycow :: (Real a, Fractional b) => a -> b
holycow = realToFrac
main = do
print (realToFrac (0/0 :: Double) :: Double)
print (holycow (0/0 :: Double) :: Double)
</pre><p>
If you compile without optimisation, both lines print -Infinity as originally described. However, if you enable optimisation, the first line prints NaN (yay!) but the second line still prints -Infinity. GHC seems to be optimising based on the incorrect assumption that realToFrac :: Double -> Double is the identity function.
</p>
TicketsimonpjWed, 10 Feb 2010 11:51:50 GMT
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:5
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:5
<p>
Yes, that's because the library <tt>GHC.Float</tt> contains the (clearly incorrect) RULE
</p>
<pre class="wiki">"realToFrac/Double->Double" realToFrac = id :: Double -> Double
</pre><p>
This ticket really really needs someone who knows and cares about numbers to figure out what the Right Thing is. It's a library issue involving only Haskell source code. Would someone like to help?
</p>
<p>
Thanks
</p>
<p>
Simon
</p>
TicketisaacdupreeWed, 10 Feb 2010 18:32:24 GMT
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:6
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:6
<p>
my gut feeling has been there should be a function or functions specifically to convert between different floating types. I think this is something that needs to be invented and proposed on libraries@ . (I'm not volunteering to do it this week because I'm going to a conference, but I'm willing to get the conversation started next week)
</p>
<p>
-Isaac Dupree
</p>
TicketsimonmarThu, 11 Feb 2010 10:02:38 GMTstatus changed; owner set
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:7
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:7
<ul>
<li><strong>owner</strong>
set to <em>simonmar</em>
</li>
<li><strong>status</strong>
changed from <em>new</em> to <em>assigned</em>
</li>
</ul>
TicketsimonmarThu, 11 Feb 2010 10:12:02 GMT
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:8
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:8
<p>
I think it's pretty straightforward. <tt>toRational</tt> should be able to convert NaN, Infinity and -Infinity correctly to Rational, and there is already code in fromRational to convert them back to the floating point types correctly (I imagine it's there in order to make things like <tt>read "NaN"</tt> work).
</p>
<p>
The Haskell report doesn't say anything about the implementation of <tt>toRational</tt>, so I think this is correct. I'm validating a patch now.
</p>
TicketdraconxThu, 11 Feb 2010 13:27:20 GMT
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:9
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:9
<p>
Don't forget negative zero.
</p>
TicketsimonmarThu, 11 Feb 2010 14:10:08 GMT
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:10
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:10
<p>
I've fixed this, but I'm not sure the fix is an improvement.
</p>
<pre class="wiki">Thu Feb 11 02:19:55 PST 2010 Simon Marlow <marlowsd@gmail.com>
* Handle NaN, -Infinity and Infinity in the toRational for Float/Double (#3676)
M ./GHC/Float.lhs -2 +9
M ./GHC/Real.lhs -1 +2
</pre><p>
and in ghc:
</p>
<pre class="wiki">Thu Feb 11 05:15:43 PST 2010 Simon Marlow <marlowsd@gmail.com>
* don't constant fold division that would result in negative zero (#3676)
</pre><p>
It has some slightly odd consequences. Someone ought to take a broader look at this in the context of refining the Haskell standard. e.g.
</p>
<pre class="wiki">Prelude> toRational (-0 :: Double) == - toRational (0 :: Double)
False
</pre><p>
whereas previously this was True.
</p>
<p>
The underlying problem is that the Rational type doesn't really include these strange values; it has been hacked in a few places so that certain unnormalised Rational values are used to represent strange floating point values, e.g. <tt>0 :% 0</tt> represents <tt>NaN</tt> and <tt>fromRational</tt> knows about it, but there's no official way to generate one of these (<tt>0/0 :: Rational</tt> gives an exception), they are only used "internally" by the libraries to make <tt>read "NaN" :: Double</tt> work.
</p>
<p>
The change I made is to make <tt>toRational (0/0 :: Double)</tt> generate <tt>0 :% 0</tt>, and similarly for the other strange values, but these had odd consequences, because these special <tt>Rational</tt> values are not normalised. However, this does fix a real bug (the subject of this ticket).
</p>
<p>
What do people think?
</p>
TicketdraconxThu, 11 Feb 2010 14:40:07 GMT
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:11
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:11
<p>
While it presumably goes against the report, I personally feel that floating point types do not belong as members of the Real class *at all*. The same is true for a number of other classes that they belong to.
</p>
<p>
In the altfloat package, where I've been experimenting with solutions to this and other problems, I used a class
</p>
<pre class="wiki">class FloatConvert a b where
toFloating :: a -> b
</pre><p>
which provides the necessary conversions: both between floating point types, and from real types to floating point types. It should be noted that the functions <tt>double2Float#</tt> and <tt>float2Double#</tt> in GHC.Exts perform correct conversion.
</p>
TicketmaederThu, 11 Feb 2010 16:42:52 GMT
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:12
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:12
<p>
I suppose arithmetic on these non-normalised rationals does also not behave like on doubles (IEEE).
I would expect toRational to fail on NaN and Infinity and to convert -0 to 0.
</p>
<p>
But rather than introducing a multi-parameter type class I would just add non-overloaded plain conversion functions between Float, Double, CDouble, etc.
</p>
TicketdraconxThu, 11 Feb 2010 18:22:50 GMT
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:13
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:13
<p>
The problem with non-overloaded functions is that there are a lot of them: 12 with the four floating types that currently exist, and 30 if long double is ever added (two types: LDouble & CLDouble). This problem could be mitigated by only providing the trivial conversions when FFI types are involved: that is, provide Double <=> CDouble, but don't bother with Float <=> CDouble.
</p>
<p>
On the other hand, we could use a single parameter version of altfloat's technique as follows:
</p>
<pre class="wiki">class FloatConvert a where
toDouble :: a -> Double
</pre><p>
since Double can represent all values of every other floating type (at least until long double exists). Then we could have a <tt>fromDouble</tt> member in the <a class="missing wiki">RealFloat?</a> class, and define <tt>toFloating = fromDouble . toDouble</tt>
</p>
TicketsimonmarTue, 23 Feb 2010 11:08:25 GMTmilestone changed
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:14
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:14
<ul>
<li><strong>milestone</strong>
changed from <em>6.12.2</em> to <em>_|_</em>
</li>
</ul>
<p>
I rolled back the patch to libraries/base:
</p>
<pre class="wiki">Tue Feb 23 10:16:03 GMT 2010 Simon Marlow <marlowsd@gmail.com>
* UNDO: Handle NaN, -Infinity and Infinity in the toRational for Float/Double (#3676)
</pre><p>
in retrospect it wasn't a good idea to shoehorn these strange values into Rational.
</p>
<p>
The position now is that <tt>toRational (1/0 :: Double)</tt> is undefined. The Haskell standard does not define the result (it should), and GHC gives unpredictable results. Hence, <tt>realToFrac</tt> cannot convert accurately between floating-point types.
</p>
<p>
Someone should:
</p>
<ul><li>propose a proer API for conversion between floating-point types, e.g.
<tt>FloatConvert</tt> above, but make a full proposal.
</li></ul><ul><li>decide whether <tt>toRational (1/0::Double)</tt> should be undefined or an
exception. If it is an exception, then we cannot optimise <tt>realToFrac</tt>
to a direct conversion, e.g. <tt>floatToDouble#</tt>.
</li></ul>
TicketdraconxTue, 23 Feb 2010 14:12:20 GMT
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:15
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:15
<p>
I agree that the "funny rationals" are a bit too strange a notion, especially when it comes to defining how arithmetic operations should behave.
</p>
<p>
After thinking about it, I'm fond of the toDouble/fromDouble approach as it seems consistent with the themes elsewhere in the library, except that it makes it harder to later add floating types that are wider than Double.
</p>
<p>
I can post this idea to the haskell-prime mailing list for discussion, if that's the right one.
</p>
TicketdraconxThu, 04 Mar 2010 18:23:36 GMT
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:16
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:16
<p>
OK, I've posted this to the haskell-prime mailing list.
</p>
<p>
<a class="ext-link" href="http://thread.gmane.org/gmane.comp.lang.haskell.prime/3146"><span class="icon"></span>http://thread.gmane.org/gmane.comp.lang.haskell.prime/3146</a>
</p>
TicketsimonmarWed, 11 Aug 2010 13:25:24 GMTowner deleted
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:17
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:17
<ul>
<li><strong>owner</strong>
<em>simonmar</em> deleted
</li>
</ul>
Ticketdaniel.is.fischerFri, 29 Jul 2011 23:39:37 GMTcc set
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:18
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:18
<ul>
<li><strong>cc</strong>
<em>daniel.is.fischer@…</em> added
</li>
</ul>
<p>
Summarising:
</p>
<p>
There is no good solution to the problem while the conversions go (potentially, without rewrite rules) via <tt>Rational</tt>. Having 'illegal' <tt>Rational</tt> values like <tt>1 :% 0</tt> generated by parsing <tt>Infinity</tt> and so on is not a big problem because those are immediately passed to <tt>fromRational</tt>, so their treatment is confined to one place and they can be handled meaningfully there. But if they could be generated by <tt>toRational</tt>, <em>every</em> operation on <tt>Rational</tt>s would have to check for these values (and it's not always clear how to handle them).
</p>
<p>
The options I'm aware of are
</p>
<ol><li>non-overloaded conversion functions between each pair of types
</li><li>a two-parameter conversion class
</li><li>adding <tt>toMaxFloat</tt> and <tt>fromMaxFloat</tt> functions to <tt>RealFloat</tt>
</li><li>adding a dedicated pseudo-numeric type for conversions
</li></ol><p>
All of them have downsides. 1. and 2. would require a number of functions resp. instances quadratic in the number of involved types - yuck. 1. would even rule out polymorphic code. 3. has the problem of what <tt>MaxFloat</tt> should be and any choice would cause breakage when a larger <tt>RealFloat</tt> type is introduced. 4: A type with the sole purpose of facilitating numeric conversions in the standard libraries is quite wartish (and new types could introduce new unrepresentable values).
</p>
<p>
I think 4. is the least undesirable option, but still it would be a lot of work to implement (it involves a change to the language report, among other things) for dubious benefit.
</p>
<p>
Unless renewed concern is expressed, I think we should mark it as low priority.
</p>
TicketcarterSat, 20 Jul 2013 05:57:54 GMTdifficulty set
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:19
http://ghc.haskell.org/trac/ghc/ticket/3676#comment:19
<ul>
<li><strong>difficulty</strong>
set to <em>Unknown</em>
</li>
</ul>
<p>
Hrm, so we need a good conversion story / semantics for mapping between different floating point types? This is near and dear to some of my current work, I'll look into coming up with a good story, and if successful, propose it to the libraries list.
</p>
Ticket