src/N.Basic.cpp

#include <symbolism/N.h>
#include <string>
#include <sstream>

using namespace std;

namespace symbolism {
namespace ring {

// newsize != size => newsize > 0
void N::reserve_small(size_t newsize)
{
    size = newsize;
    void* space = malloc((size + 1)  * sizeof(limb_t));
    data = reinterpret_cast<limb_t*>(space);
}

void N::reserve_large(size_t newsize)
{
    if(newsize)
    {
        size = newsize;
        void* space = realloc(data, (size + 1)  * sizeof(limb_t));
        data = reinterpret_cast<limb_t*>(space);
    }
    else
    {
        size = 0;
        free(data);
    }
}

void N::resize_small(size_t newsize)
{
    limb_t value = limb;
    size = newsize;
    void* space = malloc((size + 1)  * sizeof(limb_t));
    data = reinterpret_cast<limb_t*>(space);
    data[0] = value;
    for(size_t i = 1; i <= size; i++)
    {
        data[i] = 0;
    }
}

void N::resize_large(size_t newsize)
{
    size_t oldsize = size;
    size = newsize;
    void* space = realloc(data, (size + 1)  * sizeof(limb_t));
    data = reinterpret_cast<limb_t*>(space);
    for(size_t i = oldsize + 1; i <= size; i++)
    {
        data[i] = 0;
    }
}

void N::truncate_large(size_t newsize)
{
    // newsize < size > 0
    if(newsize)
    {
        size = newsize;
        void* space = realloc(data, (size + 1)  * sizeof(limb_t));
        data = reinterpret_cast<limb_t*>(space);
    }
    else
    {
        size = 0;
        limb_t value = data[0];
        free(data);
        limb = value;
    }
}

void N::normalize_large()
{
    size_t newsize = size;
    while(data[newsize] == 0 && --newsize);
    if(newsize)
    {
        size = newsize;
        void* space = realloc(data, (size + 1)  * sizeof(limb_t));
        data = reinterpret_cast<limb_t*>(space);
    }
    else
    {
        limb_t value = data[0];
        size = 0;
        free(data);
        limb = value;
    }
}

}}

---

Copyright © 2007-2008 Remco Bloemen.

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

Hosted by