1850{
1853 if (r==0)
1856 {
1857 int*
exp=(
int*)
omAlloc((src->N+1)*
sizeof(int));
1859 {
1861 {
1863 for(
int j=src->N;
j>0;
j--)
1864 {
1865 int w=(*wdegree)[
j-1];
1867 {
1868 WerrorS(
"weights must be positive");
1870 }
1872 }
1874 #ifdef PDEBUG
1876 #endif
1877 }
1878 }
1880 }
1883
1885 int *exp_q=(
int*)
omAlloc((src->N+1)*
sizeof(int));
1887 for (
int i=1;
i<r;
i++)
1888 {
1890 for(int ii=src->N;ii>0;ii--)
1893 for(
int ii=0;ii<
i;ii++)
1894 {
1896 }
1898
1899
1901 for (
int ii=
IDELEMS(J)-1;ii>=0;ii--)
1902 {
1904 {
1907 }
1908 }
1911 poly tmp;
1913 {
1914
1917
1920 {
1922 }
1924 }
1925
1927
1930
1933 }
1936
1938}
static FORCE_INLINE number n_Init(long i, const coeffs r)
a number representing i in the given coeff field/ring r
void WerrorS(const char *s)
static BOOLEAN p_Div_hi(poly p, const int *exp_q, const ring src)
static void id_DelDiv_hi(ideal id, BOOLEAN *bad, const ring r)
gmp_float exp(const gmp_float &a)
#define omFreeSize(addr, size)
poly p_Power(poly p, int i, 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_GetExpV(poly p, int *ev, const ring r)
static poly p_Init(const ring r, omBin bin)
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 (...