/* ------------------------------------ */
/*  Save as :   we_det1.h                */
/* ------------------------------------ */
nb_Z det_Z(
double **A
)
{
double **T = c_mZ(A,i_duplicate_mZ(A[R_SIZE][C0],
                                   A[C_SIZE][C0]));
nb_Z deter=i_Z(1.,0.);
int r;
int c;

    for( r=R1, c=C1; r<T[R_SIZE][C0]; r++,c+=C2)
      {
        deter = smul_Z( pivotbest_mZ(T,r,c),
                        deter);

        deter = mul_Z( deter,
                       zero_under_pivot_mZ(T,r,c));
       }

   f_mZ(T);

return(deter);
}
/* ------------------------------------ */
/* ------------------------------------ */