yecopyid.h

[[ Complex Matrices : Language c]]

/* xvcopyid.h                     free ware                [e-mail]  */

/*------------------------------- FUNCTION ---------------------------------- */
/*                                                                            */
/* -------------------------------------------------------------------------- */
void createAid_mzF(
PmzF mA,
PmzF mAID
)
{
int i;
int j;
int k = mA->rows * C4;
int n = mA->rows;

double pbID[MXR][MXC*C4];mzF mID={MXR,MXC*C4,&pbID[0][0]};

               mID.rows = n   ;
               mID.cols = n*C4;
/* ------------------------------------------------------------------ copy mA */
  for (i = 0 ; i < mA->rows ; i++)
    for ( j = 0; j < mA->cols ; j++)
          *(mAID->pb+i *mAID->cols+j) = *(mA->pb+i *mA->cols+j);
/* ----------------------------------------------------------------- copy mID */
  mid_mzF(&mID);

  for (i = 0 ; i < mA->rows ; i++)
    for ( j = 0; j < mA->cols ; j++)
          *(mAID->pb+i *mAID->cols+j+k) = *(mID.pb+i *mID.cols+j);
}

/*------------------------------- FUNCTION ---------------------------------- */
/*                                                                            */
/* -------------------------------------------------------------------------- */
void copyinv_mzF(
PmzF mAID,
PmzF minvA)
{
int i;
int j;
int k = mAID->rows * C4;

  for (i = 0 ; i < minvA->rows ; i++)
    for ( j = 0; j < minvA->cols ; j++)
          *(minvA->pb+i *minvA->cols+j) = *(mAID->pb+i *mAID->cols+k+j);
}