/* ------------------------------------ */
/*  Save as :   vgjinv.h                   */
/* ------------------------------------ */
double **invgj_T_mR(
double **Ab
)
{
double pivot = 1.;

int r= R1;

    while( pivot &&
          (r<Ab[R_SIZE  ][C0])
         )
       {
        pivot = pivotbestAId_mR(Ab,r);
        
        if(fabs(pivot)>ERROR_E) 
                 zero_under_pivot_gj1Ab_mR(Ab,r);
        r++;
       }

    while( (r>R1) )

        zero_below_pivot_gj1Ab_mR(Ab,--r);
        
return(Ab);
}
/* ------------------------------------ */
double **invgj_mR(
double **A,
double **invA
)
{
double **ID = i_duplicate_mR(A[R_SIZE][C0],A[C_SIZE][C0]);
double **AID;

int r = A[R_SIZE][C0];

  r--;

  AID = i_AbR0_mR(r,r,r);

       eye_mR(ID);
  c_A_b_Ab_mR(A,ID,AID);
   invgj_T_mR(AID);
    sort_c_mR(AID);
    sort_r_mR(AID);
   c_Inv_A_mR(AID,invA);

  f_mR(AID);
  f_mR(ID);

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