Linear Algebra and the C Language/a0kg
Install and compile this file in your working directory.
/* ------------------------------------ */
/* Save as : c00e.c */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
/* ------------------------------------ */
double **X_sort_r_mR(
double **Ab
)
{
double pivot;
int c;
int r;
int r_pivot = R1;
clrscrn();
printf(" Ab : Before the function sort_r_mR(); \n");
p_mR(Ab, S8,P4,C6);
stop();
clrscrn();
for ( c=C1; c<Ab[C_SIZE_A][C0]; c++)
{
r = r_pivot;
pivot = 0.;
while( (pivot<ERROR_E) && r<Ab[R_SIZE][C0] )
{
pivot = fabs( Ab[r][c] );
clrscrn();
printf("/* --------- DEBUG 1 ---------- */\n\n");
printf(" Ab : pivot[%d][%d] = %f \n",r,c,pivot);
p_mR(Ab, S8,P4,C8);
if(pivot>ERROR_E)
{
swapR_mR(Ab, r, r_pivot);
r_pivot++;
}
printf(" Ab");
p_mR(Ab, S8,P4,C8);
stop();
clrscrn();
/* --------- DEBUG 1 ---------- */
r++;
}
}
return(Ab);
}
/* ------------------------------------ */
/* ------------------------------------ */
double **X_gj_TP_mR(
double **Ab
)
{
gj1_T_mR(Ab);
sort_c_mR(Ab);
X_sort_r_mR(Ab);
return(Ab);
}
/* ------------------------------------ */
/* ------------------------------------ */
void fun(int r)
{
double **A = r_mR( i_mR(r,r ),999.);
double **b = r_mR( i_mR(r,C1),999.);
double **Ab = c_A_b_Ab_mR(A,b, i_Abr_Ac_bc_mR(r, r, C1));
/* i_Abr_Ac_bc_mR(RAb, CA, Cb)); */
clrscrn();
X_gj_TP_mR(Ab);
printf("Ab : After the function: sort_r_mR();\n\n"
" gj_TP_mR(Ab);");
p_mR(Ab, S10,P4,C7);
f_mR(Ab);
f_mR(b);
f_mR(A);
}
/* ------------------------------------ */
int main(void)
{
time_t t;
srand(time(&t));
do
{
fun(R5);
} while(stop_w());
return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */
Screen output example:
Ab : Before the function sort_r_mR();
+0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +4.5128
+0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +4.8739
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 -1.5244
+0.0000 -0.0000 -0.0000 +1.0000 -0.0000 -0.5550
+1.0000 +0.0000 +0.0000 +0.0000 +0.0000 -10.7817
Press return to continue.
/* --------- DEBUG 1 ---------- */
Ab : pivot[1][1] = 0.000000
+0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +4.5128
+0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +4.8739
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 -1.5244
+0.0000 -0.0000 -0.0000 +1.0000 -0.0000 -0.5550
+1.0000 +0.0000 +0.0000 +0.0000 +0.0000 -10.7817
Ab
+0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +4.5128
+0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +4.8739
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 -1.5244
+0.0000 -0.0000 -0.0000 +1.0000 -0.0000 -0.5550
+1.0000 +0.0000 +0.0000 +0.0000 +0.0000 -10.7817
Press return to continue.
/* --------- DEBUG 1 ---------- */
Ab : pivot[2][1] = 0.000000
+0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +4.5128
+0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +4.8739
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 -1.5244
+0.0000 -0.0000 -0.0000 +1.0000 -0.0000 -0.5550
+1.0000 +0.0000 +0.0000 +0.0000 +0.0000 -10.7817
Ab
+0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +4.5128
+0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +4.8739
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 -1.5244
+0.0000 -0.0000 -0.0000 +1.0000 -0.0000 -0.5550
+1.0000 +0.0000 +0.0000 +0.0000 +0.0000 -10.7817
Press return to continue.
/* --------- DEBUG 1 ---------- */
Ab : pivot[3][1] = 0.000000
+0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +4.5128
+0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +4.8739
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 -1.5244
+0.0000 -0.0000 -0.0000 +1.0000 -0.0000 -0.5550
+1.0000 +0.0000 +0.0000 +0.0000 +0.0000 -10.7817
Ab
+0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +4.5128
+0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +4.8739
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 -1.5244
+0.0000 -0.0000 -0.0000 +1.0000 -0.0000 -0.5550
+1.0000 +0.0000 +0.0000 +0.0000 +0.0000 -10.7817
Press return to continue.
/* --------- DEBUG 1 ---------- */
Ab : pivot[4][1] = 0.000000
+0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +4.5128
+0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +4.8739
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 -1.5244
+0.0000 -0.0000 -0.0000 +1.0000 -0.0000 -0.5550
+1.0000 +0.0000 +0.0000 +0.0000 +0.0000 -10.7817
Ab
+0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +4.5128
+0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +4.8739
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 -1.5244
+0.0000 -0.0000 -0.0000 +1.0000 -0.0000 -0.5550
+1.0000 +0.0000 +0.0000 +0.0000 +0.0000 -10.7817
Press return to continue.
/* --------- DEBUG 1 ---------- */
Ab : pivot[5][1] = 1.000000
+0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +4.5128
+0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +4.8739
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 -1.5244
+0.0000 -0.0000 -0.0000 +1.0000 -0.0000 -0.5550
+1.0000 +0.0000 +0.0000 +0.0000 +0.0000 -10.7817
Ab
+1.0000 +0.0000 +0.0000 +0.0000 +0.0000 -10.7817
+0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +4.8739
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 -1.5244
+0.0000 -0.0000 -0.0000 +1.0000 -0.0000 -0.5550
+0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +4.5128
Press return to continue.
.... .... ....
Ab : After the function: sort_r_mR();
gj_TP_mR(Ab);
+1.0000 +0.0000 +0.0000 +0.0000 +0.0000 -10.7817
+0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +4.5128
+0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +4.8739
+0.0000 -0.0000 -0.0000 +1.0000 -0.0000 -0.5550
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 -1.5244
Press return to continue
Press X return to stop