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 #include <linux/config.h>
  29
  30 #define _FP_W_TYPE_SIZE         32
  31 #define _FP_W_TYPE              unsigned long
  32 #define _FP_WS_TYPE             signed long
  33 #define _FP_I_TYPE              long
  34
  35 #define _FP_MUL_MEAT_S(R,X,Y)                                   \
  36   _FP_MUL_MEAT_1_wide(_FP_WFRACBITS_S,R,X,Y,umul_ppmm)
  37 #define _FP_MUL_MEAT_D(R,X,Y)                                   \
  38   _FP_MUL_MEAT_2_wide(_FP_WFRACBITS_D,R,X,Y,umul_ppmm)
  39 #define _FP_MUL_MEAT_Q(R,X,Y)                                   \
  40   _FP_MUL_MEAT_4_wide(_FP_WFRACBITS_Q,R,X,Y,umul_ppmm)
  41
  42 #define _FP_DIV_MEAT_S(R,X,Y)   _FP_DIV_MEAT_1_udiv(S,R,X,Y)
  43 #define _FP_DIV_MEAT_D(R,X,Y)   _FP_DIV_MEAT_2_udiv(D,R,X,Y)
  44 #define _FP_DIV_MEAT_Q(R,X,Y)   _FP_DIV_MEAT_4_udiv(Q,R,X,Y)
  45
  46 #define _FP_NANFRAC_S           ((_FP_QNANBIT_S << 1) - 1)
  47 #define _FP_NANFRAC_D           ((_FP_QNANBIT_D << 1) - 1), -1
  48 #define _FP_NANFRAC_Q           ((_FP_QNANBIT_Q << 1) - 1), -1, -1, -1
  49 #define _FP_NANSIGN_S           0
  50 #define _FP_NANSIGN_D           0
  51 #define _FP_NANSIGN_Q           0
  52
  53 #define _FP_KEEPNANFRACP 1
  54
  55 /*
  56  * If one NaN is signaling and the other is not,
  57  * we choose that one, otherwise we choose X.
  58  */
  59 #define _FP_CHOOSENAN(fs, wc, R, X, Y, OP)                      \
  60   do {                                                          \
  61     if ((_FP_FRAC_HIGH_RAW_##fs(X) & _FP_QNANBIT_##fs)          \
  62         && !(_FP_FRAC_HIGH_RAW_##fs(Y) & _FP_QNANBIT_##fs))     \
  63       {                                                         \
  64         R##_s = Y##_s;                                          \
  65         _FP_FRAC_COPY_##wc(R,Y);                                \
  66       }                                                         \
  67     else                                                        \
  68       {                                                         \
  69         R##_s = X##_s;                                          \
  70         _FP_FRAC_COPY_##wc(R,X);                                \
  71       }                                                         \
  72     R##_c = FP_CLS_NAN;                                         \
  73   } while (0)
  74
  75 /* Some assembly to speed things up. */
  76 #define __FP_FRAC_ADD_3(r2,r1,r0,x2,x1,x0,y2,y1,y0) ({          \
  77         unsigned int __r2 = (x2) + (y2);                        \
  78         unsigned int __r1 = (x1);                               \
  79         unsigned int __r0 = (x0);                               \
  80         __asm__ ("   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__ ("   alr %1,%2\n"                               \
  90                  "   brc 12,0f\n"                               \
  91                  "   ahi %0,1\n"                                \
  92                  "0:"                                           \
  93                  : "+&d" (__r2), "+&d" (__r1)                   \
  94                  : "d" (y1) : "cc" );                           \
  95         (r2) = __r2;                                            \
  96         (r1) = __r1;                                            \
  97         (r0) = __r0;                                            \
  98 })
  99
 100 #define __FP_FRAC_SUB_3(r2,r1,r0,x2,x1,x0,y2,y1,y0) ({          \
 101         unsigned int __r2 = (x2) - (y2);                        \
 102         unsigned int __r1 = (x1);                               \
 103         unsigned int __r0 = (x0);                               \
 104         __asm__ ("   slr %2,%3\n"                               \
 105                  "   brc 3,0f\n"                                \
 106                  "   lhi 0,1\n"                                 \
 107                  "   slr %1,0\n"                                \
 108                  "   brc 3,0f\n"                                \
 109                  "   slr %0,0\n"                                \
 110                  "0:"                                           \
 111                  : "+&d" (__r2), "+&d" (__r1), "+&d" (__r0)     \
 112                  : "d" (y0) : "cc", "0" );                      \
 113         __asm__ ("   slr %1,%2\n"                               \
 114                  "   brc 3,0f\n"                                \
 115                  "   ahi %0,-1\n"                               \
 116                  "0:"                                           \
 117                  : "+&d" (__r2), "+&d" (__r1)                   \
 118                  : "d" (y1) : "cc" );                           \
 119         (r2) = __r2;                                            \
 120         (r1) = __r1;                                            \
 121         (r0) = __r0;                                            \
 122 })
 123
 124 #define __FP_FRAC_DEC_3(x2,x1,x0,y2,y1,y0) __FP_FRAC_SUB_3(x2,x1,x0,x2,x1,x0,y2,y1,y0)
 125
 126 /* Obtain the current rounding mode. */
 127 #define FP_ROUNDMODE    mode
 128
 129 /* Exception flags. */
 130 #define FP_EX_INVALID           0x800000
 131 #define FP_EX_DIVZERO           0x400000
 132 #define FP_EX_OVERFLOW          0x200000
 133 #define FP_EX_UNDERFLOW         0x100000
 134 #define FP_EX_INEXACT           0x080000
 135
 136 /* We write the results always */
 137 #define FP_INHIBIT_RESULTS 0
 138
 139 #endif