mp_permmatrix Class Reference

Public Member Functions
	mp_permmatrix ()
	mp_permmatrix (matrix, const 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 ()
	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)
poly *	mpColAdr (int)
void	mpRowWeight (float *)
void	mpColWeight (float *)
void	mpRowSwap (int, int)
void	mpColSwap (int, int)
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
ring	_R

Detailed Description

Definition at line 193 of file linalg_from_matpol.cc.

Constructor & Destructor Documentation

mp_permmatrix::mp_permmatrix ( )

inline

Definition at line 209 of file linalg_from_matpol.cc.

: a_m(0), R(NULL) {}

mp_permmatrix::mp_permmatrix	(	matrix	A,
		const ring	r
	)

Definition at line 372 of file linalg_from_matpol.cc.

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

mp_permmatrix::mp_permmatrix

(

mp_permmatrix *

M

)

Definition at line 380 of file linalg_from_matpol.cc.

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

  a_m = M->s_m;
  a_n = M->s_n;
  sign = M->sign;
  R = M->R;

  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 407 of file linalg_from_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));
    }
  }
}

mp_permmatrix::mp_permmatrix ( )

inline

Definition at line 847 of file matpol.cc.

: a_m(0) {}

mp_permmatrix::mp_permmatrix	(	matrix	,
		ring
	)

mp_permmatrix::mp_permmatrix

(

mp_permmatrix *

)

mp_permmatrix::~mp_permmatrix

(

)

Member Function Documentation

poly * mp_permmatrix::mpColAdr ( int  c )

private

Definition at line 733 of file linalg_from_matpol.cc.

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

poly* mp_permmatrix::mpColAdr ( int  c )

inlineprivate

Definition at line 840 of file matpol.cc.

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

void mp_permmatrix::mpColReorder

(

)

Definition at line 693 of file linalg_from_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::mpColReorder

(

)

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

private

Definition at line 793 of file linalg_from_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::mpColSwap ( int  , int    )

private

void mp_permmatrix::mpColWeight ( float *  wcol )

private

Definition at line 758 of file linalg_from_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::mpColWeight ( float *  )

private

void mp_permmatrix::mpDelElem	(	int	r,
		int	c
	)

Definition at line 444 of file linalg_from_matpol.cc.

{
  p_Delete(&Xarray[a_n*qrow[r]+qcol[c]], R);
}

void mp_permmatrix::mpDelElem	(	int	,
		int
	)

void mp_permmatrix::mpElimBareiss

(

poly

div

)

Definition at line 452 of file linalg_from_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;
        }
      }
    }
  }
}

void mp_permmatrix::mpElimBareiss

(

poly

)

int mp_permmatrix::mpGetCdim

(

)

Definition at line 426 of file linalg_from_matpol.cc.

{ return s_n; }

int mp_permmatrix::mpGetCdim ( )

inline

Definition at line 854 of file matpol.cc.

{ return s_n; }

int mp_permmatrix::mpGetCol

(

)

Definition at line 813 of file linalg_from_matpol.cc.

{
  return qcol[s_n];
}

int mp_permmatrix::mpGetCol

(

)

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

Definition at line 434 of file linalg_from_matpol.cc.

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

poly mp_permmatrix::mpGetElem	(	int	,
		int
	)

int mp_permmatrix::mpGetRdim

(

)

Definition at line 424 of file linalg_from_matpol.cc.

{ return s_m; }

int mp_permmatrix::mpGetRdim ( )

inline

Definition at line 853 of file matpol.cc.

{ return s_m; }

int mp_permmatrix::mpGetRow

(

)

Definition at line 808 of file linalg_from_matpol.cc.

{
  return qrow[s_m];
}

int mp_permmatrix::mpGetRow

(

)

int mp_permmatrix::mpGetSign

(

)

Definition at line 428 of file linalg_from_matpol.cc.

{ return sign; }

int mp_permmatrix::mpGetSign ( )

inline

Definition at line 855 of file matpol.cc.

{ return sign; }

void mp_permmatrix::mpInitMat ( )

private

Definition at line 715 of file linalg_from_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;
}

void mp_permmatrix::mpInitMat ( )

private

int mp_permmatrix::mpPivotBareiss

(

row_col_weight *

C

)

Definition at line 511 of file linalg_from_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::mpPivotBareiss

(

row_col_weight *

)

int mp_permmatrix::mpPivotRow	(	row_col_weight *	C,
		int	row
	)

Definition at line 592 of file linalg_from_matpol.cc.

{
  poly p, *a;
  int 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)
  {
    p = this->mpRowAdr(row)[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 = row;
      }
      else
        p_Delete(&(this->mpRowAdr(row)[qcol[0]]), R);
    }
    if (iopt >= 0)
      mpReplace(iopt, s_m, sign, qrow);
    return 0;
  }
  this->mpRowWeight(dr);
  this->mpColWeight(dc);
  sum = 0.0;
  for(j=s_m; j>=0; j--)
    sum += dr[j];
  r = dr[row];
  a = this->mpRowAdr(row);
  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 = row;
        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 )

private

Definition at line 728 of file linalg_from_matpol.cc.

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

poly* mp_permmatrix::mpRowAdr ( int  r )

inlineprivate

Definition at line 838 of file matpol.cc.

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

void mp_permmatrix::mpRowReorder

(

)

Definition at line 672 of file linalg_from_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::mpRowReorder

(

)

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

private

Definition at line 778 of file linalg_from_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::mpRowSwap ( int  , int    )

private

void mp_permmatrix::mpRowWeight ( float *  wrow )

private

Definition at line 738 of file linalg_from_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::mpRowWeight ( float *  )

private

void mp_permmatrix::mpSaveArray

(

)

Definition at line 432 of file linalg_from_matpol.cc.

{ Xarray = NULL; }

void mp_permmatrix::mpSaveArray ( )

inline

Definition at line 857 of file matpol.cc.

{ Xarray = NULL; }

void mp_permmatrix::mpSetElem	(	poly	p,
		int	r,
		int	c
	)

Definition at line 439 of file linalg_from_matpol.cc.

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

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

void mp_permmatrix::mpSetSearch

(

int

s

)

Definition at line 430 of file linalg_from_matpol.cc.

{ piv_s = s; }

void mp_permmatrix::mpSetSearch

(

int

s

)

void mp_permmatrix::mpToIntvec

(

intvec *

v

)

Definition at line 664 of file linalg_from_matpol.cc.

{
  int i;

  for (i=v->rows()-1; i>=0; i--)
    (*v)[i] = qcol[i]+1;
}

void mp_permmatrix::mpToIntvec

(

intvec *

)

Field Documentation

ring mp_permmatrix::_R

private

Definition at line 836 of file matpol.cc.

int mp_permmatrix::a_m

private

Definition at line 196 of file linalg_from_matpol.cc.

int mp_permmatrix::a_n

private

Definition at line 196 of file linalg_from_matpol.cc.

int mp_permmatrix::piv_s

private

Definition at line 196 of file linalg_from_matpol.cc.

int * mp_permmatrix::qcol

private

Definition at line 197 of file linalg_from_matpol.cc.

int * mp_permmatrix::qrow

private

Definition at line 197 of file linalg_from_matpol.cc.

ring mp_permmatrix::R

private

Definition at line 199 of file linalg_from_matpol.cc.

int mp_permmatrix::s_m

private

Definition at line 196 of file linalg_from_matpol.cc.

int mp_permmatrix::s_n

private

Definition at line 196 of file linalg_from_matpol.cc.

int mp_permmatrix::sign

private

Definition at line 196 of file linalg_from_matpol.cc.

poly * mp_permmatrix::Xarray

private

Definition at line 198 of file linalg_from_matpol.cc.

---

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

• libpolys/polys/linalg_from_matpol.cc
• libpolys/polys/matpol.cc