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

/*------------------------------- FUNCTION ----------------------- copymAID() */
/* Debug : printmF(m);  getchar();                                            */
/* -------------------------------------------------------------------------- */
void copymAID (
pmatrix mA,
pmatrix mAID
)
{
int i;
int k;
int n = mA->rows;
int j = mA->rows;

double pbID[MXR][MXC*TWOCOL];matrix mID={MXR,MXC*TWOCOL,&pbID[0][0]};

               mID.rows = n       ;
               mID.cols = n*TWOCOL;
/* ------------------------------------------------------------------ copy mA */
    for (k = 0; k < mA->rows ; k++ )
    {
     for (i = 0; i < mA->rows ; i++ )
      {
            *(mAID->pblock+i *mAID->cols+k*TWOCOL   ) =
              *(mA->pblock+i   *mA->cols+k*TWOCOL   );
            *(mAID->pblock+i *mAID->cols+k*TWOCOL+1 ) =
              *(mA->pblock+i   *mA->cols+k*TWOCOL+1 );
      }
    }
/* ----------------------------------------------------------------- copy mID */

  mid(&mID);
    for (k = 0; k < mA->rows ; k++ )
    {
     for (i = 0; i < mA->rows ; i++ )
      {
            *(mAID->pblock+i *mAID->cols+j*TWOCOL   ) =
              *(mID.pblock+i   *mID.cols+k*TWOCOL   );
            *(mAID->pblock+i *mAID->cols+j*TWOCOL+1 ) =
              *(mID.pblock+i   *mID.cols+k*TWOCOL+1 );
      } j++;
    }
}


/*------------------------------- FUNCTION ------------------------ copyinv() */
/* Debug : printmF(m);  getchar();                                            */
/* -------------------------------------------------------------------------- */
void copyinv (
pmatrix mA,
pmatrix mB)
{
int i;
int j = mA->rows;
int k = 0;

    for (k = 0; k < mA->rows ; k++ )
    {
     for (i = 0; i < mA->rows ; i++ )
      {
            *(mB->pblock+i *mB->cols+k*TWOCOL   ) =
            *(mA->pblock+i *mA->cols+j*TWOCOL   );
            *(mB->pblock+i *mB->cols+k*TWOCOL+1 ) =
            *(mA->pblock+i *mA->cols+j*TWOCOL+1 );
      }j++;
    }
}