Mercurial > hg > octave-kai > gnulib-hg
annotate tests/test-roundl-ieee.c @ 14079:97fc9a21a8fb
maint: update almost all copyright ranges to include 2011
Run the new "make update-copyright" rule.
| author | Jim Meyering <meyering@redhat.com> |
|---|---|
| date | Sat, 01 Jan 2011 20:17:23 +0100 |
| parents | b91bc81df379 |
| children | 8250f2777afc |
| rev | line source |
|---|---|
| 13988 | 1 /* Test of rounding to nearest, breaking ties away from zero. |
|
14079
97fc9a21a8fb
maint: update almost all copyright ranges to include 2011
Jim Meyering <meyering@redhat.com>
parents:
13994
diff
changeset
|
2 Copyright (C) 2010-2011 Free Software Foundation, Inc. |
| 13988 | 3 |
| 4 This program is free software; you can redistribute it and/or modify | |
| 5 it under the terms of the GNU General Public License as published by | |
| 6 the Free Software Foundation; either version 3, or (at your option) | |
| 7 any later version. | |
| 8 | |
| 9 This program is distributed in the hope that it will be useful, | |
| 10 but WITHOUT ANY WARRANTY; without even the implied warranty of | |
| 11 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the | |
| 12 GNU General Public License for more details. | |
| 13 | |
| 14 You should have received a copy of the GNU General Public License | |
| 15 along with this program; if not, write to the Free Software Foundation, | |
| 16 Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. */ | |
| 17 | |
| 18 #include <config.h> | |
| 19 | |
| 20 #include <math.h> | |
| 21 | |
| 22 #include "fpucw.h" | |
| 23 #include "minus-zero.h" | |
| 24 #include "macros.h" | |
| 25 | |
| 26 int | |
| 27 main () | |
| 28 { | |
| 29 DECL_LONG_DOUBLE_ROUNDING | |
| 30 | |
| 31 BEGIN_LONG_DOUBLE_ROUNDING (); | |
| 32 | |
|
13994
b91bc81df379
round: Implement result sign according to IEEE 754.
Bruno Haible <bruno@clisp.org>
parents:
13988
diff
changeset
|
33 /* See IEEE 754, section 6.3: |
|
b91bc81df379
round: Implement result sign according to IEEE 754.
Bruno Haible <bruno@clisp.org>
parents:
13988
diff
changeset
|
34 "the sign of the result of the round floating-point number to |
|
b91bc81df379
round: Implement result sign according to IEEE 754.
Bruno Haible <bruno@clisp.org>
parents:
13988
diff
changeset
|
35 integral value operation is the sign of the operand. These rules |
|
b91bc81df379
round: Implement result sign according to IEEE 754.
Bruno Haible <bruno@clisp.org>
parents:
13988
diff
changeset
|
36 shall apply even when operands or results are zero or infinite." */ |
|
b91bc81df379
round: Implement result sign according to IEEE 754.
Bruno Haible <bruno@clisp.org>
parents:
13988
diff
changeset
|
37 |
| 13988 | 38 /* Zero. */ |
| 39 ASSERT (!signbit (roundl (0.0L))); | |
| 40 ASSERT (!!signbit (roundl (minus_zerol)) == !!signbit (minus_zerol)); | |
|
13994
b91bc81df379
round: Implement result sign according to IEEE 754.
Bruno Haible <bruno@clisp.org>
parents:
13988
diff
changeset
|
41 /* Positive numbers. */ |
|
b91bc81df379
round: Implement result sign according to IEEE 754.
Bruno Haible <bruno@clisp.org>
parents:
13988
diff
changeset
|
42 ASSERT (!signbit (roundl (0.3L))); |
|
b91bc81df379
round: Implement result sign according to IEEE 754.
Bruno Haible <bruno@clisp.org>
parents:
13988
diff
changeset
|
43 ASSERT (!signbit (roundl (0.7L))); |
|
b91bc81df379
round: Implement result sign according to IEEE 754.
Bruno Haible <bruno@clisp.org>
parents:
13988
diff
changeset
|
44 /* Negative numbers. */ |
|
b91bc81df379
round: Implement result sign according to IEEE 754.
Bruno Haible <bruno@clisp.org>
parents:
13988
diff
changeset
|
45 ASSERT (!!signbit (roundl (-0.3L)) == !!signbit (minus_zerol)); |
|
b91bc81df379
round: Implement result sign according to IEEE 754.
Bruno Haible <bruno@clisp.org>
parents:
13988
diff
changeset
|
46 ASSERT (!!signbit (roundl (-0.7L)) == !!signbit (minus_zerol)); |
| 13988 | 47 |
| 48 return 0; | |
| 49 } |