rintegers.cc File Reference

#include <misc/auxiliary.h>
#include <omalloc/omalloc.h>
#include <factory/factory.h>
#include <misc/mylimits.h>
#include <reporter/reporter.h>
#include "coeffs.h"
#include "numbers.h"
#include "si_gmp.h"
#include "mpr_complex.h"
#include "rintegers.h"
#include "rmodulon.h"
#include "longrat.h"
#include <string.h>

Go to the source code of this file.

Functions
number	nrzCopy (number a, const coeffs r)
int	nrzSize (number a, const coeffs r)
void	nrzDelete (number *a, const coeffs r)
BOOLEAN	nrzGreaterZero (number k, const coeffs r)
number	nrzMult (number a, number b, const coeffs r)
number	nrzInit (long i, const coeffs r)
long	nrzInt (number &n, const coeffs r)
number	nrzAdd (number a, number b, const coeffs r)
number	nrzSub (number a, number b, const coeffs r)
void	nrzPower (number a, int i, number *result, const coeffs r)
BOOLEAN	nrzIsZero (number a, const coeffs r)
BOOLEAN	nrzIsOne (number a, const coeffs r)
BOOLEAN	nrzIsMOne (number a, const coeffs r)
BOOLEAN	nrzIsUnit (number a, const coeffs r)
number	nrzGetUnit (number a, const coeffs r)
number	nrzDiv (number a, number b, const coeffs r)
number	nrzExactDiv (number a, number b, const coeffs r)
number	nrzIntMod (number a, number b, const coeffs r)
number	nrzNeg (number c, const coeffs r)
number	nrzInvers (number c, const coeffs r)
BOOLEAN	nrzGreater (number a, number b, const coeffs r)
BOOLEAN	nrzDivBy (number a, number b, const coeffs r)
int	nrzDivComp (number a, number b, const coeffs r)
BOOLEAN	nrzEqual (number a, number b, const coeffs r)
number	nrzLcm (number a, number b, const coeffs r)
number	nrzGcd (number a, number b, const coeffs r)
number	nrzExtGcd (number a, number b, number *s, number *t, const coeffs r)
nMapFunc	nrzSetMap (const coeffs src, const coeffs dst)
void	nrzWrite (number &a, const coeffs r)
const char *	nrzRead (const char *s, number *a, const coeffs r)
char *	nrzName (number n, const coeffs r)
void	nrzCoeffWrite (const coeffs r, BOOLEAN details)
BOOLEAN	nrzDBTest (number a, const char *f, const int l, const coeffs r)
void	nrzSetExp (int c, coeffs r)
void	nrzInitExp (int c, coeffs r)
coeffs	nrzQuot1 (number c, const coeffs r)
number	nrzMapQ (number from, const coeffs src, const coeffs dst)
number	nrzMapMachineInt (number from, const coeffs, const coeffs)
number	nrzMapZp (number from, const coeffs, const coeffs)
static const char *	nlEatLongC (char *s, mpz_ptr i)
static CanonicalForm	nrzConvSingNFactoryN (number n, BOOLEAN setChar, const coeffs)
static number	nrzConvFactoryNSingN (const CanonicalForm n, const coeffs r)
static char *	nrzCoeffString (const coeffs)
BOOLEAN	nrzInitChar (coeffs r, void *)

Variables
static const n_coeffType	ID = n_Z
	Our Type! More...
omBin	gmp_nrz_bin = omGetSpecBin(sizeof(mpz_t))

Function Documentation

static const char* nlEatLongC ( char *  s, mpz_ptr  i  )

static

Definition at line 400 of file rintegers.cc.

{
  const char * start=s;

  if (*s<'0' || *s>'9')
  {
    mpz_set_si(i,1);
    return s;
  }
  while (*s >= '0' && *s <= '9') s++;
  if (*s=='\0')
  {
    mpz_set_str(i,start,10);
  }
  else
  {
    char c=*s;
    *s='\0';
    mpz_set_str(i,start,10);
    *s=c;
  }
  return s;
}

number nrzAdd	(	number	a,
		number	b,
		const coeffs	r
	)

Definition at line 190 of file rintegers.cc.

{
  mpz_ptr erg = (mpz_ptr) omAllocBin(gmp_nrz_bin);
  mpz_init(erg);
  mpz_add(erg, (mpz_ptr) a, (mpz_ptr) b);
  return (number) erg;
}

static char* nrzCoeffString ( const coeffs  )

static

Definition at line 464 of file rintegers.cc.

{
  return omStrDup("integer");
}

void nrzCoeffWrite	(	const coeffs	r,
		BOOLEAN	details
	)

Definition at line 459 of file rintegers.cc.

{
  PrintS("//   coeff. ring is : Integers\n");
}

static number nrzConvFactoryNSingN ( const CanonicalForm  n, const coeffs  r  )

static

Definition at line 436 of file rintegers.cc.

{
  if (n.isImm())
    return nrzInit(n.intval(),r);
  else
  {
    mpz_ptr m = (mpz_ptr) omAllocBin(gmp_nrz_bin);
    gmp_numerator(n,m);
    return (number) m;
  }
}

static CanonicalForm nrzConvSingNFactoryN ( number  n, BOOLEAN  setChar, const coeffs    )

static

Definition at line 425 of file rintegers.cc.

{
  if (setChar) setCharacteristic( 0 );

  CanonicalForm term;
  mpz_t num;
  mpz_init_set(num, *((mpz_t*)n));
  term = make_cf(num);
  return term;
}

number nrzCopy	(	number	a,
		const coeffs	r
	)

Definition at line 161 of file rintegers.cc.

{
  if (a==NULL) return NULL;
  mpz_ptr erg = (mpz_ptr) omAllocBin(gmp_nrz_bin);
  mpz_init_set(erg, (mpz_ptr) a);
  return (number) erg;
}

BOOLEAN nrzDBTest	(	number	a,
		const char *	f,
		const int	l,
		const coeffs	r
	)

Definition at line 374 of file rintegers.cc.

{
  return TRUE;//TODO
}

void nrzDelete	(	number *	a,
		const coeffs	r
	)

Definition at line 153 of file rintegers.cc.

{
  if (*a == NULL) return;
  mpz_clear((mpz_ptr) *a);
  omFreeBin((ADDRESS) *a, gmp_nrz_bin);
  *a = NULL;
}

number nrzDiv	(	number	a,
		number	b,
		const coeffs	r
	)

Definition at line 262 of file rintegers.cc.

{
  mpz_ptr erg = (mpz_ptr) omAllocBin(gmp_nrz_bin);
  mpz_init(erg);
  mpz_ptr r = (mpz_ptr) omAllocBin(gmp_nrz_bin);
  mpz_init(r);
  mpz_tdiv_qr(erg, r, (mpz_ptr) a, (mpz_ptr) b);
  //if (!nrzIsZero((number) r, R))
  //{
  //  WerrorS("Division by non divisible element.");
  //  WerrorS("Result is without remainder.");
  //}
  mpz_clear(r);
  omFreeBin(r, gmp_nrz_bin);
  return (number) erg;
}

BOOLEAN nrzDivBy	(	number	a,
		number	b,
		const coeffs	r
	)

Definition at line 257 of file rintegers.cc.

{
  return mpz_divisible_p((mpz_ptr) a, (mpz_ptr) b) != 0;
}

int nrzDivComp	(	number	a,
		number	b,
		const coeffs	r
	)

Definition at line 246 of file rintegers.cc.

{
  if (nrzDivBy(a, b, r))
  {
    if (nrzDivBy(b, a, r)) return 2;
    return -1;
  }
  if (nrzDivBy(b, a, r)) return 1;
  return 0;
}

BOOLEAN nrzEqual	(	number	a,
		number	b,
		const coeffs	r
	)

Definition at line 231 of file rintegers.cc.

{
  return 0 == mpz_cmp((mpz_ptr) a, (mpz_ptr) b);
}

number nrzExactDiv	(	number	a,
		number	b,
		const coeffs	r
	)

Definition at line 279 of file rintegers.cc.

{
  mpz_ptr erg = (mpz_ptr) omAllocBin(gmp_nrz_bin);
  mpz_init(erg);
  mpz_tdiv_q(erg, (mpz_ptr) a, (mpz_ptr) b);
  return (number) erg;
}

number nrzExtGcd	(	number	a,
		number	b,
		number *	s,
		number *	t,
		const coeffs	r
	)

Definition at line 121 of file rintegers.cc.

{
  mpz_ptr erg = (mpz_ptr) omAllocBin(gmp_nrz_bin);
  mpz_ptr bs = (mpz_ptr) omAllocBin(gmp_nrz_bin);
  mpz_ptr bt = (mpz_ptr) omAllocBin(gmp_nrz_bin);
  mpz_init(erg);
  mpz_init(bs);
  mpz_init(bt);
  mpz_gcdext(erg, bs, bt, (mpz_ptr) a, (mpz_ptr) b);
  *s = (number) bs;
  *t = (number) bt;
  return (number) erg;
}

number nrzGcd	(	number	a,
		number	b,
		const coeffs	r
	)

Definition at line 109 of file rintegers.cc.

{
  mpz_ptr erg = (mpz_ptr) omAllocBin(gmp_nrz_bin);
  mpz_init(erg);
  mpz_gcd(erg, (mpz_ptr) a, (mpz_ptr) b);
  return (number) erg;
}

number nrzGetUnit	(	number	a,
		const coeffs	r
	)

Definition at line 206 of file rintegers.cc.

{
  return nrzInit(1, r);
}

BOOLEAN nrzGreater	(	number	a,
		number	b,
		const coeffs	r
	)

Definition at line 236 of file rintegers.cc.

{
  return 0 < mpz_cmp((mpz_ptr) a, (mpz_ptr) b);
}

BOOLEAN nrzGreaterZero	(	number	k,
		const coeffs	r
	)

Definition at line 241 of file rintegers.cc.

{
  return 0 < mpz_cmp_si((mpz_ptr) k, 0);
}

number nrzInit	(	long	i,
		const coeffs	r
	)

Definition at line 146 of file rintegers.cc.

{
  mpz_ptr erg = (mpz_ptr) omAllocBin(gmp_nrz_bin);
  mpz_init_set_si(erg, i);
  return (number) erg;
}

BOOLEAN nrzInitChar	(	coeffs	r,
		void *
	)

Definition at line 482 of file rintegers.cc.

{
  assume( getCoeffType(r) == ID );

  r->is_field=FALSE;
  r->is_domain=TRUE;
  r->rep=n_rep_gmp;

  //r->nCoeffIsEqual = ndCoeffIsEqual;
  r->cfCoeffString = nrzCoeffString;
  //r->cfKillChar = ndKillChar;
  r->cfMult  = nrzMult;
  r->cfSub   = nrzSub;
  r->cfAdd   = nrzAdd;
  r->cfDiv   = nrzDiv;
  r->cfIntMod= nrzIntMod;
  r->cfExactDiv= nrzExactDiv;
  r->cfInit = nrzInit;
  r->cfSize  = nrzSize;
  r->cfInt  = nrzInt;
  //#ifdef HAVE_RINGS
  r->cfDivComp = nrzDivComp; // only for ring stuff
  r->cfIsUnit = nrzIsUnit; // only for ring stuff
  r->cfGetUnit = nrzGetUnit; // only for ring stuff
  r->cfExtGcd = nrzExtGcd; // only for ring stuff
  r->cfDivBy = nrzDivBy; // only for ring stuff
  //#endif
  r->cfInpNeg   = nrzNeg;
  r->cfInvers= nrzInvers;
  r->cfCopy  = nrzCopy;
  r->cfWriteLong = nrzWrite;
  r->cfRead = nrzRead;
  r->cfGreater = nrzGreater;
  r->cfEqual = nrzEqual;
  r->cfIsZero = nrzIsZero;
  r->cfIsOne = nrzIsOne;
  r->cfIsMOne = nrzIsMOne;
  r->cfGreaterZero = nrzGreaterZero;
  r->cfPower = nrzPower;
  r->cfGcd  = nrzGcd;
  r->cfLcm  = nrzLcm;
  r->cfDelete= nrzDelete;
  r->cfSetMap = nrzSetMap;
  r->cfCoeffWrite = nrzCoeffWrite;
  r->cfQuot1 = nrzQuot1;
  r->convSingNFactoryN=nrzConvSingNFactoryN;
  r->convFactoryNSingN=nrzConvFactoryNSingN;
  // debug stuff

#ifdef LDEBUG
  r->cfDBTest=nrzDBTest;
#endif

  r->nNULL = 0;
  r->ch = 0;
  r->has_simple_Alloc=FALSE;
  r->has_simple_Inverse=FALSE;
  return FALSE;
}

void nrzInitExp	(	int	c,
		coeffs	r
	)

Definition at line 369 of file rintegers.cc.

{
}

long nrzInt	(	number &	n,
		const coeffs	r
	)

Definition at line 185 of file rintegers.cc.

{
  return mpz_get_si( (mpz_ptr)n);
}

number nrzIntMod	(	number	a,
		number	b,
		const coeffs	r
	)

Definition at line 287 of file rintegers.cc.

{
  mpz_ptr erg = (mpz_ptr) omAllocBin(gmp_nrz_bin);
  mpz_init(erg);
  mpz_ptr r = (mpz_ptr) omAllocBin(gmp_nrz_bin);
  mpz_init(r);
  mpz_tdiv_qr(erg, r, (mpz_ptr) a, (mpz_ptr) b);
  mpz_clear(erg);
  omFreeBin(erg, gmp_nrz_bin);
  return (number) r;
}

number nrzInvers	(	number	c,
		const coeffs	r
	)

Definition at line 299 of file rintegers.cc.

{
  if (!nrzIsUnit((number) c, r))
  {
    WerrorS("Non invertible element.");
    return (number)0; //TODO
  }
  return nrzCopy(c,r);
}

BOOLEAN nrzIsMOne	(	number	a,
		const coeffs	r
	)

Definition at line 226 of file rintegers.cc.

{
  return (a!=NULL) && (0 == mpz_cmp_si((mpz_ptr) a, -1));
}

BOOLEAN nrzIsOne	(	number	a,
		const coeffs	r
	)

Definition at line 221 of file rintegers.cc.

{
  return (a!=NULL) && (0 == mpz_cmp_si((mpz_ptr) a, 1));
}

BOOLEAN nrzIsUnit	(	number	a,
		const coeffs	r
	)

Definition at line 211 of file rintegers.cc.

{
  return 0 == mpz_cmpabs_ui((mpz_ptr) a, 1);
}

BOOLEAN nrzIsZero	(	number	a,
		const coeffs	r
	)

Definition at line 216 of file rintegers.cc.

{
  return 0 == mpz_cmpabs_ui((mpz_ptr) a, 0);
}

number nrzLcm	(	number	a,
		number	b,
		const coeffs	r
	)

Definition at line 97 of file rintegers.cc.

{
  mpz_ptr erg = (mpz_ptr) omAllocBin(gmp_nrz_bin);
  mpz_init(erg);
  mpz_lcm(erg, (mpz_ptr) a, (mpz_ptr) b);
  return (number) erg;
}

number nrzMapMachineInt	(	number	from,
		const coeffs	,
		const coeffs
	)

Definition at line 316 of file rintegers.cc.

{
  mpz_ptr erg = (mpz_ptr) omAllocBin(gmp_nrz_bin);
  mpz_init_set_ui(erg, (unsigned long) from);
  return (number) erg;
}

number nrzMapQ	(	number	from,
		const coeffs	src,
		const coeffs	dst
	)

Definition at line 330 of file rintegers.cc.

{
  mpz_ptr erg = (mpz_ptr) omAllocBin(gmp_nrz_bin);
  mpz_init(erg);
  nlGMP(from, (number) erg, src); // FIXME? TODO? // extern void   nlGMP(number &i, number n, const coeffs r); // to be replaced with n_MPZ(erg, from, src); // ?
  return (number) erg;
}

number nrzMapZp	(	number	from,
		const coeffs	,
		const coeffs
	)

Definition at line 323 of file rintegers.cc.

{
  mpz_ptr erg = (mpz_ptr) omAllocBin(gmp_nrz_bin);
  mpz_init_set_si(erg, (long) from);
  return (number) erg;
}

number nrzMult	(	number	a,
		number	b,
		const coeffs	r
	)

Definition at line 86 of file rintegers.cc.

{
  mpz_ptr erg = (mpz_ptr) omAllocBin(gmp_nrz_bin);
  mpz_init(erg);
  mpz_mul(erg, (mpz_ptr) a, (mpz_ptr) b);
  return (number) erg;
}

char* nrzName	(	number	n,
		const coeffs	r
	)

number nrzNeg	(	number	c,
		const coeffs	r
	)

Definition at line 309 of file rintegers.cc.

{
// nNeg inplace !!!
  mpz_mul_si((mpz_ptr) c, (mpz_ptr) c, -1);
  return c;
}

void nrzPower	(	number	a,
		int	i,
		number *	result,
		const coeffs	r
	)

Definition at line 135 of file rintegers.cc.

{
  mpz_ptr erg = (mpz_ptr) omAllocBin(gmp_nrz_bin);
  mpz_init(erg);
  mpz_pow_ui(erg, (mpz_ptr) a, i);
  *result = (number) erg;
}

coeffs nrzQuot1	(	number	c,
		const coeffs	r
	)

Definition at line 469 of file rintegers.cc.

{
    long ch = r->cfInt(c, r);
    mpz_ptr dummy;
    dummy = (mpz_ptr) omAlloc(sizeof(mpz_t));
    mpz_init_set_ui(dummy, ch);
    ZnmInfo info;
    info.base = dummy;
    info.exp = (unsigned long) 1;
    coeffs rr = nInitChar(n_Zn, (void*)&info);
    return(rr);
}

const char * nrzRead	(	const char *	s,
		number *	a,
		const coeffs	r
	)

Definition at line 448 of file rintegers.cc.

{
  mpz_ptr z = (mpz_ptr) omAllocBin(gmp_nrz_bin);
  {
    mpz_init(z);
    s = nlEatLongC((char *) s, z);
  }
  *a = (number) z;
  return s;
}

void nrzSetExp	(	int	c,
		coeffs	r
	)

Definition at line 365 of file rintegers.cc.

{
}

nMapFunc nrzSetMap	(	const coeffs	src,
		const coeffs	dst
	)

Definition at line 338 of file rintegers.cc.

{
  /* dst = currRing */
  if (nCoeff_is_Ring_Z(src) || nCoeff_is_Ring_ModN(src) || nCoeff_is_Ring_PtoM(src))
  {
    return ndCopyMap; //nrzCopyMap;
  }
  if (nCoeff_is_Ring_2toM(src))
  {
    return nrzMapMachineInt;
  }
  if (nCoeff_is_Zp(src))
  {
    return nrzMapZp;
  }
  if (getCoeffType(src)==n_Q /*nCoeff_is_Q(src) or coeffs_BIGINT*/)
  {
    return nrzMapQ;
  }
  return NULL;      // default
}

int nrzSize	(	number	a,
		const coeffs	r
	)

Definition at line 176 of file rintegers.cc.

{
  if (a == NULL) return 0;
  return ((mpz_ptr)a)->_mp_alloc;
}

number nrzSub	(	number	a,
		number	b,
		const coeffs	r
	)

Definition at line 198 of file rintegers.cc.

{
  mpz_ptr erg = (mpz_ptr) omAllocBin(gmp_nrz_bin);
  mpz_init(erg);
  mpz_sub(erg, (mpz_ptr) a, (mpz_ptr) b);
  return (number) erg;
}

void nrzWrite	(	number &	a,
		const coeffs	r
	)

Definition at line 380 of file rintegers.cc.

{
  char *s,*z;
  if (a==NULL)
  {
    StringAppendS("o");
  }
  else
  {
    int l=mpz_sizeinbase((mpz_ptr) a, 10) + 2;
    s=(char*)omAlloc(l);
    z=mpz_get_str(s,10,(mpz_ptr) a);
    StringAppendS(z);
    omFreeSize((ADDRESS)s,l);
  }
}

Variable Documentation

omBin gmp_nrz_bin = omGetSpecBin(sizeof(mpz_t))

Definition at line 80 of file rintegers.cc.

const n_coeffType ID = n_Z

static

Our Type!

Definition at line 32 of file rintegers.cc.