Return to t-misc.c CVS log

Up to [local] / OpenXM_contrib / gmp / mpf / tests

Annotation of OpenXM_contrib/gmp/mpf/tests/t-misc.c, Revision 1.1

1.1     ! maekawa     1: /* Exercise various mpf functions. */
        !             2:
        !             3: /*
        !             4: Copyright (C) 2000 Free Software Foundation, Inc.
        !             5:
        !             6: This file is part of the GNU MP Library.
        !             7:
        !             8: The GNU MP Library is free software; you can redistribute it and/or modify
        !             9: it under the terms of the GNU Lesser General Public License as published by
        !            10: the Free Software Foundation; either version 2.1 of the License, or (at your
        !            11: option) any later version.
        !            12:
        !            13: The GNU MP Library is distributed in the hope that it will be useful, but
        !            14: WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
        !            15: or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
        !            16: License for more details.
        !            17:
        !            18: You should have received a copy of the GNU Lesser General Public License
        !            19: along with the GNU MP Library; see the file COPYING.LIB.  If not, write to
        !            20: the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
        !            21: MA 02111-1307, USA.
        !            22: */
        !            23:
        !            24: #include <stdio.h>
        !            25: #include "gmp.h"
        !            26: #include "gmp-impl.h"
        !            27:
        !            28: #define SGN(x)       ((x) < 0 ? -1 : (x) == 0 ? 0 : 1)
        !            29:
        !            30:
        !            31: void
        !            32: check_mpf_set_si (void)
        !            33: {
        !            34:   static const struct {
        !            35:     long       x;
        !            36:     mp_size_t  want_size;
        !            37:     mp_limb_t  want_limb;
        !            38:   } data[] = {
        !            39:
        !            40:     {  0L,  0 },
        !            41:     {  1L,  1, 1 },
        !            42:     { -1L, -1, 1 },
        !            43:
        !            44:     {  LONG_MAX,  1, LONG_MAX },
        !            45:     { -LONG_MAX, -1, LONG_MAX },
        !            46:
        !            47:     { LONG_HIGHBIT, -1, ULONG_HIGHBIT },
        !            48:   };
        !            49:
        !            50:   mpf_t  x;
        !            51:   int    i;
        !            52:
        !            53:   mpf_init (x);
        !            54:   for (i = 0; i < numberof (data); i++)
        !            55:     {
        !            56:       mpf_init (x);
        !            57:       mpf_set_si (x, data[i].x);
        !            58:       if (x->_mp_size != data[i].want_size
        !            59:           || (x->_mp_size != 0 && x->_mp_d[0] != data[i].want_limb)
        !            60:           || (x->_mp_exp != (data[i].x != 0)))
        !            61:         {
        !            62:           printf ("mpf_set_si wrong on data[%d]\n", i);
        !            63:           abort();
        !            64:         }
        !            65:       mpf_clear (x);
        !            66:
        !            67:       mpf_init_set_si (x, data[i].x);
        !            68:       if (x->_mp_size != data[i].want_size
        !            69:           || (x->_mp_size != 0 && x->_mp_d[0] != data[i].want_limb)
        !            70:           || (x->_mp_exp != (data[i].x != 0)))
        !            71:         {
        !            72:           printf ("mpf_init_set_si wrong on data[%d]\n", i);
        !            73:           abort();
        !            74:         }
        !            75:       mpf_clear (x);
        !            76:     }
        !            77: }
        !            78:
        !            79:
        !            80: void
        !            81: check_mpf_cmp_si (void)
        !            82: {
        !            83:   static const struct {
        !            84:     int         a_base;
        !            85:     const char  *a;
        !            86:     const char  *b;
        !            87:     int         want;
        !            88:   } data[] = {
        !            89:     { 10, "0",  "1", -1 },
        !            90:     { 10, "0",  "0",  0 },
        !            91:     { 10, "0", "-1",  1 },
        !            92:
        !            93:     { 10, "1",  "1", 0 },
        !            94:     { 10, "1",  "0", 1 },
        !            95:     { 10, "1", "-1", 1 },
        !            96:
        !            97:     { 10, "-1",  "1", -1 },
        !            98:     { 10, "-1",  "0", -1 },
        !            99:     { 10, "-1", "-1", 0 },
        !           100:
        !           101:     { 16,         "0", "-0x80000000",  1 },
        !           102:     { 16,  "80000000", "-0x80000000",  1 },
        !           103:     { 16,  "80000001", "-0x80000000",  1 },
        !           104:     { 16, "-80000000", "-0x80000000",  0 },
        !           105:     { 16, "-80000001", "-0x80000000", -1 },
        !           106:
        !           107:     { 16,                 "0", "-0x8000000000000000",  1 },
        !           108:     { 16,  "8000000000000000", "-0x8000000000000000",  1 },
        !           109:     { 16,  "8000000000000001", "-0x8000000000000000",  1 },
        !           110:     { 16, "-8000000000000000", "-0x8000000000000000",  0 },
        !           111:     { 16, "-8000000000000001", "-0x8000000000000000", -1 },
        !           112:   };
        !           113:
        !           114:   mpf_t  a;
        !           115:   mpz_t  bz;
        !           116:   long   b;
        !           117:   int    got;
        !           118:   int    i;
        !           119:
        !           120:   mpf_init (a);
        !           121:   mpz_init (bz);
        !           122:   for (i = 0; i < numberof (data); i++)
        !           123:     {
        !           124:       mpf_set_str (a, data[i].a, data[i].a_base);
        !           125:       mpz_set_str (bz, data[i].b, 0);
        !           126:
        !           127:       if (mpz_fits_slong_p (bz))
        !           128:         {
        !           129:           b = mpz_get_si (bz);
        !           130:           got = mpf_cmp_si (a, b);
        !           131:           if (SGN (got) != data[i].want)
        !           132:             {
        !           133:               printf ("mpf_cmp_si wrong on data[%d]\n", i);
        !           134:               printf ("  a="); mpf_out_str (stdout, 10, 0, a);
        !           135:               printf (" (%s)\n", data[i].a);
        !           136:               printf ("  b=%ld (%s)\n", b, data[i].b);
        !           137:               printf ("  got=%d\n", got);
        !           138:               printf ("  want=%d\n", data[i].want);
        !           139:               abort();
        !           140:             }
        !           141:         }
        !           142:     }
        !           143:
        !           144:   mpf_clear (a);
        !           145:   mpz_clear (bz);
        !           146: }
        !           147:
        !           148:
        !           149: int
        !           150: main (void)
        !           151: {
        !           152:   check_mpf_set_si ();
        !           153:   check_mpf_cmp_si ();
        !           154:   exit (0);
        !           155: }