tests/add.dat

   1 # Data file for mpc_add.
   2 #
   3 # Copyright (C) 2008 Philippe Th\'eveny
   4 #
   5 # This file is part of the MPC Library.
   6 #
   7 # The MPC 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 MPC 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 MPC 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 # The line format respects the parameter order in function prototype as
  23 # follow:
  24 #
  25 # INEX_RE INEX_IM PREC_ROP_RE  ROP_RE  PREC_ROP_IM  ROP_IM  PREC_OP1_RE  OP1_RE  PREC_OP1_IM  OP1_IM  PREC_OP2_RE  OP2_RE  PREC_OP2_IM  OP2_IM  RND_RE  RND_IM
  26 #
  27 # where op1 = op1_re + i * op1_im, op2 = op2_re + i * op2_im, rop = rop_re + i * rop_im,
  28 # The data are read from the file and stored in variables op1, op2, rop using
  29 # rounding to nearest when needed, for instance: rop_re is ROP_RE rounded to
  30 # nearest to the precision of PREC_ROP_RE.
  31 # ROP_RE is checked against Re(op1 + op2) rounded to the precision PREC_ROP_RE
  32 #   in the direction RND_RE
  33 # ROP_IM is checked against Im(op1 + op2) rounded to the precision PREC_ROP_IM
  34 #   in the direction RND_IM
  35 # INEX_RE is the ternary value for the real part with the following notation:
  36 # "?" ternary value not checked
  37 # "+" if ROP_RE is greater than the exact mathematical result
  38 # "0" if ROP_RE is exactly the mathematical result
  39 # "-" if ROP_RE is less than the exact mathematical result
  40 # (m.m. INEX_IM)
  41 # rounding modes notation:
  42 # "N" is rounding to nearest
  43 # "Z" is rounding towards zero
  44 # "U" is rounding towards plus infinity
  45 # "D" is rounding towards minus infinity
  46 # Use prefixes "0b" for values in base two, "0x" for values in base sixteen,
  47 #   no prefix for value in base ten.
  48 # In all bases, "nan" is NaN, "inf" is infinity;
  49 # The sign of the result is checked with "+inf", "-inf", "-0", or "+0".
  50
  51 # special values (following ISO C99 standard)
  52 0 0 53 -inf 53 -inf    53 -inf 53 -inf 53 -inf 53 -inf N Z
  53 0 0 53 -inf 53  nan    53 -inf 53 +inf 53   +1 53 -inf Z U
  54 0 0 53 +inf 53 -inf    53 +inf 53 -inf 53 +inf 53   +0 U D
  55 0 0 53 +inf 53 +inf    53 +inf 53 +inf 53   -0 53   -1 D N
  56 0 0 53 -inf 53 -inf    53 -inf 53 -inf 53   -0 53   -1 N U
  57 0 0 53  nan 53  nan    53 -inf 53 +inf 53 +inf 53  nan Z D # check this is true in ISO C99
  58 0 0 53  nan 53 -inf    53 +inf 53 -inf 53  nan 53   -1 U N
  59 0 0 53 +inf 53  nan    53 +inf 53 +inf 53   -0 53  nan D Z
  60 0 0 53  nan 53  nan    53 -inf 53 -inf 53  nan 53  nan N D
  61
  62 0 0 53 +inf 53 -inf    53   -1 53 -inf 53 +inf 53   -1 N D
  63 0 0 53  nan 53   +1    53 -inf 53   +1 53 +inf 53   -0 Z N
  64 0 0 53   +1 53 -inf    53   +1 53 -inf 53   -0 53   +1 U Z
  65 0 0 53 +inf 53   +1    53 +inf 53   +1 53   -0 53   -0 D U
  66 0 0 53 +inf 53  nan    53   -1 53 -inf 53 +inf 53  nan N N
  67 0 0 53  nan 53   +2    53 -inf 53   +1 53  nan 53   +1 Z Z
  68 0 0 53   +1 53  nan    53   +1 53 -inf 53   -0 53  nan U U # check this
  69 0 0 53  nan 53  nan    53 +inf 53   +1 53  nan 53  nan D D
  70
  71 0 0 53   -0 53  nan    53   -0 53 -inf 53   +0 53 +inf D D
  72 0 0 53 -inf 53   +0    53 -inf 53   +0 53   +1 53   -0 N Z
  73 0 0 53   +0 53 -inf    53   +0 53 -inf 53   -0 53   -0 Z U
  74 0 0 53  nan 53  nan    53 +inf 53   +0 53 -inf 53  nan U D # check
  75 0 0 53  nan 53 -inf    53   -0 53 -inf 53  nan 53   -1 D N
  76 0 0 53 -inf 53  nan    53 -inf 53   +0 53   +0 53  nan N U
  77 0 0 53  nan 53  nan    53   +0 53 -inf 53  nan 53  nan Z D
  78
  79 0 0 53   +0 53   -0    53   +0 53   +1 53   -0 53   -1 Z D
  80 0 0 53   -1 53   -0    53   -1 53   -0 53   +0 53   -0 U N
  81 0 0 53  nan 53 +inf    53   -0 53   +1 53  nan 53 +inf D Z
  82 0 0 53   +0 53  nan    53   +1 53   -0 53   -1 53  nan N D
  83 0 0 53  nan 53   +1    53   +0 53   +1 53  nan 53   -0 Z N
  84 0 0 53  nan 53  nan    53   -1 53   -0 53  nan 53  nan U Z
  85
  86 0 0 53   +0 53   +0    53   -0 53   +0 53   +0 53   -0 U Z
  87 0 0 53  nan 53 -inf    53   +0 53   -0 53  nan 53 -inf D U
  88 0 0 53   -1 53  nan    53   +0 53   +0 53   -1 53  nan N N
  89 0 0 53  nan 53   -0    53   -0 53   -0 53  nan 53   -0 Z Z
  90 0 0 53  nan 53  nan    53   -0 53   +0 53  nan 53  nan U U
  91
  92 0 0 53  nan 53  nan    53  nan 53 -inf 53  nan 53 +inf U U # check
  93 0 0 53 +inf 53  nan    53 +inf 53  nan 53   -1 53  nan D D
  94 0 0 53  nan 53 -inf    53  nan 53 -inf 53  nan 53   -0 N Z
  95 0 0 53  nan 53  nan    53 -inf 53  nan 53  nan 53  nan Z U # check
  96
  97 0 0 53  nan 53  nan    53   +1 53  nan 53  nan 53   -1 Z U
  98 0 0 53  nan 53  nan    53  nan 53   +1 53   -0 53  nan U D
  99 0 0 53  nan 53  nan    53   -1 53  nan 53  nan 53  nan D N
 100
 101 0 0 53  nan 53  nan    53  nan 53   +0 53   +0 53  nan D N
 102 0 0 53  nan 53  nan    53   +0 53  nan 53  nan 53  nan N U
 103
 104 0 0 53  nan 53  nan    53  nan 53  nan 53  nan 53  nan N U
 105
 106 # pure real argument
 107 + 0 53  0x10000000000001p-52  53 -0    53 +1 53 -0 53 0x10000000000001p-105 53 -0 N N
 108 - 0 53  0x10000000000000p-52  53 -0    53 +1 53 -0 53 0x10000000000001p-105 53 -0 Z Z
 109 + 0 53  0x10000000000001p-52  53 -0    53 +1 53 -0 53 0x10000000000001p-105 53 -0 U U
 110 - 0 53  0x10000000000000p-52  53 -0    53 +1 53 -0 53 0x10000000000001p-105 53 -0 D D
 111
 112 # pure imaginary argument
 113 0 + 53 -0 53  0x10000000000001p-52     53 -0 53 0x10000000000001p-105 53 -0 53 +1 N N
 114 0 - 53 +0 53  0x10000000000000p-52     53 +0 53 0x10000000000001p-105 53 -0 53 +1 Z Z
 115 0 + 53 +0 53  0x10000000000001p-52     53 +0 53 0x10000000000001p-105 53 -0 53 +1 U U
 116 0 - 53 -0 53  0x10000000000000p-52     53 -0 53 0x10000000000001p-105 53 -0 53 +1 D D
 117