67cc32ebfd
Signed-off-by: Veres Lajos <vlajos@gmail.com> Signed-off-by: Michael Tokarev <mjt@tls.msk.ru>
434 lines
19 KiB
C
434 lines
19 KiB
C
/* Decimal context module for the decNumber C Library.
|
|
Copyright (C) 2005, 2007 Free Software Foundation, Inc.
|
|
Contributed by IBM Corporation. Author Mike Cowlishaw.
|
|
|
|
This file is part of GCC.
|
|
|
|
GCC is free software; you can redistribute it and/or modify it under
|
|
the terms of the GNU General Public License as published by the Free
|
|
Software Foundation; either version 2, or (at your option) any later
|
|
version.
|
|
|
|
In addition to the permissions in the GNU General Public License,
|
|
the Free Software Foundation gives you unlimited permission to link
|
|
the compiled version of this file into combinations with other
|
|
programs, and to distribute those combinations without any
|
|
restriction coming from the use of this file. (The General Public
|
|
License restrictions do apply in other respects; for example, they
|
|
cover modification of the file, and distribution when not linked
|
|
into a combine executable.)
|
|
|
|
GCC is distributed in the hope that it will be useful, but WITHOUT ANY
|
|
WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
|
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
|
for more details.
|
|
|
|
You should have received a copy of the GNU General Public License
|
|
along with GCC; see the file COPYING. If not, write to the Free
|
|
Software Foundation, 51 Franklin Street, Fifth Floor, Boston, MA
|
|
02110-1301, USA. */
|
|
|
|
/* ------------------------------------------------------------------ */
|
|
/* Decimal Context module */
|
|
/* ------------------------------------------------------------------ */
|
|
/* This module comprises the routines for handling arithmetic */
|
|
/* context structures. */
|
|
/* ------------------------------------------------------------------ */
|
|
|
|
#include <string.h> /* for strcmp */
|
|
#include <stdio.h> /* for printf if DECCHECK */
|
|
#include "libdecnumber/dconfig.h"
|
|
#include "libdecnumber/decContext.h"
|
|
#include "libdecnumber/decNumberLocal.h"
|
|
|
|
#if DECCHECK
|
|
/* compile-time endian tester [assumes sizeof(Int)>1] */
|
|
static const Int mfcone=1; /* constant 1 */
|
|
static const Flag *mfctop=(Flag *)&mfcone; /* -> top byte */
|
|
#define LITEND *mfctop /* named flag; 1=little-endian */
|
|
#endif
|
|
|
|
/* ------------------------------------------------------------------ */
|
|
/* round-for-reround digits */
|
|
/* ------------------------------------------------------------------ */
|
|
const uByte DECSTICKYTAB[10]={1,1,2,3,4,6,6,7,8,9}; /* used if sticky */
|
|
|
|
/* ------------------------------------------------------------------ */
|
|
/* Powers of ten (powers[n]==10**n, 0<=n<=9) */
|
|
/* ------------------------------------------------------------------ */
|
|
const uLong DECPOWERS[19] = {1, 10, 100, 1000, 10000, 100000, 1000000,
|
|
10000000, 100000000, 1000000000, 10000000000ULL, 100000000000ULL,
|
|
1000000000000ULL, 10000000000000ULL, 100000000000000ULL, 1000000000000000ULL,
|
|
10000000000000000ULL, 100000000000000000ULL, 1000000000000000000ULL, };
|
|
|
|
/* ------------------------------------------------------------------ */
|
|
/* decContextClearStatus -- clear bits in current status */
|
|
/* */
|
|
/* context is the context structure to be queried */
|
|
/* mask indicates the bits to be cleared (the status bit that */
|
|
/* corresponds to each 1 bit in the mask is cleared) */
|
|
/* returns context */
|
|
/* */
|
|
/* No error is possible. */
|
|
/* ------------------------------------------------------------------ */
|
|
decContext *decContextClearStatus(decContext *context, uInt mask) {
|
|
context->status&=~mask;
|
|
return context;
|
|
} /* decContextClearStatus */
|
|
|
|
/* ------------------------------------------------------------------ */
|
|
/* decContextDefault -- initialize a context structure */
|
|
/* */
|
|
/* context is the structure to be initialized */
|
|
/* kind selects the required set of default values, one of: */
|
|
/* DEC_INIT_BASE -- select ANSI X3-274 defaults */
|
|
/* DEC_INIT_DECIMAL32 -- select IEEE 754r defaults, 32-bit */
|
|
/* DEC_INIT_DECIMAL64 -- select IEEE 754r defaults, 64-bit */
|
|
/* DEC_INIT_DECIMAL128 -- select IEEE 754r defaults, 128-bit */
|
|
/* For any other value a valid context is returned, but with */
|
|
/* Invalid_operation set in the status field. */
|
|
/* returns a context structure with the appropriate initial values. */
|
|
/* ------------------------------------------------------------------ */
|
|
decContext * decContextDefault(decContext *context, Int kind) {
|
|
/* set defaults... */
|
|
context->digits=9; /* 9 digits */
|
|
context->emax=DEC_MAX_EMAX; /* 9-digit exponents */
|
|
context->emin=DEC_MIN_EMIN; /* .. balanced */
|
|
context->round=DEC_ROUND_HALF_UP; /* 0.5 rises */
|
|
context->traps=DEC_Errors; /* all but informational */
|
|
context->status=0; /* cleared */
|
|
context->clamp=0; /* no clamping */
|
|
#if DECSUBSET
|
|
context->extended=0; /* cleared */
|
|
#endif
|
|
switch (kind) {
|
|
case DEC_INIT_BASE:
|
|
/* [use defaults] */
|
|
break;
|
|
case DEC_INIT_DECIMAL32:
|
|
context->digits=7; /* digits */
|
|
context->emax=96; /* Emax */
|
|
context->emin=-95; /* Emin */
|
|
context->round=DEC_ROUND_HALF_EVEN; /* 0.5 to nearest even */
|
|
context->traps=0; /* no traps set */
|
|
context->clamp=1; /* clamp exponents */
|
|
#if DECSUBSET
|
|
context->extended=1; /* set */
|
|
#endif
|
|
break;
|
|
case DEC_INIT_DECIMAL64:
|
|
context->digits=16; /* digits */
|
|
context->emax=384; /* Emax */
|
|
context->emin=-383; /* Emin */
|
|
context->round=DEC_ROUND_HALF_EVEN; /* 0.5 to nearest even */
|
|
context->traps=0; /* no traps set */
|
|
context->clamp=1; /* clamp exponents */
|
|
#if DECSUBSET
|
|
context->extended=1; /* set */
|
|
#endif
|
|
break;
|
|
case DEC_INIT_DECIMAL128:
|
|
context->digits=34; /* digits */
|
|
context->emax=6144; /* Emax */
|
|
context->emin=-6143; /* Emin */
|
|
context->round=DEC_ROUND_HALF_EVEN; /* 0.5 to nearest even */
|
|
context->traps=0; /* no traps set */
|
|
context->clamp=1; /* clamp exponents */
|
|
#if DECSUBSET
|
|
context->extended=1; /* set */
|
|
#endif
|
|
break;
|
|
|
|
default: /* invalid Kind */
|
|
/* use defaults, and .. */
|
|
decContextSetStatus(context, DEC_Invalid_operation); /* trap */
|
|
}
|
|
|
|
#if DECCHECK
|
|
if (LITEND!=DECLITEND) {
|
|
const char *adj;
|
|
if (LITEND) adj="little";
|
|
else adj="big";
|
|
printf("Warning: DECLITEND is set to %d, but this computer appears to be %s-endian\n",
|
|
DECLITEND, adj);
|
|
}
|
|
#endif
|
|
return context;} /* decContextDefault */
|
|
|
|
/* ------------------------------------------------------------------ */
|
|
/* decContextGetRounding -- return current rounding mode */
|
|
/* */
|
|
/* context is the context structure to be queried */
|
|
/* returns the rounding mode */
|
|
/* */
|
|
/* No error is possible. */
|
|
/* ------------------------------------------------------------------ */
|
|
enum rounding decContextGetRounding(decContext *context) {
|
|
return context->round;
|
|
} /* decContextGetRounding */
|
|
|
|
/* ------------------------------------------------------------------ */
|
|
/* decContextGetStatus -- return current status */
|
|
/* */
|
|
/* context is the context structure to be queried */
|
|
/* returns status */
|
|
/* */
|
|
/* No error is possible. */
|
|
/* ------------------------------------------------------------------ */
|
|
uInt decContextGetStatus(decContext *context) {
|
|
return context->status;
|
|
} /* decContextGetStatus */
|
|
|
|
/* ------------------------------------------------------------------ */
|
|
/* decContextRestoreStatus -- restore bits in current status */
|
|
/* */
|
|
/* context is the context structure to be updated */
|
|
/* newstatus is the source for the bits to be restored */
|
|
/* mask indicates the bits to be restored (the status bit that */
|
|
/* corresponds to each 1 bit in the mask is set to the value of */
|
|
/* the corresponding bit in newstatus) */
|
|
/* returns context */
|
|
/* */
|
|
/* No error is possible. */
|
|
/* ------------------------------------------------------------------ */
|
|
decContext *decContextRestoreStatus(decContext *context,
|
|
uInt newstatus, uInt mask) {
|
|
context->status&=~mask; /* clear the selected bits */
|
|
context->status|=(mask&newstatus); /* or in the new bits */
|
|
return context;
|
|
} /* decContextRestoreStatus */
|
|
|
|
/* ------------------------------------------------------------------ */
|
|
/* decContextSaveStatus -- save bits in current status */
|
|
/* */
|
|
/* context is the context structure to be queried */
|
|
/* mask indicates the bits to be saved (the status bits that */
|
|
/* correspond to each 1 bit in the mask are saved) */
|
|
/* returns the AND of the mask and the current status */
|
|
/* */
|
|
/* No error is possible. */
|
|
/* ------------------------------------------------------------------ */
|
|
uInt decContextSaveStatus(decContext *context, uInt mask) {
|
|
return context->status&mask;
|
|
} /* decContextSaveStatus */
|
|
|
|
/* ------------------------------------------------------------------ */
|
|
/* decContextSetRounding -- set current rounding mode */
|
|
/* */
|
|
/* context is the context structure to be updated */
|
|
/* newround is the value which will replace the current mode */
|
|
/* returns context */
|
|
/* */
|
|
/* No error is possible. */
|
|
/* ------------------------------------------------------------------ */
|
|
decContext *decContextSetRounding(decContext *context,
|
|
enum rounding newround) {
|
|
context->round=newround;
|
|
return context;
|
|
} /* decContextSetRounding */
|
|
|
|
/* ------------------------------------------------------------------ */
|
|
/* decContextSetStatus -- set status and raise trap if appropriate */
|
|
/* */
|
|
/* context is the context structure to be updated */
|
|
/* status is the DEC_ exception code */
|
|
/* returns the context structure */
|
|
/* */
|
|
/* Control may never return from this routine, if there is a signal */
|
|
/* handler and it takes a long jump. */
|
|
/* ------------------------------------------------------------------ */
|
|
decContext * decContextSetStatus(decContext *context, uInt status) {
|
|
context->status|=status;
|
|
if (status & context->traps) raise(SIGFPE);
|
|
return context;} /* decContextSetStatus */
|
|
|
|
/* ------------------------------------------------------------------ */
|
|
/* decContextSetStatusFromString -- set status from a string + trap */
|
|
/* */
|
|
/* context is the context structure to be updated */
|
|
/* string is a string exactly equal to one that might be returned */
|
|
/* by decContextStatusToString */
|
|
/* */
|
|
/* The status bit corresponding to the string is set, and a trap */
|
|
/* is raised if appropriate. */
|
|
/* */
|
|
/* returns the context structure, unless the string is equal to */
|
|
/* DEC_Condition_MU or is not recognized. In these cases NULL is */
|
|
/* returned. */
|
|
/* ------------------------------------------------------------------ */
|
|
decContext * decContextSetStatusFromString(decContext *context,
|
|
const char *string) {
|
|
if (strcmp(string, DEC_Condition_CS)==0)
|
|
return decContextSetStatus(context, DEC_Conversion_syntax);
|
|
if (strcmp(string, DEC_Condition_DZ)==0)
|
|
return decContextSetStatus(context, DEC_Division_by_zero);
|
|
if (strcmp(string, DEC_Condition_DI)==0)
|
|
return decContextSetStatus(context, DEC_Division_impossible);
|
|
if (strcmp(string, DEC_Condition_DU)==0)
|
|
return decContextSetStatus(context, DEC_Division_undefined);
|
|
if (strcmp(string, DEC_Condition_IE)==0)
|
|
return decContextSetStatus(context, DEC_Inexact);
|
|
if (strcmp(string, DEC_Condition_IS)==0)
|
|
return decContextSetStatus(context, DEC_Insufficient_storage);
|
|
if (strcmp(string, DEC_Condition_IC)==0)
|
|
return decContextSetStatus(context, DEC_Invalid_context);
|
|
if (strcmp(string, DEC_Condition_IO)==0)
|
|
return decContextSetStatus(context, DEC_Invalid_operation);
|
|
#if DECSUBSET
|
|
if (strcmp(string, DEC_Condition_LD)==0)
|
|
return decContextSetStatus(context, DEC_Lost_digits);
|
|
#endif
|
|
if (strcmp(string, DEC_Condition_OV)==0)
|
|
return decContextSetStatus(context, DEC_Overflow);
|
|
if (strcmp(string, DEC_Condition_PA)==0)
|
|
return decContextSetStatus(context, DEC_Clamped);
|
|
if (strcmp(string, DEC_Condition_RO)==0)
|
|
return decContextSetStatus(context, DEC_Rounded);
|
|
if (strcmp(string, DEC_Condition_SU)==0)
|
|
return decContextSetStatus(context, DEC_Subnormal);
|
|
if (strcmp(string, DEC_Condition_UN)==0)
|
|
return decContextSetStatus(context, DEC_Underflow);
|
|
if (strcmp(string, DEC_Condition_ZE)==0)
|
|
return context;
|
|
return NULL; /* Multiple status, or unknown */
|
|
} /* decContextSetStatusFromString */
|
|
|
|
/* ------------------------------------------------------------------ */
|
|
/* decContextSetStatusFromStringQuiet -- set status from a string */
|
|
/* */
|
|
/* context is the context structure to be updated */
|
|
/* string is a string exactly equal to one that might be returned */
|
|
/* by decContextStatusToString */
|
|
/* */
|
|
/* The status bit corresponding to the string is set; no trap is */
|
|
/* raised. */
|
|
/* */
|
|
/* returns the context structure, unless the string is equal to */
|
|
/* DEC_Condition_MU or is not recognized. In these cases NULL is */
|
|
/* returned. */
|
|
/* ------------------------------------------------------------------ */
|
|
decContext * decContextSetStatusFromStringQuiet(decContext *context,
|
|
const char *string) {
|
|
if (strcmp(string, DEC_Condition_CS)==0)
|
|
return decContextSetStatusQuiet(context, DEC_Conversion_syntax);
|
|
if (strcmp(string, DEC_Condition_DZ)==0)
|
|
return decContextSetStatusQuiet(context, DEC_Division_by_zero);
|
|
if (strcmp(string, DEC_Condition_DI)==0)
|
|
return decContextSetStatusQuiet(context, DEC_Division_impossible);
|
|
if (strcmp(string, DEC_Condition_DU)==0)
|
|
return decContextSetStatusQuiet(context, DEC_Division_undefined);
|
|
if (strcmp(string, DEC_Condition_IE)==0)
|
|
return decContextSetStatusQuiet(context, DEC_Inexact);
|
|
if (strcmp(string, DEC_Condition_IS)==0)
|
|
return decContextSetStatusQuiet(context, DEC_Insufficient_storage);
|
|
if (strcmp(string, DEC_Condition_IC)==0)
|
|
return decContextSetStatusQuiet(context, DEC_Invalid_context);
|
|
if (strcmp(string, DEC_Condition_IO)==0)
|
|
return decContextSetStatusQuiet(context, DEC_Invalid_operation);
|
|
#if DECSUBSET
|
|
if (strcmp(string, DEC_Condition_LD)==0)
|
|
return decContextSetStatusQuiet(context, DEC_Lost_digits);
|
|
#endif
|
|
if (strcmp(string, DEC_Condition_OV)==0)
|
|
return decContextSetStatusQuiet(context, DEC_Overflow);
|
|
if (strcmp(string, DEC_Condition_PA)==0)
|
|
return decContextSetStatusQuiet(context, DEC_Clamped);
|
|
if (strcmp(string, DEC_Condition_RO)==0)
|
|
return decContextSetStatusQuiet(context, DEC_Rounded);
|
|
if (strcmp(string, DEC_Condition_SU)==0)
|
|
return decContextSetStatusQuiet(context, DEC_Subnormal);
|
|
if (strcmp(string, DEC_Condition_UN)==0)
|
|
return decContextSetStatusQuiet(context, DEC_Underflow);
|
|
if (strcmp(string, DEC_Condition_ZE)==0)
|
|
return context;
|
|
return NULL; /* Multiple status, or unknown */
|
|
} /* decContextSetStatusFromStringQuiet */
|
|
|
|
/* ------------------------------------------------------------------ */
|
|
/* decContextSetStatusQuiet -- set status without trap */
|
|
/* */
|
|
/* context is the context structure to be updated */
|
|
/* status is the DEC_ exception code */
|
|
/* returns the context structure */
|
|
/* */
|
|
/* No error is possible. */
|
|
/* ------------------------------------------------------------------ */
|
|
decContext * decContextSetStatusQuiet(decContext *context, uInt status) {
|
|
context->status|=status;
|
|
return context;} /* decContextSetStatusQuiet */
|
|
|
|
/* ------------------------------------------------------------------ */
|
|
/* decContextStatusToString -- convert status flags to a string */
|
|
/* */
|
|
/* context is a context with valid status field */
|
|
/* */
|
|
/* returns a constant string describing the condition. If multiple */
|
|
/* (or no) flags are set, a generic constant message is returned. */
|
|
/* ------------------------------------------------------------------ */
|
|
const char *decContextStatusToString(const decContext *context) {
|
|
Int status=context->status;
|
|
|
|
/* test the five IEEE first, as some of the others are ambiguous when */
|
|
/* DECEXTFLAG=0 */
|
|
if (status==DEC_Invalid_operation ) return DEC_Condition_IO;
|
|
if (status==DEC_Division_by_zero ) return DEC_Condition_DZ;
|
|
if (status==DEC_Overflow ) return DEC_Condition_OV;
|
|
if (status==DEC_Underflow ) return DEC_Condition_UN;
|
|
if (status==DEC_Inexact ) return DEC_Condition_IE;
|
|
|
|
if (status==DEC_Division_impossible ) return DEC_Condition_DI;
|
|
if (status==DEC_Division_undefined ) return DEC_Condition_DU;
|
|
if (status==DEC_Rounded ) return DEC_Condition_RO;
|
|
if (status==DEC_Clamped ) return DEC_Condition_PA;
|
|
if (status==DEC_Subnormal ) return DEC_Condition_SU;
|
|
if (status==DEC_Conversion_syntax ) return DEC_Condition_CS;
|
|
if (status==DEC_Insufficient_storage ) return DEC_Condition_IS;
|
|
if (status==DEC_Invalid_context ) return DEC_Condition_IC;
|
|
#if DECSUBSET
|
|
if (status==DEC_Lost_digits ) return DEC_Condition_LD;
|
|
#endif
|
|
if (status==0 ) return DEC_Condition_ZE;
|
|
return DEC_Condition_MU; /* Multiple errors */
|
|
} /* decContextStatusToString */
|
|
|
|
/* ------------------------------------------------------------------ */
|
|
/* decContextTestSavedStatus -- test bits in saved status */
|
|
/* */
|
|
/* oldstatus is the status word to be tested */
|
|
/* mask indicates the bits to be tested (the oldstatus bits that */
|
|
/* correspond to each 1 bit in the mask are tested) */
|
|
/* returns 1 if any of the tested bits are 1, or 0 otherwise */
|
|
/* */
|
|
/* No error is possible. */
|
|
/* ------------------------------------------------------------------ */
|
|
uInt decContextTestSavedStatus(uInt oldstatus, uInt mask) {
|
|
return (oldstatus&mask)!=0;
|
|
} /* decContextTestSavedStatus */
|
|
|
|
/* ------------------------------------------------------------------ */
|
|
/* decContextTestStatus -- test bits in current status */
|
|
/* */
|
|
/* context is the context structure to be updated */
|
|
/* mask indicates the bits to be tested (the status bits that */
|
|
/* correspond to each 1 bit in the mask are tested) */
|
|
/* returns 1 if any of the tested bits are 1, or 0 otherwise */
|
|
/* */
|
|
/* No error is possible. */
|
|
/* ------------------------------------------------------------------ */
|
|
uInt decContextTestStatus(decContext *context, uInt mask) {
|
|
return (context->status&mask)!=0;
|
|
} /* decContextTestStatus */
|
|
|
|
/* ------------------------------------------------------------------ */
|
|
/* decContextZeroStatus -- clear all status bits */
|
|
/* */
|
|
/* context is the context structure to be updated */
|
|
/* returns context */
|
|
/* */
|
|
/* No error is possible. */
|
|
/* ------------------------------------------------------------------ */
|
|
decContext *decContextZeroStatus(decContext *context) {
|
|
context->status=0;
|
|
return context;
|
|
} /* decContextZeroStatus */
|