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


/* --------------------------------- FUNCTION ------------------------------  */
/* Do    :                                                                    */
/*                                                                            */
/* Call  :                                                                    */
/* Debug :                                                                    */
/* -------------------------------------------------------------------------- */
void mdiag0F(
pmatrix m,
int maxnumer,
int maxdenom)
{
int i;
int j;
int row0;
int row;

     row0 = rndp0I(m->rows);

  if (m->rows * TWOCOL != m->cols)
    {
     printf("\n mdiag0F() error - matrix must be square");
     printf("\n Press return to continue");
     getchar();
     exit(1);
    }
  for (   i = 0 ; i < m->rows ; i++)
    {
    for ( j = 0 ; j < m->cols ; j++,j++)
       {
       *(m->pblock+i *m->cols+j  ) = 0;
       *(m->pblock+i *m->cols+j+1) = 1;

       if( (i*TWOCOL) == j)
         {
          *(m->pblock+i *m->cols+j)     = rndI(maxnumer);
          *(m->pblock+i *m->cols+j+1) = rndpI(maxdenom);

          if( i == row0)
           {
            *(m->pblock+i *m->cols+j)   = 0;
            *(m->pblock+i *m->cols+j+1) = 1;
           }
         }
       }
    }

    for ( row = 0; row < m->rows ; row++)
     {
      frowminiF(m,row);                                 /* a) fraction mini   */
     }
}


/* --------------------------------- FUNCTION ------------------------------  */
/* Do    :                                                                    */
/*                                                                            */
/* Call  :                                                                    */
/* Debug :                                                                    */
/* -------------------------------------------------------------------------- */
void mdiag0I(
pmatrix m,
int maxnumer)
{
int i;
int j;
int row0;

     row0 = rndp0I(m->rows);

  if (m->rows * TWOCOL != m->cols)
    {
     printf("\n mid() error - matrix must be square");
     printf("\n Press return to continue");
     getchar();
     exit(1);
    }
  for (   i = 0 ; i < m->rows ; i++)
    {
    for ( j = 0 ; j < m->cols ; j++,j++)
       {
       *(m->pblock+i *m->cols+j  ) = 0;
       *(m->pblock+i *m->cols+j+1) = 1;

       if( (i*TWOCOL) == j)
         {
          *(m->pblock+i *m->cols+j) = rndI(maxnumer);
          if(i==row0)
            {
             *(m->pblock+i *m->cols+j) = 0;
            }
         }
       }
    }
}