Return to t-scan.c CVS log

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

Annotation of OpenXM_contrib/gmp/tests/mpz/t-scan.c, Revision 1.1

1.1     ! ohara       1: /* Tests of mpz_scan0 and mpz_scan1.
        !             2:
        !             3: Copyright 2000, 2001, 2002 Free Software Foundation, Inc.
        !             4:
        !             5: This file is part of the GNU MP Library.
        !             6:
        !             7: The GNU MP Library is free software; you can redistribute it and/or modify
        !             8: it under the terms of the GNU Lesser General Public License as published by
        !             9: the Free Software Foundation; either version 2.1 of the License, or (at your
        !            10: option) any later version.
        !            11:
        !            12: The GNU MP Library is distributed in the hope that it will be useful, but
        !            13: WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
        !            14: or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
        !            15: License for more details.
        !            16:
        !            17: You should have received a copy of the GNU Lesser General Public License
        !            18: along with the GNU MP Library; see the file COPYING.LIB.  If not, write to
        !            19: the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
        !            20: MA 02111-1307, USA. */
        !            21:
        !            22: #include <stdio.h>
        !            23: #include <stdlib.h>
        !            24: #include "gmp.h"
        !            25: #include "gmp-impl.h"
        !            26: #include "tests.h"
        !            27:
        !            28:
        !            29: unsigned long
        !            30: refmpz_scan (mpz_srcptr z, unsigned long i, int sought)
        !            31: {
        !            32:   unsigned long  z_bits = (unsigned long) ABSIZ(z) * GMP_NUMB_BITS;
        !            33:
        !            34:   do
        !            35:     {
        !            36:       if (mpz_tstbit (z, i) == sought)
        !            37:         return i;
        !            38:       i++;
        !            39:     }
        !            40:   while (i <= z_bits);
        !            41:
        !            42:   return ULONG_MAX;
        !            43: }
        !            44:
        !            45: unsigned long
        !            46: refmpz_scan0 (mpz_srcptr z, unsigned long starting_bit)
        !            47: {
        !            48:   return refmpz_scan (z, starting_bit, 0);
        !            49: }
        !            50:
        !            51: unsigned long
        !            52: refmpz_scan1 (mpz_srcptr z, unsigned long starting_bit)
        !            53: {
        !            54:   return refmpz_scan (z, starting_bit, 1);
        !            55: }
        !            56:
        !            57:
        !            58: void
        !            59: check_ref (void)
        !            60: {
        !            61:   mpz_t          z;
        !            62:   int            test, size, neg, sought;
        !            63:   unsigned long  i, got, want;
        !            64:
        !            65:   mpz_init (z);
        !            66:   for (test = 0; test < 10; test++)
        !            67:     {
        !            68:       for (size = 0; size < 5; size++)
        !            69:         {
        !            70:           mpz_random2 (z, size);
        !            71:
        !            72:           for (neg = 0; neg <= 1; neg++)
        !            73:             {
        !            74:               if (neg)
        !            75:                 mpz_neg (z, z);
        !            76:
        !            77:               for (i = 0; i < size * GMP_NUMB_BITS + 8; i++)
        !            78:                 {
        !            79:                   for (sought = 0; sought <= 1; sought++)
        !            80:                     {
        !            81:                       if (sought == 0)
        !            82:                         {
        !            83:                           got = mpz_scan0 (z, i);
        !            84:                           want = refmpz_scan0 (z, i);
        !            85:                         }
        !            86:                       else
        !            87:                         {
        !            88:                           got = mpz_scan1 (z, i);
        !            89:                           want = refmpz_scan1 (z, i);
        !            90:                         }
        !            91:
        !            92:                       if (got != want)
        !            93:                         {
        !            94:                           printf ("wrong at test=%d, size=%d, neg=%d, i=%lu, sought=%d\n",
        !            95:                                   test, size, neg, i, sought);
        !            96:                           printf ("   z 0x");
        !            97:                           mpz_out_str (stdout, -16, z);
        !            98:                           printf ("\n");
        !            99:                           printf ("   got=%lu, want=%lu\n", got, want);
        !           100:                           exit (1);
        !           101:                         }
        !           102:                     }
        !           103:                 }
        !           104:             }
        !           105:         }
        !           106:     }
        !           107:   mpz_clear (z);
        !           108: }
        !           109:
        !           110:
        !           111: int
        !           112: main (int argc, char *argv[])
        !           113: {
        !           114:   tests_start ();
        !           115:
        !           116:   check_ref ();
        !           117:
        !           118:   tests_end ();
        !           119:   exit (0);
        !           120: }