Linear Algebra and the C Language/a09d


Code study: GJ_TPfreeV_mR(Ab,new_Ab);

  • The algorithm consists to copy Ab into new_Ab.
  • Put the free variables.
  • Solve the system.
/* ------------------------------------ */
/* ------------------------------------ */
double **GJ_TPfreeV_mR(
double **Ab,
double **new_Ab)
{
int r;
	            
  c_mR(Ab,new_Ab);
                                                         
  put_freeV_TP_mR(new_Ab,Ab[R_SIZE][C0]);
  
  r = rsize_R(new_Ab);
  
  while(r>R1)
    
        zero_above_pivot_gj1Ab_mR(new_Ab,r--);	
	
return(new_Ab);
}
/* ------------------------------------ */
/* ------------------------------------ */