/* tfprintf.c -- test file for mpfr_fprintf and mpfr_vfprintf
Copyright 2008-2017 Free Software Foundation, Inc.
Contributed by the AriC and Caramba projects, INRIA.
This file is part of the GNU MPFR Library.
The GNU MPFR 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 3 of the License, or (at your
option) any later version.
The GNU MPFR 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 MPFR Library; see the file COPYING.LESSER. If not, see
http://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc.,
51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA. */
#ifdef HAVE_STDARG
#include <stdarg.h>
#include <stdio.h>
#include <stdlib.h>
#include <float.h>
#include <stddef.h>
#include "mpfr-intmax.h"
#include "mpfr-test.h"
#if MPFR_VERSION >= MPFR_VERSION_NUM(2,4,0)
#define QUOTE(X) NAME(X)
#define NAME(X) #X
#define check_length(num_test, var, value, var_spec) \
if ((var) != (value)) \
{ \
printf ("Error in test #%d: mpfr_vfprintf printed %"QUOTE(var_spec) \
" characters instead of %d\n", (num_test), (var), (value)); \
exit (1); \
}
#define check_length_with_cmp(num_test, var, value, cmp, var_spec) \
if (cmp != 0) \
{ \
mpfr_printf ("Error in test #%d, mpfr_vfprintf printed %" \
QUOTE(var_spec)" characters instead of %d\n", \
(num_test), (var), (value)); \
exit (1); \
}
/* limit for random precision in random() */
const int prec_max_printf = 5000;
static void
check (FILE *fout, const char *fmt, mpfr_t x)
{
if (mpfr_fprintf (fout, fmt, x) == -1)
{
mpfr_printf ("Error in mpfr_fprintf(fout, \"%s\", %Re)\n",
fmt, x);
exit (1);
}
fputc ('\n', fout);
}
static void
check_vfprintf (FILE *fout, const char *fmt, ...)
{
va_list ap;
va_start (ap, fmt);
if (mpfr_vfprintf (fout, fmt, ap) == -1)
{
mpfr_printf ("Error in mpfr_vfprintf(fout, \"%s\", ...)\n", fmt);
va_end (ap);
exit (1);
}
va_end (ap);
fputc ('\n', fout);
}
static void
check_special (FILE *fout)
{
mpfr_t x;
mpfr_init (x);
mpfr_set_inf (x, 1);
check (fout, "%Ra", x);
check (fout, "%Rb", x);
check (fout, "%Re", x);
check (fout, "%Rf", x);
check (fout, "%Rg", x);
check_vfprintf (fout, "%Ra", x);
check_vfprintf (fout, "%Rb", x);
check_vfprintf (fout, "%Re", x);
check_vfprintf (fout, "%Rf", x);
check_vfprintf (fout, "%Rg", x);
mpfr_set_inf (x, -1);
check (fout, "%Ra", x);
check (fout, "%Rb", x);
check (fout, "%Re", x);
check (fout, "%Rf", x);
check (fout, "%Rg", x);
check_vfprintf (fout, "%Ra", x);
check_vfprintf (fout, "%Rb", x);
check_vfprintf (fout, "%Re", x);
check_vfprintf (fout, "%Rf", x);
check_vfprintf (fout, "%Rg", x);
mpfr_set_nan (x);
check (fout, "%Ra", x);
check (fout, "%Rb", x);
check (fout, "%Re", x);
check (fout, "%Rf", x);
check (fout, "%Rg", x);
check_vfprintf (fout, "%Ra", x);
check_vfprintf (fout, "%Rb", x);
check_vfprintf (fout, "%Re", x);
check_vfprintf (fout, "%Rf", x);
check_vfprintf (fout, "%Rg", x);
mpfr_clear (x);
}
static void
check_mixed (FILE *fout)
{
int ch = 'a';
#ifndef NPRINTF_HH
signed char sch = -1;
unsigned char uch = 1;
#endif
short sh = -1;
unsigned short ush = 1;
int i = -1;
int j = 1;
unsigned int ui = 1;
long lo = -1;
unsigned long ulo = 1;
float f = -1.25;
double d = -1.25;
#if !defined(NPRINTF_T) || !defined(NPRINTF_L)
long double ld = -1.25;
#endif
#ifndef NPRINTF_T
ptrdiff_t p = 1, saved_p;
#endif
size_t sz = 1;
mpz_t mpz;
mpq_t mpq;
mpf_t mpf;
mpfr_rnd_t rnd = MPFR_RNDN;
mp_size_t limb_size = 3;
mp_limb_t limb[3];
mpfr_t mpfr;
mpfr_prec_t prec = 53;
mpz_init (mpz);
mpz_set_ui (mpz, ulo);
mpq_init (mpq);
mpq_set_si (mpq, lo, ulo);
mpf_init (mpf);
mpf_set_q (mpf, mpq);
mpfr_init2 (mpfr, prec);
mpfr_set_f (mpfr, mpf, MPFR_RNDN);
limb[0] = limb[1] = limb[2] = ~ (mp_limb_t) 0;
check_vfprintf (fout, "a. %Ra, b. %u, c. %lx%n", mpfr, ui, ulo, &j);
check_length (1, j, 22, d);
check_vfprintf (fout, "a. %c, b. %Rb, c. %u, d. %li%ln", i, mpfr, i,
lo, &ulo);
check_length (2, ulo, 36, lu);
check_vfprintf (fout, "a. %hi, b. %*f, c. %Re%hn", ush, 3, f, mpfr, &ush);
check_length (3, ush, 29, hu);
check_vfprintf (fout, "a. %hi, b. %f, c. %#.2Rf%n", sh, d, mpfr, &i);
check_length (4, i, 29, d);
check_vfprintf (fout, "a. %R*A, b. %Fe, c. %i%zn", rnd, mpfr, mpf, sz,
&sz);
check_length (5, (unsigned long) sz, 34, lu); /* no format specifier "%zu" in C89 */
check_vfprintf (fout, "a. %Pu, b. %c, c. %Zi%Zn", prec, ch, mpz, &mpz);
check_length_with_cmp (6, mpz, 17, mpz_cmp_ui (mpz, 17), Zi);
check_vfprintf (fout, "%% a. %#.0RNg, b. %Qx%Rn, c. %p", mpfr, mpq, &mpfr,
(void *) &i);
check_length_with_cmp (7, mpfr, 15, mpfr_cmp_ui (mpfr, 15), Rg);
#ifndef NPRINTF_T
saved_p = p;
check_vfprintf (fout, "%% a. %RNg, b. %Qx, c. %td%tn", mpfr, mpq, p, &p);
if (p != 20)
mpfr_fprintf (stderr, "Error in test 8, got '%% a. %RNg, b. %Qx, c. %td'\n", mpfr, mpq, saved_p);
check_length (8, (long) p, 20, ld); /* no format specifier "%td" in C89 */
#endif
#ifndef NPRINTF_L
check_vfprintf (fout, "a. %RA, b. %Lf, c. %QX%zn", mpfr, ld, mpq, &sz);
check_length (9, (unsigned long) sz, 30, lu); /* no format specifier "%zu" in C89 */
#endif
#ifndef NPRINTF_HH
check_vfprintf (fout, "a. %hhi, b. %RA, c. %hhu%hhn", sch, mpfr, uch, &uch);
check_length (10, (unsigned int) uch, 22, u); /* no format specifier "%hhu" in C89 */
#endif
#if (__GNU_MP_VERSION * 10 + __GNU_MP_VERSION_MINOR) >= 42
/* The 'M' specifier was added in gmp 4.2.0 */
check_vfprintf (fout, "a. %Mx b. %Re%Mn", limb[0], mpfr, &limb[0]);
if (limb[0] != 14 + GMP_NUMB_BITS / 4 || limb[1] != ~ (mp_limb_t) 0
|| limb[2] != ~ (mp_limb_t) 0)
{
printf ("Error in test #11: mpfr_vfprintf did not print %d characters"
" as expected\n", 14 + (int) GMP_NUMB_BITS / 4);
exit (1);
}
limb[0] = ~ (mp_limb_t) 0;
/* we tell vfprintf that limb array is 2 cells wide
and check it doesn't go through */
check_vfprintf (fout, "a. %Re .b %Nx%Nn", mpfr, limb, limb_size, limb,
limb_size - 1);
if (limb[0] != 14 + 3 * GMP_NUMB_BITS / 4 || limb[1] != (mp_limb_t) 0
|| limb[2] != ~ (mp_limb_t) 0)
{
printf ("Error in test #12: mpfr_vfprintf did not print %d characters"
" as expected\n", 14 + (int) GMP_NUMB_BITS / 4);
exit (1);
}
#endif
#if defined(HAVE_LONG_LONG) && !defined(NPRINTF_LL)
{
long long llo = -1;
unsigned long long ullo = 1;
check_vfprintf (fout, "a. %Re, b. %llx%Qn", mpfr, ullo, &mpq);
check_length_with_cmp (21, mpq, 16, mpq_cmp_ui (mpq, 16, 1), Qu);
check_vfprintf (fout, "a. %lli, b. %Rf%Fn", llo, mpfr, &mpf);
check_length_with_cmp (22, mpf, 19, mpf_cmp_ui (mpf, 19), Fg);
}
#endif
#if defined(_MPFR_H_HAVE_INTMAX_T) && !defined(NPRINTF_J)
{
intmax_t im = -1;
uintmax_t uim = 1;
check_vfprintf (fout, "a. %*RA, b. %ji%Qn", 10, mpfr, im, &mpq);
check_length_with_cmp (31, mpq, 20, mpq_cmp_ui (mpq, 20, 1), Qu);
check_vfprintf (fout, "a. %.*Re, b. %jx%Fn", 10, mpfr, uim, &mpf);
check_length_with_cmp (32, mpf, 25, mpf_cmp_ui (mpf, 25), Fg);
}
#endif
mpfr_clear (mpfr);
mpf_clear (mpf);
mpq_clear (mpq);
mpz_clear (mpz);
}
static void
check_random (FILE *fout, int nb_tests)
{
int i;
mpfr_t x;
mpfr_rnd_t rnd;
char flag[] =
{
'-',
'+',
' ',
'#',
'0', /* no ambiguity: first zeros are flag zero*/
'\''
};
char specifier[] =
{
'a',
'b',
'e',
'f',
'g'
};
mpfr_exp_t old_emin, old_emax;
old_emin = mpfr_get_emin ();
old_emax = mpfr_get_emax ();
mpfr_init (x);
for (i = 0; i < nb_tests; ++i)
{
int ret;
int j, jmax;
int spec, prec;
#define FMT_SIZE 13
char fmt[FMT_SIZE]; /* at most something like "%-+ #0'.*R*f" */
char *ptr = fmt;
tests_default_random (x, 256, MPFR_EMIN_MIN, MPFR_EMAX_MAX, 0);
rnd = RND_RAND ();
spec = (int) (randlimb () % 5);
jmax = (spec == 3 || spec == 4) ? 6 : 5; /* ' flag only with %f or %g */
/* advantage small precision */
prec = (int) (randlimb () % ((randlimb () % 2) ? 10 : prec_max_printf));
if (spec == 3
&& (mpfr_get_exp (x) > prec_max_printf
|| mpfr_get_exp (x) < -prec_max_printf))
/* change style 'f' to style 'e' when number x is large */
--spec;
*ptr++ = '%';
for (j = 0; j < jmax; j++)
{
if (randlimb () % 3 == 0)
*ptr++ = flag[j];
}
*ptr++ = '.';
*ptr++ = '*';
*ptr++ = 'R';
*ptr++ = '*';
*ptr++ = specifier[spec];
*ptr = '\0';
MPFR_ASSERTD (ptr - fmt < FMT_SIZE);
mpfr_fprintf (fout, "mpfr_fprintf(fout, \"%s\", %d, %s, %Re)\n",
fmt, prec, mpfr_print_rnd_mode (rnd), x);
ret = mpfr_fprintf (fout, fmt, prec, rnd, x);
if (ret == -1)
{
if (spec == 3
&& (MPFR_GET_EXP (x) > INT_MAX || MPFR_GET_EXP (x) < -INT_MAX))
/* normal failure: x is too large to be output with full precision */
{
mpfr_fprintf (fout, "too large !");
}
else
{
mpfr_printf ("Error in mpfr_fprintf(fout, \"%s\", %d, %s, %Re)\n",
fmt, prec, mpfr_print_rnd_mode (rnd), x);
exit (1);
}
}
mpfr_fprintf (fout, "\n");
}
mpfr_set_emin (old_emin);
mpfr_set_emax (old_emax);
mpfr_clear (x);
}
static void
bug_20090316 (FILE *fout)
{
mpfr_t x;
mpfr_init2 (x, 53);
/* bug 20090316: fixed in r6112 */
mpfr_set_ui_2exp (x, 0x60fa2916, -30, MPFR_RNDN);
check (fout, "%-#.4095RDg\n", x);
mpfr_clear (x);
}
int
main (int argc, char *argv[])
{
FILE *fout;
int N;
tests_start_mpfr ();
/* with no argument: prints to /dev/null,
tfprintf N: prints N tests to stdout */
if (argc == 1)
{
N = 1000;
fout = fopen ("/dev/null", "w");
/* If we failed to open this device, try with a dummy file */
if (fout == NULL)
{
fout = fopen ("tfprintf_out.txt", "w");
if (fout == NULL)
{
printf ("Can't open /dev/null or a temporary file\n");
exit (1);
}
}
}
else
{
fout = stdout;
N = atoi (argv[1]);
}
check_special (fout);
check_mixed (fout);
check_random (fout, N);
bug_20090316 (fout);
fclose (fout);
tests_end_mpfr ();
return 0;
}
#else /* MPFR_VERSION */
int
main (void)
{
printf ("Warning! Test disabled for this MPFR version.\n");
return 0;
}
#endif /* MPFR_VERSION */
#else /* HAVE_STDARG */
int
main (void)
{
/* We have nothing to test. */
return 77;
}
#endif /* HAVE_STDARG */