Linear Algebra and the C Language/a0m2

__NOTOC__ Install and compile this file in your working directory. /* ------------------------------------ */ /* Save as : c02c.c */ /* ------------------------------------ */ #include "v_a.h" /* ------------------------------------ */ #define RCA RC4 /* ------------------------------------ */ #define EV R1 /* ------------------------------------ */ int main(void) { double a[RCA*RCA]={ +5.701637292731, +0.991969965586, +0.254041088748, +1.125039107310, +0.991969965586, +6.402440296173, +0.359161539264, +1.590572531025, +0.254041088748, +0.359161539264, +5.091980394202, +0.407341745750, +1.125039107310, +1.590572531025, +0.407341745750, +6.803942016894 }; double **A = ca_A_mR(a, i_mR(RCA,RCA)); double **EValue = eigs_mR(A, i_mR(RCA,RCA)); double **Ide = eye_mR( i_mR(RCA,RCA)); double **sIde = smul_mR(EValue[EV][C1],Ide, i_mR(RCA,RCA)); double **AmnssIde = sub_mR(A,sIde, i_mR(RCA,RCA)); double **b = m0_mR(i_mR(RCA,C1)); double **Ab = c_A_b_Ab_mR(AmnssIde,b, i_Abr_Ac_bc_mR(RCA,RCA,C1)); double **Ab_free = i_Abr_Ac_bc_mR(RCA,RCA,RCA); double **b_free = i_mR(RCA,RCA); double **V = i_mR(RCA,RCA); clrscrn(); printf(" An eigenvector associated with an eigenvalue\n" " with single multiplicity is unique, up to a\n" " change in the sign. \n\n\n" " A:"); p_mR(A,S10,P5,C10); gj_PP_mR(Ab,NO); put_zeroR_mR(Ab,Ab_free); put_freeV_mR(Ab_free); gj_PP_mR(Ab_free,YES); c_Ab_b_mR(Ab_free,b_free); c_c_mR(b_free,C2,V,EV); Normalize_mR(V); printf("\n\n\n" " Copy V%d into the last file\n\n" " V%d:",EV,EV); P_mR(V,S10,P12,C10); stop(); f_mR(A); f_mR(EValue); f_mR(Ide); f_mR(sIde); f_mR(AmnssIde); f_mR(Ab); f_mR(b); f_mR(Ab_free); f_mR(b_free); f_mR(V); return 0; } /* ------------------------------------ */ /* ------------------------------------ */ Screen output example: An eigenvector associated with an eigenvalue with single multiplicity is unique, up to a change in the sign. A: +5.70164 +0.99197 +0.25404 +1.12504 +0.99197 +6.40244 +0.35916 +1.59057 +0.25404 +0.35916 +5.09198 +0.40734 +1.12504 +1.59057 +0.40734 +6.80394 Copy V1 into the last file V1: +0.418818962301, -nan, -nan, -nan, +0.592123360495, -nan, -nan, -nan, +0.151641348419, -nan, -nan, -nan, +0.671554543000, -nan, -nan, -nan Press return to continue.