weight.h File Reference

#include "polys/monomials/ring.h"

Go to the source code of this file.

Functions
void	kEcartWeights (poly *s, int sl, short *eweight, const ring R)
long	maxdegreeWecart (poly p, int *l, ring r)
long	totaldegreeWecart (poly p, ring r)
long	totaldegreeWecart_IV (poly p, ring r, const int *w)
int *	iv2array (intvec *iv, const ring R)
double	wFunctionalBuch (int *degw, int *lpol, int npol, double *rel, double wx, double wNsqr)
void	wCall (poly *s, int sl, int *x, double wNsqr, const ring R)

Variables
EXTERN_VAR short *	ecartWeights
THREAD_VAR double(*	wFunctional )(int *degw, int *lpol, int npol, double *rel, double wx, double wNsqr)

Function Documentation

iv2array()

int * iv2array	(	intvec *	iv,
		const ring	R
	)

Definition at line 200 of file weight.cc.

{
  int *s=(int *)omAlloc0((rVar(R)+1)*sizeof(int));
  int len=0;
  if(iv!=NULL)
    len=si_min(iv->length(),rVar(R)); // usually: rVar(R)==length()
  int i;
  //for(i=pVariables;i>len;i--) s[i]=1;
  for(i=len;i>0;i--)               s[i]=(*iv)[i-1];
  return s;
}

kEcartWeights()

void kEcartWeights	(	poly *	s,
		int	sl,
		short *	eweight,
		const ring	R
	)

Definition at line 182 of file weight.cc.

{
  int  n, i;
  int  *x;
 
  *eweight = 0;
  n = rVar(R);
  if (rHasLocalOrMixedOrdering(R))
    wFunctional = wFunctionalMora;
  else
    wFunctional = wFunctionalBuch;
  x = (int * )omAlloc(2 * (n + 1) * sizeof(int));
  wCall(s, sl, x, (double)2.0 / (double)n, R);
  for (i = n; i!=0; i--)
    eweight[i] = x[i + n + 1];
  omFreeSize((ADDRESS)x, 2 * (n + 1) * sizeof(int));
}

maxdegreeWecart()

long maxdegreeWecart	(	poly	p,
		int *	l,
		ring	r
	)

Definition at line 247 of file weight.cc.

{
  short k=p_GetComp(p, r);
  int ll=1;
  long t,max;
 
  max=totaldegreeWecart(p, r);
  pIter(p);
  while ((p!=NULL) && (p_GetComp(p, r)==k))
  {
    t=totaldegreeWecart(p, r);
    if (t>max) max=t;
    ll++;
    pIter(p);
  }
  *l=ll;
  return max;
}

totaldegreeWecart()

long totaldegreeWecart	(	poly	p,
		ring	r
	)

Definition at line 217 of file weight.cc.

{
  int i;
  long j =0;
 
  for (i=rVar(r); i>0; i--)
    j += (int)(p_GetExp(p,i,r) * ecartWeights[i]);
  return  j;
}

totaldegreeWecart_IV()

long totaldegreeWecart_IV	(	poly	p,
		ring	r,
		const int *	w
	)

Definition at line 231 of file weight.cc.

{
  int i;
  long j =0;
 
  for (i=rVar(r); i>0; i--)
    j += (long)((int)(p_GetExp(p,i,r) * w[i]));
  return  j;
}

wCall()

void wCall	(	poly *	s,
		int	sl,
		int *	x,
		double	wNsqr,
		const ring	R
	)

Definition at line 108 of file weight.cc.

{
  int  n, q, npol, mons, i;
  int  *A, *xopt, *lpol, *degw;
  double  f1, fx, eps, *rel;
  void *adr;
 
  n = rVar(R);
  lpol = (int * )omAlloc((sl + 1) * sizeof(int));
  wDimensions(s, sl, lpol, &npol, &mons);
  xopt = x + (n + 1);
  for (i = n; i!=0; i--)
    xopt[i] = 1;
  if (mons==0)
  {
    omFreeSize((ADDRESS)lpol, (sl + 1) * sizeof(int));
    return;
  }
  adr = (void * )omAllocAligned(npol * sizeof(double));
  rel = (double*)adr;
  q = (n + 1) * mons * sizeof(int);
  A = (int * )omAlloc(q);
  wInit(s, sl, mons, A, R);
  degw = A + (n * mons);
  memset(degw, 0, mons * sizeof(int));
  for (i = n; i!=0; i--)
    wAdd(A, mons, i, 1, rVar(R));
  wNorm(degw, lpol, npol, rel);
  f1 = (*wFunctional)(degw, lpol, npol, rel, (double)1.0, wNsqr);
  if (TEST_OPT_PROT) Print("// %e\n",f1);
  eps = f1;
  fx = (double)2.0 * eps;
  memset(x, 0, (n + 1) * sizeof(int));
  wFirstSearch(A, x, mons, lpol, npol, rel, &fx, wNsqr, rVar(R));
  if (TEST_OPT_PROT) Print("// %e\n",fx);
  memcpy(x + 1, xopt + 1, n * sizeof(int));
  memset(degw, 0, mons * sizeof(int));
  for (i = n; i!=0; i--)
  {
    x[i] *= 16;
    wAdd(A, mons, i, x[i], rVar(R));
  }
  wSecondSearch(A, x, lpol, npol, mons, rel, &fx, wNsqr, rVar(R));
  if (TEST_OPT_PROT) Print("// %e\n",fx);
  if (fx >= eps)
  {
    for (i = n; i!=0; i--)
      xopt[i] = 1;
  }
  else
  {
    wGcd(xopt, n);
//    if (BTEST1(22))
//    {
//      f1 = fx + (double)0.1 * (f1 - fx);
//      wSimple(x, n);
//      memset(degw, 0, mons * sizeof(int));
//      for (i = n; i!=0; i--)
//        wAdd(A, mons, i, x[i], rVar(R));
//      eps = wPrWeight(x, n);
//      fx = (*wFunctional)(degw, lpol, npol, rel, eps);
//      if (fx < f1)
//      {
//        if (TEST_OPT_PROT) Print("// %e\n",fx);
//        memcpy(xopt + 1, x + 1, n * sizeof(int));
//      }
//    }
  }
  omFreeSize((ADDRESS)A, q);
  omFreeSize((ADDRESS)lpol, (sl + 1) * sizeof(int));
  omFreeSize((ADDRESS)adr, npol * sizeof(double));
}

wFunctionalBuch()

double wFunctionalBuch	(	int *	degw,
		int *	lpol,
		int	npol,
		double *	rel,
		double	wx,
		double	wNsqr
	)

Definition at line 78 of file weight0.cc.

{
  int  i, j, ecl, ecu, ec;
  int  *ex;
  double gfmax, ghom, pfmax;
  double *r;
 
  ex = degw;
  r = rel;
  gfmax = (double)0.0;
  ghom = (double)1.0;
  for (i = 0; i < npol; i++)
  {
    ecu = ecl = *ex++;
    for (j = lpol[i] - 1; j!=0 ; j--)
    {
      ec = *ex++;
      if (ec < ecl)
        ecl = ec;
      else if (ec > ecu)
        ecu = ec;
    }
    pfmax = (double)ecl / (double)ecu;
    if (pfmax < ghom)
      ghom = pfmax;
    gfmax += (double)(ecu * ecu) * (*r++);
  }
  if (ghom > (double)0.5)
    gfmax *= ((double)1.0 - (ghom * ghom)) / (double)0.75;
  return gfmax / pow(wx, wNsqr);
}

Variable Documentation

ecartWeights

EXTERN_VAR short* ecartWeights

Definition at line 12 of file weight.h.

wFunctional

THREAD_VAR double(* wFunctional) (int *degw, int *lpol, int npol, double *rel, double wx, double wNsqr)	(	int *	degw,
		int *	lpol,
		int	npol,
		double *	rel,
		double	wx,
		double	wNsqr
	)

Definition at line 26 of file weight.h.