Mathc matrices/a10
Installer et compiler ces fichiers dans votre répertoire de travail.
|
e05.c |
|---|
/* ------------------------------------ */
/* Save as : f01.c */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
/* ------------------------------------ */
void XX_p_freeV(
double **b,
int e,
int d)
{
int m;
int r;
int c;
char v;
m = 0;
for(r=R1; r<b[R_SIZE][C0]; r++)
{
v = 'y' - (b[R_SIZE][C0]-C1);
printf("\n");
m++;
printf(" x%d = ",m);
for(c=C1; c<b[C_SIZE][C0]; c++)
{
if(c==C1) printf(" %+*.*f",e,d,b[r][c]);
else
{
if(fabs(b[r][c])>ERROR_E) printf(" %+*.*f*%c",e,d,b[r][c],v);
else printf(" %+*.*f*%c",e,d,0.,v);
v++;
}
}
}
printf("\n");
}
/* ------------------------------------ */
/* ------------------------------------ */
int main(void)
{
double ab[R4*C7]={
+1, -3, +4, -2, +5, +4, 0,
+2, -6, +9, -1, +8, +2, 0,
+2, -6, +9, -1, +9, +7, 0,
-1, 3, -4, 2, -5, -4, 0,
};
double **Ab = ca_A_mR(ab,i_Abr_Ac_bc_mR(R4,C6,C1));
double **A = c_Ab_A_mR(Ab,i_mR(R4,C6));
double **b = c_Ab_b_mR(Ab,i_mR(R4,C1));
double **Ab_free = i_Abr_Ac_bc_mR(csize_A_R(Ab),csize_A_R(Ab),C1+C3);
double **b_free = i_mR(rsize_R(Ab_free),csize_b_R(Ab_free));
int r;
clrscrn();
printf(" Ab :");
p_mR(Ab,S6,P1,C10);
printf(" Ab : gj_PP_mR(Ab,NO) :");
gj_PP_mR(Ab,NO);
p_mR(Ab,S6,P1,C10);
stop();
clrscrn();
put_zeroR_mR(Ab,Ab_free);
printf(" Ab_free : put_zeroR_mR(Ab,Ab_free);");
p_mR(Ab_free,S6,P1,C10);
stop();
clrscrn();
put_freeV_mR(Ab_free);
printf(" Ab_free : put_freeV_mR(Ab_free);");
p_mR(Ab_free,S6,P1,C10);
stop();
clrscrn();
r = rsize_R(Ab_free);
while(r>R1)
zero_below_pivot_gj1Ab_mR(Ab_free,r--);
printf(" Ab_free : zero_below_pivot_gj1Ab_mR(Ab_free,r--);");
p_mR(Ab_free,S6,P1,C10);
stop();
clrscrn();
c_Ab_b_mR(Ab_free,b_free);
printf(" b_free :");
p_mR(b_free,S10,P0,C7);
printf(" b_free : free variables");
XX_p_freeV(b_free,S8,P0);
stop();
f_mR(Ab);
f_mR(Ab_free);
f_mR(b_free);
f_mR(b);
f_mR(A);
return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */
Exemple de sortie écran :
------------------------------------
Ab :
+1.0 -3.0 +4.0 -2.0 +5.0 +4.0 +0.0
+2.0 -6.0 +9.0 -1.0 +8.0 +2.0 +0.0
+2.0 -6.0 +9.0 -1.0 +9.0 +7.0 +0.0
-1.0 +3.0 -4.0 +2.0 -5.0 -4.0 +0.0
Ab : gj_PP_mR(Ab,NO) :
+1.0 -3.0 +4.5 -0.5 +4.0 +1.0 +0.0
-0.0 -0.0 +1.0 +3.0 -2.0 -6.0 -0.0
+0.0 +0.0 +0.0 +0.0 +1.0 +5.0 +0.0
+0.0 +0.0 +0.0 +0.0 +0.0 +0.0 +0.0
Press return to continue.
------------------------------------
Ab_free : put_zeroR_mR(Ab,Ab_free);
+1.0 -3.0 +4.5 -0.5 +4.0 +1.0 +0.0 +0.0 +0.0 +0.0
+0.0 +0.0 +0.0 +0.0 +0.0 +0.0 +0.0 +0.0 +0.0 +0.0
-0.0 -0.0 +1.0 +3.0 -2.0 -6.0 -0.0 +0.0 +0.0 +0.0
+0.0 +0.0 +0.0 +0.0 +0.0 +0.0 +0.0 +0.0 +0.0 +0.0
+0.0 +0.0 +0.0 +0.0 +1.0 +5.0 +0.0 +0.0 +0.0 +0.0
+0.0 +0.0 +0.0 +0.0 +0.0 +0.0 +0.0 +0.0 +0.0 +0.0
Press return to continue.
------------------------------------
Ab_free : put_freeV_mR(Ab_free);
+1.0 -3.0 +4.5 -0.5 +4.0 +1.0 +0.0 +0.0 +0.0 +0.0
+0.0 +1.0 +0.0 +0.0 +0.0 +0.0 +0.0 +1.0 +0.0 +0.0
-0.0 -0.0 +1.0 +3.0 -2.0 -6.0 -0.0 +0.0 +0.0 +0.0
+0.0 +0.0 +0.0 +1.0 +0.0 +0.0 +0.0 +0.0 +1.0 +0.0
+0.0 +0.0 +0.0 +0.0 +1.0 +5.0 +0.0 +0.0 +0.0 +0.0
+0.0 +0.0 +0.0 +0.0 +0.0 +1.0 +0.0 +0.0 +0.0 +1.0
Press return to continue.
------------------------------------
Ab_free : zero_below_pivot_gj1Ab_mR(Ab_free,r--);
+1.0 +0.0 +0.0 +0.0 +0.0 +0.0 +0.0 +3.0 +14.0 +37.0
+0.0 +1.0 +0.0 +0.0 +0.0 +0.0 +0.0 +1.0 +0.0 +0.0
+0.0 +0.0 +1.0 +0.0 +0.0 +0.0 +0.0 +0.0 -3.0 -4.0
+0.0 +0.0 +0.0 +1.0 +0.0 +0.0 +0.0 +0.0 +1.0 +0.0
+0.0 +0.0 +0.0 +0.0 +1.0 +0.0 +0.0 +0.0 +0.0 -5.0
+0.0 +0.0 +0.0 +0.0 +0.0 +1.0 +0.0 +0.0 +0.0 +1.0
Press return to continue.
------------------------------------
b_free :
+0 +3 +14 +37
+0 +1 +0 +0
+0 +0 -3 -4
+0 +0 +1 +0
+0 +0 +0 -5
+0 +0 +0 +1
b_free : free variables
x1 = +0 +3*s +14*t +37*u
x2 = +0 +1*s +0*t +0*u
x3 = +0 +0*s -3*t -4*u
x4 = +0 +0*s +1*t +0*u
x5 = +0 +0*s +0*t -5*u
x6 = +0 +0*s +0*t +1*u
Press return to continue.