factory's class for integers More...

#include <int_int.h>

Public Member Functions
void *	operator new (size_t)

void	operator delete (void *addr, size_t)

	InternalInteger ()

	InternalInteger (const InternalCF &)

	InternalInteger (const int i)

	InternalInteger (const long i)

	InternalInteger (const char *str, const int base=10)

	InternalInteger (const mpz_ptr mpi)

	~InternalInteger ()

InternalCF *	deepCopyObject () const

const char *	classname () const

void	print (OSTREAM &, char *)

InternalCF *	genZero ()

InternalCF *	genOne ()

bool	is_imm () const

int	levelcoeff () const

InternalCF *	neg ()
	InternalCF * InternalInteger::neg () More...

int	comparesame (InternalCF *)

InternalCF *	addsame (InternalCF *)

InternalCF *	subsame (InternalCF *)

InternalCF *	mulsame (InternalCF *)

InternalCF *	dividesame (InternalCF *)

InternalCF *	modulosame (InternalCF *)

InternalCF *	divsame (InternalCF *)

InternalCF *	modsame (InternalCF *)

void	divremsame (InternalCF , InternalCF &, InternalCF *&)

bool	divremsamet (InternalCF , InternalCF &, InternalCF *&)

int	comparecoeff (InternalCF *)

InternalCF *	addcoeff (InternalCF *)

InternalCF *	subcoeff (InternalCF *, bool)

InternalCF *	mulcoeff (InternalCF *)

InternalCF *	dividecoeff (InternalCF *, bool)

InternalCF *	modulocoeff (InternalCF *, bool)

InternalCF *	divcoeff (InternalCF *, bool)

InternalCF *	modcoeff (InternalCF *, bool)

void	divremcoeff (InternalCF , InternalCF &, InternalCF *&, bool)

bool	divremcoefft (InternalCF , InternalCF &, InternalCF *&, bool)

InternalCF *	bgcdsame (const InternalCF *const) const

InternalCF *	bgcdcoeff (const InternalCF *const)

InternalCF *	bextgcdsame (InternalCF *, CanonicalForm &, CanonicalForm &)

InternalCF *	bextgcdcoeff (InternalCF *, CanonicalForm &, CanonicalForm &)

long	intval () const

int	intmod (int p) const

int	sign () const
	int InternalInteger::sign () const More...

InternalCF *	sqrt ()
	InternalCF * InternalInteger::sqrt () More...

int	ilog2 ()
	int InternalInteger::ilog2 () More...

Public Member Functions inherited from InternalCF
	InternalCF ()

	InternalCF (const InternalCF &)

virtual	~InternalCF ()

int	deleteObject ()

InternalCF *	copyObject ()

virtual InternalCF *	deepCopyObject () const PVIRT_INTCF("deepCopyObject")

virtual const char *	classname () const PVIRT_CHARCC("classname")

virtual InternalCF *	genZero () PVIRT_INTCF("genZero")

virtual InternalCF *	genOne () PVIRT_INTCF("genOne")

virtual int	level () const

virtual int	levelcoeff () const

virtual int	type () const

virtual Variable	variable () const

virtual void	print (OSTREAM &, char *) PVIRT_VOID("print")

virtual bool	inBaseDomain () const

virtual bool	inExtension () const

virtual bool	inCoeffDomain () const

virtual bool	inPolyDomain () const

virtual bool	inQuotDomain () const

virtual bool	isZero () const

virtual bool	isOne () const
	bool InternalCF::isOne, isZero () const More...

virtual bool	isUnivariate () const

virtual long	intval () const

virtual int	intmod (int) const

virtual int	sign () const PVIRT_INT("sign")

virtual InternalCF *	num ()
	InternalCF * InternalCF::num (), den () More...

virtual InternalCF *	den ()

virtual InternalCF *	neg () PVIRT_INTCF("neg")

virtual InternalCF *	invert ()

virtual InternalCF *	tryInvert (const CanonicalForm &, bool &)

virtual int	comparesame (InternalCF *) PVIRT_INT("comparesame")

virtual int	comparecoeff (InternalCF *) PVIRT_INT("comparecoeff")

virtual InternalCF *	addsame (InternalCF *) PVIRT_INTCF("addsame")

virtual InternalCF *	subsame (InternalCF *) PVIRT_INTCF("subsame")

virtual InternalCF *	mulsame (InternalCF *) PVIRT_INTCF("mulsame")

virtual InternalCF *	tryMulsame (InternalCF *, const CanonicalForm &)

virtual InternalCF *	dividesame (InternalCF *) PVIRT_INTCF("dividesame")

virtual InternalCF *	modulosame (InternalCF *) PVIRT_INTCF("modulosame")

virtual InternalCF *	divsame (InternalCF *) PVIRT_INTCF("divsame")

virtual InternalCF *	tryDivsame (InternalCF *, const CanonicalForm &, bool &)

virtual InternalCF *	modsame (InternalCF *) PVIRT_INTCF("modsame")

virtual void	divremsame (InternalCF , InternalCF &, InternalCF *&) PVIRT_VOID("divremsame")

virtual bool	divremsamet (InternalCF , InternalCF &, InternalCF *&) PVIRT_BOOL("divremsamet")

virtual bool	tryDivremsamet (InternalCF , InternalCF &, InternalCF *&, const CanonicalForm &, bool &)

virtual InternalCF *	addcoeff (InternalCF *) PVIRT_INTCF("addcoeff")

virtual InternalCF *	subcoeff (InternalCF *, bool) PVIRT_INTCF("subcoeff")

virtual InternalCF *	mulcoeff (InternalCF *) PVIRT_INTCF("mulcoeff")

virtual InternalCF *	dividecoeff (InternalCF *, bool) PVIRT_INTCF("dividecoeff")

virtual InternalCF *	tryDividecoeff (InternalCF *, bool, const CanonicalForm &, bool &)

virtual InternalCF *	modulocoeff (InternalCF *, bool) PVIRT_INTCF("dividecoeff")

virtual InternalCF *	divcoeff (InternalCF *, bool) PVIRT_INTCF("divcoeff")

virtual InternalCF *	tryDivcoeff (InternalCF *, bool, const CanonicalForm &, bool &)

virtual InternalCF *	modcoeff (InternalCF *, bool) PVIRT_INTCF("modcoeff")

virtual void	divremcoeff (InternalCF , InternalCF &, InternalCF *&, bool) PVIRT_VOID("divremcoeff")

virtual bool	divremcoefft (InternalCF , InternalCF &, InternalCF *&, bool) PVIRT_BOOL("divremcoefft")

virtual bool	tryDivremcoefft (InternalCF , InternalCF &, InternalCF *&, bool, const CanonicalForm &, bool &)

virtual InternalCF *	bgcdsame (const InternalCF *const) const
	InternalCF * InternalCF::bgcdsame, bgcdcoeff ( const InternalCF * const ) More...

virtual InternalCF *	bgcdcoeff (const InternalCF *const)

virtual InternalCF *	bextgcdsame (InternalCF *, CanonicalForm &, CanonicalForm &)
	InternalCF * InternalCF::bextgcdsame ( InternalCF *, CanonicalForm & a, CanonicalForm & b ) More...

virtual InternalCF *	bextgcdcoeff (InternalCF *, CanonicalForm &, CanonicalForm &)

virtual InternalCF *	sqrt ()
	InternalCF * InternalCF::sqrt () More...

virtual int	ilog2 ()
	int InternalCF::ilog2 () More...

virtual CanonicalForm	lc ()
	CanonicalForm InternalCF::lc (), Lc (), LC () More...

virtual CanonicalForm	Lc ()

virtual CanonicalForm	LC ()

virtual CanonicalForm	coeff (int i)
	CanonicalForm InternalCF::coeff ( int i ) More...

virtual int	degree ()
	int InternalCF::degree () More...

virtual int	taildegree ()

virtual CanonicalForm	tailcoeff ()
	CanonicalForm InternalCF::tailcoeff (), int InternalCF::taildegree () More...

Public Member Functions inherited from omallocClass
void *	operator new (size_t size) throw (std::bad_alloc)

void	operator delete (void *block) throw ()

void *	operator new[] (size_t size) throw (std::bad_alloc)

void	operator delete[] (void *block) throw ()

void *	operator new (size_t size, const std::nothrow_t &) throw ()

void *	operator new[] (size_t size, const std::nothrow_t &) throw ()

Private Member Functions
InternalCF *	normalizeMyself ()
	normalizeMyself(), uiNormalizeMyself() - normalize CO. More...

InternalCF *	uiNormalizeMyself ()
	‘uiNormalizeMyself()’ is the same as ‘normalizeMyself()’ except that CO is expected to be non-negative. More...

Static Private Member Functions
static InternalCF *	normalizeMPI (mpz_ptr)
	normalizeMPI(), uiNormalizeMPI() - normalize a mpi. More...

static InternalCF *	uiNormalizeMPI (mpz_ptr)
	‘uiNormalizeMPI()’ is the same as ‘normalizeMPI()’ except that ‘aMpi’ is expected to be non-begative. More...

static mpz_ptr	MPI (const InternalCF *const c)
	MPI() - return underlying mpz_t of ‘c’. More...

Private Attributes
mpz_t	thempi

Static Private Attributes
static const omBin	InternalInteger_bin = omGetSpecBin(sizeof(InternalInteger))

Friends
class	InternalRational

void	gmp_numerator (const CanonicalForm &f, mpz_ptr result)

void	gmp_denominator (const CanonicalForm &f, mpz_ptr result)

void	getmpi (InternalCF *value, mpz_t mpi)

Additional Inherited Members
Protected Member Functions inherited from InternalCF
int	getRefCount ()

void	incRefCount ()

int	decRefCount ()

Detailed Description

factory's class for integers

an integer is represented as an mpz_t thempi

See also: InternalRational

Definition at line 40 of file int_int.h.

Constructor & Destructor Documentation

◆ InternalInteger() [1/6]

InternalInteger::InternalInteger ( )

inline

Definition at line 70 of file int_int.h.

70{ mpz_init( thempi ); }

InternalInteger::thempi

mpz_t thempi

Definition: int_int.h:43

◆ InternalInteger() [2/6]

InternalInteger::InternalInteger ( const InternalCF & )

inline

Definition at line 71 of file int_int.h.

    {
        ASSERT( 0, "ups there is something wrong in your code" );
    }

◆ InternalInteger() [3/6]

InternalInteger::InternalInteger ( const int i )

inline

Definition at line 75 of file int_int.h.

75{ mpz_init_set_si( thempi, (long)i );}

i

int i

Definition: cfEzgcd.cc:132

◆ InternalInteger() [4/6]

InternalInteger::InternalInteger ( const long i )

inline

Definition at line 76 of file int_int.h.

76{ mpz_init_set_si( thempi, i );}

◆ InternalInteger() [5/6]

InternalInteger::InternalInteger	(	const char *	str,
		const int	base = `10`
	)

inline

Definition at line 77 of file int_int.h.

78 { mpz_init_set_str( thempi, str, base ); }

CxxTest::base

char N base

Definition: ValueTraits.h:144

LibThread::str

char * str(leftv arg)

Definition: shared.cc:704

◆ InternalInteger() [6/6]

InternalInteger::InternalInteger ( const mpz_ptr mpi )

inline

Definition at line 79 of file int_int.h.

79{thempi[0]=*mpi;}

◆ ~InternalInteger()

InternalInteger::~InternalInteger ( )

inline

Definition at line 80 of file int_int.h.

80{ mpz_clear( thempi ); }

Member Function Documentation

◆ addcoeff()

InternalCF * InternalInteger::addcoeff ( InternalCF * c )

virtual

Reimplemented from InternalCF.

Definition at line 205 of file int_int.cc.

{
    ASSERT( ::is_imm( c ) == INTMARK, "incompatible base coefficients" );
    long cc = imm2int( c );
    if ( getRefCount() > 1 )
    {
        decRefCount();
        mpz_t dummy;
        mpz_init( dummy );
        if ( cc < 0 )
            mpz_sub_ui( dummy, thempi, -cc );
        else
            mpz_add_ui( dummy, thempi, cc );
        if ( mpz_is_imm( dummy ) )
        {
            InternalCF * res = int2imm( mpz_get_si( dummy ) );
            mpz_clear( dummy );
            return res;
        }
        else
            return new InternalInteger( dummy );
    }
    else
    {
        if ( cc < 0 )
            mpz_sub_ui( thempi, thempi, -cc );
        else
            mpz_add_ui( thempi, thempi, cc );
        if ( mpz_is_imm( thempi ) )
        {
            InternalCF * res = int2imm( mpz_get_si( thempi ) );
            delete this;
            return res;
        }
        else
            return this;
    }
}

◆ addsame()

InternalCF * InternalInteger::addsame ( InternalCF * c )

virtual

Reimplemented from InternalCF.

Definition at line 85 of file int_int.cc.

{
    if ( getRefCount() > 1 )
    {
        decRefCount();
        mpz_t dummy;
        mpz_init( dummy );
        mpz_add( dummy, thempi, MPI( c ) );
        if ( mpz_is_imm( dummy ) )
        {
            InternalCF * res = int2imm( mpz_get_si( dummy ) );
            mpz_clear( dummy );
            return res;
        }
        else
            return new InternalInteger( dummy );
    }
    else
    {
        mpz_add( thempi, thempi, MPI( c ) );
        if ( mpz_is_imm( thempi ) )
        {
            InternalCF * res = int2imm( mpz_get_si( thempi ) );
            delete this;
            return res;
        }
        else
            return this;
    }
}

◆ bextgcdcoeff()

InternalCF * InternalInteger::bextgcdcoeff	(	InternalCF *	c,
		CanonicalForm &	a,
		CanonicalForm &	b
	)

virtual

See also: CanonicalForm::bextgcd(), InternalInteger::bextgcdsame()

Reimplemented from InternalCF.

Definition at line 464 of file int_int.cc.

{
    ASSERT( ::is_imm( c ) == INTMARK, "incompatible base coefficients" );
 
    // simply return 1 if we are calculating over the rationals
    if ( cf_glob_switches.isOn( SW_RATIONAL ) )
    {
        a = 1/CanonicalForm( copyObject() ); b = 0;
        return int2imm( 1 );
    }
 
    long cInt = imm2int( c );
 
    // trivial cases
    if ( cInt == 1 || cInt == -1 )
    {
        a = 0; b = cInt;
        return int2imm( 1 );
    }
    else if ( cInt == 0 )
    {
        a = 1; b = 0;
        return copyObject();
    }
 
    // calculate q and r such that CO = q*cInt + r
    InternalCF * q = 0, * r = 0;
    divremcoeff( c, q, r, false );
 
    // we do not repeat all the code to calculate the gcd of two
    // immediates.  Note that r is an immediate since c != 0, so
    // we do not have to destroy it.  q is destroyed by the
    // CanonicalForm destructor, hence we do not need to worry
    // about it, either.
    CanonicalForm aPrime, bPrime;
    CanonicalForm result = bextgcd( c, r, aPrime, bPrime );
    a = bPrime;
    b = aPrime - CanonicalForm( q ) * bPrime;
 
    return result.getval();
}

◆ bextgcdsame()

InternalCF * InternalInteger::bextgcdsame	(	InternalCF *	c,
		CanonicalForm &	a,
		CanonicalForm &	b
	)

virtual

See also: CanonicalForm::bextgcd(), InternalInteger::bextgcdcoeff()

Reimplemented from InternalCF.

Definition at line 409 of file int_int.cc.

{
    ASSERT( ! ::is_imm( c ) && c->levelcoeff() == IntegerDomain, "incompatible base coefficients" );
 
    // simply return 1 if we are calculating over the rationals
    if ( cf_glob_switches.isOn( SW_RATIONAL ) )
    {
        a = 1/CanonicalForm( copyObject() ); b = 0;
        return int2imm( 1 );
    }
 
    // calculate extended gcd
    mpz_t result, aMPI, bMPI;
    mpz_init( result );
    mpz_init( aMPI );
    mpz_init( bMPI );
    mpz_gcdext( result, aMPI, bMPI, thempi, MPI( c ) );
 
    // check and modify signs
    if ( mpz_sgn( result ) < 0 )
    {
        mpz_neg( result, result );
        mpz_neg( aMPI, aMPI );
        mpz_neg( bMPI, bMPI );
    }
 
    // postconditioning of result
    if ( mpz_is_imm( aMPI ) )
    {
        a = CanonicalForm( int2imm( mpz_get_si( aMPI ) ) );
        mpz_clear( aMPI );
    }
    else
        a = CanonicalForm( new InternalInteger( aMPI ) );
    if ( mpz_is_imm( bMPI ) )
    {
        b = CanonicalForm( int2imm( mpz_get_si( bMPI ) ) );
        mpz_clear( bMPI );
    }
    else
        b = CanonicalForm( new InternalInteger( bMPI ) );
    if ( mpz_is_imm( result ) )
    {
        InternalCF * res = int2imm( mpz_get_si( result ) );
        mpz_clear( result );
        return res;
    }
    else
        return new InternalInteger( result );
}

◆ bgcdcoeff()

InternalCF * InternalInteger::bgcdcoeff ( const InternalCF * const c )

virtual

See also: CanonicalForm::bgcd(), InternalInteger::bgcdsame()

Reimplemented from InternalCF.

Definition at line 377 of file int_int.cc.

{
    ASSERT( ::is_imm( c ) == INTMARK, "incompatible base coefficients" );
 
    // simply return 1 if we are calculating over the rationals
    if ( cf_glob_switches.isOn( SW_RATIONAL ) )
         return int2imm( 1 );
 
    long cInt = imm2int( c );
 
    // trivial cases
    if ( cInt == 1 || cInt == -1 )
        return int2imm( 1 );
    else if ( cInt == 0 )
        return copyObject();
 
    // calculate gcd.  We need a positive operand since
    // `mpz_gcd_ui()' operates an unsigned int's only.
    if ( cInt < 0 ) cInt = -cInt;
    mpz_t dummy;
    mpz_init( dummy );
    // we do not need dummy since we know that cInt != 0
    cInt = mpz_gcd_ui( dummy, thempi, cInt );
    mpz_clear( dummy );
    if ( cInt < 0 ) cInt = -cInt;
    return int2imm( cInt );
}

◆ bgcdsame()

InternalCF * InternalInteger::bgcdsame ( const InternalCF * const c ) const

virtual

See also: CanonicalForm::bgcd(), InternalInteger::bgcdcoeff()

Reimplemented from InternalCF.

Definition at line 348 of file int_int.cc.

{
    ASSERT( ! ::is_imm( c ) && c->levelcoeff() == IntegerDomain, "incompatible base coefficients" );
 
    // simply return 1 if we are calculating over the rationals
    if ( cf_glob_switches.isOn( SW_RATIONAL ) )
         return int2imm( 1 );
 
    // calculate gcd
    mpz_t result;
    mpz_init( result );
    mpz_gcd( result, thempi, MPI( c ) );
    mpz_abs( result, result );
 
    // check for immediate result
    if ( mpz_is_imm( result ) )
    {
        InternalCF * res = int2imm( mpz_get_si( result ) );
        mpz_clear( result );
        return res;
    }
    else
        return new InternalInteger( result );
}

◆ classname()

const char * InternalInteger::classname ( ) const

inlinevirtual

Reimplemented from InternalCF.

Definition at line 82 of file int_int.h.

82{ return "InternalInteger"; }

◆ comparecoeff()

int InternalInteger::comparecoeff ( InternalCF * c )

virtual

See also: CanonicalForm::operator <(), CanonicalForm::operator ==(), InternalInteger::comparesame()

Reimplemented from InternalCF.

Definition at line 198 of file int_int.cc.

{
    ASSERT( ::is_imm( c ) == INTMARK, "incompatible base coefficients" );
    return mpz_cmp_si( thempi, imm2int( c ) );
}

◆ comparesame()

int InternalInteger::comparesame ( InternalCF * c )

virtual

See also: CanonicalForm::operator <(), CanonicalForm::operator ==(), InternalInteger::comparecoeff()

Reimplemented from InternalCF.

Definition at line 188 of file int_int.cc.

{
    ASSERT( ! ::is_imm( c ) && c->levelcoeff() == IntegerDomain, "incompatible base coefficients" );
    return mpz_cmp( thempi, MPI( c ) );
}

◆ deepCopyObject()

InternalCF * InternalInteger::deepCopyObject ( ) const

virtual

Reimplemented from InternalCF.

Definition at line 18 of file int_int.cc.

{
    mpz_t dummy;
    mpz_init_set( dummy, thempi );
    return new InternalInteger( dummy );
}

◆ divcoeff()

InternalCF * InternalInteger::divcoeff	(	InternalCF *	c,
		bool	invert
	)

virtual

See also: CanonicalForm::div(), InternalInteger::divsame()

Reimplemented from InternalCF.

Definition at line 151 of file int_intdiv.cc.

{
    ASSERT( ::is_imm( c ) == INTMARK,
            "type error: immediate integer expected" );
    ASSERT( invert || imm2int( c ) != 0,
            "math error: divide by zero" );
    ASSERT( ! invert || imm2int( c ) == 0,
            "math error: c does not divide CO" );
 
    if ( invert ) {
        if ( deleteObject() ) delete this;
        // this may happen iff `c' == 0
        return int2imm( 0 );
    } else if ( getRefCount() > 1 ) {
        deleteObject();
        mpz_t mpiC;
        mpz_t mpiResult;
        mpz_init_set_si( mpiC, imm2int( c ) );
        mpz_init( mpiResult );
        mpz_divexact( mpiResult, thempi, mpiC );
        mpz_clear( mpiC );
        return normalizeMPI( mpiResult );
    } else {
        mpz_t mpiC;
        mpz_init_set_si( mpiC, imm2int( c ) );
        mpz_divexact( thempi, thempi, mpiC );
        mpz_clear( mpiC );
        return normalizeMyself();
    }
}

◆ dividecoeff()

InternalCF * InternalInteger::dividecoeff	(	InternalCF *	c,
		bool	invert
	)

virtual

See also: CanonicalForm::operator /(), InternalInteger::dividesame()

Reimplemented from InternalCF.

Definition at line 69 of file int_intdiv.cc.

{
    ASSERT( ::is_imm( c ) == INTMARK,
            "type error: immediate integer expected" );
    ASSERT( invert || imm2int( c ) != 0,
            "math error: divide by zero" );
 
    long intC = imm2int( c );
 
    if ( cf_glob_switches.isOn( SW_RATIONAL ) ) {
        mpz_t n, d;
        if ( invert ) {
            mpz_init_set_si( n, intC );
            mpz_init_set( d, thempi );
        } else {
            mpz_init_set( n, thempi );
            mpz_init_set_si( d, intC );
        }
        if ( deleteObject() ) delete this;
        InternalRational * result = new InternalRational( n, d );
        return result->normalize_myself();
    }
 
    if ( invert ) {
        int mpiSign = mpz_sgn( thempi );
        if ( deleteObject() ) delete this;
        if ( intC >= 0 )
            return int2imm( 0 );
        else
            return int2imm( -mpiSign );
    } else if ( getRefCount() > 1 ) {
        decRefCount();
        mpz_t mpiResult;
        mpz_init( mpiResult );
        if ( intC > 0 )
            mpz_fdiv_q_ui( mpiResult, thempi, intC );
        else {
            mpz_fdiv_q_ui( mpiResult, thempi, -intC );
            mpz_neg( mpiResult, mpiResult );
        }
        return normalizeMPI( mpiResult );
    } else {
        if ( intC > 0 )
            mpz_fdiv_q_ui( thempi, thempi, intC );
        else {
            mpz_fdiv_q_ui( thempi, thempi, -intC );
            mpz_neg( thempi, thempi );
        }
        return normalizeMyself();
    }
}

◆ dividesame()

InternalCF * InternalInteger::dividesame ( InternalCF * c )

virtual

See also: CanonicalForm::operator /(), InternalInteger::dividecoeff()

Reimplemented from InternalCF.

Definition at line 28 of file int_intdiv.cc.

{
    ASSERT( ! ::is_imm( c ) && c->levelcoeff() == IntegerDomain,
            "type error: InternalInteger expected" );
 
    if ( c == this ) {
        if ( deleteObject() ) delete this;
        return int2imm( 1 );
    }
 
    if ( cf_glob_switches.isOn( SW_RATIONAL ) ) {
        mpz_t n, d;
        mpz_init_set( n, thempi );
        mpz_init_set( d, MPI( c ) );
        if ( deleteObject() ) delete this;
        InternalRational * result = new InternalRational( n, d );
        return result->normalize_myself();
    }
 
    if ( getRefCount() > 1 ) {
        decRefCount();
        mpz_t mpiResult;
        mpz_init( mpiResult );
        if ( mpz_sgn( MPI( c ) ) > 0 )
            mpz_fdiv_q( mpiResult, thempi, MPI( c ) );
        else
            mpz_cdiv_q( mpiResult, thempi, MPI( c ) );
        return normalizeMPI( mpiResult );
    } else {
        if ( mpz_sgn( MPI( c ) ) > 0 )
            mpz_fdiv_q( thempi, thempi, MPI( c ) );
        else
            mpz_cdiv_q( thempi, thempi, MPI( c ) );
        return normalizeMyself();
    }
}

◆ divremcoeff()

void InternalInteger::divremcoeff	(	InternalCF *	c,
		InternalCF *&	quot,
		InternalCF *&	rem,
		bool	invert
	)

virtual

See also: CanonicalForm::divrem(), InternalInteger::divremsame()

Reimplemented from InternalCF.

Definition at line 308 of file int_intdiv.cc.

{
    ASSERT( ::is_imm( c ) == INTMARK,
            "type error: immediate integer expected" );
    ASSERT( invert || imm2int( c ) != 0,
            "math error: divide by zero" );
 
    long intC = imm2int( c );
 
    if ( cf_glob_switches.isOn( SW_RATIONAL ) ) {
        mpz_t n, d;
        if ( invert ) {
            mpz_init_set_si( n, intC );
            mpz_init_set( d, thempi );
        } else {
            mpz_init_set( n, thempi );
            mpz_init_set_si( d, intC );
        }
        InternalRational * result = new InternalRational( n, d );
        quot = result->normalize_myself();
        rem = int2imm( 0 );
        return;
    }
 
    if ( invert ) {
        if ( intC >= 0 ) {
            rem = c;
            quot = int2imm( 0 );
        } else {
            mpz_t mpiResult;
            mpz_init_set( mpiResult, thempi );
            mpz_abs( mpiResult, mpiResult );
            mpz_sub_ui( mpiResult, mpiResult, -intC );
            rem = uiNormalizeMPI( mpiResult );
            quot = int2imm( -mpz_sgn( thempi ) );
        }
    } else {
        mpz_t q;
        mpz_t dummy;
        mpz_init( q ); mpz_init( dummy );
        if ( intC > 0 ) {
            rem = int2imm( mpz_fdiv_qr_ui( q, dummy, thempi, intC ) );
            quot = normalizeMPI( q );
        } else {
            rem = int2imm( mpz_fdiv_qr_ui( q, dummy, thempi, -intC ) );
            mpz_neg( q, q );
            quot = normalizeMPI( q );
        }
        mpz_clear( dummy );
    }
}

◆ divremcoefft()

bool InternalInteger::divremcoefft	(	InternalCF *	c,
		InternalCF *&	quot,
		InternalCF *&	rem,
		bool	invert
	)

virtual

See also: CanonicalForm::divremt(), InternalInteger::divremsamet()

Reimplemented from InternalCF.

Definition at line 374 of file int_intdiv.cc.

{
    divremcoeff( c, quot, rem, invert );
    return true;
}

◆ divremsame()

void InternalInteger::divremsame	(	InternalCF *	c,
		InternalCF *&	quot,
		InternalCF *&	rem
	)

virtual

See also: CanonicalForm::divrem(), InternalInteger::divremcoeff()

Reimplemented from InternalCF.

Definition at line 271 of file int_intdiv.cc.

{
    ASSERT( ! ::is_imm( c ) && c->levelcoeff() == IntegerDomain,
            "type error: InternalInteger expected" );
 
    if ( c == this ) {
        quot = int2imm( 1 );
        rem = int2imm( 0 );
        return;
    }
 
    if ( cf_glob_switches.isOn( SW_RATIONAL ) ) {
        mpz_t n, d;
        mpz_init_set( n, thempi );
        mpz_init_set( d, MPI( c ) );
        InternalRational * result = new InternalRational( n, d );
        quot = result->normalize_myself();
        rem = int2imm( 0 );
        return;
    }
 
    mpz_t q;
    mpz_t r;
    mpz_init( q ); mpz_init( r );
    if ( mpz_sgn( MPI( c ) ) > 0 )
        mpz_fdiv_qr( q, r, thempi, MPI( c ) );
    else
        mpz_cdiv_qr( q, r, thempi, MPI( c ) );
 
    quot = normalizeMPI( q );
    rem = uiNormalizeMPI( r );
}

◆ divremsamet()

bool InternalInteger::divremsamet	(	InternalCF *	c,
		InternalCF *&	quot,
		InternalCF *&	rem
	)

virtual

See also: CanonicalForm::divremt(), InternalInteger::divremcoefft()

Reimplemented from InternalCF.

Definition at line 364 of file int_intdiv.cc.

{
    divremsame( c, quot, rem );
    return true;
}

◆ divsame()

InternalCF * InternalInteger::divsame ( InternalCF * c )

virtual

See also: CanonicalForm::div(), InternalInteger::divcoeff()

Reimplemented from InternalCF.

Definition at line 125 of file int_intdiv.cc.

{
    ASSERT( ! ::is_imm( c ) && c->levelcoeff() == IntegerDomain,
            "type error: InternalInteger expected" );
 
    if ( c == this ) {
        if ( deleteObject() ) delete this;
        return int2imm( 1 );
    }
 
    if ( getRefCount() > 1 ) {
        deleteObject();
        mpz_t mpiResult;
        mpz_init( mpiResult );
        mpz_divexact( mpiResult, thempi, MPI( c ) );
        return normalizeMPI( mpiResult );
    } else {
        mpz_divexact( thempi, thempi, MPI( c ) );
        return normalizeMyself();
    }
}

◆ genOne()

InternalCF * InternalInteger::genOne ( )

virtual

Reimplemented from InternalCF.

Definition at line 54 of file int_int.cc.

{
    if ( isOne() )
        return copyObject();
    else
        return new InternalInteger( 1 );
}

◆ genZero()

InternalCF * InternalInteger::genZero ( )

virtual

Reimplemented from InternalCF.

Definition at line 46 of file int_int.cc.

{
    if ( isZero() )
        return copyObject();
    else
        return new InternalInteger();
}

◆ ilog2()

int InternalInteger::ilog2 ( )

virtual

int InternalInteger::ilog2 ()

See also: CanonicalForm::ilog2()

Reimplemented from InternalCF.

Definition at line 549 of file int_int.cc.

{
    ASSERT( mpz_cmp_si( thempi, 0 ) > 0, "log() argument <= 0" );
    return mpz_sizeinbase( thempi, 2 ) - 1;
}

◆ intmod()

int InternalInteger::intmod ( int p ) const

virtual

Reimplemented from InternalCF.

Definition at line 511 of file int_int.cc.

{
  return (int)mpz_fdiv_ui( thempi, (unsigned long)p );
}

◆ intval()

long InternalInteger::intval ( ) const

virtual

Reimplemented from InternalCF.

Definition at line 506 of file int_int.cc.

{
  return mpz_get_si( thempi );
}

◆ is_imm()

bool InternalInteger::is_imm ( ) const

Definition at line 41 of file int_int.cc.

{
    return mpz_is_imm( thempi );
}

◆ levelcoeff()

int InternalInteger::levelcoeff ( ) const

inlinevirtual

Reimplemented from InternalCF.

Definition at line 91 of file int_int.h.

91{ return IntegerDomain; }

◆ modcoeff()

InternalCF * InternalInteger::modcoeff	(	InternalCF *	c,
		bool	invert
	)

virtual

See also: see CanonicalForm::mod(), InternalInteger::modsame()

Reimplemented from InternalCF.

Definition at line 262 of file int_intdiv.cc.

{
    return modulocoeff( c, invert );
}

◆ modsame()

InternalCF * InternalInteger::modsame ( InternalCF * c )

virtual

See also: see CanonicalForm::mod(), InternalInteger::modcoeff()

Reimplemented from InternalCF.

Definition at line 253 of file int_intdiv.cc.

{
    return modulosame( c );
}

◆ modulocoeff()

InternalCF * InternalInteger::modulocoeff	(	InternalCF *	c,
		bool	invert
	)

virtual

See also: CanonicalForm::operator %(), InternalInteger::modulosame()

Reimplemented from InternalCF.

Definition at line 212 of file int_intdiv.cc.

{
    ASSERT( ::is_imm( c ) == INTMARK,
            "type error: immediate integer expected" );
    ASSERT( invert || imm2int( c ) != 0,
            "math error: divide by zero" );
 
    if ( cf_glob_switches.isOn( SW_RATIONAL ) ) {
        if ( deleteObject() ) delete this;
        return int2imm( 0 );
    }
 
    long intC = imm2int( c );
 
    if ( invert ) {
        if ( intC >= 0 ) {
            if ( deleteObject() ) delete this;
            return c;
        } else {
            // no checks for refCount == 1 are done.  It is not worth ...
            mpz_t mpiResult;
            mpz_init_set( mpiResult, thempi );
            mpz_abs( mpiResult, mpiResult );
            mpz_sub_ui( mpiResult, mpiResult, -intC );
            if ( deleteObject() ) delete this;
            return uiNormalizeMPI( mpiResult );
        }
    } else {
        mpz_t dummy;
        mpz_init( dummy );
        InternalCF * result = int2imm( mpz_mod_ui( dummy, thempi, tabs( intC ) ) );
        mpz_clear( dummy );
        if ( deleteObject() ) delete this;
        return result;
    }
}

◆ modulosame()

InternalCF * InternalInteger::modulosame ( InternalCF * c )

virtual

See also: CanonicalForm::operator %(), InternalInteger::modulocoeff()

Reimplemented from InternalCF.

Definition at line 186 of file int_intdiv.cc.

{
    ASSERT( ! ::is_imm( c ) && c->levelcoeff() == IntegerDomain,
            "type error: InternalInteger expected" );
 
    if ( (c == this) || cf_glob_switches.isOn( SW_RATIONAL ) ) {
        if ( deleteObject() ) delete this;
        return int2imm( 0 );
    }
 
    if ( getRefCount() > 1 ) {
        decRefCount();
        mpz_t mpiResult;
        mpz_init( mpiResult );
        mpz_mod( mpiResult, thempi, MPI( c ) );
        return uiNormalizeMPI( mpiResult );
    } else {
        mpz_mod( thempi, thempi, MPI( c ) );
        return uiNormalizeMyself();
    }
}

◆ MPI()

mpz_ptr InternalInteger::MPI ( const InternalCF *const c )

inlinestaticprivate

MPI() - return underlying mpz_t of ‘c’.

‘c’ is expected to be an ‘InternalInteger *’. `c's underlying mpz_t is returned.

Definition at line 232 of file int_int.h.

{
    return (((InternalInteger*)c)->thempi);
}

◆ mulcoeff()

InternalCF * InternalInteger::mulcoeff ( InternalCF * c )

virtual

Reimplemented from InternalCF.

Definition at line 299 of file int_int.cc.

{
    ASSERT( ::is_imm( c ) == INTMARK, "incompatible base coefficients" );
    long cc = imm2int( c );
    if ( getRefCount() > 1 )
    {
        decRefCount();
        mpz_t dummy;
        mpz_init( dummy );
        if ( cc < 0 )
        {
            mpz_mul_ui( dummy, thempi, -cc );
            mpz_neg( dummy, dummy );
        }
        else
            mpz_mul_ui( dummy, thempi, cc );
        if ( mpz_is_imm( dummy ) )
        {
            InternalCF * res = int2imm( mpz_get_si( dummy ) );
            mpz_clear( dummy );
            return res;
        }
        else
            return new InternalInteger( dummy );
    }
    else
    {
        if ( cc < 0 )
        {
            mpz_mul_ui( thempi, thempi, -cc );
            mpz_neg( thempi, thempi );
        }
        else
            mpz_mul_ui( thempi, thempi, cc );
        if ( mpz_is_imm( thempi ) )
        {
            InternalCF * res = int2imm( mpz_get_si( thempi ) );
            delete this;
            return res;
        }
        else
            return this;
    }
}

◆ mulsame()

InternalCF * InternalInteger::mulsame ( InternalCF * c )

virtual

Reimplemented from InternalCF.

Definition at line 147 of file int_int.cc.

{
    if ( getRefCount() > 1 )
    {
        decRefCount();
        mpz_t dummy;
        mpz_init( dummy );
        mpz_mul( dummy, thempi, MPI( c ) );
        #if 0
        if ( mpz_is_imm( dummy ) )
        {
        // can this happen ???
            InternalCF * res = int2imm( mpz_get_si( dummy ) );
            mpz_clear( dummy );
            return res;
        }
        else
        #endif
            return new InternalInteger( dummy );
    }
    else
    {
        mpz_mul( thempi, thempi, MPI( c ) );
        #if 0
        if ( mpz_is_imm( &thempi ) )
        {
        // can this happen ???
            InternalCF * res = int2imm( mpz_get_si( &thempi ) );
            delete this;
            return res;
        }
        else
        #endif
            return this;
    }
}

◆ neg()

InternalCF * InternalInteger::neg ( )

virtual

InternalCF * InternalInteger::neg ()

See also: CanonicalForm::operator -()

Reimplemented from InternalCF.

Definition at line 66 of file int_int.cc.

{
    if ( getRefCount() > 1 )
    {
        decRefCount();
        mpz_t dummy;
        mpz_init_set( dummy, thempi );
        mpz_neg( dummy, dummy );
        return new InternalInteger( dummy );
    }
    else
    {
        mpz_neg( thempi, thempi );
        return this;
    }
}

◆ normalizeMPI()

InternalCF * InternalInteger::normalizeMPI ( mpz_ptr aMpi )

inlinestaticprivate

normalizeMPI(), uiNormalizeMPI() - normalize a mpi.

If ‘aMpi’ fits into an immediate integer, clear ‘aMpi’ and return the immediate. Otherwise, return a new ‘InternalInteger’ with ‘aMpi’ as underlying mpi.

Definition at line 196 of file int_int.h.

{
    if ( mpz_is_imm( aMpi ) ) {
        InternalCF * result = int2imm( mpz_get_si( aMpi ) );
        mpz_clear( aMpi );
        return result;
    } else
        return new InternalInteger( aMpi );
}

◆ normalizeMyself()

InternalCF * InternalInteger::normalizeMyself ( )

inlineprivate

normalizeMyself(), uiNormalizeMyself() - normalize CO.

If CO fits into an immediate integer, delete CO and return the immediate. Otherwise, return a pointer to CO.

Note: We do not mind reference counting at this point! CO is deleted unconditionally!

Definition at line 152 of file int_int.h.

{
    ASSERT( getRefCount() == 1, "internal error: must not delete CO" );
 
    if ( mpz_is_imm( thempi ) ) {
        InternalCF * result = int2imm( mpz_get_si( thempi ) );
        delete this;
        return result;
    } else
        return this;
}

◆ operator delete()

void InternalInteger::operator delete	(	void *	addr,
		size_t
	)

inline

Definition at line 64 of file int_int.h.

    {
      omFreeBin(addr, InternalInteger_bin);
    }

◆ operator new()

void * InternalInteger::operator new ( size_t )

inline

Definition at line 58 of file int_int.h.

    {
      void* addr;
      omTypeAllocBin(void*, addr, InternalInteger_bin);
      return addr;
    }

◆ print()

void InternalInteger::print	(	OSTREAM &	os,
		char *	c
	)

virtual

Reimplemented from InternalCF.

Definition at line 26 of file int_int.cc.

{
    if ( *c == '*' && mpz_cmp_si( thempi, 1 ) == 0 )
        os << c+1;
    else if ( *c == '*' && mpz_cmp_si( thempi, -1 ) == 0 )
        os << '-' << c+1;
    else {
        char * str = new char[mpz_sizeinbase( thempi, 10 ) + 2];
        str = mpz_get_str( str, 10, thempi );
        os << str << c;
        delete [] str;
    }
}

◆ sign()

int InternalInteger::sign ( ) const

virtual

int InternalInteger::sign () const

See also: CanonicalForm::sign()

Reimplemented from InternalCF.

Definition at line 520 of file int_int.cc.

{
    return mpz_sgn( thempi );
}

◆ sqrt()

InternalCF * InternalInteger::sqrt ( )

virtual

InternalCF * InternalInteger::sqrt ()

See also: CanonicalForm::sqrt()

Reimplemented from InternalCF.

Definition at line 529 of file int_int.cc.

{
    ASSERT( mpz_cmp_si( thempi, 0 ) >= 0, "sqrt() argument < 0" );
    mpz_t result;
    mpz_init( result );
    mpz_sqrt( result, thempi );
    if ( mpz_is_imm( result ) )
    {
        InternalCF * res = int2imm( mpz_get_si( result ) );
        mpz_clear( result );
        return res;
    }
    else
        return new InternalInteger( result );
}

◆ subcoeff()

InternalCF * InternalInteger::subcoeff	(	InternalCF *	c,
		bool	negate
	)

virtual

Reimplemented from InternalCF.

Definition at line 244 of file int_int.cc.

{
    ASSERT( ::is_imm( c ) == INTMARK, "incompatible base coefficients" );
    long cc = imm2int( c );
    if ( getRefCount() > 1 )
    {
        decRefCount();
        mpz_t dummy;
        if ( negate )
        {
            mpz_init_set_si( dummy, cc );
            mpz_sub( dummy, dummy, thempi );
        }
        else
        {
            mpz_init( dummy );
            if ( cc < 0 )
                mpz_add_ui( dummy, thempi, -cc );
            else
                mpz_sub_ui( dummy, thempi, cc );
        }
        if ( mpz_is_imm( dummy ) )
        {
            InternalCF * res = int2imm( mpz_get_si( dummy ) );
            mpz_clear( dummy );
            return res;
        }
        else
            return new InternalInteger( dummy );
    }
    else
    {
        if ( negate )
        {
            mpz_t dummy;
            mpz_init_set_si( dummy, cc );
            mpz_sub( thempi, dummy, thempi );
            mpz_clear( dummy );
        }
        else
            if ( cc < 0 )
                mpz_add_ui( thempi, thempi, -cc );
            else
                mpz_sub_ui( thempi, thempi, cc );
        if ( mpz_is_imm( thempi ) )
        {
            InternalCF * res = int2imm( mpz_get_si( thempi ) );
            delete this;
            return res;
        }
        else
            return this;
    }
}

◆ subsame()

InternalCF * InternalInteger::subsame ( InternalCF * c )

virtual

Reimplemented from InternalCF.

Definition at line 116 of file int_int.cc.

{
    if ( getRefCount() > 1 )
    {
        decRefCount();
        mpz_t dummy;
        mpz_init( dummy );
        mpz_sub( dummy, thempi, MPI( c ) );
        if ( mpz_is_imm( dummy ) )
        {
            InternalCF * res = int2imm( mpz_get_si( dummy ) );
            mpz_clear( dummy );
            return res;
        }
        else
            return new InternalInteger( dummy );
    }
    else
    {
        mpz_sub( thempi, thempi, MPI( c ) );
        if ( mpz_is_imm( thempi ) )
        {
            InternalCF * res = int2imm( mpz_get_si( thempi ) );
            delete this;
            return res;
        }
        else
            return this;
    }
}

◆ uiNormalizeMPI()

InternalCF * InternalInteger::uiNormalizeMPI ( mpz_ptr aMpi )

inlinestaticprivate

‘uiNormalizeMPI()’ is the same as ‘normalizeMPI()’ except that ‘aMpi’ is expected to be non-begative.

In this case, we may use ‘mpz_get_ui()’ to convert ‘aMpi’ into an immediate which is slightly faster than the signed variant.

Definition at line 213 of file int_int.h.

{
    if ( mpz_is_imm( aMpi ) ) {
        InternalCF * result = int2imm( mpz_get_ui( aMpi ) );
        mpz_clear( aMpi );
        return result;
    } else
        return new InternalInteger( aMpi );
}

◆ uiNormalizeMyself()

InternalCF * InternalInteger::uiNormalizeMyself ( )

inlineprivate

‘uiNormalizeMyself()’ is the same as ‘normalizeMyself()’ except that CO is expected to be non-negative.

In this case, we may use ‘mpz_get_ui()’ to convert the underlying mpi into an immediate which is slightly faster than the signed variant.

Note: We do not mind reference counting at this point! CO is deleted unconditionally!

Definition at line 174 of file int_int.h.

{
    ASSERT( getRefCount() == 1, "internal error: must not delete CO" );
 
    if ( mpz_is_imm( thempi ) ) {
        InternalCF * result = int2imm( mpz_get_ui( thempi ) );
        delete this;
        return result;
    } else
        return this;
}

Friends And Related Function Documentation

◆ getmpi

void getmpi	(	InternalCF *	value,
		mpz_t	mpi
	)

friend

Definition at line 303 of file cf_factory.cc.

{
    ASSERT( ! is_imm( value ) && (value->levelcoeff() == IntegerDomain || value->levelcoeff() == PrimePowerDomain), "illegal operation" );
    mpz_init_set (mpi, ((InternalInteger*)value)->thempi);
}

◆ gmp_denominator

void gmp_denominator	(	const CanonicalForm &	f,
		mpz_ptr	result
	)

friend

Definition at line 40 of file singext.cc.

{
    InternalCF * ff = f.getval();
    ASSERT( ! is_imm( ff ), "illegal type" );
    if ( ff->levelcoeff() == IntegerDomain )
    {
        mpz_init_set_si( result, 1 );
        ff->deleteObject();
    }
    else  if ( ff->levelcoeff() == RationalDomain )
    {
        mpz_init_set( result, (InternalRational::MPQDEN( ff )) );
        ff->deleteObject();
    }
    else
    {
        ASSERT( 0, "illegal type" );
    }
}

◆ gmp_numerator

void gmp_numerator	(	const CanonicalForm &	f,
		mpz_ptr	result
	)

friend

Definition at line 20 of file singext.cc.

{
    InternalCF * ff = f.getval();
    ASSERT( ! is_imm( ff ), "illegal type" );
    if ( ff->levelcoeff() == IntegerDomain )
    {
        mpz_init_set( result, (InternalInteger::MPI( ff )) );
        ff->deleteObject();
    }
    else  if ( ff->levelcoeff() == RationalDomain )
    {
        mpz_init_set( result, (InternalRational::MPQNUM( ff )) );
        ff->deleteObject();
    }
    else
    {
        ASSERT( 0, "illegal type" );
    }
}

◆ InternalRational

friend class InternalRational

friend

Definition at line 134 of file int_int.h.

Field Documentation

◆ InternalInteger_bin

const omBin InternalInteger::InternalInteger_bin = omGetSpecBin(sizeof(InternalInteger))

staticprivate

Definition at line 54 of file int_int.h.

◆ thempi

mpz_t InternalInteger::thempi

private

Definition at line 43 of file int_int.h.

The documentation for this class was generated from the following files:

factory/int_int.h
factory/int_int.cc
factory/int_intdiv.cc

Public Member Functions

Private Member Functions

Static Private Member Functions

Private Attributes

Static Private Attributes

Friends

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ InternalInteger() [1/6]

◆ InternalInteger() [2/6]

◆ InternalInteger() [3/6]

◆ InternalInteger() [4/6]

◆ InternalInteger() [5/6]

◆ InternalInteger() [6/6]

◆ ~InternalInteger()

Member Function Documentation

◆ addcoeff()

◆ addsame()

◆ bextgcdcoeff()

◆ bextgcdsame()

◆ bgcdcoeff()

◆ bgcdsame()

◆ classname()

◆ comparecoeff()

◆ comparesame()

◆ deepCopyObject()

◆ divcoeff()

◆ dividecoeff()

◆ dividesame()

◆ divremcoeff()

◆ divremcoefft()

◆ divremsame()

◆ divremsamet()

◆ divsame()

◆ genOne()

◆ genZero()

◆ ilog2()

◆ intmod()

◆ intval()

◆ is_imm()

◆ levelcoeff()

◆ modcoeff()

◆ modsame()

◆ modulocoeff()

◆ modulosame()

◆ MPI()

◆ mulcoeff()

◆ mulsame()

◆ neg()

◆ normalizeMPI()

◆ normalizeMyself()

◆ operator delete()

◆ operator new()

◆ print()

◆ sign()

◆ sqrt()

◆ subcoeff()

◆ subsame()

◆ uiNormalizeMPI()

◆ uiNormalizeMyself()

Friends And Related Function Documentation

◆ getmpi

◆ gmp_denominator

◆ gmp_numerator

◆ InternalRational

Field Documentation

◆ InternalInteger_bin

◆ thempi