Subject: Re: cygwin-b20 gcc: "A - floor(A)" returns a negative value
From: [email-masked] (Earnie Boyd)

11 Nov 1998 11:48:04 -0800
Message-ID: <[email-masked]>
Reply-To: [email-masked]
Mime-Version: 1.0
Content-Type: text/plain; charset=us-ascii
To: Peter Kabal <[email-masked]>, [email-masked]

Hi Peter,

Although your example indeed shows a potential bug in the optimizer
IMHO your example is a little buggy itself.  I believe your intent
with the printf is to show the values that gave you the result and to
do this you need two variables.  When I modified your example to use
two variables, AV1 and AV2, the result was the same for all values of
optimization.

If this was not your intent then please forgive my claims. 

---Peter Kabal <[email-masked]> wrote:
>
> Switching from cygwin-b19 to cygwin-b20 caused a C-routine to give
> incorrect results.  This was traced to a calculation (in double) of
> "A - float(A)" which returns a negative value.  By definition it
should
> always be positive.  A short test program which demonstrates this
> problem is included.  If the program is further simplified, the
problem
> goes away.  To exercise the bug, the program must be compiled with
> optimization level -O3.
> 
> % gcc -O3 tfloor.c -o tfloor
> % ./tfloor
> AV = 20, AV - floor(AV) = -1.11022e-15
> 
> ----------------
> #include <math.h>
> #include <stdio.h>
> 
> int
> main (argc, argv)
>      int argc;
>      char *argv[];
> {
>   double AV;
> 
>   AV = 0.05;
>   AV = 1.0 / (AV - floor (AV));
>   printf ("AV = %g, AV - floor(AV) = %g\n", AV, AV - floor(AV));
> 
>   return 0;
> };
> ----------------
> Peter Kabal    [email-masked]
> Dept. Electrical & Computer Engineering
> McGill University
> -
> For help on using this list (especially unsubscribing), send a
message to
> "[email-masked]" with one line of text: "help".
> 

==
-                  \\||//
-------------o0O0--Earnie--0O0o--------------
--          [email-masked]          --
-- http://www.freeyellow.com/members5/gw32 --
----------------ooo0O--O0ooo-----------------

PS: Newbie's, you should visit my page.
_________________________________________________________
DO YOU YAHOO!?
Get your free @yahoo.com address at http://mail.yahoo.com

-
For help on using this list (especially unsubscribing), send a message to
"[email-masked]" with one line of text: "help".