Mathc complexes/08p
Installer et compiler ce fichier dans votre répertoire de travail.
|
c00c.c |
|---|
/* ------------------------------------ */
/* Save as : c00c.c */
/* ------------------------------------ */
#include "w_a.h"
/* ------------------------------------ */
#define ARRAY A7
#define RCA R4
/* ------------------------------------ */
void fun(void)
{
double **A[ARRAY];
double **P[ARRAY];
double **InvP[ARRAY];
double **T = i_mZ(RCA,RCA);
double **EigsValue = i_mZ(RCA,A1);
int c;
for(c=A0; c<ARRAY; c++)
{
A[c] = i_mZ(RCA,RCA);
P[c] = r_mZ( i_mZ(RCA,RCA),9);
InvP[c] = inv_mZ(P[c], i_mZ(RCA,RCA));
}
rcsymmetric_mZ(A[A0],9);
for(c=A0; c<(ARRAY-A1); c++)
{
mul_mZ(InvP[c],A[c],T);
mul_mZ(T,P[c],A[c+A1]);
}
clrscrn();
for(c=A0; c<ARRAY; c++)
{
printf(" A[%d] : ",c);
p_mZ(A[c], S8,P1, S8,P1,C6);
printf(" EigsValue of A[%d] : ",c);
eigs_mZ(A[c],EigsValue);
p_mZ(EigsValue,S10,P2,S10,P2,C6);
printf("\n\n");
}
for(c=A0; c<ARRAY; c++)
{
f_mZ(A[c]);
f_mZ(P[c]);
f_mZ(InvP[c]);
}
f_mZ(T);
f_mZ(EigsValue);
}
/* ------------------------------------ */
int main(void)
{
time_t t;
srand(time(&t));
do{
fun();
}while(stop_w());
return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */
Des matrices semblable ont les mêmes valeurs propres.
Exemple de sortie écran :
A[0] :
+262.0 +0.0i +105.0 -26.0i +166.0 +33.0i -40.0 +7.0i
+105.0 +26.0i +147.0 +0.0i -11.0 +40.0i +20.0 +26.0i
+166.0 -33.0i -11.0 -40.0i +207.0 +0.0i +50.0 +11.0i
-40.0 -7.0i +20.0 -26.0i +50.0 -11.0i +270.0 +0.0i
EigsValue of A[0] :
+439.74 +0.00i
+286.09 -0.00i
+160.03 +0.00i
+0.14 +0.00i
A[1] :
+16.6 -177.8i +2.4 -228.9i +8.7 -68.1i -187.0 -48.5i
+77.3 +214.7i +244.7 +343.8i +131.1 +90.7i +126.5 +97.9i
+48.3 -296.5i +428.4 -367.2i +289.6 -248.5i -146.7 -151.7i
+60.1 +59.9i -16.9 -195.7i -134.9 -105.9i +335.0 +82.6i
EigsValue of A[1] :
+439.74 +0.00i
+286.09 +0.00i
+160.03 +0.00i
+0.14 +0.00i
A[2] :
+628.9 -22.8i +49.5 -239.3i +142.9 -9.6i -4.1 -31.3i
+27.9 -122.9i +244.4 -73.4i +22.6 -56.2i -11.8 -80.0i
-725.7 +75.6i -355.4 +675.9i -129.7 +52.7i +34.1 +324.8i
-144.5 -74.0i -97.8 +71.8i -56.6 -22.6i +142.4 +43.5i
EigsValue of A[2] :
+439.74 -0.00i
+286.09 +0.00i
+160.03 +0.00i
+0.14 +0.00i
A[3] :
+374.4 -329.7i -306.4 +776.5i -272.9 +190.2i -975.5 -225.2i
+372.5 +236.1i -497.5 -521.3i -72.0 -390.0i +477.1 -877.5i
-452.7 +245.6i +1029.4 -242.9i +690.2 +213.6i +454.8 +1183.3i
-373.0 +33.8i +493.5 -56.4i +154.6 +130.7i +318.9 +637.5i
EigsValue of A[3] :
+439.74 +0.00i
+286.09 +0.00i
+160.03 +0.00i
+0.14 -0.00i
A[4] :
+1352.8 -372.2i +150.0 +376.1i +830.7 -65.3i -119.1 +203.3i
+1295.1 +732.0i +209.0 +388.4i +863.4 +692.4i -362.6 +155.4i
-854.7 +572.5i -149.5 -189.8i -451.2 +257.0i -4.4 -205.5i
+365.9 +2600.4i -772.2 +145.7i -156.5 +1787.3i -224.6 -273.2i
EigsValue of A[4] :
+439.74 +0.00i
+286.09 +0.00i
+160.03 +0.00i
+0.14 -0.00i
A[5] :
+1854.6 +398.5i -357.9 -2110.1i +241.1 -896.9i +20.0 +2081.2i
-214.9 -1455.0i -1474.0 +1064.1i -809.0 +104.6i +1793.7 -650.8i
-575.7 -243.4i +65.8 +990.1i +51.2 +296.6i +278.2 -882.3i
-1289.2 -309.5i +298.1 +1964.8i -165.7 +677.9i +454.3 -1759.2i
EigsValue of A[5] :
+439.74 +0.00i
+286.09 -0.00i
+160.03 +0.00i
+0.14 +0.00i
A[6] :
-94.1 -3888.2i +360.1 +56.4i -536.7 -1002.5i -938.3 -1865.5i
+1698.4 -3987.4i +1020.5 +622.1i -203.6 -961.3i +385.8 -2664.0i
+1257.2 +6295.7i -740.7 +149.4i +1162.9 +1462.6i +2025.5 +2982.8i
-3567.6 +2810.9i +157.8 -699.5i -430.7 +1344.5i -1203.3 +1803.5i
EigsValue of A[6] :
+439.74 -0.00i
+286.09 +0.00i
+160.03 -0.00i
+0.14 -0.00i