On Tue, Sep 18, 2012 at 7:23 AM, Tijl Coosemans <t...@coosemans.org> wrote:

> On 15-09-2012 17:39, Mehmet Erol Sanliturk wrote:
> > On Sat, Sep 15, 2012 at 7:30 AM, Tijl Coosemans <t...@coosemans.org>
> wrote:
> >> On 15-09-2012 16:09, Roman Divacky wrote:
> >>> Is this correct?
> >>>
> >>> lev ~$ ./cos 1.23456789e20
> >>> 6.031937e-01
> >>> -9.629173e-02
> >>> 2.814722e-01
> >>
> >> Yes, that's what the libm call returns.
> >
> > Linux z 3.5.3-1.fc17.x86_64 #1 SMP Wed Aug 29 18:46:34 UTC 2012 x86_64
> > x86_64 x86_64 GNU/Linux
> >
> > clang version 3.0 (tags/RELEASE_30/final)
> > Target: x86_64-redhat-linux-gnu
> > Thread model: posix
> >
> >
> > Output of the initial program is the following :
> >
> > #include <math.h>
> > #include <stdio.h>
> > #include <stdlib.h>
> >
> > int
> > main( int argc, char **argv ) {
> >         double d = strtod( argv[ 1 ], NULL );
> >
> >         printf( " cos : %e\n", ( double ) cos( d ));
> >         printf( "cosf : %e\n", ( double ) cosf( d ));
> >         printf( "cosl : %e\n", ( double ) cosl( d ));
> >         return( 0 );
> > }
> >
> >
> > cos : 2.814722e-01
> > cosf : -9.629173e-02
> > cosl : 7.738403e-01
>
> This is probably because SSE instructions are used on amd64.
>
> > Output of the following program is different :
>
> The reason is that...
>
> > #include <math.h>
> > #include <stdio.h>
> > #include <stdlib.h>
> >
> > int
> > main( int argc, char **argv ) {
> >         double d ;
> >         double two_pi ;
> >         double f ;
> >         double v ;
> >
> >     two_pi = 2 * 3.14159265358979323846 ;
> >         d = strtod( argv[ 1 ], NULL );
> >
> >     f = floor ( d / two_pi ) ;
> >         v = d - f * two_pi ;
>
> ...this is a poor way to compute a remainder. Try to use fmod() or
> remainder() instead.
>


My C knowledge is NOT very well . Thanks .



>
> >         printf( "  given : %e\n", ( double ) d );
> >         printf( "  multiplier : %e\n", ( double ) f );
> >         printf( "reduced : %e\n", ( double ) v );
> >
> >
> >         printf( " cos ( %e ) : %e\n", d , ( double ) cos( d ));
> >         printf( "cosf ( %e ) : %e\n", d , ( double ) cosf( d ));
> >         printf( "cosl ( %e ) : %e\n", d , ( double ) cosl( d ));
> >
> >
> >         printf( " cos ( %e ) : %e\n", v , ( double ) cos( v ));
> >         printf( "cosf ( %e ) : %e\n", v , ( double ) cosf( v ));
> >         printf( "cosl ( %e ) : %e\n", v , ( double ) cosl( v ));
> >
> >
> >     return( 0 );
> > }
> >
> >
> >   given : 1.234568e+20
> >   multiplier : 1.964876e+19
> > reduced : 1.638400e+04
> >
> >
> >  cos ( 1.234568e+20 ) : 2.814722e-01
> > cosf ( 1.234568e+20 ) : -9.629173e-02
> > cosl ( 1.234568e+20 ) : 7.738403e-01
> >
> >  cos ( 1.638400e+04 ) : -8.285342e-01
> > cosf ( 1.638400e+04 ) : -8.285342e-01
> > cosl ( 1.638400e+04 ) : -8.285342e-01
>
>
My intention was to check whether there is a difference between Clang
compiled programs in different operating systems .


The GCC output is as follows :


Linux z 3.5.3-1.fc17.x86_64 #1 SMP Wed Aug 29 18:46:34 UTC 2012 x86_64
x86_64 x86_64 GNU/Linux

cc (GCC) 4.7.0 20120507 (Red Hat 4.7.0-5)
Copyright (C) 2012 Free Software Foundation, Inc.
This is free software; see the source for copying conditions.  There is NO
warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.

     given : 1.234568e+20
 cos ( 1.234568e+20 ) : 2.814722e-01
cosf ( 1.234568e+20 ) : -9.629173e-02
cosl ( 1.234568e+20 ) : 7.738403e-01

multiplier : 1.964876e+19
   reduced : 1.638400e+04
 cos ( 1.638400e+04 ) : -8.285342e-01
cosf ( 1.638400e+04 ) : -8.285342e-01
cosl ( 1.638400e+04 ) : -8.285342e-01

multiplier : 2.607000e+03
   reduced : 3.735904e+00
 cos ( 3.735904e+00 ) : -8.285342e-01
cosf ( 3.735904e+00 ) : -8.285342e-01
cosl ( 3.735904e+00 ) : -8.285342e-01


This shows that GCC is NOT better than Clang .

Thank you very much .

Mehmet Erol Sanliturk
_______________________________________________
freebsd-current@freebsd.org mailing list
http://lists.freebsd.org/mailman/listinfo/freebsd-current
To unsubscribe, send any mail to "freebsd-current-unsubscr...@freebsd.org"

Reply via email to