/* ------------------------------------ */
/*  Save as :   wgj4.h               */
/* ------------------------------------ */
void gj4_T_mZ(
double **Ab
)
{
int r = R1;
int c = C1;

    while((r<Ab[R_SIZE][C0]) &&
          (c<Ab[C_SIZE_A][C0])
         )
       {
        if(pivotbest4Ab_r0_mZ(Ab,r,c)>ERROR_E)  
           zero_under_pivot_gjAb_mZ(Ab,r,c);
        r++;
        c+=C2;
       }
}
/* ------------------------------------ */
void gj4_mZ(
double **Ab
)
{
double **Ab_T = c_mZ(Ab,i_duplicateAbR0_Ab_c_A_c_mZ(

                             Ab[R_SIZE][C0],
                             Ab[C_SIZE][C0],
                             Ab[C_SIZE_A][C0]));
double **new_Ab;
int r_not_0;

  gj4_T_mZ(Ab_T);
  sort4_r_coef_r0_mZ(Ab_T);
  r_not_0 =  is_system_consistent_mZ(Ab,Ab_T);

  if(r_not_0>2)
     {
       new_Ab = i_duplicateAbR0_Ab_c_A_c_mZ(
                             r_not_0,
                             Ab[C_SIZE][C0],
                             Ab[C_SIZE_A][C0]);
       c_n_r_mZ(Ab,r_not_0,new_Ab);
       gj2_mZ(new_Ab);
       m0_mZ(Ab);
       c_mZ(new_Ab,Ab);
       f_mZ(new_Ab);
     }
   else printf(" The system has only an equation !!\n\n");

f_mZ(Ab_T);
}
/* ------------------------------------ */
/* ------------------------------------ */