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 825 of file matpol.cc.

Constructor & Destructor Documentation

mp_permmatrix() [1/3]

mp_permmatrix::mp_permmatrix ( )

inline

Definition at line 842 of file matpol.cc.

: a_m(0) {}

mp_permmatrix() [2/3]

mp_permmatrix::mp_permmatrix	(	matrix	A,
		ring	R
	)

Definition at line 863 of file matpol.cc.

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

mp_permmatrix() [3/3]

mp_permmatrix::mp_permmatrix

(

mp_permmatrix *

M

)

Definition at line 872 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::~mp_permmatrix

(

)

Definition at line 898 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

mpColAdr()

poly* mp_permmatrix::mpColAdr ( int  c )

inlineprivate

Definition at line 835 of file matpol.cc.

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

mpColReorder()

void mp_permmatrix::mpColReorder

(

)

Definition at line 998 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;
    }
  }
}

mpColSwap()

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

private

Definition at line 971 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;
   }
}

mpColWeight()

void mp_permmatrix::mpColWeight ( float *  wcol )

private

Definition at line 917 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;
  }
}

mpDelElem()

void mp_permmatrix::mpDelElem	(	int	,
		int
	)

mpElimBareiss()

void mp_permmatrix::mpElimBareiss

(

poly

div

)

Definition at line 1144 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;
        }
      }
    }
  }
}

mpGetCdim()

int mp_permmatrix::mpGetCdim ( )

inline

Definition at line 849 of file matpol.cc.

{ return s_n; }

mpGetCol()

int mp_permmatrix::mpGetCol

(

)

mpGetElem()

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

Definition at line 1136 of file matpol.cc.

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

mpGetRdim()

int mp_permmatrix::mpGetRdim ( )

inline

Definition at line 848 of file matpol.cc.

{ return s_m; }

mpGetRow()

int mp_permmatrix::mpGetRow

(

)

mpGetSign()

int mp_permmatrix::mpGetSign ( )

inline

Definition at line 850 of file matpol.cc.

{ return sign; }

mpInitMat()

void mp_permmatrix::mpInitMat ( )

private

Definition at line 985 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;
}

mpPivotBareiss()

int mp_permmatrix::mpPivotBareiss

(

row_col_weight *

C

)

Definition at line 1059 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;
}

mpPivotRow()

int mp_permmatrix::mpPivotRow	(	row_col_weight *	,
		int
	)

mpRowAdr()

poly* mp_permmatrix::mpRowAdr ( int  r )

inlineprivate

Definition at line 833 of file matpol.cc.

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

mpRowReorder()

void mp_permmatrix::mpRowReorder

(

)

Definition at line 1019 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;
    }
  }
}

mpRowSwap()

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

private

Definition at line 956 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;
   }
}

mpRowWeight()

void mp_permmatrix::mpRowWeight ( float *  wrow )

private

Definition at line 936 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;
  }
}

mpSaveArray()

void mp_permmatrix::mpSaveArray ( )

inline

Definition at line 852 of file matpol.cc.

{ Xarray = NULL; }

mpSetElem()

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

mpSetSearch()

void mp_permmatrix::mpSetSearch

(

int

s

)

mpToIntvec()

void mp_permmatrix::mpToIntvec

(

intvec *

)

Field Documentation

_R

ring mp_permmatrix::_R

private

Definition at line 831 of file matpol.cc.

a_m

int mp_permmatrix::a_m

private

Definition at line 828 of file matpol.cc.

a_n

int mp_permmatrix::a_n

private

Definition at line 828 of file matpol.cc.

piv_s

int mp_permmatrix::piv_s

private

Definition at line 828 of file matpol.cc.

qcol

int * mp_permmatrix::qcol

private

Definition at line 829 of file matpol.cc.

qrow

int* mp_permmatrix::qrow

private

Definition at line 829 of file matpol.cc.

s_m

int mp_permmatrix::s_m

private

Definition at line 828 of file matpol.cc.

s_n

int mp_permmatrix::s_n

private

Definition at line 828 of file matpol.cc.

sign

int mp_permmatrix::sign

private

Definition at line 828 of file matpol.cc.

Xarray

poly* mp_permmatrix::Xarray

private

Definition at line 830 of file matpol.cc.

---

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

• libpolys/polys/matpol.cc