Mathc complexes/038
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-4
#define RCA RC6
/* ------------------------------------ */
/* ------------------------------------ */
void fun(void)
{
double a[RCA*(RCA*C2)] ={
+23.037247782,-0.000000000, -9.054210143,-4.760094043, +7.055780246,-8.116625940, +5.837328527,+3.181102856, -5.497293954,+1.987290486, -3.853659006,-3.688738932,
-9.054210143,+4.760094043, +21.661298903,+0.000000000, -0.063767253,+2.540467712, -0.356099425,+12.951777333, -0.578233130,+5.883310886, +7.278478707,-7.059870759,
+7.055780246,+8.116625940, -0.063767253,-2.540467712, +41.444943135,-0.000000000, -11.011924524,-22.763984410, +4.447732552,-27.866895323, +1.844442206,-1.924521124,
+5.837328527,-3.181102856, -0.356099425,-12.951777333, -11.011924524,+22.763984410, +42.110215120,+0.000000000, +20.713304065,+9.756936697, -11.927657780,-5.975712843,
-5.497293954,-1.987290486, -0.578233130,-5.883310886, +4.447732552,+27.866895323, +20.713304065,-9.756936697, +42.740823156,+0.000000000, +0.112009750,+7.491403896,
-3.853659006,+3.688738932, +7.278478707,+7.059870759, +1.844442206,+1.924521124, -11.927657780,+5.975712843, +0.112009750,-7.491403896, +34.005471906,-0.000000000
};
double **A = ca_A_mZ(a, i_mZ(RCA,RCA));
double **V = i_mZ(RCA,RCA);
double **EigsValue = i_mZ(RCA,C1);
clrscrn();
printf(" Copy/Past into the octave windows \n\n");
p_Octave_mZ(A,"a",P9,P9);
printf(" format short e\n"
" [V, E] = eigs (a,%d)\n\n", RCA);
stop();
clrscrn();
eigs_V_mZ(A,V,FACTOR_E);
printf(" V :");
p_mZ(V, S10,P4, S8,P4, C4);
eigs_mZ(A,EigsValue);
printf(" EigsValue :");
p_mZ(EigsValue, S10,P4, S8,P4, C4);
stop();
f_mZ(A);
f_mZ(V);
f_mZ(EigsValue);
}
/* ------------------------------------ */
int main(void)
{
fun();
return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */
Exemple de référence : Il permet de vérifier les calculs.
Exemple de sortie écran :
Copy/Past into the octave windows
a=[
+23.037247782-0.000000000*i,-9.054210143-4.760094043*i,+7.055780246-8.116625940*i,+5.837328527+3.181102856*i,-5.497293954+1.987290486*i,-3.853659006-3.688738932*i;
-9.054210143+4.760094043*i,+21.661298903+0.000000000*i,-0.063767253+2.540467712*i,-0.356099425+12.951777333*i,-0.578233130+5.883310886*i,+7.278478707-7.059870759*i;
+7.055780246+8.116625940*i,-0.063767253-2.540467712*i,+41.444943135-0.000000000*i,-11.011924524-22.763984410*i,+4.447732552-27.866895323*i,+1.844442206-1.924521124*i;
+5.837328527-3.181102856*i,-0.356099425-12.951777333*i,-11.011924524+22.763984410*i,+42.110215120+0.000000000*i,+20.713304065+9.756936697*i,-11.927657780-5.975712843*i;
-5.497293954-1.987290486*i,-0.578233130-5.883310886*i,+4.447732552+27.866895323*i,+20.713304065-9.756936697*i,+42.740823156+0.000000000*i,+0.112009750+7.491403896*i;
-3.853659006+3.688738932*i,+7.278478707+7.059870759*i,+1.844442206+1.924521124*i,-11.927657780+5.975712843*i,+0.112009750-7.491403896*i,+34.005471906-0.000000000*i]
format short e
[V, E] = eigs (a,6)
Press return to continue.
V :
+0.0472 -0.0305i -0.3753 -0.2176i +0.5549 +0.2953i -0.0601 -0.3678i
-0.0358 -0.1350i +0.3886 -0.1975i -0.2947 -0.1076i -0.4218 +0.0222i
+0.2788 +0.4857i -0.0064 -0.3123i -0.0172 +0.1230i -0.4521 -0.3000i
-0.5611 +0.0519i -0.2410 -0.2193i +0.2069 -0.1780i +0.1740 +0.1319i
-0.4598 +0.3367i +0.1809 +0.0620i +0.2151 +0.1181i -0.4732 +0.2943i
+0.1461 +0.0000i +0.6176 +0.0000i +0.5970 +0.0000i +0.1737 -0.0000i
-0.0601 +0.1285i -0.4652 -0.1951i
+0.1630 +0.5438i -0.3752 -0.2251i
-0.0201 +0.1228i +0.4052 +0.3152i
+0.0323 +0.5181i +0.4202 -0.0900i
-0.1936 -0.3835i -0.0450 -0.2851i
+0.4338 +0.0000i +0.1501 -0.0000i
EigsValue :
+96.0000 +0.0000i
+53.0000 -0.0000i
+25.0000 +0.0000i
+19.0000 +0.0000i
+12.0000 -0.0000i
+0.0000 +0.0000i
Press return to continue.