Mathc matrices/Fichiers h : vrm sy3

Installer ce fichier dans votre répertoire de travail.

	`vrm_com4.h`

/* ------------------------------------ */
/* ------------------------------------ */
/*  Save as :   vrm_com4.h              */
/* ------------------------------------ */

/* ------------------------------------ */
void r_commute4_mR(
double **A,
double **B,
double **C,
double **E,
int      n
)
{
double **P; 
double **Pinv; 
double **D;
double **T;
double **Dt;

int r=rsize_R(A);

 P    = i_mR(r,r);
 Pinv = i_mR(r,r);
 D    = i_mR(r,r);
 T    = i_mR(r,r);
 Dt   = i_mR(r,r);
 
 while(!det_R(r_mR(P,n)));
 rdiag_mR(D,n);
 
 inv_mR(P,Pinv);
 mul_mR(P,Diag_sin_mR(D,Dt),T);
 mul_mR(T,Pinv,A);

 mul_mR(P,Diag_cos_mR(D,Dt),T);
 mul_mR(T,Pinv,B);
 
 mul_mR(P,Diag_sinh_mR(D,Dt),T);
 mul_mR(T,Pinv,C);

 mul_mR(P,Diag_cosh_mR(D,Dt),T);
 mul_mR(T,Pinv,E);
 
 
 f_mR(P);
 f_mR(Pinv);
 f_mR(D);
 f_mR(T);
 f_mR(Dt);
}
/* ------------------------------------ */
/* ------------------------------------ */

Quatre matrices commutatives