Linear Algebra and the C Language/a0ac

__NOTOC__ Install and compile this file in your working directory. /* ------------------------------------ */ /* Save as : c00i.c */ /* ------------------------------------ */ #include "v_a.h" /* ------------------------------------ */ void fun(int r) { double **U = r_mR(i_mR(r,r),9.); double **Orth = i_mR(r,r); double **Orth_Normal = i_mR(r,r); double **Orth_Normal_T = i_mR(r,r); clrscrn(); printf(" U :"); p_mR(U, S8,P4, C6); printf(" Orth_Normal : orth_mR(U,Orth,YES);"); orth_mR(U,Orth,YES); p_mR(Orth, S8,P4, C6); stop(); clrscrn(); printf(" Orth_Normal :"); p_mR(Orth, S8,P4, C6); printf(" Orth_Normal_T : transpose_mR(Orth,Orth_Normal_T);"); transpose_mR(Orth,Orth_Normal_T); p_mR(Orth_Normal_T, S8,P4, C6); printf(" Orth_Normal_T Orth_Normal :"); mul_mR(Orth_Normal_T,Orth, U); p_mR(U, S8,P4, C6); f_mR(U); f_mR(Orth); f_mR(Orth_Normal); f_mR(Orth_Normal_T); } /* ------------------------------------ */ int main(void) { time_t t; srand(time(&t)); do { fun(rp_I(R2)+R2); } while(stop_w()); return 0; } /* ------------------------------------ */ /* ------------------------------------ */ Screen output example: U : +8.0000 -3.0000 -1.0000 +1.0000 +8.0000 +7.0000 +1.0000 +5.0000 +1.0000 -3.0000 -6.0000 +9.0000 +6.0000 -6.0000 +9.0000 -2.0000 Orth_Normal : orth_mR(U,Orth,YES); +0.6228 -0.2625 -0.3959 -0.6217 +0.6228 +0.7242 +0.0386 +0.2935 +0.0778 -0.2918 -0.7000 +0.6471 +0.4671 -0.5669 +0.5930 +0.3297 Press return to continue. Orth_Normal : +0.6228 -0.2625 -0.3959 -0.6217 +0.6228 +0.7242 +0.0386 +0.2935 +0.0778 -0.2918 -0.7000 +0.6471 +0.4671 -0.5669 +0.5930 +0.3297 Orth_Normal_T : transpose_mR(Orth,Orth_Normal_T); +0.6228 +0.6228 +0.0778 +0.4671 -0.2625 +0.7242 -0.2918 -0.5669 -0.3959 +0.0386 -0.7000 +0.5930 -0.6217 +0.2935 +0.6471 +0.3297 Orth_Normal_T Orth_Normal : +1.0000 +0.0000 +0.0000 +0.0000 +0.0000 +1.0000 -0.0000 -0.0000 +0.0000 -0.0000 +1.0000 +0.0000 +0.0000 -0.0000 +0.0000 +1.0000 Press return to continue Press X return to stop