Linear Algebra and the C Language/a01x

__NOTOC__ Install this file in your working directory. /* ------------------------------------ */ /* Save as : vgc_ab.h */ /* ------------------------------------ */ double **c_A_b_Ab_mR( double **A, double **b, double **Ab ) { int r; int c; int c_A; for (r=R1; r<A[R_SIZE][C0]; r++) for (c=C1; c<A[C_SIZE][C0]; c++) Ab[r][c] = A[r][c]; c_A = --c; for (r=R1; r<A[R_SIZE][C0]; r++) for (c=C1; c<b[C_SIZE][C0]; c++) Ab[r][c+c_A] = b[r][c]; return(Ab); } /* ------------------------------------ */ double **c_Ab_A_mR( double **Ab, double **A ) { int r; int c; for (r=R1; r<Ab[R_SIZE][C0]; r++) for (c=C1; c<Ab[C_SIZE_A][C0]; c++) A[r][c] = Ab[r][c]; return(A); } /* ------------------------------------ */ double **c_Ab_b_mR( double **Ab, double **b ) { int r; int c; int b_c; for (r=R1; r<Ab[R_SIZE][C0]; r++) for (c=Ab[C_SIZE_A][C0],b_c=C1; c<Ab[C_SIZE][C0]; c++,b_c++) b[r][b_c] = Ab[r][c]; return(b); } /* ------------------------------------ */ double **c_Ab_subArxr_mR( double **Ab, double **subA ) { int r; int c; for (r=R1; r<Ab[R_SIZE][C0]; r++) for (c=C1; c<Ab[R_SIZE][C0]; c++) subA[r][c] = Ab[r][c]; return(subA); } /* ------------------------------------ */ /* ------------------------------------ */