Mathc complexes/027

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          RC3 
#define EVALUE       RC1  
/* ------------------------------------ */       
/* ------------------------------------ */
void fun(void)
{                          
double a[RCA*(RCA*C2)] ={   
 +26191,+1471,   -3444,+9191,   +18535, -3803,
  -3444,-9191,  +27183,+1470,    -3591, -9206,
 +18535,+3803,   -3591,+9206,   +29984, +5471};
                       
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, C3);  
  
         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("     zIde : EVALUE%d * Ide",EVALUE);                       p_mZ(zIde, S9,P5, S8,P5, C4);  

//  printf(" sAmnszIde : A - (EVALUE%d*Ide)",EVALUE);                p_mZ(sAmnszIde, S9,P5, S8,P5, C4); 
  
//  printf("       Ab : [A-(EVALUE%d*Ide)] V%d = 0 ",EVALUE,EVALUE); p_mZ(Ab, S9,P5, S8,P5, C4);  

//  printf("   b_free :");                                           p_mZ(b_free, S9,P5, S8,P5, C4); 
        
//  printf("      V%d :",EVALUE);                                    p_mZ(V, S9,P5, S8,P5, C4); 
      
  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;
}
/* ------------------------------------ */
/* ------------------------------------ */

Étudier les calculs intermédiaires :  En supprimant les "//" vous pouvez observer les calculs intermédiaires

Exemple de sortie écran :

        A :
  +26191 +1471i    -3444 +9191i   +18535 -3803i 
   -3444 -9191i   +27183 +1470i    -3591 -9206i 
  +18535 +3803i    -3591 +9206i   +29984 +5471i 

     zIde : EVALUE1 * Ide
 +0.54119+0.03188i  +0.00000+0.00000i  +0.00000+0.00000i 
 +0.00000+0.00000i  +0.54119+0.03188i  +0.00000+0.00000i 
 +0.00000+0.00000i  +0.00000+0.00000i  +0.54119+0.03188i 

 Press return to continue.