Linear Algebra and the C Language/a0j7

__NOTOC__ Install and compile this file in your working directory. /* ------------------------------------ */ /* Save as : c00a.c */ /* ------------------------------------ */ #include "v_a.h" /* ------------------------------------ */ /* ------------------------------------ */ double **X_c_withR0_mR( double **A, double **B ) { int r; int c; for (r=R0; r<A[R_SIZE][C0]; r++) for (c=C1; c<A[C_SIZE][C0]; c++) B[r][c] = A[r][c]; return(B); } /* ------------------------------------ */ /* ------------------------------------ */ void fun(int r,int c) { double **A = r_mR(i_mR(r,c),99.); double **B = i_mR(r,c); clrscrn(); printf(" A:\n"); pall_mR(swapC_withR0_mR(A,C1,C2), S5,P0,C10); printf(" B: c_withR0_mR(A,B);\n"); pall_mR(c_withR0_mR(A,B), S5,P0,C10); f_mR(A); f_mR(B); } /* ------------------------------------ */ int main(void) { time_t t; srand(time(&t)); int i; do { i = rp_I(R3)+R1; fun(i,i+C1); } while(stop_w()); return 0; } /* ------------------------------------ */ /* ------------------------------------ */ If the index of the columns of matrix A has been modified, and if we copy A into B, this function allows us to preserve the position of the indices of A in B. Screen output example: A: +5 +2 +1 +3 +4 +5 +6 -58 -71 +71 -28 -96 +0 -97 +59 -33 +87 -20 +0 -88 -14 +10 -35 -46 +0 +57 -26 +4 -40 +31 B: c_withR0_mR(A,B); +5 +2 +1 +3 +4 +5 +6 -58 -71 +71 -28 -96 +0 -97 +59 -33 +87 -20 +0 -88 -14 +10 -35 -46 +0 +57 -26 +4 -40 +31 Press return to continue Press X return to stop