10 #if SI_INTEGER_VARIANT == 2
17 static number nrzMult (number a, number
b,
const coeffs)
21 mpz_mul(erg, (mpz_ptr) a, (mpz_ptr)
b);
28 static number nrzLcm (number a,number
b,
const coeffs)
32 mpz_lcm(erg, (mpz_ptr) a, (mpz_ptr)
b);
40 static number nrzGcd (number a,number
b,
const coeffs)
44 mpz_gcd(erg, (mpz_ptr) a, (mpz_ptr)
b);
52 static number nrzExtGcd (number a, number
b, number *
s, number *t,
const coeffs)
60 mpz_gcdext(erg, bs, bt, (mpz_ptr) a, (mpz_ptr)
b);
66 static number nrzXExtGcd (number a, number
b, number *
s, number *t, number *u, number *
v,
const coeffs )
75 mpz_gcdext(erg, bs, bt, (mpz_ptr)a, (mpz_ptr)
b);
80 mpz_init_set(bu, (mpz_ptr)
b);
81 mpz_init_set(bv, (mpz_ptr) a);
83 assume(mpz_cmp_si(erg, 0));
97 static void nrzPower (number a,
int i, number *
result,
const coeffs)
101 mpz_pow_ui(erg, (mpz_ptr) a,
i);
111 mpz_init_set_si(erg,
i);
119 mpz_clear((mpz_ptr) *a);
125 static number nrzCopy(number a,
const coeffs)
129 mpz_init_set(erg, (mpz_ptr) a);
134 number nrzCopyMap(number a,
const coeffs ,
const coeffs dst)
136 return nrzCopy(a,dst);
142 mpz_ptr
p=(mpz_ptr)a;
144 if (
s==1)
s=(mpz_cmp_ui(
p,0)!=0);
152 static long nrzInt(number &n,
const coeffs)
154 return mpz_get_si( (mpz_ptr)n);
157 static number nrzAdd (number a, number
b,
const coeffs)
161 mpz_add(erg, (mpz_ptr) a, (mpz_ptr)
b);
165 static number nrzSub (number a, number
b,
const coeffs)
169 mpz_sub(erg, (mpz_ptr) a, (mpz_ptr)
b);
173 static number nrzGetUnit (number,
const coeffs r)
180 return 0 == mpz_cmpabs_ui((mpz_ptr) a, 1);
185 return 0 == mpz_cmpabs_ui((mpz_ptr) a, 0);
190 return (0 == mpz_cmp_ui((mpz_ptr) a, 1));
195 return (0 == mpz_cmp_si((mpz_ptr) a, -1));
200 return 0 == mpz_cmp((mpz_ptr) a, (mpz_ptr)
b);
205 return 0 < mpz_cmp((mpz_ptr) a, (mpz_ptr)
b);
215 return mpz_divisible_p((mpz_ptr) a, (mpz_ptr)
b) != 0;
218 static int nrzDivComp(number a, number
b,
const coeffs r)
220 if (nrzDivBy(a,
b, r))
222 if (nrzDivBy(
b, a, r))
return 2;
225 if (nrzDivBy(
b, a, r))
return 1;
229 static number nrzDiv (number a,number
b,
const coeffs r)
241 mpz_tdiv_qr(erg, r, (mpz_ptr) a, (mpz_ptr)
b);
248 static number nrzExactDiv (number a,number
b,
const coeffs r)
258 mpz_tdiv_q(erg, (mpz_ptr) a, (mpz_ptr)
b);
263 static number nrzEucNorm (number a,
const coeffs )
267 mpz_abs(
abs, (mpz_ptr)a);
272 static number nrzSmallestQuotRem (number a, number
b, number * r,
const coeffs )
278 int gsign = mpz_sgn((mpz_ptr)
b);
282 mpz_abs(gg, (mpz_ptr)
b);
283 mpz_fdiv_qr(qq, rr, (mpz_ptr) a, gg);
284 mpz_tdiv_q_2exp(ghalf, gg, 1);
285 if (mpz_cmp(rr, ghalf) > 0)
288 mpz_add_ui(qq, qq, 1);
290 if (gsign < 0) mpz_neg(qq, qq);
306 static number nrzQuotRem (number a, number
b, number * r,
const coeffs )
312 mpz_tdiv_qr(qq, rr, (mpz_ptr) a, (mpz_ptr)
b);
325 static number nrzIntMod (number a,number
b,
const coeffs)
331 mpz_tdiv_qr(erg, r, (mpz_ptr) a, (mpz_ptr)
b);
337 static number nrzInvers (number c,
const coeffs r)
339 if (!nrzIsUnit((number) c, r))
341 WerrorS(
"Non invertible element.");
347 static number nrzNeg (number c,
const coeffs)
354 static number nrzMapMachineInt(number from,
const coeffs ,
const coeffs )
357 mpz_init_set_ui(erg, (
unsigned long) from);
361 static number nrzMapZp(number from,
const coeffs ,
const coeffs )
364 mpz_init_set_si(erg, (
long) from);
368 static number nrzMapQ(number from,
const coeffs src,
const coeffs )
372 nlMPZ(erg, from, src);
391 return nrzMapMachineInt;
408 void nrzSetExp(
int,
coeffs)
412 void nrzInitExp(
int,
coeffs)
417 static BOOLEAN nrzDBTest (number,
const char *,
const int,
const coeffs)
432 int l=mpz_sizeinbase((mpz_ptr) a, 10) + 2;
434 z=mpz_get_str(
s,10,(mpz_ptr) a);
443 static const char * nlEatLongC(
char *
s, mpz_ptr
i)
445 const char * start=
s;
447 if (*s<'0' || *s>
'9')
452 while (*
s >=
'0' && *
s <=
'9')
s++;
455 mpz_set_str(
i,start,10);
461 mpz_set_str(
i,start,10);
474 mpz_init_set(
num, *((mpz_t*)n));
489 WarnS(
"denominator is not 1 in factory");
495 static const char * nrzRead (
const char *
s, number *a,
const coeffs)
500 s = nlEatLongC((
char *)
s, z);
508 long ch = r->cfInt(c, r);
510 mpz_init_set_ui(dummy, ch);
513 info.exp = (
unsigned long) 1;
519 static number nrzInitMPZ(mpz_t
m,
const coeffs)
526 static void nrzMPZ(mpz_t
res, number &a,
const coeffs)
528 mpz_init_set(
res, (mpz_ptr) a);
531 static number nrzFarey(number r, number
N,
const coeffs R)
533 number a0 = nrzCopy(
N,
R);
535 number a1 = nrzCopy(r,
R);
539 PrintS(
"Farey start with ");
547 number as = nrzMult(a1, a1,
R);
549 if (nrzGreater(
N, as,
R))
555 number q = nrzDiv(a0, a1,
R);
556 number t = nrzMult(a1, q,
R),
557 s = nrzSub(a0, t,
R);
563 t = nrzMult(b1, q,
R);
564 s = nrzSub(b0, t,
R);
571 number as = nrzMult(b1, b1,
R);
574 if (nrzGreater(as,
N,
R))
631 r->cfIntMod= nrzIntMod;
632 r->cfExactDiv= nrzExactDiv;
634 r->cfInitMPZ = nrzInitMPZ;
638 r->cfDivComp = nrzDivComp;
639 r->cfIsUnit = nrzIsUnit;
640 r->cfGetUnit = nrzGetUnit;
641 r->cfExtGcd = nrzExtGcd;
642 r->cfXExtGcd = nrzXExtGcd;
643 r->cfDivBy = nrzDivBy;
644 r->cfEucNorm = nrzEucNorm;
645 r->cfQuotRem = nrzSmallestQuotRem;
646 r->cfInpNeg = nrzNeg;
647 r->cfInvers= nrzInvers;
651 r->cfGreater = nrzGreater;
652 r->cfEqual = nrzEqual;
653 r->cfIsZero = nrzIsZero;
654 r->cfIsOne = nrzIsOne;
655 r->cfIsMOne = nrzIsMOne;
656 r->cfGreaterZero = nrzGreaterZero;
657 r->cfPower = nrzPower;
661 r->cfSetMap = nrzSetMap;
662 r->cfQuot1 = nrzQuot1;
663 r->convSingNFactoryN=nrzConvSingNFactoryN;
664 r->convFactoryNSingN=nrzConvFactoryNSingN;
672 r->cfDBTest=nrzDBTest;
676 r->has_simple_Alloc=
FALSE;
677 r->has_simple_Inverse=
FALSE;
Rational abs(const Rational &a)
const CanonicalForm CFMap CFMap & N
static FORCE_INLINE BOOLEAN nCoeff_is_Z(const coeffs r)
number ndCopyMap(number a, const coeffs src, const coeffs dst)
static FORCE_INLINE BOOLEAN nCoeff_is_Ring_PtoM(const coeffs r)
@ n_Q
rational (GMP) numbers
@ n_Zn
only used if HAVE_RINGS is defined
@ n_Z
only used if HAVE_RINGS is defined
void n_Print(number &a, const coeffs r)
print a number (BEWARE of string buffers!) mostly for debugging
static FORCE_INLINE nMapFunc n_SetMap(const coeffs src, const coeffs dst)
set the mapping function pointers for translating numbers from src to dst
static FORCE_INLINE number n_Div(number a, number b, const coeffs r)
return the quotient of 'a' and 'b', i.e., a/b; raises an error if 'b' is not invertible in r exceptio...
coeffs nInitChar(n_coeffType t, void *parameter)
one-time initialisations for new coeffs in case of an error return NULL
static FORCE_INLINE n_coeffType getCoeffType(const coeffs r)
Returns the type of coeffs domain.
static FORCE_INLINE void n_Delete(number *p, const coeffs r)
delete 'p'
static FORCE_INLINE BOOLEAN nCoeff_is_Zn(const coeffs r)
static FORCE_INLINE BOOLEAN nCoeff_is_Zp(const coeffs r)
static FORCE_INLINE BOOLEAN nCoeff_is_Ring_2toM(const coeffs r)
static FORCE_INLINE void n_InpMult(number &a, number b, const coeffs r)
multiplication of 'a' and 'b'; replacement of 'a' by the product a*b
@ n_rep_gap_gmp
(), see rinteger.h, new impl.
@ n_rep_gmp
(mpz_ptr), see rmodulon,h
number(* nMapFunc)(number a, const coeffs src, const coeffs dst)
maps "a", which lives in src, into dst
void nKillChar(coeffs r)
undo all initialisations
const CanonicalForm int s
const Variable & v
< [in] a sqrfree bivariate poly
const ExtensionInfo & info
< [in] sqrfree poly
CanonicalForm FACTORY_PUBLIC make_cf(const mpz_ptr n)
void FACTORY_PUBLIC gmp_numerator(const CanonicalForm &f, mpz_ptr result)
void WerrorS(const char *s)
void mpz_mul_si(mpz_ptr r, mpz_srcptr s, long int si)
void nlMPZ(mpz_t m, number &n, const coeffs r)
number nlChineseRemainderSym(number *x, number *q, int rl, BOOLEAN sym, CFArray &inv_cache, const coeffs CF)
The main handler for Singular numbers which are suitable for Singular polynomials.
const char *const nDivBy0
#define omFreeSize(addr, size)
#define omFreeBin(addr, bin)
void StringAppendS(const char *st)
void PrintS(const char *s)
static char * nrzCoeffName(const coeffs)
number nrzReadFd(const ssiInfo *d, const coeffs)
void nrzWrite(number a, const coeffs r)
void nrzDelete(number *a, const coeffs)
int nrzSize(number a, const coeffs)
BOOLEAN nrzInitChar(coeffs r, void *parameter)
number nrzInit(long i, const coeffs r)
void nrzWriteFd(number n, const ssiInfo *d, const coeffs)
void s_readmpz_base(s_buff F, mpz_ptr a, int base)