bnum.h File Reference

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
struct	bnum

Macros
#define	EXP_BITS   11
#define	MANT_BITS   53
#define	LIMB_DIGITS   9 /* each DWORD stores up to 9 digits */
#define	LIMB_MAX   1000000000 /* 10^9 */
#define	BNUM_PREC64   128 /* data size needed to store 64-bit double */
#define	BNUM_PREC80   2048 /* data size needed to store 80-bit double */

Functions
static int	bnum_idx (struct bnum *b, int idx)
static BOOL	bnum_lshift (struct bnum *b, int shift)
static BOOL	bnum_rshift (struct bnum *b, int shift)
static void	bnum_mult (struct bnum *b, int mult)

Variables
static const int	p10s [] = { 1, 10, 100, 1000, 10000, 100000, 1000000, 10000000, 100000000, 1000000000 }

Macro Definition Documentation

BNUM_PREC64

#define BNUM_PREC64   128 /* data size needed to store 64-bit double */

Definition at line 30 of file bnum.h.

BNUM_PREC80

#define BNUM_PREC80   2048 /* data size needed to store 80-bit double */

Definition at line 31 of file bnum.h.

EXP_BITS

#define EXP_BITS   11

Definition at line 22 of file bnum.h.

LIMB_DIGITS

#define LIMB_DIGITS   9 /* each DWORD stores up to 9 digits */

Definition at line 27 of file bnum.h.

LIMB_MAX

#define LIMB_MAX   1000000000 /* 10^9 */

Definition at line 28 of file bnum.h.

MANT_BITS

#define MANT_BITS   53

Definition at line 23 of file bnum.h.

Function Documentation

bnum_idx()

static int bnum_idx ( struct bnum *  b, int  idx  )

inlinestatic

Definition at line 41 of file bnum.h.

{
    return idx & (b->size - 1);
}

Referenced by bnum_lshift(), bnum_mult(), bnum_rshift(), bnum_to_mant(), fpnum_parse_bnum(), and pf_output_fp().

bnum_lshift()

static BOOL bnum_lshift ( struct bnum *  b, int  shift  )

inlinestatic

Definition at line 47 of file bnum.h.

{
    DWORD rest = 0;
    ULONGLONG tmp;
    int i;
 
    /* The limbs number can change by up to 1 so shift <= 29 */
    assert(shift <= 29);
 
    for(i=b->b; i<b->e; i++) {
        tmp = ((ULONGLONG)b->data[bnum_idx(b, i)] << shift) + rest;
        rest = tmp / LIMB_MAX;
        b->data[bnum_idx(b, i)] = tmp % LIMB_MAX;
 
        if(i == b->b && !b->data[bnum_idx(b, i)])
            b->b++;
    }
 
    if(rest) {
        b->data[bnum_idx(b, b->e)] = rest;
        b->e++;
 
        if(bnum_idx(b, b->b) == bnum_idx(b, b->e)) {
            if(b->data[bnum_idx(b, b->b)]) b->data[bnum_idx(b, b->b+1)] |= 1;
            b->b++;
        }
        return TRUE;
    }
    return FALSE;
}

Referenced by fpnum_parse_bnum(), and pf_output_fp().

bnum_mult()

static void bnum_mult ( struct bnum *  b, int  mult  )

inlinestatic

Definition at line 109 of file bnum.h.

{
    DWORD rest = 0;
    ULONGLONG tmp;
    int i;
 
    assert(mult <= LIMB_MAX);
 
    for(i=b->b; i<b->e; i++) {
        tmp = ((ULONGLONG)b->data[bnum_idx(b, i)] * mult) + rest;
        rest = tmp / LIMB_MAX;
        b->data[bnum_idx(b, i)] = tmp % LIMB_MAX;
 
        if(i == b->b && !b->data[bnum_idx(b, i)])
            b->b++;
    }
 
    if(rest) {
        b->data[bnum_idx(b, b->e)] = rest;
        b->e++;
 
        if(bnum_idx(b, b->b) == bnum_idx(b, b->e)) {
            if(b->data[bnum_idx(b, b->b)]) b->data[bnum_idx(b, b->b+1)] |= 1;
            b->b++;
        }
    }
}

Referenced by fpnum_parse_bnum().

bnum_rshift()

static BOOL bnum_rshift ( struct bnum *  b, int  shift  )

inlinestatic

Definition at line 79 of file bnum.h.

{
    DWORD tmp, rest = 0;
    BOOL ret = FALSE;
    int i;
 
    /* Compute LIMB_MAX << shift without accuracy loss */
    assert(shift <= 9);
 
    for(i=b->e-1; i>=b->b; i--) {
        tmp = b->data[bnum_idx(b, i)] & ((1<<shift)-1);
        b->data[bnum_idx(b, i)] = (b->data[bnum_idx(b, i)] >> shift) + rest;
        rest = (LIMB_MAX >> shift) * tmp;
        if(i==b->e-1 && !b->data[bnum_idx(b, i)]) {
            b->e--;
            ret = TRUE;
        }
    }
 
    if(rest) {
        if(bnum_idx(b, b->b-1) == bnum_idx(b, b->e)) {
            if(rest) b->data[bnum_idx(b, b->b)] |= 1;
        } else {
            b->b--;
            b->data[bnum_idx(b, b->b)] = rest;
        }
    }
    return ret;
}

Referenced by fpnum_parse_bnum(), and pf_output_fp().

Variable Documentation

p10s

const int p10s[] = { 1, 10, 100, 1000, 10000, 100000, 1000000, 10000000, 100000000, 1000000000 }

static

Definition at line 25 of file bnum.h.

Referenced by fpnum_parse_bnum(), log10i(), and pf_output_fp().