adjoint


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


c00c.c
/* ------------------------------------ */
/*  Save as :   c00c.c                  */
/* ------------------------------------ */
#include "w_a.h"
/* ------------------------------------ */
void fun(int r)
{
double **A                =       r_mZ(  i_mZ(r,r),9.);
double **Aadjoint         = adjoint_mZ(A,i_mZ(r,r)   );

  clrscrn();
  printf(" A :");
  p_mZ(A,S4,P0, S3,P0, C6);
  printf(" Adjoint(A) :");
  p_mZ(Aadjoint, S6,P0, S5,P0, C6);

  printf(" det(adjoint(A)) = ");
  p_Z(det_Z(Aadjoint),S0,P0,S0,P0);printf("\n");
  printf(" det(A)**(r-R1)  = ");
  p_Z(pow_Z((r-R1),det_Z(A)),S0,P0,S0,P0);printf("\n"); 
  
  f_mZ(A);
  f_mZ(Aadjoint);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));
do
{
  fun(rp_I(RC3)+RC1);

} while(stop_w());

  return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */


det(adjoint(A)) = det(A)**(r-R1)


Exemple de sortie écran :
 A :
  -9 -5i   +8 -2i   -1 +7i   -5 +2i 
  -6 -2i   -2 -2i   +7 -9i   +8 -7i 
  -9 -9i   -7 +3i   +5 -3i   -3 -6i 
  -7 +1i   -2 +3i   -8 +4i   +6 -2i 

 Adjoint(A) :
   -72-1466i   +256 -440i   -168-1090i   +874 -734i 
 -2574 +834i  -1448+1544i  +2034-1566i  +1044 -596i 
  +495 -243i  -1594 -344i    +43+1127i  +1926 -520i 
    -8 -790i  -2042 +606i   +686-1080i  -1106+1434i 

 det(adjoint(A)) = +62673255583040+15409522815104i 
 det(A)**(r-R1)  = +62673255583040+15409522815104i 

 Press   return to continue
 Press X return to stop