LCOV - code coverage report
Current view: top level - gcc - dfp.c (source / functions) Hit Total Coverage
Test: gcc.info Lines: 302 352 85.8 %
Date: 2020-04-04 11:58:09 Functions: 22 23 95.7 %
Legend: Lines: hit not hit | Branches: + taken - not taken # not executed Branches: 0 0 -

           Branch data     Line data    Source code
       1                 :            : /* Decimal floating point support.
       2                 :            :    Copyright (C) 2005-2020 Free Software Foundation, Inc.
       3                 :            : 
       4                 :            : This file is part of GCC.
       5                 :            : 
       6                 :            : GCC is free software; you can redistribute it and/or modify it under
       7                 :            : the terms of the GNU General Public License as published by the Free
       8                 :            : Software Foundation; either version 3, or (at your option) any later
       9                 :            : version.
      10                 :            : 
      11                 :            : GCC is distributed in the hope that it will be useful, but WITHOUT ANY
      12                 :            : WARRANTY; without even the implied warranty of MERCHANTABILITY or
      13                 :            : FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
      14                 :            : for more details.
      15                 :            : 
      16                 :            : You should have received a copy of the GNU General Public License
      17                 :            : along with GCC; see the file COPYING3.  If not see
      18                 :            : <http://www.gnu.org/licenses/>.  */
      19                 :            : 
      20                 :            : #include "config.h"
      21                 :            : #include "system.h"
      22                 :            : #include "coretypes.h"
      23                 :            : #include "tm.h"
      24                 :            : #include "tree.h"
      25                 :            : #include "dfp.h"
      26                 :            : 
      27                 :            : /* The order of the following headers is important for making sure
      28                 :            :    decNumber structure is large enough to hold decimal128 digits.  */
      29                 :            : 
      30                 :            : #include "decimal128.h"
      31                 :            : #include "decimal64.h"
      32                 :            : #include "decimal32.h"
      33                 :            : 
      34                 :            : #ifndef WORDS_BIGENDIAN
      35                 :            : #define WORDS_BIGENDIAN 0
      36                 :            : #endif
      37                 :            : 
      38                 :            : /* Initialize R (a real with the decimal flag set) from DN.  Can
      39                 :            :    utilize status passed in via CONTEXT, if a previous operation had
      40                 :            :    interesting status.  */
      41                 :            : 
      42                 :            : static void
      43                 :     859762 : decimal_from_decnumber (REAL_VALUE_TYPE *r, decNumber *dn, decContext *context)
      44                 :            : {
      45                 :     859762 :   memset (r, 0, sizeof (REAL_VALUE_TYPE));
      46                 :            : 
      47                 :     859762 :   r->cl = rvc_normal;
      48                 :     859762 :   if (decNumberIsNaN (dn))
      49                 :        227 :     r->cl = rvc_nan;
      50                 :     859762 :   if (decNumberIsInfinite (dn))
      51                 :       1177 :     r->cl = rvc_inf;
      52                 :     859762 :   if (context->status & DEC_Overflow)
      53                 :        850 :     r->cl = rvc_inf;
      54                 :     859762 :   if (decNumberIsNegative (dn))
      55                 :     232191 :     r->sign = 1;
      56                 :     859762 :   r->decimal = 1;
      57                 :            : 
      58                 :     859762 :   if (r->cl != rvc_normal)
      59                 :            :     return;
      60                 :            : 
      61                 :     858358 :   decContextDefault (context, DEC_INIT_DECIMAL128);
      62                 :     858358 :   context->traps = 0;
      63                 :            : 
      64                 :     858358 :   decimal128FromNumber ((decimal128 *) r->sig, dn, context);
      65                 :            : }
      66                 :            : 
      67                 :            : /* Create decimal encoded R from string S.  */
      68                 :            : 
      69                 :            : void
      70                 :     124005 : decimal_real_from_string (REAL_VALUE_TYPE *r, const char *s)
      71                 :            : {
      72                 :     124005 :   decNumber dn;
      73                 :     124005 :   decContext set;
      74                 :     124005 :   decContextDefault (&set, DEC_INIT_DECIMAL128);
      75                 :     124005 :   set.traps = 0;
      76                 :            : 
      77                 :     124005 :   decNumberFromString (&dn, s, &set);
      78                 :            : 
      79                 :            :   /* It would be more efficient to store directly in decNumber format,
      80                 :            :      but that is impractical from current data structure size.
      81                 :            :      Encoding as a decimal128 is much more compact.  */
      82                 :     124005 :   decimal_from_decnumber (r, &dn, &set);
      83                 :     124005 : }
      84                 :            : 
      85                 :            : /* Initialize a decNumber from a REAL_VALUE_TYPE.  */
      86                 :            : 
      87                 :            : static void
      88                 :    1046600 : decimal_to_decnumber (const REAL_VALUE_TYPE *r, decNumber *dn)
      89                 :            : {
      90                 :    1046600 :   decContext set;
      91                 :    1046600 :   decContextDefault (&set, DEC_INIT_DECIMAL128);
      92                 :    1046600 :   set.traps = 0;
      93                 :            : 
      94                 :    1046600 :   switch (r->cl)
      95                 :            :     {
      96                 :          9 :     case rvc_zero:
      97                 :          9 :       decNumberZero (dn);
      98                 :          9 :       break;
      99                 :        385 :     case rvc_inf:
     100                 :        385 :       decNumberFromString (dn, "Infinity", &set);
     101                 :        385 :       break;
     102                 :        105 :     case rvc_nan:
     103                 :        105 :       if (r->signalling)
     104                 :          0 :         decNumberFromString (dn, "snan", &set);
     105                 :            :       else
     106                 :        105 :         decNumberFromString (dn, "nan", &set);
     107                 :            :       break;
     108                 :    1046100 :     case rvc_normal:
     109                 :    1046100 :       if (!r->decimal)
     110                 :            :         {
     111                 :            :           /* dconst{1,2,m1,half} are used in various places in
     112                 :            :              the middle-end and optimizers, allow them here
     113                 :            :              as an exception by converting them to decimal.  */
     114                 :        481 :           if (memcmp (r, &dconst1, sizeof (*r)) == 0)
     115                 :            :             {
     116                 :          0 :               decNumberFromString (dn, "1", &set);
     117                 :          0 :               break;
     118                 :            :             }
     119                 :        481 :           if (memcmp (r, &dconst2, sizeof (*r)) == 0)
     120                 :            :             {
     121                 :          1 :               decNumberFromString (dn, "2", &set);
     122                 :          1 :               break;
     123                 :            :             }
     124                 :        480 :           if (memcmp (r, &dconstm1, sizeof (*r)) == 0)
     125                 :            :             {
     126                 :        480 :               decNumberFromString (dn, "-1", &set);
     127                 :        480 :               break;
     128                 :            :             }
     129                 :          0 :           if (memcmp (r, &dconsthalf, sizeof (*r)) == 0)
     130                 :            :             {
     131                 :          0 :               decNumberFromString (dn, "0.5", &set);
     132                 :          0 :               break;
     133                 :            :             }
     134                 :          0 :           gcc_unreachable ();
     135                 :            :         }
     136                 :    1045620 :       decimal128ToNumber ((const decimal128 *) r->sig, dn);
     137                 :    1045620 :       break;
     138                 :          0 :     default:
     139                 :          0 :       gcc_unreachable ();
     140                 :            :     }
     141                 :            : 
     142                 :            :   /* Fix up sign bit.  */
     143                 :    1046600 :   if (r->sign != decNumberIsNegative (dn))
     144                 :        130 :     dn->bits ^= DECNEG;
     145                 :    1046600 : }
     146                 :            : 
     147                 :            : /* Encode a real into an IEEE 754 decimal32 type.  */
     148                 :            : 
     149                 :            : void
     150                 :       9256 : encode_decimal32 (const struct real_format *fmt ATTRIBUTE_UNUSED,
     151                 :            :                   long *buf, const REAL_VALUE_TYPE *r)
     152                 :            : {
     153                 :       9256 :   decNumber dn;
     154                 :       9256 :   decimal32 d32;
     155                 :       9256 :   decContext set;
     156                 :       9256 :   int32_t image;
     157                 :            : 
     158                 :       9256 :   decContextDefault (&set, DEC_INIT_DECIMAL128);
     159                 :       9256 :   set.traps = 0;
     160                 :            : 
     161                 :       9256 :   decimal_to_decnumber (r, &dn);
     162                 :       9256 :   decimal32FromNumber (&d32, &dn, &set);
     163                 :            : 
     164                 :       9256 :   memcpy (&image, d32.bytes, sizeof (int32_t));
     165                 :       9256 :   buf[0] = image;
     166                 :       9256 : }
     167                 :            : 
     168                 :            : /* Decode an IEEE 754 decimal32 type into a real.  */
     169                 :            : 
     170                 :            : void
     171                 :        967 : decode_decimal32 (const struct real_format *fmt ATTRIBUTE_UNUSED,
     172                 :            :                   REAL_VALUE_TYPE *r, const long *buf)
     173                 :            : {
     174                 :        967 :   decNumber dn;
     175                 :        967 :   decimal32 d32;
     176                 :        967 :   decContext set;
     177                 :        967 :   int32_t image;
     178                 :            : 
     179                 :        967 :   decContextDefault (&set, DEC_INIT_DECIMAL128);
     180                 :        967 :   set.traps = 0;
     181                 :            : 
     182                 :        967 :   image = buf[0];
     183                 :        967 :   memcpy (&d32.bytes, &image, sizeof (int32_t));
     184                 :            : 
     185                 :        967 :   decimal32ToNumber (&d32, &dn);
     186                 :        967 :   decimal_from_decnumber (r, &dn, &set);
     187                 :        967 : }
     188                 :            : 
     189                 :            : /* Encode a real into an IEEE 754 decimal64 type.  */
     190                 :            : 
     191                 :            : void
     192                 :       8807 : encode_decimal64 (const struct real_format *fmt ATTRIBUTE_UNUSED,
     193                 :            :                   long *buf, const REAL_VALUE_TYPE *r)
     194                 :            : {
     195                 :       8807 :   decNumber dn;
     196                 :       8807 :   decimal64 d64;
     197                 :       8807 :   decContext set;
     198                 :       8807 :   int32_t image;
     199                 :            : 
     200                 :       8807 :   decContextDefault (&set, DEC_INIT_DECIMAL128);
     201                 :       8807 :   set.traps = 0;
     202                 :            : 
     203                 :       8807 :   decimal_to_decnumber (r, &dn);
     204                 :       8807 :   decimal64FromNumber (&d64, &dn, &set);
     205                 :            : 
     206                 :       8807 :   if (WORDS_BIGENDIAN == FLOAT_WORDS_BIG_ENDIAN)
     207                 :            :     {
     208                 :       8807 :       memcpy (&image, &d64.bytes[0], sizeof (int32_t));
     209                 :       8807 :       buf[0] = image;
     210                 :       8807 :       memcpy (&image, &d64.bytes[4], sizeof (int32_t));
     211                 :       8807 :       buf[1] = image;
     212                 :            :     }
     213                 :            :   else
     214                 :            :     {
     215                 :            :       memcpy (&image, &d64.bytes[4], sizeof (int32_t));
     216                 :            :       buf[0] = image;
     217                 :            :       memcpy (&image, &d64.bytes[0], sizeof (int32_t));
     218                 :            :       buf[1] = image;
     219                 :            :     }
     220                 :       8807 : }
     221                 :            : 
     222                 :            : /* Decode an IEEE 754 decimal64 type into a real.  */
     223                 :            : 
     224                 :            : void
     225                 :       3716 : decode_decimal64 (const struct real_format *fmt ATTRIBUTE_UNUSED,
     226                 :            :                   REAL_VALUE_TYPE *r, const long *buf)
     227                 :            : {
     228                 :       3716 :   decNumber dn;
     229                 :       3716 :   decimal64 d64;
     230                 :       3716 :   decContext set;
     231                 :       3716 :   int32_t image;
     232                 :            : 
     233                 :       3716 :   decContextDefault (&set, DEC_INIT_DECIMAL128);
     234                 :       3716 :   set.traps = 0;
     235                 :            : 
     236                 :       3716 :   if (WORDS_BIGENDIAN == FLOAT_WORDS_BIG_ENDIAN)
     237                 :            :     {
     238                 :       3716 :       image = buf[0];
     239                 :       3716 :       memcpy (&d64.bytes[0], &image, sizeof (int32_t));
     240                 :       3716 :       image = buf[1];
     241                 :       3716 :       memcpy (&d64.bytes[4], &image, sizeof (int32_t));
     242                 :            :     }
     243                 :            :   else
     244                 :            :     {
     245                 :            :       image = buf[1];
     246                 :            :       memcpy (&d64.bytes[0], &image, sizeof (int32_t));
     247                 :            :       image = buf[0];
     248                 :            :       memcpy (&d64.bytes[4], &image, sizeof (int32_t));
     249                 :            :     }
     250                 :            : 
     251                 :       3716 :   decimal64ToNumber (&d64, &dn);
     252                 :       3716 :   decimal_from_decnumber (r, &dn, &set);
     253                 :       3716 : }
     254                 :            : 
     255                 :            : /* Encode a real into an IEEE 754 decimal128 type.  */
     256                 :            : 
     257                 :            : void
     258                 :       8755 : encode_decimal128 (const struct real_format *fmt ATTRIBUTE_UNUSED,
     259                 :            :                    long *buf, const REAL_VALUE_TYPE *r)
     260                 :            : {
     261                 :       8755 :   decNumber dn;
     262                 :       8755 :   decContext set;
     263                 :       8755 :   decimal128 d128;
     264                 :       8755 :   int32_t image;
     265                 :            : 
     266                 :       8755 :   decContextDefault (&set, DEC_INIT_DECIMAL128);
     267                 :       8755 :   set.traps = 0;
     268                 :            : 
     269                 :       8755 :   decimal_to_decnumber (r, &dn);
     270                 :       8755 :   decimal128FromNumber (&d128, &dn, &set);
     271                 :            : 
     272                 :       8755 :   if (WORDS_BIGENDIAN == FLOAT_WORDS_BIG_ENDIAN)
     273                 :            :     {
     274                 :       8755 :       memcpy (&image, &d128.bytes[0], sizeof (int32_t));
     275                 :       8755 :       buf[0] = image;
     276                 :       8755 :       memcpy (&image, &d128.bytes[4], sizeof (int32_t));
     277                 :       8755 :       buf[1] = image;
     278                 :       8755 :       memcpy (&image, &d128.bytes[8], sizeof (int32_t));
     279                 :       8755 :       buf[2] = image;
     280                 :       8755 :       memcpy (&image, &d128.bytes[12], sizeof (int32_t));
     281                 :       8755 :       buf[3] = image;
     282                 :            :     }
     283                 :            :   else
     284                 :            :     {
     285                 :            :       memcpy (&image, &d128.bytes[12], sizeof (int32_t));
     286                 :            :       buf[0] = image;
     287                 :            :       memcpy (&image, &d128.bytes[8], sizeof (int32_t));
     288                 :            :       buf[1] = image;
     289                 :            :       memcpy (&image, &d128.bytes[4], sizeof (int32_t));
     290                 :            :       buf[2] = image;
     291                 :            :       memcpy (&image, &d128.bytes[0], sizeof (int32_t));
     292                 :            :       buf[3] = image;
     293                 :            :     }
     294                 :       8755 : }
     295                 :            : 
     296                 :            : /* Decode an IEEE 754 decimal128 type into a real.  */
     297                 :            : 
     298                 :            : void
     299                 :       3458 : decode_decimal128 (const struct real_format *fmt ATTRIBUTE_UNUSED,
     300                 :            :                    REAL_VALUE_TYPE *r, const long *buf)
     301                 :            : {
     302                 :       3458 :   decNumber dn;
     303                 :       3458 :   decimal128 d128;
     304                 :       3458 :   decContext set;
     305                 :       3458 :   int32_t image;
     306                 :            : 
     307                 :       3458 :   decContextDefault (&set, DEC_INIT_DECIMAL128);
     308                 :       3458 :   set.traps = 0;
     309                 :            : 
     310                 :       3458 :   if (WORDS_BIGENDIAN == FLOAT_WORDS_BIG_ENDIAN)
     311                 :            :     {
     312                 :       3458 :       image = buf[0];
     313                 :       3458 :       memcpy (&d128.bytes[0],  &image, sizeof (int32_t));
     314                 :       3458 :       image = buf[1];
     315                 :       3458 :       memcpy (&d128.bytes[4],  &image, sizeof (int32_t));
     316                 :       3458 :       image = buf[2];
     317                 :       3458 :       memcpy (&d128.bytes[8],  &image, sizeof (int32_t));
     318                 :       3458 :       image = buf[3];
     319                 :       3458 :       memcpy (&d128.bytes[12], &image, sizeof (int32_t));
     320                 :            :     }
     321                 :            :   else
     322                 :            :     {
     323                 :            :       image = buf[3];
     324                 :            :       memcpy (&d128.bytes[0],  &image, sizeof (int32_t));
     325                 :            :       image = buf[2];
     326                 :            :       memcpy (&d128.bytes[4],  &image, sizeof (int32_t));
     327                 :            :       image = buf[1];
     328                 :            :       memcpy (&d128.bytes[8],  &image, sizeof (int32_t));
     329                 :            :       image = buf[0];
     330                 :            :       memcpy (&d128.bytes[12], &image, sizeof (int32_t));
     331                 :            :     }
     332                 :            : 
     333                 :       3458 :   decimal128ToNumber (&d128, &dn);
     334                 :       3458 :   decimal_from_decnumber (r, &dn, &set);
     335                 :       3458 : }
     336                 :            : 
     337                 :            : /* Helper function to convert from a binary real internal
     338                 :            :    representation.  */
     339                 :            : 
     340                 :            : static void
     341                 :        110 : decimal_to_binary (REAL_VALUE_TYPE *to, const REAL_VALUE_TYPE *from,
     342                 :            :                    const real_format *fmt)
     343                 :            : {
     344                 :        110 :   char string[256];
     345                 :        110 :   if (from->cl == rvc_normal)
     346                 :            :     {
     347                 :        109 :       const decimal128 *const d128 = (const decimal128 *) from->sig;
     348                 :        109 :       decimal128ToString (d128, string);
     349                 :            :     }
     350                 :            :   else
     351                 :          1 :     real_to_decimal (string, from, sizeof (string), 0, 1);
     352                 :        110 :   real_from_string3 (to, string, fmt);
     353                 :        110 : }
     354                 :            : 
     355                 :            : 
     356                 :            : /* Helper function to convert from a binary real internal
     357                 :            :    representation.  */
     358                 :            : 
     359                 :            : static void
     360                 :      91812 : decimal_from_binary (REAL_VALUE_TYPE *to, const REAL_VALUE_TYPE *from)
     361                 :            : {
     362                 :      91812 :   char string[256];
     363                 :            : 
     364                 :            :   /* We convert to string, then to decNumber then to decimal128.  */
     365                 :      91812 :   real_to_decimal (string, from, sizeof (string), 0, 1);
     366                 :      91812 :   decimal_real_from_string (to, string);
     367                 :      91812 : }
     368                 :            : 
     369                 :            : /* Helper function to real.c:do_compare() to handle decimal internal
     370                 :            :    representation including when one of the operands is still in the
     371                 :            :    binary internal representation.  */
     372                 :            : 
     373                 :            : int
     374                 :     336941 : decimal_do_compare (const REAL_VALUE_TYPE *a, const REAL_VALUE_TYPE *b,
     375                 :            :                     int nan_result)
     376                 :            : {
     377                 :     336941 :   decContext set;
     378                 :     336941 :   decNumber dn, dn2, dn3;
     379                 :     336941 :   REAL_VALUE_TYPE a1, b1;
     380                 :            : 
     381                 :            :   /* If either operand is non-decimal, create temporary versions.  */
     382                 :     336941 :   if (!a->decimal)
     383                 :            :     {
     384                 :       5208 :       decimal_from_binary (&a1, a);
     385                 :       5208 :       a = &a1;
     386                 :            :     }
     387                 :     336941 :   if (!b->decimal)
     388                 :            :     {
     389                 :      80813 :       decimal_from_binary (&b1, b);
     390                 :      80813 :       b = &b1;
     391                 :            :     }
     392                 :            : 
     393                 :            :   /* Convert into decNumber form for comparison operation.  */
     394                 :     336941 :   decContextDefault (&set, DEC_INIT_DECIMAL128);
     395                 :     336941 :   set.traps = 0;
     396                 :     336941 :   decimal128ToNumber ((const decimal128 *) a->sig, &dn2);
     397                 :     336941 :   decimal128ToNumber ((const decimal128 *) b->sig, &dn3);
     398                 :            : 
     399                 :            :   /* Finally, do the comparison.  */
     400                 :     336941 :   decNumberCompare (&dn, &dn2, &dn3, &set);
     401                 :            : 
     402                 :            :   /* Return the comparison result.  */
     403                 :     336941 :   if (decNumberIsNaN (&dn))
     404                 :            :     return nan_result;
     405                 :     336941 :   else if (decNumberIsZero (&dn))
     406                 :            :     return 0;
     407                 :     331023 :   else if (decNumberIsNegative (&dn))
     408                 :            :     return -1;
     409                 :            :   else
     410                 :     246066 :     return 1;
     411                 :            : }
     412                 :            : 
     413                 :            : /* Helper to round_for_format, handling decimal float types.  */
     414                 :            : 
     415                 :            : void
     416                 :     574922 : decimal_round_for_format (const struct real_format *fmt, REAL_VALUE_TYPE *r)
     417                 :            : {
     418                 :     574922 :   decNumber dn;
     419                 :     574922 :   decContext set;
     420                 :            : 
     421                 :            :   /* Real encoding occurs later.  */
     422                 :     574922 :   if (r->cl != rvc_normal)
     423                 :     357520 :     return;
     424                 :            : 
     425                 :     574394 :   decContextDefault (&set, DEC_INIT_DECIMAL128);
     426                 :     574394 :   set.traps = 0;
     427                 :     574394 :   decimal128ToNumber ((decimal128 *) r->sig, &dn);
     428                 :            : 
     429                 :     574394 :   if (fmt == &decimal_quad_format)
     430                 :            :     {
     431                 :            :       /* The internal format is already in this format.  */
     432                 :            :       return;
     433                 :            :     }
     434                 :     217402 :   else if (fmt == &decimal_single_format)
     435                 :            :     {
     436                 :      62298 :       decimal32 d32;
     437                 :      62298 :       decContextDefault (&set, DEC_INIT_DECIMAL32);
     438                 :      62298 :       set.traps = 0;
     439                 :            : 
     440                 :      62298 :       decimal32FromNumber (&d32, &dn, &set);
     441                 :      62298 :       decimal32ToNumber (&d32, &dn);
     442                 :            :     }
     443                 :     155104 :   else if (fmt == &decimal_double_format)
     444                 :            :     {
     445                 :     155104 :       decimal64 d64;
     446                 :     155104 :       decContextDefault (&set, DEC_INIT_DECIMAL64);
     447                 :     155104 :       set.traps = 0;
     448                 :            : 
     449                 :     155104 :       decimal64FromNumber (&d64, &dn, &set);
     450                 :     155104 :       decimal64ToNumber (&d64, &dn);
     451                 :            :     }
     452                 :            :   else
     453                 :          0 :     gcc_unreachable ();
     454                 :            : 
     455                 :     217402 :   decimal_from_decnumber (r, &dn, &set);
     456                 :            : }
     457                 :            : 
     458                 :            : /* Extend or truncate to a new mode.  Handles conversions between
     459                 :            :    binary and decimal types.  */
     460                 :            : 
     461                 :            : void
     462                 :     549028 : decimal_real_convert (REAL_VALUE_TYPE *r, const real_format *fmt,
     463                 :            :                       const REAL_VALUE_TYPE *a)
     464                 :            : {
     465                 :     549028 :   if (a->decimal && fmt->b == 10)
     466                 :            :     return;
     467                 :       5756 :   if (a->decimal)
     468                 :        110 :       decimal_to_binary (r, a, fmt);
     469                 :            :   else
     470                 :       5646 :       decimal_from_binary (r, a);
     471                 :            : }
     472                 :            : 
     473                 :            : /* Render R_ORIG as a decimal floating point constant.  Emit DIGITS
     474                 :            :    significant digits in the result, bounded by BUF_SIZE.  If DIGITS
     475                 :            :    is 0, choose the maximum for the representation.  If
     476                 :            :    CROP_TRAILING_ZEROS, strip trailing zeros.  Currently, not honoring
     477                 :            :    DIGITS or CROP_TRAILING_ZEROS.  */
     478                 :            : 
     479                 :            : void
     480                 :         46 : decimal_real_to_decimal (char *str, const REAL_VALUE_TYPE *r_orig,
     481                 :            :                          size_t buf_size,
     482                 :            :                          size_t digits ATTRIBUTE_UNUSED,
     483                 :            :                          int crop_trailing_zeros ATTRIBUTE_UNUSED)
     484                 :            : {
     485                 :         46 :   const decimal128 *const d128 = (const decimal128*) r_orig->sig;
     486                 :            : 
     487                 :            :   /* decimal128ToString requires space for at least 24 characters;
     488                 :            :      Require two more for suffix.  */
     489                 :         46 :   gcc_assert (buf_size >= 24);
     490                 :         46 :   decimal128ToString (d128, str);
     491                 :         46 : }
     492                 :            : 
     493                 :            : static bool
     494                 :      11891 : decimal_do_add (REAL_VALUE_TYPE *r, const REAL_VALUE_TYPE *op0,
     495                 :            :                 const REAL_VALUE_TYPE *op1, int subtract_p)
     496                 :            : {
     497                 :      11891 :   decNumber dn;
     498                 :      11891 :   decContext set;
     499                 :      11891 :   decNumber dn2, dn3;
     500                 :            : 
     501                 :      11891 :   decimal_to_decnumber (op0, &dn2);
     502                 :      11891 :   decimal_to_decnumber (op1, &dn3);
     503                 :            : 
     504                 :      11891 :   decContextDefault (&set, DEC_INIT_DECIMAL128);
     505                 :      11891 :   set.traps = 0;
     506                 :            : 
     507                 :      11891 :   if (subtract_p)
     508                 :       3931 :     decNumberSubtract (&dn, &dn2, &dn3, &set);
     509                 :            :   else
     510                 :       7960 :     decNumberAdd (&dn, &dn2, &dn3, &set);
     511                 :            : 
     512                 :      11891 :   decimal_from_decnumber (r, &dn, &set);
     513                 :            : 
     514                 :            :   /* Return true, if inexact.  */
     515                 :      11891 :   return (set.status & DEC_Inexact);
     516                 :            : }
     517                 :            : 
     518                 :            : /* Compute R = OP0 * OP1.  */
     519                 :            : 
     520                 :            : static bool
     521                 :     497844 : decimal_do_multiply (REAL_VALUE_TYPE *r, const REAL_VALUE_TYPE *op0,
     522                 :            :                      const REAL_VALUE_TYPE *op1)
     523                 :            : {
     524                 :     497844 :   decContext set;
     525                 :     497844 :   decNumber dn, dn2, dn3;
     526                 :            : 
     527                 :     497844 :   decimal_to_decnumber (op0, &dn2);
     528                 :     497844 :   decimal_to_decnumber (op1, &dn3);
     529                 :            : 
     530                 :     497844 :   decContextDefault (&set, DEC_INIT_DECIMAL128);
     531                 :     497844 :   set.traps = 0;
     532                 :            : 
     533                 :     497844 :   decNumberMultiply (&dn, &dn2, &dn3, &set);
     534                 :     497844 :   decimal_from_decnumber (r, &dn, &set);
     535                 :            : 
     536                 :            :   /* Return true, if inexact.  */
     537                 :     497844 :   return (set.status & DEC_Inexact);
     538                 :            : }
     539                 :            : 
     540                 :            : /* Compute R = OP0 / OP1.  */
     541                 :            : 
     542                 :            : static bool
     543                 :        155 : decimal_do_divide (REAL_VALUE_TYPE *r, const REAL_VALUE_TYPE *op0,
     544                 :            :                    const REAL_VALUE_TYPE *op1)
     545                 :            : {
     546                 :        155 :   decContext set;
     547                 :        155 :   decNumber dn, dn2, dn3;
     548                 :            : 
     549                 :        155 :   decimal_to_decnumber (op0, &dn2);
     550                 :        155 :   decimal_to_decnumber (op1, &dn3);
     551                 :            : 
     552                 :        155 :   decContextDefault (&set, DEC_INIT_DECIMAL128);
     553                 :        155 :   set.traps = 0;
     554                 :            : 
     555                 :        155 :   decNumberDivide (&dn, &dn2, &dn3, &set);
     556                 :        155 :   decimal_from_decnumber (r, &dn, &set);
     557                 :            : 
     558                 :            :   /* Return true, if inexact.  */
     559                 :        155 :   return (set.status & DEC_Inexact);
     560                 :            : }
     561                 :            : 
     562                 :            : /* Set R to A truncated to an integral value toward zero (decimal
     563                 :            :    floating point).  */
     564                 :            : 
     565                 :            : void
     566                 :        324 : decimal_do_fix_trunc (REAL_VALUE_TYPE *r, const REAL_VALUE_TYPE *a)
     567                 :            : {
     568                 :        324 :   decNumber dn, dn2;
     569                 :        324 :   decContext set;
     570                 :            : 
     571                 :        324 :   decContextDefault (&set, DEC_INIT_DECIMAL128);
     572                 :        324 :   set.traps = 0;
     573                 :        324 :   set.round = DEC_ROUND_DOWN;
     574                 :        324 :   decimal128ToNumber ((const decimal128 *) a->sig, &dn2);
     575                 :            : 
     576                 :        324 :   decNumberToIntegralValue (&dn, &dn2, &set);
     577                 :        324 :   decimal_from_decnumber (r, &dn, &set);
     578                 :        324 : }
     579                 :            : 
     580                 :            : /* Render decimal float value R as an integer.  */
     581                 :            : 
     582                 :            : HOST_WIDE_INT
     583                 :          0 : decimal_real_to_integer (const REAL_VALUE_TYPE *r)
     584                 :            : {
     585                 :          0 :   decContext set;
     586                 :          0 :   decNumber dn, dn2, dn3;
     587                 :          0 :   REAL_VALUE_TYPE to;
     588                 :          0 :   char string[256];
     589                 :            : 
     590                 :          0 :   decContextDefault (&set, DEC_INIT_DECIMAL128);
     591                 :          0 :   set.traps = 0;
     592                 :          0 :   set.round = DEC_ROUND_DOWN;
     593                 :          0 :   decimal128ToNumber ((const decimal128 *) r->sig, &dn);
     594                 :            : 
     595                 :          0 :   decNumberToIntegralValue (&dn2, &dn, &set);
     596                 :          0 :   decNumberZero (&dn3);
     597                 :          0 :   decNumberRescale (&dn, &dn2, &dn3, &set);
     598                 :            : 
     599                 :            :   /* Convert to REAL_VALUE_TYPE and call appropriate conversion
     600                 :            :      function.  */
     601                 :          0 :   decNumberToString (&dn, string);
     602                 :          0 :   real_from_string (&to, string);
     603                 :          0 :   return real_to_integer (&to);
     604                 :            : }
     605                 :            : 
     606                 :            : /* Likewise, but returns a wide_int with PRECISION.  *FAIL is set if the
     607                 :            :    value does not fit.  */
     608                 :            : 
     609                 :            : wide_int
     610                 :        238 : decimal_real_to_integer (const REAL_VALUE_TYPE *r, bool *fail, int precision)
     611                 :            : {
     612                 :        238 :   decContext set;
     613                 :        238 :   decNumber dn, dn2, dn3;
     614                 :        238 :   REAL_VALUE_TYPE to;
     615                 :        238 :   char string[256];
     616                 :            : 
     617                 :        238 :   decContextDefault (&set, DEC_INIT_DECIMAL128);
     618                 :        238 :   set.traps = 0;
     619                 :        238 :   set.round = DEC_ROUND_DOWN;
     620                 :        238 :   decimal128ToNumber ((const decimal128 *) r->sig, &dn);
     621                 :            : 
     622                 :        238 :   decNumberToIntegralValue (&dn2, &dn, &set);
     623                 :        238 :   decNumberZero (&dn3);
     624                 :        238 :   decNumberRescale (&dn, &dn2, &dn3, &set);
     625                 :            : 
     626                 :            :   /* Convert to REAL_VALUE_TYPE and call appropriate conversion
     627                 :            :      function.  */
     628                 :        238 :   decNumberToString (&dn, string);
     629                 :        238 :   real_from_string (&to, string);
     630                 :        238 :   return real_to_integer (&to, fail, precision);
     631                 :            : }
     632                 :            : 
     633                 :            : /* Perform the decimal floating point operation described by CODE.
     634                 :            :    For a unary operation, OP1 will be NULL.  This function returns
     635                 :            :    true if the result may be inexact due to loss of precision.  */
     636                 :            : 
     637                 :            : bool
     638                 :     516484 : decimal_real_arithmetic (REAL_VALUE_TYPE *r, enum tree_code code,
     639                 :            :                          const REAL_VALUE_TYPE *op0,
     640                 :            :                          const REAL_VALUE_TYPE *op1)
     641                 :            : {
     642                 :     516484 :   REAL_VALUE_TYPE a, b;
     643                 :            : 
     644                 :            :   /* If either operand is non-decimal, create temporaries.  */
     645                 :     516484 :   if (!op0->decimal)
     646                 :            :     {
     647                 :         67 :       decimal_from_binary (&a, op0);
     648                 :         67 :       op0 = &a;
     649                 :            :     }
     650                 :     516484 :   if (op1 && !op1->decimal)
     651                 :            :     {
     652                 :         78 :       decimal_from_binary (&b, op1);
     653                 :         78 :       op1 = &b;
     654                 :            :     }
     655                 :            : 
     656                 :     516484 :   switch (code)
     657                 :            :     {
     658                 :       7960 :     case PLUS_EXPR:
     659                 :       7960 :       return decimal_do_add (r, op0, op1, 0);
     660                 :            : 
     661                 :       3931 :     case MINUS_EXPR:
     662                 :       3931 :       return decimal_do_add (r, op0, op1, 1);
     663                 :            : 
     664                 :     497844 :     case MULT_EXPR:
     665                 :     497844 :       return decimal_do_multiply (r, op0, op1);
     666                 :            : 
     667                 :        155 :     case RDIV_EXPR:
     668                 :        155 :       return decimal_do_divide (r, op0, op1);
     669                 :            : 
     670                 :          0 :     case MIN_EXPR:
     671                 :          0 :       if (op1->cl == rvc_nan)
     672                 :          0 :         *r = *op1;
     673                 :          0 :       else if (real_compare (UNLT_EXPR, op0, op1))
     674                 :          0 :         *r = *op0;
     675                 :            :       else
     676                 :          0 :         *r = *op1;
     677                 :            :       return false;
     678                 :            : 
     679                 :          0 :     case MAX_EXPR:
     680                 :          0 :       if (op1->cl == rvc_nan)
     681                 :          0 :         *r = *op1;
     682                 :          0 :       else if (real_compare (LT_EXPR, op0, op1))
     683                 :          0 :         *r = *op1;
     684                 :            :       else
     685                 :          0 :         *r = *op0;
     686                 :            :       return false;
     687                 :            : 
     688                 :       6594 :     case NEGATE_EXPR:
     689                 :       6594 :       {
     690                 :       6594 :         *r = *op0;
     691                 :            :         /* Flip sign bit.  */
     692                 :       6594 :         decimal128FlipSign ((decimal128 *) r->sig);
     693                 :            :         /* Keep sign field in sync.  */
     694                 :       6594 :         r->sign ^= 1;
     695                 :            :       }
     696                 :       6594 :       return false;
     697                 :            : 
     698                 :          0 :     case ABS_EXPR:
     699                 :          0 :       {
     700                 :          0 :         *r = *op0;
     701                 :            :         /* Clear sign bit.  */
     702                 :          0 :         decimal128ClearSign ((decimal128 *) r->sig);
     703                 :            :         /* Keep sign field in sync.  */
     704                 :          0 :         r->sign = 0;
     705                 :            :       }
     706                 :          0 :       return false;
     707                 :            : 
     708                 :          0 :     case FIX_TRUNC_EXPR:
     709                 :          0 :       decimal_do_fix_trunc (r, op0);
     710                 :          0 :       return false;
     711                 :            : 
     712                 :          0 :     default:
     713                 :          0 :       gcc_unreachable ();
     714                 :            :     }
     715                 :            : }
     716                 :            : 
     717                 :            : /* Fills R with the largest finite value representable in mode MODE.
     718                 :            :    If SIGN is nonzero, R is set to the most negative finite value.  */
     719                 :            : 
     720                 :            : void
     721                 :         82 : decimal_real_maxval (REAL_VALUE_TYPE *r, int sign, machine_mode mode)
     722                 :            : {
     723                 :         82 :   const char *max;
     724                 :            : 
     725                 :         82 :   switch (mode)
     726                 :            :     {
     727                 :            :     case E_SDmode:
     728                 :            :       max = "9.999999E96";
     729                 :            :       break;
     730                 :         20 :     case E_DDmode:
     731                 :         20 :       max = "9.999999999999999E384";
     732                 :         20 :       break;
     733                 :         16 :     case E_TDmode:
     734                 :         16 :       max = "9.999999999999999999999999999999999E6144";
     735                 :         16 :       break;
     736                 :          0 :     default:
     737                 :          0 :       gcc_unreachable ();
     738                 :            :     }
     739                 :            : 
     740                 :         82 :   decimal_real_from_string (r, max);
     741                 :         82 :   if (sign)
     742                 :         35 :     decimal128SetSign ((decimal128 *) r->sig, 1);
     743                 :         82 : }

Generated by: LCOV version 1.0

LCOV profile is generated on x86_64 machine using following configure options: configure --disable-bootstrap --enable-coverage=opt --enable-languages=c,c++,fortran,go,jit,lto --enable-host-shared. GCC test suite is run with the built compiler.