monomial.h

Go to the documentation of this file.

// Copyright (C) 2009 Sebastian Pancratz                                     //
//                                                                           //
// Distributed under the terms of the GNU General Public License as          //
// published by the Free Software Foundation; either version 2 of the        //
// License, or (at your option) any later version.                           //
//                                                                           //
// http://www.gnu.org/licenses/                                              //

#ifndef _MONOMIAL_H_
#define _MONOMIAL_H_

#include <inttypes.h>


// Data types                                                                //

typedef uint64_t monomial_t;

typedef uint8_t exponent_t;

// Object-like macros                                                        //

#define MON_BITS_PER_EXP  8

#define MON_BITMASK_BLOCK  0x00000000000000FFULL

#define MON_BITMASK_ALL  0xFFFFFFFFFFFFFFFFULL

#define MON_VARS(n)  (n == 1 ? "XY" : n == 2 ? "XYZ" : n == 3 ? "WXYZ" : n == 4 ? "VWXYZ" : n == 5 ? "UVWXYZ" : n == 6 ? "TUVWXYZ" : "STUVWXYZ")

// Function-like macros                                                      //

// Initialization

#define MON_INIT(x)  ((x) = 0)

#define MON_CLEAR(x)

#define MON_SET(x, y)  ((x) = (y))

#define MON_SWAP(x, y)  { monomial_t _t_ = (x); (x) = (y); (y) = _t_; }

#define MON_ONE(x)  ((x) = 0)


// Access

#define MON_GET_EXP(x, i)  ((exponent_t) (((x) >> ((i) * MON_BITS_PER_EXP)) & MON_BITMASK_BLOCK))

#define MON_SET_EXP(x, i, e)  ((x) = (((x) & (MON_BITMASK_ALL - (MON_BITMASK_BLOCK << ((i) * MON_BITS_PER_EXP)))) | (((monomial_t)(e)) << ((i) * MON_BITS_PER_EXP))))

#define MON_INC_EXP(x, i, e)  ((x) = ((x) + (((monomial_t)(e)) << ((i) * MON_BITS_PER_EXP))))

#define MON_DEC_EXP(x, i, e)  ((x) = ((x) - (((monomial_t)(e)) << ((i) * MON_BITS_PER_EXP))))


// Comparison

#define MON_CMP_INVLEX(x, y)  (((x) < (y)) ? -1 : ((x) > (y) ? 1 : 0))

#define MON_IS_ONE(x)  ((x) == 0)


// Multiplication and division

#define MON_MUL(x, y, z)  ((x) = (y) + (z))

#define MON_DIV(x, y, z)  ((x) = (y) - (z))


// Function declarations                                                     //

uint32_t mon_degree(monomial_t op);

int32_t mon_divides(monomial_t x, monomial_t y);

char * mon_to_string(monomial_t x, const uint8_t n);
char * mon_to_string_pretty(monomial_t x, const uint8_t n, const char * vars);
monomial_t mon_from_string(char * str);

monomial_t * mon_generate_by_degree(uint32_t * len, const uint8_t n, 
                                    const uint32_t d);
monomial_t * mon_generate_by_degree_invlex(uint32_t * len, const uint8_t n, 
                                           const uint32_t d);

#endif

---