src/N.Addition.cpp

#include <symbolism/N.h>

namespace symbolism {
namespace ring {

void N::add_one_large()
{
    limb_t carry = 1;
    for(size_t i=0; i <= size && carry; i++)
    {
        limb_t old = data[i];
        data[i] += carry;
        carry = (limb < old) ? 1 : 0;
    }
    add_carry_large(carry);
}

void N::add_large_small(const N& a, const limb_t b)
{
    reserve(a.size);
    N::limb_t carry = mpn_add_1(data, a.data, a.size + 1, b);
    add_carry_large(carry);
}

void N::add_large_large(const N& a, const N& b)
{
    reserve(max(a.size, b.size));
    N::limb_t carry;
    if(a.size > b.size)
    {
        carry = mpn_add(data, a.data, a.size + 1, b.data, b.size + 1);
    }
    else if(a.size == b.size)
    {
        carry = mpn_add_n(data, a.data, b.data, a.size + 1);
    }
    else
    {
        carry = mpn_add(data, b.data, b.size + 1, a.data, a.size + 1);
    }
    add_carry_large(carry);
}

void N::add_large_small(const limb_t a)
{
    limb_t carry;
    carry = mpn_add_1(data, data, size + 1, a);
    add_carry_large(carry);
}

void N::add_small_large(const N& a)
{
    limb_t old = limb;
    limb_t carry;
    reserve(a.size);
    carry = mpn_add_1(data, a.data, a.size + 1, old);
    add_carry_large(carry);
}

void N::add_large_large(const N& a)
{
    if(size < a.size)
    {
        resize(a.size);
    }
    if(size == a.size)
    {
        limb_t carry = mpn_add_n(data, data, a.data, size + 1);
        add_carry_large(carry);
    }
    else
    {
        limb_t carry = mpn_add(data, data, size + 1, a.data, a.size + 1);
        add_carry_large(carry);
    }
}

}}

---

Copyright © 2007-2008 Remco Bloemen.

Generated on Tue Jan 22 17:35:31 2008 for symbolism by doxygen 1.5.4

Hosted by