Aller au contenu

Mathc matrices/040

Un livre de Wikilivres.


Étude du code


Installer et compiler ces fichiers dans votre répertoire de travail.


c00a.c
/* ------------------------------------ */
/*  Save as :   c00a.c                  */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
#define   RCA    RC3
/* ------------------------------------ */
#define    EV    R2  
/* ------------------------------------ */
int main(void)
{
double a[RCA*RCA]={
+8.146673651126, +0.924044002095,  +1.085385018334, 
+0.924044002095, +13.821477213201, +6.837925615505, 
+1.085385018334, +6.837925615506,  +16.031849135673   
};

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     =            i_mR(RCA,RCA);
double **AmnssIde =            i_mR(RCA,RCA);

double **Ab       =            i_Abr_Ac_bc_mR(RCA,RCA,C1);
double **b        =            i_mR(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(" A :");
  p_mR(A,S10,P5,C10);
  
         smul_mR(EValue[EV][C1],Ide,sIde);          
        MmnsD_mR(A,sIde,AmnssIde);
        
     c_A_b_Ab_mR(AmnssIde,b,Ab); 
        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); 
          c_c_mR(b_free,C3,V,EV+C1); 
                    
    Normalize_mR(V);
           
          printf(" V%d : Copy V%d into the last file",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;
}
/* ------------------------------------ */
/* ------------------------------------ */


Calculer les deux premiers vecteurs propres


Exemple de sortie écran :
 A :
  +8.14667   +0.92404   +1.08539 
  +0.92404  +13.82148   +6.83793 
  +1.08539   +6.83793  +16.03185 

 V2 : Copy V2 into the last file
      -nan, -0.987635528464, -0.990992430410, 
      -nan, +0.156767544201, +0.000000000000, 
      -nan, +0.000000000000, +0.133917896001  

 Press return to continue.