diff options
author | John Hauser <jhauser@eecs.berkeley.edu> | 2016-07-22 18:03:04 -0700 |
---|---|---|
committer | John Hauser <jhauser@eecs.berkeley.edu> | 2016-07-22 18:03:04 -0700 |
commit | cb5087cd7403acf31ac24ac4be8e019a51904895 (patch) | |
tree | 3eeb55d6ad63e33dc8e3be33614e94bbe8a8cac5 /source/s_remStepMBy32.c | |
parent | 45fdcf1c6583e4af380b147ac568f5aa721b7ba8 (diff) | |
download | berkeley-softfloat-3-cb5087cd7403acf31ac24ac4be8e019a51904895.zip berkeley-softfloat-3-cb5087cd7403acf31ac24ac4be8e019a51904895.tar.gz berkeley-softfloat-3-cb5087cd7403acf31ac24ac4be8e019a51904895.tar.bz2 |
Release 3b. See "doc/SoftFloat-history.html".
Diffstat (limited to 'source/s_remStepMBy32.c')
-rw-r--r-- | source/s_remStepMBy32.c | 17 |
1 files changed, 8 insertions, 9 deletions
diff --git a/source/s_remStepMBy32.c b/source/s_remStepMBy32.c index 429d3dc..2849d94 100644 --- a/source/s_remStepMBy32.c +++ b/source/s_remStepMBy32.c @@ -2,9 +2,9 @@ /*============================================================================ This C source file is part of the SoftFloat IEEE Floating-Point Arithmetic -Package, Release 3a, by John R. Hauser. +Package, Release 3b, by John R. Hauser. -Copyright 2011, 2012, 2013, 2014, 2015 The Regents of the University of +Copyright 2011, 2012, 2013, 2014, 2015, 2016 The Regents of the University of California. All rights reserved. Redistribution and use in source and binary forms, with or without @@ -44,34 +44,33 @@ void softfloat_remStepMBy32( uint_fast8_t size_words, const uint32_t *remPtr, - uint_fast8_t count, + uint_fast8_t dist, const uint32_t *bPtr, uint32_t q, uint32_t *zPtr ) { - uint_fast8_t negCount; unsigned int index, lastIndex; uint64_t dwordProd; uint32_t wordRem, wordShiftedRem, wordProd; - uint_fast8_t borrow; + uint_fast8_t uNegDist, borrow; - negCount = -count; index = indexWordLo( size_words ); lastIndex = indexWordHi( size_words ); dwordProd = (uint64_t) bPtr[index] * q; wordRem = remPtr[index]; - wordShiftedRem = wordRem<<count; + wordShiftedRem = wordRem<<dist; wordProd = dwordProd; zPtr[index] = wordShiftedRem - wordProd; if ( index != lastIndex ) { + uNegDist = -dist; borrow = (wordShiftedRem < wordProd); for (;;) { - wordShiftedRem = wordRem>>(negCount & 31); + wordShiftedRem = wordRem>>(uNegDist & 31); index += wordIncr; dwordProd = (uint64_t) bPtr[index] * q + (dwordProd>>32); wordRem = remPtr[index]; - wordShiftedRem |= wordRem<<count; + wordShiftedRem |= wordRem<<dist; wordProd = dwordProd; zPtr[index] = wordShiftedRem - wordProd - borrow; if ( index == lastIndex ) break; |