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

Constructor & Destructor Documentation

mp_permmatrix() [1/3]

mp_permmatrix::mp_permmatrix ( )

inline

Definition at line 935 of file matpol.cc.

: a_m(0) {}

mp_permmatrix() [2/3]

mp_permmatrix::mp_permmatrix	(	matrix	A,
		ring	R
	)

Definition at line 956 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 965 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 991 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 928 of file matpol.cc.

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

mpColReorder()

void mp_permmatrix::mpColReorder

(

)

Definition at line 1091 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 1064 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 1010 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 1237 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 942 of file matpol.cc.

{ return s_n; }

mpGetCol()

int mp_permmatrix::mpGetCol

(

)

mpGetElem()

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

Definition at line 1229 of file matpol.cc.

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

mpGetRdim()

int mp_permmatrix::mpGetRdim ( )

inline

Definition at line 941 of file matpol.cc.

{ return s_m; }

mpGetRow()

int mp_permmatrix::mpGetRow

(

)

mpGetSign()

int mp_permmatrix::mpGetSign ( )

inline

Definition at line 943 of file matpol.cc.

{ return sign; }

mpInitMat()

void mp_permmatrix::mpInitMat ( )

private

Definition at line 1078 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 1152 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 926 of file matpol.cc.

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

mpRowReorder()

void mp_permmatrix::mpRowReorder

(

)

Definition at line 1112 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 1049 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 1029 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 945 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 924 of file matpol.cc.

a_m

int mp_permmatrix::a_m

private

Definition at line 921 of file matpol.cc.

a_n

int mp_permmatrix::a_n

private

Definition at line 921 of file matpol.cc.

piv_s

int mp_permmatrix::piv_s

private

Definition at line 921 of file matpol.cc.

qcol

int * mp_permmatrix::qcol

private

Definition at line 922 of file matpol.cc.

qrow

int* mp_permmatrix::qrow

private

Definition at line 922 of file matpol.cc.

s_m

int mp_permmatrix::s_m

private

Definition at line 921 of file matpol.cc.

s_n

int mp_permmatrix::s_n

private

Definition at line 921 of file matpol.cc.

sign

int mp_permmatrix::sign

private

Definition at line 921 of file matpol.cc.

Xarray

poly* mp_permmatrix::Xarray

private

Definition at line 923 of file matpol.cc.

---

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

• libpolys/polys/matpol.cc