26 extern "C" double (*
wFunctional)(
int *degw,
int *lpol,
int npol,
27 double *rel,
double wx,
double wNsqr);
29 double *rel,
double wx,
double wNsqr);
31 double *rel,
double wx,
double wNsqr);
32 extern "C" void wAdd(
int *
A,
int mons,
int kn,
int xx,
int rvar);
33 extern "C" void wNorm(
int *degw,
int *lpol,
int npol,
double *rel);
35 int *lpol,
int npol,
double *rel,
double *fopt,
double wNsqr,
int rvar);
37 int npol,
int mons,
double *rel,
double *fk,
double wNsqr,
int rvar);
38 extern "C" void wGcd(
int *
x,
int n);
40 static void wDimensions(poly*
s,
int sl,
int *lpol,
int *npol,
int *mons)
46 for (i = 0; i <= sl; i++)
70 static void wInit(poly*
s,
int sl,
int mons,
int *
A,
const ring
R)
72 int n, a,
i,
j, *
B, *C;
78 a = (n + 1) *
sizeof(
int);
80 for (i = 0; i <= sl; i++)
91 for (j = 0; j < n; j++)
102 for (j = 0; j < n; j++)
114 void wCall(poly*
s,
int sl,
int *
x,
double wNsqr,
const ring
R)
116 int n, q, npol, mons,
i;
117 int *
A, *xopt, *lpol, *degw;
118 double f1, fx, eps, *rel;
122 lpol = (
int * )
omAlloc((sl + 1) *
sizeof(int));
125 for (i = n; i!=0; i--)
134 q = (n + 1) * mons *
sizeof(
int);
136 wInit(s, sl, mons, A, R);
137 degw = A + (n * mons);
138 memset(degw, 0, mons *
sizeof(
int));
139 for (i = n; i!=0; i--)
141 wNorm(degw, lpol, npol, rel);
142 f1 = (*wFunctional)(degw, lpol, npol, rel, (double)1.0, wNsqr);
145 fx = (double)2.0 * eps;
146 memset(x, 0, (n + 1) *
sizeof(
int));
149 memcpy(x + 1, xopt + 1, n *
sizeof(
int));
150 memset(degw, 0, mons *
sizeof(
int));
151 for (i = n; i!=0; i--)
160 for (i = n; i!=0; i--)
199 x = (
int * )
omAlloc(2 * (n + 1) *
sizeof(int));
200 wCall(s, sl, x, (
double)2.0 / (
double)n, R);
201 for (i = n; i!=0; i--)
202 eweight[i] = x[i + n + 1];
214 for(i=len;i>0;i--) s[i]=(*iv)[i-1];
228 for (i=
rVar(r); i>0; i--)
242 for (i=
rVar(r); i>0; i--)
243 j += (
long)((int)(
p_GetExp(p,i,r) * w[
i]));
BOOLEAN rHasLocalOrMixedOrdering(const ring r)
const CanonicalForm int s
static int si_min(const int a, const int b)
void kEcartWeights(poly *s, int sl, short *eweight, const ring R)
long totaldegreeWecart(poly p, ring r)
short * iv2array(intvec *iv, const ring R)
static void p_GetExpV(poly p, int *ev, const ring r)
#define omFreeSize(addr, size)
static short rVar(const ring r)
#define rVar(r) (r->N)
void wAdd(int *A, int mons, int kn, int xx, int rvar)
long totaldegreeWecart_IV(poly p, ring r, const short *w)
static long p_GetExp(const poly p, const unsigned long iBitmask, const int VarOffset)
get a single variable exponent : the integer VarOffset encodes:
static int max(int a, int b)
void wSecondSearch(int *A, int *x, int *lpol, int npol, int mons, double *rel, double *fk, double wNsqr, int rvar)
double(* wFunctional)(int *degw, int *lpol, int npol, double *rel, double wx, double wNsqr)
double wFunctionalMora(int *degw, int *lpol, int npol, double *rel, double wx, double wNsqr)
static void wDimensions(poly *s, int sl, int *lpol, int *npol, int *mons)
static void wInit(poly *s, int sl, int mons, int *A, const ring R)
void wCall(poly *s, int sl, int *x, double wNsqr, const ring R)
long maxdegreeWecart(poly p, int *l, ring r)
void wFirstSearch(int *A, int *x, int mons, int *lpol, int npol, double *rel, double *fopt, double wNsqr, int rvar)
void wNorm(int *degw, int *lpol, int npol, double *rel)
double wFunctionalBuch(int *degw, int *lpol, int npol, double *rel, double wx, double wNsqr)