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

/* --------------------------------- FUNCTION ------------------------------  */
/* Do    :                                                                    */
/*                                                                            */
/* Call  :                                                                    */
/* Debug :                                                                    */
/* -------------------------------------------------------------------------- */
void msymetricF(
pmatrix m,
int maxnumer,
int maxdenom)
{
int i;
int j;
int row;
double Tnumer;
double Tdenom;


  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++)
       {
        Tnumer = rndI (maxnumer);
        Tdenom = rndpI(maxdenom);

       *(m->pblock+i     *m->cols+j)         = Tnumer;
       *(m->pblock+i     *m->cols+j+1)        = Tdenom;
       *(m->pblock+(j/2) *m->cols+i*TWOCOL)   = Tnumer;
       *(m->pblock+(j/2) *m->cols+i*TWOCOL+1) = Tdenom;

       }
    }

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

/* ------------------------------ FUNCTION ---------------------------- mid() */
/*                                                                            */
/* -------------------------------------------------------------------------- */
void msymetricI(
pmatrix m,
int maxnumer)
{
int i;
int j;
double Tnumer;


  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++)
       {
        Tnumer = rndI (maxnumer);

       *(m->pblock+i     *m->cols+j)          = Tnumer;
       *(m->pblock+i     *m->cols+j+1)        = 1;
       *(m->pblock+(j/2) *m->cols+i*TWOCOL)   = Tnumer;
       *(m->pblock+(j/2) *m->cols+i*TWOCOL+1) = 1;

       }
    }
}