/* xcpivrow.h                   freeware                   [[Email Removed]]  */


/* ------------------------------------- FUNCTION --------------------------  */
/* Do    :                                                                    */
/*                                                                            */
/* Call  :                                                                    */
/* Debug :  printf("%.0f/%.0f,",x.numer,x.denom);                             */
/* -------------------------------------------------------------------------- */
pivotgoodplaceF(
pmatrix  m
)
{
     int i;
     int j;
     int k;
fraction f;

 for ( i = 0 ; i < m->rows ; i++)
  {
   for ( j = i*TWOCOL; j < m->cols; j++,j++)
    {
     if( !( *(m->pblock+i *m->cols+j) ) )
       {
        k = i;
         do
          {
           if( *(m->pblock+k *m->cols+j) )
             {
              swaprowF(m,i,k);
              k=m->rows;
             }
           k++;
          }
         while(k < m->rows);
       }
    }
  }
}



/* ------------------------------------- FUNCTION --------------------------  */
/* Do    : All pivot = 1.                                                     */
/*                                                                            */
/* Call  :                                                                    */
/* Debug :  printf("%.0f/%.0f,",x.numer,x.denom);                             */
/* -------------------------------------------------------------------------- */
pivot1F(
pmatrix  m
)
{
     int i;
     int j;
fraction f;

 for (i = 0 ; i < m->rows ; i++)
  {
   for ( j = 0 ; j < m->cols; j++,j++)
    {
     if( (*(m->pblock+i *m->cols+j)) )
       {
        f.numer = *(m->pblock+i *m->cols+j);
        f.denom = *(m->pblock+i *m->cols+j+1);

         f = invF(f);

          mulrowF(m,i,f);
        frowminiF(m,i  );

        j = m->cols;
       }
    }
  }
}