Linear Algebra and the C Language/a09n

__NOTOC__ Install and compile this file in your working directory. /* ------------------------------------ */ /* Save as : c00b.c */ /* ------------------------------------ */ #include "v_a.h" /* ------------------------------------ */ /* ------------------------------------ */ #define RAb R3 #define CA C3 #define Cb C1 /* ------------------------------------ */ # define FREEV C1 /* ------------------------------------ */ /* ------------------------------------ */ int main(void) { double ab[RAb*(CA+C1)]={ +346, +262, +43, -578, -560, +831, -984, -731, +346, +262, +43, -578 }; double **Ab = ca_A_mR(ab,i_Abr_Ac_bc_mR(RAb, CA, Cb)); double **Ab_free = i_Abr_Ac_bc_mR(CA, CA, Cb + FREEV); double **b_free = i_mR(CA, Cb + FREEV); clrscrn(); printf(" Ab :"); p_mR(Ab, S10,P4,C10); printf("\n Ab : gj_PP_mR(Ab,NO) :"); gj_PP_mR(Ab,NO); p_mR(Ab, S10,P4,C10); stop(); clrscrn(); printf(" Ab :"); p_mR(Ab, S10,P4,C10); put_zeroR_mR(Ab,Ab_free); printf(" Ab_free : put_zeroR_mR(Ab,Ab_free);"); p_mR(Ab_free, S10,P4,C10); put_freeV_mR(Ab_free); printf(" Ab_free : put_freeV_mR(Ab_free);"); p_mR(Ab_free, S10,P4,C10); printf("\n Ab_free : gj_PP_mR(Ab_free,YES);"); gj_PP_mR(Ab_free,YES); p_mR(Ab_free, S10,P4,C10); stop(); clrscrn(); c_Ab_b_mR(Ab_free,b_free); printf(" b_free :"); p_mR(b_free, S10,P4,C7); printf(" b FREEV"); p_freeV(b_free, S8,P4); stop(); f_mR(Ab); f_mR(Ab_free); f_mR(b_free); return 0; } /* ------------------------------------ */ /* ------------------------------------ */ Screen output example: Ab : +346.0000 +262.0000 +43.0000 -578.0000 -560.0000 +831.0000 -984.0000 -731.0000 +346.0000 +262.0000 +43.0000 -578.0000 Ab : gj_PP_mR(Ab,NO) : +1.0000 -1.4839 +1.7571 +1.3054 +0.0000 +1.0000 -0.7286 -1.3278 +0.0000 +0.0000 +0.0000 +0.0000 Press return to continue. Ab : +1.0000 -1.4839 +1.7571 +1.3054 +0.0000 +1.0000 -0.7286 -1.3278 +0.0000 +0.0000 +0.0000 +0.0000 Ab_free : put_zeroR_mR(Ab,Ab_free); +1.0000 -1.4839 +1.7571 +1.3054 +0.0000 +0.0000 +1.0000 -0.7286 -1.3278 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 Ab_free : put_freeV_mR(Ab_free); +1.0000 -1.4839 +1.7571 +1.3054 +0.0000 +0.0000 +1.0000 -0.7286 -1.3278 +0.0000 +0.0000 +0.0000 +1.0000 +0.0000 +1.0000 Ab_free : gj_PP_mR(Ab_free,YES); +1.0000 +0.0000 +0.0000 -0.6651 -0.6760 +0.0000 +1.0000 +0.0000 -1.3278 +0.7286 +0.0000 +0.0000 +1.0000 +0.0000 +1.0000 Press return to continue. b_free : -0.6651 -0.6760 -1.3278 +0.7286 +0.0000 +1.0000 b FREEV x1 = -0.6651 -0.6760*v x2 = -1.3278 +0.7286*v x3 = +0.0000 +1.0000*v Press return to continue.