arch/s390/include/asm/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 int
  31 #define _FP_WS_TYPE             signed int
  32 #define _FP_I_TYPE              int
  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 volatile(                                           \
  80                 "       alr     %2,%3\n"                        \
  81                 "       brc     12,0f\n"                        \
  82                 "       lhi     0,1\n"                          \
  83                 "       alr     %1,0\n"                         \
  84                 "       brc     12,0f\n"                        \
  85                 "       alr     %0,0\n"                         \
  86                 "0:"                                            \
  87                 : "+&d" (__r2), "+&d" (__r1), "+&d" (__r0)      \
  88                 : "d" (y0), "i" (1) : "cc", "0" );              \
  89         asm volatile(                                           \
  90                 "       alr     %1,%2\n"                        \
  91                 "       brc     12,0f\n"                        \
  92                 "       ahi     %0,1\n"                         \
  93                 "0:"                                            \
  94                 : "+&d" (__r2), "+&d" (__r1)                    \
  95                 : "d" (y1) : "cc");                             \
  96         (r2) = __r2;                                            \
  97         (r1) = __r1;                                            \
  98         (r0) = __r0;                                            \
  99 })
 100
 101 #define __FP_FRAC_SUB_3(r2,r1,r0,x2,x1,x0,y2,y1,y0) ({          \
 102         unsigned int __r2 = (x2) - (y2);                        \
 103         unsigned int __r1 = (x1);                               \
 104         unsigned int __r0 = (x0);                               \
 105         asm volatile(                                           \
 106                 "       slr   %2,%3\n"                          \
 107                 "       brc     3,0f\n"                         \
 108                 "       lhi     0,1\n"                          \
 109                 "       slr     %1,0\n"                         \
 110                 "       brc     3,0f\n"                         \
 111                 "       slr     %0,0\n"                         \
 112                 "0:"                                            \
 113                 : "+&d" (__r2), "+&d" (__r1), "+&d" (__r0)      \
 114                 : "d" (y0) : "cc", "0");                        \
 115         asm volatile(                                           \
 116                 "       slr     %1,%2\n"                        \
 117                 "       brc     3,0f\n"                         \
 118                 "       ahi     %0,-1\n"                        \
 119                 "0:"                                            \
 120                 : "+&d" (__r2), "+&d" (__r1)                    \
 121                 : "d" (y1) : "cc");                             \
 122         (r2) = __r2;                                            \
 123         (r1) = __r1;                                            \
 124         (r0) = __r0;                                            \
 125 })
 126
 127 #define __FP_FRAC_DEC_3(x2,x1,x0,y2,y1,y0) __FP_FRAC_SUB_3(x2,x1,x0,x2,x1,x0,y2,y1,y0)
 128
 129 /* Obtain the current rounding mode. */
 130 #define FP_ROUNDMODE    mode
 131
 132 /* Exception flags. */
 133 #define FP_EX_INVALID           0x800000
 134 #define FP_EX_DIVZERO           0x400000
 135 #define FP_EX_OVERFLOW          0x200000
 136 #define FP_EX_UNDERFLOW         0x100000
 137 #define FP_EX_INEXACT           0x080000
 138
 139 /* We write the results always */
 140 #define FP_INHIBIT_RESULTS 0
 141
 142 #endif