Aller au contenu

Mathc complexes/029

Un livre de Wikilivres.


Application


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


c00a.c
/* ------------------------------------ */
/*  Save as :   c00a.c                  */
/* ------------------------------------ */
#include "w_a.h"
/* ------------------------------------ */
#define FACTOR_E    +1.E-5         
#define RCA          RC4 
#define EVALUE       RC1  
/* ------------------------------------ */       
/* ------------------------------------ */
void fun(void)
{                          
double a[RCA*(RCA*C2)] ={   
  +25072,   +0,   +21293,-4811,    +5386,-6531,   +11975,-5847, 
  +21293,+4811,   +30848,   +0,    -1526,-1774,   +13478,-6658, 
   +5386,+6531,    -1526,+1774,   +20007,   +0,    +7542,-1322, 
  +11975,+5847,   +13478,+6658,    +7542,+1322,   +21250,   +0  };
                       
double **A        = ca_A_mZ(a,           i_mZ(RCA,RCA));
double **sA       = smul_mZ(FACTOR_E, A, i_mZ(RCA,RCA));
double **EValue   = eigs_mZ(         sA, i_mZ(RCA,RCA));
double **V        =                      i_mZ(RCA,RCA);

double **Ide      =  eye_mZ(             i_mZ(RCA,RCA));
double **zIde     =                      i_mZ(RCA,RCA);
double **sAmnszIde =                     i_mZ(RCA,RCA);

double **Ab       =            i_Abr_Ac_bc_mZ(RCA,RCA,C1);
double **b        =                      i_mZ(    RCA,C1);
double **Ab_free  =            i_Abr_Ac_bc_mZ(RCA,RCA,C2);
double  **b_free  =                      i_mZ(    RCA,C2);

  clrscrn(); 
  printf(" A :");
  p_mZ(A, S8,P0, S6,P0, C4);  
  
         zmul_mZ(i_Z(EValue[EVALUE][C1],EValue[EVALUE][C2]),Ide,zIde);  
          sub_mZ(sA,zIde,sAmnszIde);
      
     c_A_b_Ab_mZ(sAmnszIde,b,Ab); 

  GJ_PP_FreeV_mZ(Ab,Ab_free,b_free);
    
          c_c_mZ(b_free,C2,V,EVALUE); 

    Normalize_mZ(V);
  
  printf(" V%d :",EVALUE);
  p_mZ(V, S9,P5, S8,P5, C4); 
  printf("\n\n\n\n\n\n\n\n");  
  stop();
              
  f_mZ(A);
  f_mZ(V);  
  f_mZ(EValue);
  
  f_mZ(Ide);
  f_mZ(zIde);
  f_mZ(sAmnszIde);
  
  f_mZ(sA);
  
  f_mZ(Ab);
  f_mZ(b);
  
  f_mZ(Ab_free);
  f_mZ(b_free);    
}
/* ------------------------------------ */
int main(void)
{
    fun();
    
  return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */


Calculer un vecteur propre :  Permet de calculer le vecteur propre pour la valeur propre "EVALUE".


Exemple de sortie écran :
 A :
  +25072    +0i   +21293 -4811i    +5386 -6531i   +11975 -5847i 
  +21293 +4811i   +30848    +0i    -1526 -1774i   +13478 -6658i 
   +5386 +6531i    -1526 +1774i   +20007    +0i    +7542 -1322i 
  +11975 +5847i   +13478 +6658i    +7542 +1322i   +21250    +0i 

 V1 :
 +0.48972-0.32491i      -nan    -nani      -nan    -nani      -nan    -nani 
 +0.58126-0.25956i      -nan    -nani      -nan    -nani      -nan    -nani 
 +0.18661+0.05423i      -nan    -nani      -nan    -nani      -nan    -nani 
 +0.46001-0.00000i      -nan    -nani      -nan    -nani      -nan    -nani 

 Les nan correspondent à la normalisation d'un vecteur nul.







 Press return to continue.