/* Test mp*_class functions.

Copyright 2002, 2003 Free Software Foundation, Inc.

This file is part of the GNU MP Library.

The GNU MP Library is free software; you can redistribute it and/or modify
it under the terms of the GNU Lesser General Public License as published by
the Free Software Foundation; either version 2.1 of the License, or (at your
option) any later version.

The GNU MP Library is distributed in the hope that it will be useful, but
WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
License for more details.

You should have received a copy of the GNU Lesser General Public License
along with the GNU MP Library; see the file COPYING.LIB.  If not, write to
the Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,
MA 02110-1301, USA. */


/* Note that we don't use <climits> for LONG_MIN, but instead our own
   definitions in gmp-impl.h.  In g++ 2.95.4 (debian 3.0) under
   -mcpu=ultrasparc, limits.h sees __sparc_v9__ defined and assumes that
   means long is 64-bit long, but it's only 32-bits, causing fatal compile
   errors.  */

#include "config.h"

#include <string>

#include "gmp.h"
#ifdef WANT_MPFR
#include "mpfr.h"
#endif
#include "gmpxx.h"
#include "gmp-impl.h"
#include "tests.h"

using namespace std;


void
check_mpz (void)
{
  // mpz_class::fits_sint_p
  {
    int        fits;
    mpz_class  z;
    z = INT_MIN; fits = z.fits_sint_p(); ASSERT_ALWAYS (fits);
    z--;         fits = z.fits_sint_p(); ASSERT_ALWAYS (! fits);
    z = INT_MAX; fits = z.fits_sint_p(); ASSERT_ALWAYS (fits);
    z++;         fits = z.fits_sint_p(); ASSERT_ALWAYS (! fits);
  }

  // mpz_class::fits_uint_p
  {
    int        fits;
    mpz_class  z;
    z = 0;        fits = z.fits_uint_p(); ASSERT_ALWAYS (fits);
    z--;          fits = z.fits_uint_p(); ASSERT_ALWAYS (! fits);
    z = UINT_MAX; fits = z.fits_uint_p(); ASSERT_ALWAYS (fits);
    z++;          fits = z.fits_uint_p(); ASSERT_ALWAYS (! fits);
  }

  // mpz_class::fits_slong_p
  {
    int        fits;
    mpz_class  z;
    z = LONG_MIN; fits = z.fits_slong_p(); ASSERT_ALWAYS (fits);
    z--;          fits = z.fits_slong_p(); ASSERT_ALWAYS (! fits);
    z = LONG_MAX; fits = z.fits_slong_p(); ASSERT_ALWAYS (fits);
    z++;          fits = z.fits_slong_p(); ASSERT_ALWAYS (! fits);
  }

  // mpz_class::fits_ulong_p
  {
    int        fits;
    mpz_class  z;
    z = 0;         fits = z.fits_ulong_p(); ASSERT_ALWAYS (fits);
    z--;           fits = z.fits_ulong_p(); ASSERT_ALWAYS (! fits);
    z = ULONG_MAX; fits = z.fits_ulong_p(); ASSERT_ALWAYS (fits);
    z++;           fits = z.fits_ulong_p(); ASSERT_ALWAYS (! fits);
  }

  // mpz_class::fits_sshort_p
  {
    int        fits;
    mpz_class  z;
    z = SHRT_MIN; fits = z.fits_sshort_p(); ASSERT_ALWAYS (fits);
    z--;          fits = z.fits_sshort_p(); ASSERT_ALWAYS (! fits);
    z = SHRT_MAX; fits = z.fits_sshort_p(); ASSERT_ALWAYS (fits);
    z++;          fits = z.fits_sshort_p(); ASSERT_ALWAYS (! fits);
  }

  // mpz_class::fits_ushort_p
  {
    int        fits;
    mpz_class  z;
    z = 0;         fits = z.fits_ushort_p(); ASSERT_ALWAYS (fits);
    z--;           fits = z.fits_ushort_p(); ASSERT_ALWAYS (! fits);
    z = USHRT_MAX; fits = z.fits_ushort_p(); ASSERT_ALWAYS (fits);
    z++;           fits = z.fits_ushort_p(); ASSERT_ALWAYS (! fits);
  }

  // mpz_class::get_mpz_t
  {
    mpz_class  z(0);
    mpz_ptr    p = z.get_mpz_t();
    ASSERT_ALWAYS (mpz_cmp_ui (p, 0) == 0);
  }
  {
    mpz_class  z(0);
    mpz_srcptr p = z.get_mpz_t();
    ASSERT_ALWAYS (mpz_cmp_ui (p, 0) == 0);
  }

  // mpz_class::get_d
  // mpz_class::get_si
  // mpz_class::get_ui
  {
    mpz_class  z(123);
    { double d = z.get_d();  ASSERT_ALWAYS (d == 123.0); }
    { long   l = z.get_si(); ASSERT_ALWAYS (l == 123L); }
    { long   u = z.get_ui(); ASSERT_ALWAYS (u == 123L); }
  }
  {
    mpz_class  z(-123);
    { double d = z.get_d();  ASSERT_ALWAYS (d == -123.0); }
    { long   l = z.get_si(); ASSERT_ALWAYS (l == -123L); }
  }

  // mpz_class::get_str
  {
    mpz_class  z(123);
    string     s;
    s = z.get_str(); ASSERT_ALWAYS (s == "123");
    s = z.get_str(16); ASSERT_ALWAYS (s == "7b");
    s = z.get_str(-16); ASSERT_ALWAYS (s == "7B");
  }

  // mpz_class::set_str
  {
    mpz_class  z;
    int        ret;
    ret = z.set_str ("123", 10);  ASSERT_ALWAYS (ret == 0 && z == 123);
    ret = z.set_str ("7b",  16);  ASSERT_ALWAYS (ret == 0 && z == 123);
    ret = z.set_str ("7B",  16);  ASSERT_ALWAYS (ret == 0 && z == 123);
    ret = z.set_str ("0x7B", 0);  ASSERT_ALWAYS (ret == 0 && z == 123);

    ret = z.set_str (string("123"), 10);  ASSERT_ALWAYS (ret == 0 && z == 123);
    ret = z.set_str (string("7b"),  16);  ASSERT_ALWAYS (ret == 0 && z == 123);
    ret = z.set_str (string("7B"),  16);  ASSERT_ALWAYS (ret == 0 && z == 123);
    ret = z.set_str (string("0x7B"), 0);  ASSERT_ALWAYS (ret == 0 && z == 123);
  }
}

void
check_mpq (void)
{
  // mpq_class::canonicalize
  {
    mpq_class  q(12,9);
    q.canonicalize();
    ASSERT_ALWAYS (q.get_num() == 4);
    ASSERT_ALWAYS (q.get_den() == 3);
  }

  // mpq_class::get_d
  {
    mpq_class  q(123);
    { double d = q.get_d();  ASSERT_ALWAYS (d == 123.0); }
  }
  {
    mpq_class  q(-123);
    { double d = q.get_d();  ASSERT_ALWAYS (d == -123.0); }
  }

  // mpq_class::get_mpq_t
  {
    mpq_class  q(0);
    mpq_ptr    p = q.get_mpq_t();
    ASSERT_ALWAYS (mpq_cmp_ui (p, 0, 1) == 0);
  }
  {
    mpq_class  q(0);
    mpq_srcptr p = q.get_mpq_t();
    ASSERT_ALWAYS (mpq_cmp_ui (p, 0, 1) == 0);
  }

  // mpq_class::get_num, mpq_class::get_den
  {
    mpq_class  q(4,5);
    mpz_class  z;
    z = q.get_num(); ASSERT_ALWAYS (z == 4);
    z = q.get_den(); ASSERT_ALWAYS (z == 5);
  }

  // mpq_class::get_num_mpz_t, mpq_class::get_den_mpz_t
  {
    mpq_class  q(4,5);
    mpz_ptr    p;
    p = q.get_num_mpz_t(); ASSERT_ALWAYS (mpz_cmp_ui (p, 4) == 0);
    p = q.get_den_mpz_t(); ASSERT_ALWAYS (mpz_cmp_ui (p, 5) == 0);
  }
  {
    mpq_class  q(4,5);
    mpz_srcptr p;
    p = q.get_num_mpz_t(); ASSERT_ALWAYS (mpz_cmp_ui (p, 4) == 0);
    p = q.get_den_mpz_t(); ASSERT_ALWAYS (mpz_cmp_ui (p, 5) == 0);
  }

  // mpq_class::get_str
  {
    mpq_class  q(17,11);
    string     s;
    s = q.get_str();    ASSERT_ALWAYS (s == "17/11");
    s = q.get_str(10);  ASSERT_ALWAYS (s == "17/11");
    s = q.get_str(16);  ASSERT_ALWAYS (s == "11/b");
    s = q.get_str(-16); ASSERT_ALWAYS (s == "11/B");
  }

  // mpq_class::set_str
  {
    mpq_class  q;
    int        ret;
    ret = q.set_str ("123", 10);     ASSERT_ALWAYS (ret == 0 && q == 123);
    ret = q.set_str ("4/5", 10);     ASSERT_ALWAYS (ret == 0 && q == mpq_class(4,5));
    ret = q.set_str ("7b",  16);     ASSERT_ALWAYS (ret == 0 && q == 123);
    ret = q.set_str ("7B",  16);     ASSERT_ALWAYS (ret == 0 && q == 123);
    ret = q.set_str ("0x7B", 0);     ASSERT_ALWAYS (ret == 0 && q == 123);
    ret = q.set_str ("0x10/17", 0);  ASSERT_ALWAYS (ret == 0 && q == mpq_class(16,17));

    ret = q.set_str (string("4/5"), 10);  ASSERT_ALWAYS (ret == 0 && q == mpq_class(4,5));
    ret = q.set_str (string("123"), 10);  ASSERT_ALWAYS (ret == 0 && q == 123);
    ret = q.set_str (string("7b"),  16);  ASSERT_ALWAYS (ret == 0 && q == 123);
    ret = q.set_str (string("7B"),  16);  ASSERT_ALWAYS (ret == 0 && q == 123);
    ret = q.set_str (string("0x7B"), 0);  ASSERT_ALWAYS (ret == 0 && q == 123);
    ret = q.set_str (string("0x10/17"), 0);  ASSERT_ALWAYS (ret == 0 && q == mpq_class(16,17));
  }
}

void
check_mpf (void)
{
  // mpf_class::fits_sint_p
  {
    int        fits;
    mpf_class  f (0, 2*8*sizeof(int));
    f = INT_MIN; fits = f.fits_sint_p(); ASSERT_ALWAYS (fits);
    f--;         fits = f.fits_sint_p(); ASSERT_ALWAYS (! fits);
    f = INT_MAX; fits = f.fits_sint_p(); ASSERT_ALWAYS (fits);
    f++;         fits = f.fits_sint_p(); ASSERT_ALWAYS (! fits);
  }

  // mpf_class::fits_uint_p
  {
    int        fits;
    mpf_class  f (0, 2*8*sizeof(int));
    f = 0;        fits = f.fits_uint_p(); ASSERT_ALWAYS (fits);
    f--;          fits = f.fits_uint_p(); ASSERT_ALWAYS (! fits);
    f = UINT_MAX; fits = f.fits_uint_p(); ASSERT_ALWAYS (fits);
    f++;          fits = f.fits_uint_p(); ASSERT_ALWAYS (! fits);
  }

  // mpf_class::fits_slong_p
  {
    int        fits;
    mpf_class  f (0, 2*8*sizeof(long));
    f = LONG_MIN; fits = f.fits_slong_p(); ASSERT_ALWAYS (fits);
    f--;          fits = f.fits_slong_p(); ASSERT_ALWAYS (! fits);
    f = LONG_MAX; fits = f.fits_slong_p(); ASSERT_ALWAYS (fits);
    f++;          fits = f.fits_slong_p(); ASSERT_ALWAYS (! fits);
  }

  // mpf_class::fits_ulong_p
  {
    int        fits;
    mpf_class  f (0, 2*8*sizeof(long));
    f = 0;         fits = f.fits_ulong_p(); ASSERT_ALWAYS (fits);
    f--;           fits = f.fits_ulong_p(); ASSERT_ALWAYS (! fits);
    f = ULONG_MAX; fits = f.fits_ulong_p(); ASSERT_ALWAYS (fits);
    f++;           fits = f.fits_ulong_p(); ASSERT_ALWAYS (! fits);
  }

  // mpf_class::fits_sshort_p
  {
    int        fits;
    mpf_class  f (0, 2*8*sizeof(short));
    f = SHRT_MIN; fits = f.fits_sshort_p(); ASSERT_ALWAYS (fits);
    f--;          fits = f.fits_sshort_p(); ASSERT_ALWAYS (! fits);
    f = SHRT_MAX; fits = f.fits_sshort_p(); ASSERT_ALWAYS (fits);
    f++;          fits = f.fits_sshort_p(); ASSERT_ALWAYS (! fits);
  }

  // mpf_class::fits_ushort_p
  {
    int        fits;
    mpf_class  f (0, 2*8*sizeof(short));
    f = 0;         fits = f.fits_ushort_p(); ASSERT_ALWAYS (fits);
    f--;           fits = f.fits_ushort_p(); ASSERT_ALWAYS (! fits);
    f = USHRT_MAX; fits = f.fits_ushort_p(); ASSERT_ALWAYS (fits);
    f++;           fits = f.fits_ushort_p(); ASSERT_ALWAYS (! fits);
  }

  // mpf_class::get_d
  // mpf_class::get_si
  // mpf_class::get_ui
  {
    mpf_class  f(123);
    { double d = f.get_d();  ASSERT_ALWAYS (d == 123.0); }
    { long   l = f.get_si(); ASSERT_ALWAYS (l == 123L); }
    { long   u = f.get_ui(); ASSERT_ALWAYS (u == 123L); }
  }
  {
    mpf_class  f(-123);
    { double d = f.get_d();  ASSERT_ALWAYS (d == -123.0); }
    { long   l = f.get_si(); ASSERT_ALWAYS (l == -123L); }
  }

  // mpf_class::get_prec
  {
    mpf_class  f;
    ASSERT_ALWAYS (f.get_prec() == mpf_get_default_prec());
  }

  // mpf_class::get_str
  {
    mpf_class  f(123);
    string     s;
    mp_exp_t   e;
    s = f.get_str(e);        ASSERT_ALWAYS (s == "123" && e == 3);
    s = f.get_str(e,  16);   ASSERT_ALWAYS (s == "7b"  && e == 2);
    s = f.get_str(e, -16);   ASSERT_ALWAYS (s == "7B"  && e == 2);
    s = f.get_str(e, 10, 2); ASSERT_ALWAYS (s == "12"  && e == 3);
    s = f.get_str(e, 10, 1); ASSERT_ALWAYS (s == "1"   && e == 3);
  }

  // mpf_class::set_str
  {
    mpf_class  f;
    int        ret;
    ret = f.set_str ("123",     10);  ASSERT_ALWAYS (ret == 0 && f == 123);
    ret = f.set_str ("123e1",   10);  ASSERT_ALWAYS (ret == 0 && f == 1230);
    ret = f.set_str ("1230e-1", 10);  ASSERT_ALWAYS (ret == 0 && f == 123);
    ret = f.set_str ("7b",      16);  ASSERT_ALWAYS (ret == 0 && f == 123);
    ret = f.set_str ("7B",      16);  ASSERT_ALWAYS (ret == 0 && f == 123);
    ret = f.set_str ("7B@1",    16);  ASSERT_ALWAYS (ret == 0 && f == 1968);
    ret = f.set_str ("7B0@-1",  16);  ASSERT_ALWAYS (ret == 0 && f == 123);

    ret = f.set_str (string("123"),     10);  ASSERT_ALWAYS (ret == 0 && f == 123);
    ret = f.set_str (string("123e1"),   10);  ASSERT_ALWAYS (ret == 0 && f == 1230);
    ret = f.set_str (string("1230e-1"), 10);  ASSERT_ALWAYS (ret == 0 && f == 123);
    ret = f.set_str (string("7b"),      16);  ASSERT_ALWAYS (ret == 0 && f == 123);
    ret = f.set_str (string("7B"),      16);  ASSERT_ALWAYS (ret == 0 && f == 123);
    ret = f.set_str (string("7B@1"),    16);  ASSERT_ALWAYS (ret == 0 && f == 1968);
    ret = f.set_str (string("7B0@-1"),  16);  ASSERT_ALWAYS (ret == 0 && f == 123);
  }

  // mpf_class::set_prec
  {
    mpf_class  f;
    f.set_prec (256);
    ASSERT_ALWAYS (f.get_prec () >= 256);
  }

  // mpf_class::set_prec_raw
  {
    mpf_class  f (0, 100 * GMP_NUMB_BITS);
    f.set_prec_raw (5 * GMP_NUMB_BITS);
    ASSERT_ALWAYS (f.get_prec () >= 5 * GMP_NUMB_BITS);
    ASSERT_ALWAYS (f.get_prec () < 100 * GMP_NUMB_BITS);
    f.set_prec_raw (100 * GMP_NUMB_BITS);
  }
}


static void
check_mpfr_hypot (void)
{
#if WANT_MPFR
  // hypot(mpfr_t,mpfr_t)
  {
    mpfr_class  x(3), y(4), got;
    got = hypot (x, y);
    ASSERT_ALWAYS (got == 5);
  }

  // hypot(mpfr_t,long)
  {
    mpfr_class  x(0), got;
    long        y = LONG_MIN;
    unsigned long  want = - (unsigned long) LONG_MIN;
    got = hypot (x, y);
    ASSERT_ALWAYS (got == want);
  }
  // hypot(long,mpfr_t)
  {
    mpfr_class  y(0), got;
    long        x = LONG_MIN;
    unsigned long  want = - (unsigned long) LONG_MIN;
    got = hypot (x, y);
    ASSERT_ALWAYS (got == want);
  }

  // hypot(mpfr_t,double)
  {
    mpfr_class  x(3), got;
    got = hypot (x, -4.0);
    ASSERT_ALWAYS (got == 5);
  }
  // hypot(double,mpfr_t)
  {
    mpfr_class  y(4), got;
    got = hypot (3.0, y);
    ASSERT_ALWAYS (got == 5);
  }
#endif
}


int
main (void)
{
  tests_start();

  check_mpz();
  check_mpq();
  check_mpf();
  check_mpfr_hypot ();

  tests_end();
  return 0;
}