include/asm-s390/sfp-machine.h

   1 /* Machine-dependent software floating-point definitions.
   2    S/390 kernel version.
   3    Copyright (C) 1997,1998,1999 Free Software Foundation, Inc.
   4    This file is part of the GNU C Library.
   5    Contributed by Richard Henderson (rth@cygnus.com),
   6                   Jakub Jelinek (jj@ultra.linux.cz),
   7                   David S. Miller (davem@redhat.com) and
   8                   Peter Maydell (pmaydell@chiark.greenend.org.uk).
   9
  10    The GNU C Library is free software; you can redistribute it and/or
  11    modify it under the terms of the GNU Library General Public License as
  12    published by the Free Software Foundation; either version 2 of the
  13    License, or (at your option) any later version.
  14
  15    The GNU C Library is distributed in the hope that it will be useful,
  16    but WITHOUT ANY WARRANTY; without even the implied warranty of
  17    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  18    Library General Public License for more details.
  19
  20    You should have received a copy of the GNU Library General Public
  21    License along with the GNU C Library; see the file COPYING.LIB.  If
  22    not, write to the Free Software Foundation, Inc.,
  23    59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.  */
  24
  25 #ifndef _SFP_MACHINE_H
  26 #define _SFP_MACHINE_H
  27
  28
  29 #define _FP_W_TYPE_SIZE         32
  30 #define _FP_W_TYPE              unsigned long
  31 #define _FP_WS_TYPE             signed long
  32 #define _FP_I_TYPE              long
  33
  34 #define _FP_MUL_MEAT_S(R,X,Y)                                   \
  35   _FP_MUL_MEAT_1_wide(_FP_WFRACBITS_S,R,X,Y,umul_ppmm)
  36 #define _FP_MUL_MEAT_D(R,X,Y)                                   \
  37   _FP_MUL_MEAT_2_wide(_FP_WFRACBITS_D,R,X,Y,umul_ppmm)
  38 #define _FP_MUL_MEAT_Q(R,X,Y)                                   \
  39   _FP_MUL_MEAT_4_wide(_FP_WFRACBITS_Q,R,X,Y,umul_ppmm)
  40
  41 #define _FP_DIV_MEAT_S(R,X,Y)   _FP_DIV_MEAT_1_udiv(S,R,X,Y)
  42 #define _FP_DIV_MEAT_D(R,X,Y)   _FP_DIV_MEAT_2_udiv(D,R,X,Y)
  43 #define _FP_DIV_MEAT_Q(R,X,Y)   _FP_DIV_MEAT_4_udiv(Q,R,X,Y)
  44
  45 #define _FP_NANFRAC_S           ((_FP_QNANBIT_S << 1) - 1)
  46 #define _FP_NANFRAC_D           ((_FP_QNANBIT_D << 1) - 1), -1
  47 #define _FP_NANFRAC_Q           ((_FP_QNANBIT_Q << 1) - 1), -1, -1, -1
  48 #define _FP_NANSIGN_S           0
  49 #define _FP_NANSIGN_D           0
  50 #define _FP_NANSIGN_Q           0
  51
  52 #define _FP_KEEPNANFRACP 1
  53
  54 /*
  55  * If one NaN is signaling and the other is not,
  56  * we choose that one, otherwise we choose X.
  57  */
  58 #define _FP_CHOOSENAN(fs, wc, R, X, Y, OP)                      \
  59   do {                                                          \
  60     if ((_FP_FRAC_HIGH_RAW_##fs(X) & _FP_QNANBIT_##fs)          \
  61         && !(_FP_FRAC_HIGH_RAW_##fs(Y) & _FP_QNANBIT_##fs))     \
  62       {                                                         \
  63         R##_s = Y##_s;                                          \
  64         _FP_FRAC_COPY_##wc(R,Y);                                \
  65       }                                                         \
  66     else                                                        \
  67       {                                                         \
  68         R##_s = X##_s;                                          \
  69         _FP_FRAC_COPY_##wc(R,X);                                \
  70       }                                                         \
  71     R##_c = FP_CLS_NAN;                                         \
  72   } while (0)
  73
  74 /* Some assembly to speed things up. */
  75 #define __FP_FRAC_ADD_3(r2,r1,r0,x2,x1,x0,y2,y1,y0) ({          \
  76         unsigned int __r2 = (x2) + (y2);                        \
  77         unsigned int __r1 = (x1);                               \
  78         unsigned int __r0 = (x0);                               \
  79         __asm__ ("   alr %2,%3\n"                               \
  80                  "   brc 12,0f\n"                               \
  81                  "   lhi 0,1\n"                                 \
  82                  "   alr %1,0\n"                                \
  83                  "   brc 12,0f\n"                               \
  84                  "   alr %0,0\n"                                \
  85                  "0:"                                           \
  86                  : "+&d" (__r2), "+&d" (__r1), "+&d" (__r0)     \
  87                  : "d" (y0), "i" (1) : "cc", "0" );             \
  88         __asm__ ("   alr %1,%2\n"                               \
  89                  "   brc 12,0f\n"                               \
  90                  "   ahi %0,1\n"                                \
  91                  "0:"                                           \
  92                  : "+&d" (__r2), "+&d" (__r1)                   \
  93                  : "d" (y1) : "cc" );                           \
  94         (r2) = __r2;                                            \
  95         (r1) = __r1;                                            \
  96         (r0) = __r0;                                            \
  97 })
  98
  99 #define __FP_FRAC_SUB_3(r2,r1,r0,x2,x1,x0,y2,y1,y0) ({          \
 100         unsigned int __r2 = (x2) - (y2);                        \
 101         unsigned int __r1 = (x1);                               \
 102         unsigned int __r0 = (x0);                               \
 103         __asm__ ("   slr %2,%3\n"                               \
 104                  "   brc 3,0f\n"                                \
 105                  "   lhi 0,1\n"                                 \
 106                  "   slr %1,0\n"                                \
 107                  "   brc 3,0f\n"                                \
 108                  "   slr %0,0\n"                                \
 109                  "0:"                                           \
 110                  : "+&d" (__r2), "+&d" (__r1), "+&d" (__r0)     \
 111                  : "d" (y0) : "cc", "0" );                      \
 112         __asm__ ("   slr %1,%2\n"                               \
 113                  "   brc 3,0f\n"                                \
 114                  "   ahi %0,-1\n"                               \
 115                  "0:"                                           \
 116                  : "+&d" (__r2), "+&d" (__r1)                   \
 117                  : "d" (y1) : "cc" );                           \
 118         (r2) = __r2;                                            \
 119         (r1) = __r1;                                            \
 120         (r0) = __r0;                                            \
 121 })
 122
 123 #define __FP_FRAC_DEC_3(x2,x1,x0,y2,y1,y0) __FP_FRAC_SUB_3(x2,x1,x0,x2,x1,x0,y2,y1,y0)
 124
 125 /* Obtain the current rounding mode. */
 126 #define FP_ROUNDMODE    mode
 127
 128 /* Exception flags. */
 129 #define FP_EX_INVALID           0x800000
 130 #define FP_EX_DIVZERO           0x400000
 131 #define FP_EX_OVERFLOW          0x200000
 132 #define FP_EX_UNDERFLOW         0x100000
 133 #define FP_EX_INEXACT           0x080000
 134
 135 /* We write the results always */
 136 #define FP_INHIBIT_RESULTS 0
 137
 138 #endif