xref: /netbsd-src/external/lgpl3/mpc/dist/tests/check_data.c (revision 367b82799ab709709d3c3b541df56a2a14644d3e)
1 /* check_data.c -- Check computed data against reference result.
2 
3 Copyright (C) 2012, 2013, 2014, 2022 INRIA
4 
5 This file is part of GNU MPC.
6 
7 GNU MPC is free software; you can redistribute it and/or modify it under
8 the terms of the GNU Lesser General Public License as published by the
9 Free Software Foundation; either version 3 of the License, or (at your
10 option) any later version.
11 
12 GNU MPC is distributed in the hope that it will be useful, but WITHOUT ANY
13 WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
14 FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for
15 more details.
16 
17 You should have received a copy of the GNU Lesser General Public License
18 along with this program. If not, see http://www.gnu.org/licenses/ .
19 */
20 
21 #include "mpc-tests.h"
22 
23 #define MPC_INEX_CMP(c, r, i)                                 \
24   (((r) == TERNARY_NOT_CHECKED || (r) == MPC_INEX_RE(c))      \
25    && ((i) == TERNARY_NOT_CHECKED || (i) == MPC_INEX_IM (c)))
26 
27 static int
check_param(mpc_operand_t * got,mpc_operand_t * expected,mpc_param_t t)28 check_param  (mpc_operand_t* got, mpc_operand_t* expected, mpc_param_t t)
29 {
30   switch (t)
31     {
32     case NATIVE_INT:
33       return got->i == expected->i;
34     case NATIVE_UL:
35       return got->ui == expected->ui;
36     case NATIVE_L:
37       return got->si == expected->si;
38     case NATIVE_D:
39       return got->d == expected->d;
40 
41     case GMP_Z:
42       return mpz_cmp (got->mpz, expected->mpz);
43     case GMP_Q:
44       return mpq_cmp (got->mpq, expected->mpq);
45     case GMP_F:
46       return mpf_cmp (got->mpf, expected->mpf);
47 
48     case MPFR_INEX:
49       return expected->mpfr_inex == TERNARY_NOT_CHECKED
50         || got->mpfr_inex == expected->mpfr_inex;
51 
52     case MPFR:
53       return tpl_check_mpfr_data (got->mpfr,
54                                   expected->mpfr_data);
55 
56     case MPC_INEX:
57       return MPC_INEX_CMP (got->mpc_inex,
58                            expected->mpc_inex_data.real,
59                            expected->mpc_inex_data.imag);
60 
61     case MPC:
62       return tpl_check_mpc_data (got->mpc,
63                                  expected->mpc_data);
64 
65     case MPCC_INEX:
66       return got->mpcc_inex == expected->mpcc_inex;
67 
68     default:
69       fprintf (stderr, "check_data: unsupported type.\n");
70       exit (1);
71     }
72 }
73 
74 void
check_data(mpc_datafile_context_t * dc,mpc_fun_param_t * params,int reused_op)75 check_data (mpc_datafile_context_t* dc, mpc_fun_param_t* params, int reused_op)
76 {
77   int out, i;
78   int total = params->nbout + params->nbin;
79 
80   for (out = 0; out < params->nbout; out++)
81     {
82       if (!check_param (&(params->P[out]), &(params->P[total + out]),
83                         params->T[out]))
84         {
85           printf ("%s() failed", params->name);
86           if (dc != NULL)
87             printf (" (line %lu, file %s)",
88                     dc->test_line_number, dc->pathname);
89           else
90             printf (" with random parameter%c",
91                     params->nbin > 2 ? 's' : '\0');
92 
93           if (reused_op != 0)
94             printf (" when reusing input parameter op%d as output parameter",
95                     reused_op - params->nbout);
96           printf ("\n");
97 
98           for (i = 0; i < params->nbin; i++)
99             {
100               printf ("op%d", i + 1);
101               print_parameter (params, params->nbout + i);
102             }
103 
104           for (i = 0; i < params->nbout; i++)
105             {
106               if ((params->T[i] == MPFR_INEX && params->T[out] != MPFR_INEX)
107                   || (params->T[i] == MPC_INEX && params->T[out] != MPC_INEX))
108                 continue; /* don't print inexact flag if it is correct */
109 
110               printf ("     got%c",
111                       (total + i > params->nbout ? ' ' : i + '0'));
112               print_parameter (params, i);
113               printf ("expected%c",
114                       (total + i > params->nbout ? ' ' : i + '0'));
115               print_parameter (params, total + i);
116             }
117           printf ("\n");
118           exit (1);
119         }
120     }
121 }
122