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);
}
/* ------------------------------------ */
/* ------------------------------------ */