1759{
1762 if (r==0)
1765 {
1766 int*
exp=(
int*)
omAlloc((src->N+1)*
sizeof(int));
1768 {
1770 {
1772 for(
int j=src->N;
j>0;
j--)
1775 #ifdef PDEBUG
1777 #endif
1778 }
1779 }
1781 }
1787 int *exp_q=(
int*)
omAlloc((src->N+1)*
sizeof(int));
1788 for (
int i=1;
i<r;
i++)
1789 {
1790
1791
1792
1794 for(int ii=1;ii<=src->N;ii++)
1796 for(
int ii=0;ii<
i;ii++)
p_Div_hi(J->m[ii],exp_q,src);
1798
1800 for (
int ii=
IDELEMS(J)-1;ii>=0;ii--)
1801 {
1803 {
1806 }
1807 }
1810 poly tmp;
1812 {
1813
1820 {
1822 }
1824 }
1825
1827
1834 }
1836
1838}
static void p_Div_hi(poly p, const int *exp_q, const ring src)
gmp_float exp(const gmp_float &a)
#define omFreeSize(addr, size)
poly p_Power(poly p, int i, const ring r)
static poly p_Neg(poly p, const ring r)
static poly p_Add_q(poly p, poly q, const ring r)
static poly p_Mult_q(poly p, poly q, const ring r)
static void p_SetExpV(poly p, int *ev, const ring r)
static unsigned long p_SetExp(poly p, const unsigned long e, const unsigned long iBitmask, const int VarOffset)
set a single variable exponent @Note: VarOffset encodes the position in p->exp
static void p_Setm(poly p, const ring r)
static long p_GetExp(const poly p, const unsigned long iBitmask, const int VarOffset)
get a single variable exponent @Note: the integer VarOffset encodes:
static void p_Delete(poly *p, const ring r)
static void p_GetExpV(poly p, int *ev, const ring r)
ideal id_CopyFirstK(const ideal ide, const int k, const ring r)
copies the first k (>= 1) entries of the given ideal/module and returns these as a new ideal/module (...