s_remStepMBy32.c File Reference

#include <stdint.h>
#include "platform.h"
#include "primitiveTypes.h"

Include dependency graph for s_remStepMBy32.c:

Go to the source code of this file.

Functions
void	softfloat_remStepMBy32 (uint_fast8_t size_words, const uint32_t *remPtr, uint_fast8_t dist, const uint32_t *bPtr, uint32_t q, uint32_t *zPtr)

Function Documentation

softfloat_remStepMBy32()

void softfloat_remStepMBy32	(	uint_fast8_t	size_words,
		const uint32_t *	remPtr,
		uint_fast8_t	dist,
		const uint32_t *	bPtr,
		uint32_t	q,
		uint32_t *	zPtr )

Definition at line 44 of file s_remStepMBy32.c.

{
    unsigned int index, lastIndex;
    uint64_t dwordProd;
    uint32_t wordRem, wordShiftedRem, wordProd;
    uint_fast8_t uNegDist, borrow;
 
    index = indexWordLo( size_words );
    lastIndex = indexWordHi( size_words );
    dwordProd = (uint64_t) bPtr[index] * q;
    wordRem = remPtr[index];
    wordShiftedRem = wordRem<<dist;
    wordProd = dwordProd;
    zPtr[index] = wordShiftedRem - wordProd;
    if ( index != lastIndex ) {
        uNegDist = -dist;
        borrow = (wordShiftedRem < wordProd);
        for (;;) {
            wordShiftedRem = wordRem>>(uNegDist & 31);
            index += wordIncr;
            dwordProd = (uint64_t) bPtr[index] * q + (dwordProd>>32);
            wordRem = remPtr[index];
            wordShiftedRem |= wordRem<<dist;
            wordProd = dwordProd;
            zPtr[index] = wordShiftedRem - wordProd - borrow;
            if ( index == lastIndex ) break;
            borrow =
                borrow ? (wordShiftedRem <= wordProd)
                    : (wordShiftedRem < wordProd);
        }
    }
 
}