Linear Algebra and the C Language/a0l0

__NOTOC__ Install and compile this file in your working directory. /* ------------------------------------ */ /* Save as: c00c.c */ /* ------------------------------------ */ #include "v_a.h" /* ------------------------------------ */ /* ------------------------------------ */ void X_gj3_T_mR( double **Ab, int above ) { int r= R1; int c= C1; while( (r<Ab[R_SIZE][C0]) && (c<Ab[C_SIZE_A][C0]) ) { printf(" Ab:"); p_mR(Ab,S12,P4,C8); if(pivotbest_gj3Ab_mR(Ab,r,c)>ERROR_E) { printf(" Ab: pivot[%d][[%d]" " pivotbest_gj3Ab_mR_mR(Ab,r);",r,r); p_mR(Ab,S12,P4,C8); zero_under_pivot_gj3Ab_mR(Ab,r,c); printf(" Ab: pivot[%d][[%d]" " zero_under_pivot_gj1Ab_mR();",r,r); p_mR(Ab,S12,P4,C8); stop(); clrscrn(); } else r--; r++; c++; } r = Ab[R_SIZE][C0]; /* if(above) while(r>R1) { r--; c=C1; while(fabs(Ab[r][c])<ERROR_E && c<(Ab[C_SIZE_A][C0]-C1) ) c++; if(fabs(Ab[r][c])>ERROR_E) zero_above_pivot_gj3Ab_mR(Ab,r,c); }*/ } /* ------------------------------------ */ double **X_gj_PP_mR( double **Ab, int above ) { X_gj3_T_mR(Ab,above); 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_PP_mR(Ab,NO); printf(" gj_PP_mR(Ab,NO);"); 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: +327.0000 -155.0000 -511.0000 -507.0000 +473.0000 +633.0000 -514.0000 -180.0000 -176.0000 -453.0000 -618.0000 -466.0000 -903.0000 +335.0000 +302.0000 +746.0000 -508.0000 -251.0000 +116.0000 -566.0000 +925.0000 -528.0000 +940.0000 -986.0000 +325.0000 +108.0000 -762.0000 -383.0000 -867.0000 +982.0000 Ab: pivot[1][[1] pivotbest_gj3Ab_mR_mR(Ab,r); -903.0000 +335.0000 +302.0000 +746.0000 -508.0000 -251.0000 -514.0000 -180.0000 -176.0000 -453.0000 -618.0000 -466.0000 +327.0000 -155.0000 -511.0000 -507.0000 +473.0000 +633.0000 +116.0000 -566.0000 +925.0000 -528.0000 +940.0000 -986.0000 +325.0000 +108.0000 -762.0000 -383.0000 -867.0000 +982.0000 Ab: pivot[1][[1] zero_under_pivot_gj1Ab_mR(); +1.0000 -0.3710 -0.3344 -0.8261 +0.5626 +0.2780 +0.0000 -370.6866 -347.9025 -877.6334 -328.8394 -323.1274 +0.0000 -33.6877 -401.6379 -236.8538 +289.0399 +542.1063 +0.0000 -522.9657 +963.7951 -432.1683 +874.7420 -1018.2436 +0.0000 +228.5703 -653.3068 -114.5061 -1049.8350 +891.6622 Press return to continue. Ab: +1.0000 -0.3710 -0.3344 -0.8261 +0.5626 +0.2780 +0.0000 -370.6866 -347.9025 -877.6334 -328.8394 -323.1274 +0.0000 -33.6877 -401.6379 -236.8538 +289.0399 +542.1063 +0.0000 -522.9657 +963.7951 -432.1683 +874.7420 -1018.2436 +0.0000 +228.5703 -653.3068 -114.5061 -1049.8350 +891.6622 Ab: pivot[2][[2] pivotbest_gj3Ab_mR_mR(Ab,r); +1.0000 -0.3710 -0.3344 -0.8261 +0.5626 +0.2780 +0.0000 -522.9657 +963.7951 -432.1683 +874.7420 -1018.2436 +0.0000 -33.6877 -401.6379 -236.8538 +289.0399 +542.1063 +0.0000 -370.6866 -347.9025 -877.6334 -328.8394 -323.1274 +0.0000 +228.5703 -653.3068 -114.5061 -1049.8350 +891.6622 Ab: pivot[2][[2] zero_under_pivot_gj1Ab_mR(); +1.0000 -0.3710 -0.3344 -0.8261 +0.5626 +0.2780 -0.0000 +1.0000 -1.8429 +0.8264 -1.6727 +1.9471 +0.0000 +0.0000 -463.7223 -209.0150 +232.6919 +607.6982 +0.0000 +0.0000 -1031.0563 -571.3055 -948.8708 +398.6203 +0.0000 +0.0000 -232.0650 -303.3920 -667.5153 +446.6229 Press return to continue. Ab: +1.0000 -0.3710 -0.3344 -0.8261 +0.5626 +0.2780 -0.0000 +1.0000 -1.8429 +0.8264 -1.6727 +1.9471 +0.0000 +0.0000 -463.7223 -209.0150 +232.6919 +607.6982 +0.0000 +0.0000 -1031.0563 -571.3055 -948.8708 +398.6203 +0.0000 +0.0000 -232.0650 -303.3920 -667.5153 +446.6229 Ab: pivot[3][[3] pivotbest_gj3Ab_mR_mR(Ab,r); +1.0000 -0.3710 -0.3344 -0.8261 +0.5626 +0.2780 -0.0000 +1.0000 -1.8429 +0.8264 -1.6727 +1.9471 +0.0000 +0.0000 -1031.0563 -571.3055 -948.8708 +398.6203 +0.0000 +0.0000 -463.7223 -209.0150 +232.6919 +607.6982 +0.0000 +0.0000 -232.0650 -303.3920 -667.5153 +446.6229 Ab: pivot[3][[3] zero_under_pivot_gj1Ab_mR(); +1.0000 -0.3710 -0.3344 -0.8261 +0.5626 +0.2780 -0.0000 +1.0000 -1.8429 +0.8264 -1.6727 +1.9471 -0.0000 -0.0000 +1.0000 +0.5541 +0.9203 -0.3866 +0.0000 +0.0000 +0.0000 +47.9323 +659.4510 +428.4168 +0.0000 +0.0000 +0.0000 -174.8054 -453.9482 +356.9034 Press return to continue. Ab: +1.0000 -0.3710 -0.3344 -0.8261 +0.5626 +0.2780 -0.0000 +1.0000 -1.8429 +0.8264 -1.6727 +1.9471 -0.0000 -0.0000 +1.0000 +0.5541 +0.9203 -0.3866 +0.0000 +0.0000 +0.0000 +47.9323 +659.4510 +428.4168 +0.0000 +0.0000 +0.0000 -174.8054 -453.9482 +356.9034 Ab: pivot[4][[4] pivotbest_gj3Ab_mR_mR(Ab,r); +1.0000 -0.3710 -0.3344 -0.8261 +0.5626 +0.2780 -0.0000 +1.0000 -1.8429 +0.8264 -1.6727 +1.9471 -0.0000 -0.0000 +1.0000 +0.5541 +0.9203 -0.3866 +0.0000 +0.0000 +0.0000 -174.8054 -453.9482 +356.9034 +0.0000 +0.0000 +0.0000 +47.9323 +659.4510 +428.4168 Ab: pivot[4][[4] zero_under_pivot_gj1Ab_mR(); +1.0000 -0.3710 -0.3344 -0.8261 +0.5626 +0.2780 -0.0000 +1.0000 -1.8429 +0.8264 -1.6727 +1.9471 -0.0000 -0.0000 +1.0000 +0.5541 +0.9203 -0.3866 -0.0000 -0.0000 -0.0000 +1.0000 +2.5969 -2.0417 +0.0000 +0.0000 +0.0000 +0.0000 +534.9766 +526.2811 Press return to continue. Ab: +1.0000 -0.3710 -0.3344 -0.8261 +0.5626 +0.2780 -0.0000 +1.0000 -1.8429 +0.8264 -1.6727 +1.9471 -0.0000 -0.0000 +1.0000 +0.5541 +0.9203 -0.3866 -0.0000 -0.0000 -0.0000 +1.0000 +2.5969 -2.0417 +0.0000 +0.0000 +0.0000 +0.0000 +534.9766 +526.2811 Ab: pivot[5][[5] pivotbest_gj3Ab_mR_mR(Ab,r); +1.0000 -0.3710 -0.3344 -0.8261 +0.5626 +0.2780 -0.0000 +1.0000 -1.8429 +0.8264 -1.6727 +1.9471 -0.0000 -0.0000 +1.0000 +0.5541 +0.9203 -0.3866 -0.0000 -0.0000 -0.0000 +1.0000 +2.5969 -2.0417 +0.0000 +0.0000 +0.0000 +0.0000 +534.9766 +526.2811 Ab: pivot[5][[5] zero_under_pivot_gj1Ab_mR(); +1.0000 -0.3710 -0.3344 -0.8261 +0.5626 +0.2780 -0.0000 +1.0000 -1.8429 +0.8264 -1.6727 +1.9471 -0.0000 -0.0000 +1.0000 +0.5541 +0.9203 -0.3866 -0.0000 -0.0000 -0.0000 +1.0000 +2.5969 -2.0417 +0.0000 +0.0000 +0.0000 +0.0000 +1.0000 +0.9837 Press return to continue. gj_PP_mR(Ab,NO); +1.0000 -0.3710 -0.3344 -0.8261 +0.5626 +0.2780 -0.0000 +1.0000 -1.8429 +0.8264 -1.6727 +1.9471 -0.0000 -0.0000 +1.0000 +0.5541 +0.9203 -0.3866 -0.0000 -0.0000 -0.0000 +1.0000 +2.5969 -2.0417 +0.0000 +0.0000 +0.0000 +0.0000 +1.0000 +0.9837 Press return to continue Press X return to stop