yddeter1.h

[[ Complex Matrices : Language c]]

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


/* --------------------------------- FUNCTION ------------------------------  */
/* Do    :                                                                    */
/*                                                                            */
/* -------------------------------------------------------------------------- */
complexF det_zF(
PmzF mA)
{
int      row1;
int      row2;
complexF z_deter;
complexF z_T;

double pbAT[MXR][MXC*C4];mzF mAT ={MXR,MXC*C4,&pbAT[0][0]};

            mAT.rows=(mA->rows);   mAT.cols=(mA->rows)*C4;
/* -------------------------------------------------------------------------- */

  copy_mzF(mA, &mAT);

  z_deter = initD_zF(1,1,0,1);

   for(row1 = 0; row1<(mAT.rows-1) ; row1++)
      {
       z_T      = pivotnt0_mzF(&mAT,row1);
       z_deter  = mul_zF(z_deter,z_T);
       for(row2 = row1+1; row2<mAT.rows ; row2++)
            z_deter = pivotwork_mzF(&mAT,row1,row2,z_deter);
      }

/*------------------------------------------------------------------ Diagonal */
    for (row1=0 ; row1<mAT.rows ; row1++)
     {
       z_T.r.n = *(mAT.pb+row1 *mAT.cols+row1*C4+0);
       z_T.r.d = *(mAT.pb+row1 *mAT.cols+row1*C4+1);
       z_T.i.n = *(mAT.pb+row1 *mAT.cols+row1*C4+2);
       z_T.i.d = *(mAT.pb+row1 *mAT.cols+row1*C4+3);

       z_deter   = mul_zF(z_deter,z_T);
      }

 return(z_deter);
}