Linear Algebra and the C Language/a01u

/* ------------------------------------ */
/*  Save as :   vgj1.h                  */
/* ------------------------------------ */
double **gj1_T_mR(
double **Ab
)
{
int r= R1;

    while(r<Ab[R_SIZE][C0])
       {
        if(fabs(pivotbest_gj1Ab_mR(Ab,r))>ERROR_E)  
         
          zero_under_pivot_gj1Ab_mR(Ab,r);
          
        r++;
       } 
                
    while(r>R1)
    
        zero_above_pivot_gj1Ab_mR(Ab,--r);
        
return(Ab);
}
/* ------------------------------------ */
double **gj_TP_mR(
double **Ab
)
{
   gj1_T_mR(Ab);
  sort_c_mR(Ab);
  sort_r_mR(Ab);
  
return(Ab);
}
/* ------------------------------------ */
double **gj1_mR(
double **Ab
)
{
   gj_TP_mR(Ab);
  
return(Ab);
}
/* ------------------------------------ */
/* ------------------------------------ */