mp_permmatrix Class Reference

Public Member Functions
	mp_permmatrix ()
	mp_permmatrix (matrix, ring)
	mp_permmatrix (mp_permmatrix *)
	~mp_permmatrix ()
int	mpGetRow ()
int	mpGetCol ()
int	mpGetRdim ()
int	mpGetCdim ()
int	mpGetSign ()
void	mpSetSearch (int s)
void	mpSaveArray ()
poly	mpGetElem (int, int)
void	mpSetElem (poly, int, int)
void	mpDelElem (int, int)
void	mpElimBareiss (poly)
int	mpPivotBareiss (row_col_weight *)
int	mpPivotRow (row_col_weight *, int)
void	mpToIntvec (intvec *)
void	mpRowReorder ()
void	mpColReorder ()

Private Member Functions
void	mpInitMat ()
poly *	mpRowAdr (int r)
poly *	mpColAdr (int c)
void	mpRowWeight (float *)
void	mpColWeight (float *)
void	mpRowSwap (int, int)
void	mpColSwap (int, int)

Private Attributes
int	a_m
int	a_n
int	s_m
int	s_n
int	sign
int	piv_s
int *	qrow
int *	qcol
poly *	Xarray
ring	_R

Detailed Description

Definition at line 829 of file matpol.cc.

Constructor & Destructor Documentation

mp_permmatrix::mp_permmatrix ( )

inline

Definition at line 846 of file matpol.cc.

: a_m(0) {}

mp_permmatrix::mp_permmatrix	(	matrix	A,
		ring	R
	)

Definition at line 867 of file matpol.cc.

                                             : sign(1)
{
  a_m = A->nrows;
  a_n = A->ncols;
  this->mpInitMat();
  Xarray = A->m;
  _R=R;
}

mp_permmatrix::mp_permmatrix

(

mp_permmatrix *

M

)

Definition at line 876 of file matpol.cc.

{
  poly p, *athis, *aM;
  int i, j;

  _R=M->_R;
  a_m = M->s_m;
  a_n = M->s_n;
  sign = M->sign;
  this->mpInitMat();
  Xarray = (poly *)omAlloc0(a_m*a_n*sizeof(poly));
  for (i=a_m-1; i>=0; i--)
  {
    athis = this->mpRowAdr(i);
    aM = M->mpRowAdr(i);
    for (j=a_n-1; j>=0; j--)
    {
      p = aM[M->qcol[j]];
      if (p)
      {
        athis[j] = p_Copy(p,_R);
      }
    }
  }
}

mp_permmatrix::~mp_permmatrix

(

)

Definition at line 902 of file matpol.cc.

{
  int k;

  if (a_m != 0)
  {
    omFreeSize((ADDRESS)qrow,a_m*sizeof(int));
    omFreeSize((ADDRESS)qcol,a_n*sizeof(int));
    if (Xarray != NULL)
    {
      for (k=a_m*a_n-1; k>=0; k--)
        p_Delete(&Xarray[k],_R);
      omFreeSize((ADDRESS)Xarray,a_m*a_n*sizeof(poly));
    }
  }
}

Member Function Documentation

poly* mp_permmatrix::mpColAdr ( int  c )

inlineprivate

Definition at line 839 of file matpol.cc.

  { return &(Xarray[qcol[c]]); }

void mp_permmatrix::mpColReorder

(

)

Definition at line 1002 of file matpol.cc.

{
  int k, j, j1, j2;

  if (a_n > a_m)
    k = a_n - a_m;
  else
    k = 0;
  for (j=a_n-1; j>=k; j--)
  {
    j1 = qcol[j];
    if (j1 != j)
    {
      this->mpColSwap(j1, j);
      j2 = 0;
      while (qcol[j2] != j) j2++;
      qcol[j2] = j1;
    }
  }
}

void mp_permmatrix::mpColSwap ( int  j1, int  j2  )

private

Definition at line 975 of file matpol.cc.

{
   poly p, *a1, *a2;
   int i, k = a_n*a_m;

   a1 = &(Xarray[j1]);
   a2 = &(Xarray[j2]);
   for (i=0; i< k; i+=a_n)
   {
     p = a1[i];
     a1[i] = a2[i];
     a2[i] = p;
   }
}

void mp_permmatrix::mpColWeight ( float *  wcol )

private

Definition at line 921 of file matpol.cc.

{
  poly p, *a;
  int i, j;
  float count;

  for (j=s_n; j>=0; j--)
  {
    a = this->mpColAdr(j);
    count = 0.0;
    for(i=s_m; i>=0; i--)
    {
      p = a[a_n*qrow[i]];
      if (p)
        count += mp_PolyWeight(p,_R);
    }
    wcol[j] = count;
  }
}

void mp_permmatrix::mpDelElem	(	int	,
		int
	)

void mp_permmatrix::mpElimBareiss

(

poly

div

)

Definition at line 1148 of file matpol.cc.

{
  poly piv, elim, q1, q2, *ap, *a;
  int i, j, jj;

  ap = this->mpRowAdr(s_m);
  piv = ap[qcol[s_n]];
  for(i=s_m-1; i>=0; i--)
  {
    a = this->mpRowAdr(i);
    elim = a[qcol[s_n]];
    if (elim != NULL)
    {
      elim = p_Neg(elim,_R);
      for (j=s_n-1; j>=0; j--)
      {
        q2 = NULL;
        jj = qcol[j];
        if (ap[jj] != NULL)
        {
          q2 = SM_MULT(ap[jj], elim, div,_R);
          if (a[jj] != NULL)
          {
            q1 = SM_MULT(a[jj], piv, div,_R);
            p_Delete(&a[jj],_R);
            q2 = p_Add_q(q2, q1, _R);
          }
        }
        else if (a[jj] != NULL)
        {
          q2 = SM_MULT(a[jj], piv, div, _R);
        }
        if ((q2!=NULL) && div)
          SM_DIV(q2, div, _R);
        a[jj] = q2;
      }
      p_Delete(&a[qcol[s_n]], _R);
    }
    else
    {
      for (j=s_n-1; j>=0; j--)
      {
        jj = qcol[j];
        if (a[jj] != NULL)
        {
          q2 = SM_MULT(a[jj], piv, div, _R);
          p_Delete(&a[jj], _R);
          if (div)
            SM_DIV(q2, div, _R);
          a[jj] = q2;
        }
      }
    }
  }
}

int mp_permmatrix::mpGetCdim ( )

inline

Definition at line 853 of file matpol.cc.

{ return s_n; }

int mp_permmatrix::mpGetCol

(

)

poly mp_permmatrix::mpGetElem	(	int	r,
		int	c
	)

Definition at line 1140 of file matpol.cc.

{
  return Xarray[a_n*qrow[r]+qcol[c]];
}

int mp_permmatrix::mpGetRdim ( )

inline

Definition at line 852 of file matpol.cc.

{ return s_m; }

int mp_permmatrix::mpGetRow

(

)

int mp_permmatrix::mpGetSign ( )

inline

Definition at line 854 of file matpol.cc.

{ return sign; }

void mp_permmatrix::mpInitMat ( )

private

Definition at line 989 of file matpol.cc.

{
  int k;

  s_m = a_m;
  s_n = a_n;
  piv_s = 0;
  qrow = (int *)omAlloc(a_m*sizeof(int));
  qcol = (int *)omAlloc(a_n*sizeof(int));
  for (k=a_m-1; k>=0; k--) qrow[k] = k;
  for (k=a_n-1; k>=0; k--) qcol[k] = k;
}

int mp_permmatrix::mpPivotBareiss

(

row_col_weight *

C

)

Definition at line 1063 of file matpol.cc.

{
  poly p, *a;
  int i, j, iopt, jopt;
  float sum, f1, f2, fo, r, ro, lp;
  float *dr = C->wrow, *dc = C->wcol;

  fo = 1.0e20;
  ro = 0.0;
  iopt = jopt = -1;

  s_n--;
  s_m--;
  if (s_m == 0)
    return 0;
  if (s_n == 0)
  {
    for(i=s_m; i>=0; i--)
    {
      p = this->mpRowAdr(i)[qcol[0]];
      if (p)
      {
        f1 = mp_PolyWeight(p,_R);
        if (f1 < fo)
        {
          fo = f1;
          if (iopt >= 0)
            p_Delete(&(this->mpRowAdr(iopt)[qcol[0]]),_R);
          iopt = i;
        }
        else
          p_Delete(&(this->mpRowAdr(i)[qcol[0]]),_R);
      }
    }
    if (iopt >= 0)
      mpReplace(iopt, s_m, sign, qrow);
    return 0;
  }
  this->mpRowWeight(dr);
  this->mpColWeight(dc);
  sum = 0.0;
  for(i=s_m; i>=0; i--)
    sum += dr[i];
  for(i=s_m; i>=0; i--)
  {
    r = dr[i];
    a = this->mpRowAdr(i);
    for(j=s_n; j>=0; j--)
    {
      p = a[qcol[j]];
      if (p)
      {
        lp = mp_PolyWeight(p,_R);
        ro = r - lp;
        f1 = ro * (dc[j]-lp);
        if (f1 != 0.0)
        {
          f2 = lp * (sum - ro - dc[j]);
          f2 += f1;
        }
        else
          f2 = lp-r-dc[j];
        if (f2 < fo)
        {
          fo = f2;
          iopt = i;
          jopt = j;
        }
      }
    }
  }
  if (iopt < 0)
    return 0;
  mpReplace(iopt, s_m, sign, qrow);
  mpReplace(jopt, s_n, sign, qcol);
  return 1;
}

int mp_permmatrix::mpPivotRow	(	row_col_weight *	,
		int
	)

poly* mp_permmatrix::mpRowAdr ( int  r )

inlineprivate

Definition at line 837 of file matpol.cc.

  { return &(Xarray[a_n*qrow[r]]); }

void mp_permmatrix::mpRowReorder

(

)

Definition at line 1023 of file matpol.cc.

{
  int k, i, i1, i2;

  if (a_m > a_n)
    k = a_m - a_n;
  else
    k = 0;
  for (i=a_m-1; i>=k; i--)
  {
    i1 = qrow[i];
    if (i1 != i)
    {
      this->mpRowSwap(i1, i);
      i2 = 0;
      while (qrow[i2] != i) i2++;
      qrow[i2] = i1;
    }
  }
}

void mp_permmatrix::mpRowSwap ( int  i1, int  i2  )

private

Definition at line 960 of file matpol.cc.

{
   poly p, *a1, *a2;
   int j;

   a1 = &(Xarray[a_n*i1]);
   a2 = &(Xarray[a_n*i2]);
   for (j=a_n-1; j>= 0; j--)
   {
     p = a1[j];
     a1[j] = a2[j];
     a2[j] = p;
   }
}

void mp_permmatrix::mpRowWeight ( float *  wrow )

private

Definition at line 940 of file matpol.cc.

{
  poly p, *a;
  int i, j;
  float count;

  for (i=s_m; i>=0; i--)
  {
    a = this->mpRowAdr(i);
    count = 0.0;
    for(j=s_n; j>=0; j--)
    {
      p = a[qcol[j]];
      if (p)
        count += mp_PolyWeight(p,_R);
    }
    wrow[i] = count;
  }
}

void mp_permmatrix::mpSaveArray ( )

inline

Definition at line 856 of file matpol.cc.

{ Xarray = NULL; }

void mp_permmatrix::mpSetElem	(	poly	,
		int	,
		int
	)

void mp_permmatrix::mpSetSearch

(

int

s

)

void mp_permmatrix::mpToIntvec

(

intvec *

)

Field Documentation

ring mp_permmatrix::_R

private

Definition at line 835 of file matpol.cc.

int mp_permmatrix::a_m

private

Definition at line 832 of file matpol.cc.

int mp_permmatrix::a_n

private

Definition at line 832 of file matpol.cc.

int mp_permmatrix::piv_s

private

Definition at line 832 of file matpol.cc.

int * mp_permmatrix::qcol

private

Definition at line 833 of file matpol.cc.

int* mp_permmatrix::qrow

private

Definition at line 833 of file matpol.cc.

int mp_permmatrix::s_m

private

Definition at line 832 of file matpol.cc.

int mp_permmatrix::s_n

private

Definition at line 832 of file matpol.cc.

int mp_permmatrix::sign

private

Definition at line 832 of file matpol.cc.

poly* mp_permmatrix::Xarray

private

Definition at line 834 of file matpol.cc.

---

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

• libpolys/polys/matpol.cc