arch/mips/math-emu/dp_fmax.c

   1 // SPDX-License-Identifier: GPL-2.0-only
   2 /*
   3  * IEEE754 floating point arithmetic
   4  * double precision: MIN{,A}.f
   5  * MIN : Scalar Floating-Point Minimum
   6  * MINA: Scalar Floating-Point argument with Minimum Absolute Value
   7  *
   8  * MIN.D : FPR[fd] = minNum(FPR[fs],FPR[ft])
   9  * MINA.D: FPR[fd] = maxNumMag(FPR[fs],FPR[ft])
  10  *
  11  * MIPS floating point support
  12  * Copyright (C) 2015 Imagination Technologies, Ltd.
  13  * Author: Markos Chandras <markos.chandras@imgtec.com>
  14  */
  15
  16 #include "ieee754dp.h"
  17
  18 union ieee754dp ieee754dp_fmax(union ieee754dp x, union ieee754dp y)
  19 {
  20         COMPXDP;
  21         COMPYDP;
  22
  23         EXPLODEXDP;
  24         EXPLODEYDP;
  25
  26         FLUSHXDP;
  27         FLUSHYDP;
  28
  29         ieee754_clearcx();
  30
  31         switch (CLPAIR(xc, yc)) {
  32         case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_SNAN):
  33         case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_SNAN):
  34         case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_SNAN):
  35         case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_SNAN):
  36         case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_SNAN):
  37                 return ieee754dp_nanxcpt(y);
  38
  39         case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_SNAN):
  40         case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_QNAN):
  41         case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_ZERO):
  42         case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_NORM):
  43         case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_DNORM):
  44         case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_INF):
  45                 return ieee754dp_nanxcpt(x);
  46
  47         /*
  48          * Quiet NaN handling
  49          */
  50
  51         /*
  52          *    The case of both inputs quiet NaNs
  53          */
  54         case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_QNAN):
  55                 return x;
  56
  57         /*
  58          *    The cases of exactly one input quiet NaN (numbers
  59          *    are here preferred as returned values to NaNs)
  60          */
  61         case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_QNAN):
  62         case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_QNAN):
  63         case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_QNAN):
  64         case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_QNAN):
  65                 return x;
  66
  67         case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_ZERO):
  68         case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_NORM):
  69         case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_DNORM):
  70         case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_INF):
  71                 return y;
  72
  73         /*
  74          * Infinity and zero handling
  75          */
  76         case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_ZERO):
  77         case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_NORM):
  78         case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_DNORM):
  79         case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_ZERO):
  80         case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_ZERO):
  81                 return xs ? y : x;
  82
  83         case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_INF):
  84         case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_INF):
  85         case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_INF):
  86         case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_INF):
  87         case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_NORM):
  88         case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_DNORM):
  89                 return ys ? x : y;
  90
  91         case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_ZERO):
  92                 return ieee754dp_zero(xs & ys);
  93
  94         case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_DNORM):
  95                 DPDNORMX;
  96                 fallthrough;
  97         case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_DNORM):
  98                 DPDNORMY;
  99                 break;
 100
 101         case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_NORM):
 102                 DPDNORMX;
 103         }
 104
 105         /* Finally get to do some computation */
 106
 107         assert(xm & DP_HIDDEN_BIT);
 108         assert(ym & DP_HIDDEN_BIT);
 109
 110         /* Compare signs */
 111         if (xs > ys)
 112                 return y;
 113         else if (xs < ys)
 114                 return x;
 115
 116         /* Signs of inputs are equal, let's compare exponents */
 117         if (xs == 0) {
 118                 /* Inputs are both positive */
 119                 if (xe > ye)
 120                         return x;
 121                 else if (xe < ye)
 122                         return y;
 123         } else {
 124                 /* Inputs are both negative */
 125                 if (xe > ye)
 126                         return y;
 127                 else if (xe < ye)
 128                         return x;
 129         }
 130
 131         /* Signs and exponents of inputs are equal, let's compare mantissas */
 132         if (xs == 0) {
 133                 /* Inputs are both positive, with equal signs and exponents */
 134                 if (xm <= ym)
 135                         return y;
 136                 return x;
 137         }
 138         /* Inputs are both negative, with equal signs and exponents */
 139         if (xm <= ym)
 140                 return x;
 141         return y;
 142 }
 143
 144 union ieee754dp ieee754dp_fmaxa(union ieee754dp x, union ieee754dp y)
 145 {
 146         COMPXDP;
 147         COMPYDP;
 148
 149         EXPLODEXDP;
 150         EXPLODEYDP;
 151
 152         FLUSHXDP;
 153         FLUSHYDP;
 154
 155         ieee754_clearcx();
 156
 157         switch (CLPAIR(xc, yc)) {
 158         case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_SNAN):
 159         case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_SNAN):
 160         case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_SNAN):
 161         case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_SNAN):
 162         case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_SNAN):
 163                 return ieee754dp_nanxcpt(y);
 164
 165         case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_SNAN):
 166         case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_QNAN):
 167         case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_ZERO):
 168         case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_NORM):
 169         case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_DNORM):
 170         case CLPAIR(IEEE754_CLASS_SNAN, IEEE754_CLASS_INF):
 171                 return ieee754dp_nanxcpt(x);
 172
 173         /*
 174          * Quiet NaN handling
 175          */
 176
 177         /*
 178          *    The case of both inputs quiet NaNs
 179          */
 180         case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_QNAN):
 181                 return x;
 182
 183         /*
 184          *    The cases of exactly one input quiet NaN (numbers
 185          *    are here preferred as returned values to NaNs)
 186          */
 187         case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_QNAN):
 188         case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_QNAN):
 189         case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_QNAN):
 190         case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_QNAN):
 191                 return x;
 192
 193         case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_ZERO):
 194         case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_NORM):
 195         case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_DNORM):
 196         case CLPAIR(IEEE754_CLASS_QNAN, IEEE754_CLASS_INF):
 197                 return y;
 198
 199         /*
 200          * Infinity and zero handling
 201          */
 202         case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_INF):
 203                 return ieee754dp_inf(xs & ys);
 204
 205         case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_ZERO):
 206         case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_NORM):
 207         case CLPAIR(IEEE754_CLASS_INF, IEEE754_CLASS_DNORM):
 208         case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_ZERO):
 209         case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_ZERO):
 210                 return x;
 211
 212         case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_INF):
 213         case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_INF):
 214         case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_INF):
 215         case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_NORM):
 216         case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_DNORM):
 217                 return y;
 218
 219         case CLPAIR(IEEE754_CLASS_ZERO, IEEE754_CLASS_ZERO):
 220                 return ieee754dp_zero(xs & ys);
 221
 222         case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_DNORM):
 223                 DPDNORMX;
 224                 fallthrough;
 225         case CLPAIR(IEEE754_CLASS_NORM, IEEE754_CLASS_DNORM):
 226                 DPDNORMY;
 227                 break;
 228
 229         case CLPAIR(IEEE754_CLASS_DNORM, IEEE754_CLASS_NORM):
 230                 DPDNORMX;
 231         }
 232
 233         /* Finally get to do some computation */
 234
 235         assert(xm & DP_HIDDEN_BIT);
 236         assert(ym & DP_HIDDEN_BIT);
 237
 238         /* Compare exponent */
 239         if (xe > ye)
 240                 return x;
 241         else if (xe < ye)
 242                 return y;
 243
 244         /* Compare mantissa */
 245         if (xm < ym)
 246                 return y;
 247         else if (xm > ym)
 248                 return x;
 249         else if (xs == 0)
 250                 return x;
 251         return y;
 252 }