Mathc matrices/e1503
Installer et compiler ces fichiers dans votre répertoire de travail.
|
gj03.c |
|---|
/* ------------------------------------ */
/* Save as : gj03.c */
/* ------------------------------------ */
#include "v_a.h"
/* ------------------------------------ */
/* ------------------------------------ */
double **xx_sort_c_mR(
double **Ab
)
{
double **T = i_RC_mR(Ab[R_SIZE][C0],Ab[C_SIZE][C0]);
int c;
for(c=C1; c<Ab[C_SIZE][C0]; c++)
{
clrscrn();
printf("/* --------- DEBUG 1 ---------- */\n\n");
printf(" Copy C[%d] of Ab into C[%.0f] of T\n",c,Ab[R0][c]);
pall_mR(Ab,S8,P4,C8);
c_c_withR0_mR( Ab, c,
T, Ab[R0][c]);
printf(" T[R0][%.0f]\n",Ab[R0][c]);
pall_mR(T,S8,P4,C8);
stop();
clrscrn();
/* --------- DEBUG 1 ---------- */
}
printf(" Ab : before\n");
pall_mR(Ab,S8,P4,C8);
c_withR0_mR(T,Ab);
printf(" Ab : after \n");
pall_mR(Ab,S8,P4,C8);
stop();
clrscrn();
f_mR(T);
return(Ab);
}
/* ------------------------------------ */
/* ------------------------------------ */
double **XX_gj_TP_mR(
double **Ab
)
{
gj1_T_mR(Ab);
xx_sort_c_mR(Ab);
sort_r_mR(Ab);
return(Ab);
}
/* ------------------------------------ */
/* ------------------------------------ */
void fun(int r)
{
double **A = i_mR(r,r);
double **b = i_mR(r,C1);
double **Ab = i_Abr_Ac_bc_mR(r,r,C1);
/* i_Abr_Ac_bc_mR(R_Ab,C_A,C_b)); */
do
{
r_mR(A,999.);
printf(".");
}while(!det_R(A));
r_mR(b,999.);
c_A_b_Ab_mR(A,b,Ab);
clrscrn();
printf(" Copy/Past into the octave window.\n\n");
p_Octave_mR(Ab,"Ab",P0);
printf("\n rref(Ab,.00000000001)\n\n");
stop();
clrscrn();
XX_gj_TP_mR(Ab);
printf(" gj_TP_mR(Ab) :\n");
p_mR(Ab,S10,P4,C7);
f_mR(Ab);
f_mR(b);
f_mR(A);
}
/* -----------fun(rp_I(R2)+R2);------------------------- */
int main(void)
{
time_t t;
srand(time(&t));
do
{
fun(R5);
} while(stop_w());
return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */
Nous utilisons la fonction pall_mR(); au lieu de p_mR(); pour voir l'index des colonnes. Nous pouvons voir comment les colonnes sont triées. On copie la matrice Ab colonne par colonne dans la matrice T. On met simplement l'index de la matrice Ab dans le bon ordre.
Exemple de sortie écran :
------------------------------------
Copy/Past into the octave window.
Ab=[
+422,+390,-715,+340,+868,-505;
+30,+188,-867,-723,-319,+406;
-807,+220,+390,+254,+104,-325;
-29,+828,+514,+704,-421,-193;
-399,+140,+298,-47,+750,-415]
rref(Ab,.00000000001)
Press return to continue.
/* --------- DEBUG 1 ---------- */
Copy C[1] of Ab into C[5] of T
+6.0000 +5.0000 +3.0000 +2.0000 +4.0000 +1.0000 +6.0000
+7.0000 +1.0000 +0.0000 +0.0000 +0.0000 +0.0000 -0.4564
+6.0000 +0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +0.0064
+0.0000 +0.0000 +0.0000 +1.0000 +0.0000 +0.0000 -0.1281
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 +0.0000 -0.3933
+0.0000 -0.0000 -0.0000 -0.0000 -0.0000 +1.0000 +0.1883
T[R0][5]
+6.0000 +1.0000 +2.0000 +3.0000 +4.0000 +5.0000 +6.0000
+7.0000 +0.0000 +0.0000 +0.0000 +0.0000 +1.0000 +0.0000
+0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000
+0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000
+0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000
+0.0000 +0.0000 +0.0000 +0.0000 +0.0000 -0.0000 +0.0000
Press return to continue.
/* --------- DEBUG 1 ---------- */
Copy C[2] of Ab into C[3] of T
+6.0000 +5.0000 +3.0000 +2.0000 +4.0000 +1.0000 +6.0000
+7.0000 +1.0000 +0.0000 +0.0000 +0.0000 +0.0000 -0.4564
+6.0000 +0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +0.0064
+0.0000 +0.0000 +0.0000 +1.0000 +0.0000 +0.0000 -0.1281
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 +0.0000 -0.3933
+0.0000 -0.0000 -0.0000 -0.0000 -0.0000 +1.0000 +0.1883
T[R0][3]
+6.0000 +1.0000 +2.0000 +3.0000 +4.0000 +5.0000 +6.0000
+7.0000 +0.0000 +0.0000 +0.0000 +0.0000 +1.0000 +0.0000
+0.0000 +0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +0.0000
+0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000
+0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000
+0.0000 +0.0000 +0.0000 -0.0000 +0.0000 -0.0000 +0.0000
Press return to continue.
/* --------- DEBUG 1 ---------- */
Copy C[3] of Ab into C[2] of T
+6.0000 +5.0000 +3.0000 +2.0000 +4.0000 +1.0000 +6.0000
+7.0000 +1.0000 +0.0000 +0.0000 +0.0000 +0.0000 -0.4564
+6.0000 +0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +0.0064
+0.0000 +0.0000 +0.0000 +1.0000 +0.0000 +0.0000 -0.1281
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 +0.0000 -0.3933
+0.0000 -0.0000 -0.0000 -0.0000 -0.0000 +1.0000 +0.1883
T[R0][2]
+6.0000 +1.0000 +2.0000 +3.0000 +4.0000 +5.0000 +6.0000
+7.0000 +0.0000 +0.0000 +0.0000 +0.0000 +1.0000 +0.0000
+0.0000 +0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +0.0000
+0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +0.0000
+0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000
+0.0000 +0.0000 -0.0000 -0.0000 +0.0000 -0.0000 +0.0000
Press return to continue.
/* --------- DEBUG 1 ---------- */
Copy C[4] of Ab into C[4] of T
+6.0000 +5.0000 +3.0000 +2.0000 +4.0000 +1.0000 +6.0000
+7.0000 +1.0000 +0.0000 +0.0000 +0.0000 +0.0000 -0.4564
+6.0000 +0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +0.0064
+0.0000 +0.0000 +0.0000 +1.0000 +0.0000 +0.0000 -0.1281
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 +0.0000 -0.3933
+0.0000 -0.0000 -0.0000 -0.0000 -0.0000 +1.0000 +0.1883
T[R0][4]
+6.0000 +1.0000 +2.0000 +3.0000 +4.0000 +5.0000 +6.0000
+7.0000 +0.0000 +0.0000 +0.0000 +0.0000 +1.0000 +0.0000
+0.0000 +0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +0.0000
+0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +0.0000
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 +0.0000 +0.0000
+0.0000 +0.0000 -0.0000 -0.0000 -0.0000 -0.0000 +0.0000
Press return to continue.
/* --------- DEBUG 1 ---------- */
Copy C[5] of Ab into C[1] of T
+6.0000 +5.0000 +3.0000 +2.0000 +4.0000 +1.0000 +6.0000
+7.0000 +1.0000 +0.0000 +0.0000 +0.0000 +0.0000 -0.4564
+6.0000 +0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +0.0064
+0.0000 +0.0000 +0.0000 +1.0000 +0.0000 +0.0000 -0.1281
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 +0.0000 -0.3933
+0.0000 -0.0000 -0.0000 -0.0000 -0.0000 +1.0000 +0.1883
T[R0][1]
+6.0000 +1.0000 +2.0000 +3.0000 +4.0000 +5.0000 +6.0000
+7.0000 +0.0000 +0.0000 +0.0000 +0.0000 +1.0000 +0.0000
+0.0000 +0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +0.0000
+0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +0.0000
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 +0.0000 +0.0000
+0.0000 +1.0000 -0.0000 -0.0000 -0.0000 -0.0000 +0.0000
Press return to continue.
/* --------- DEBUG 1 ---------- */
Copy C[6] of Ab into C[6] of T
+6.0000 +5.0000 +3.0000 +2.0000 +4.0000 +1.0000 +6.0000
+7.0000 +1.0000 +0.0000 +0.0000 +0.0000 +0.0000 -0.4564
+6.0000 +0.0000 +1.0000 +0.0000 +0.0000 +0.0000 +0.0064
+0.0000 +0.0000 +0.0000 +1.0000 +0.0000 +0.0000 -0.1281
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 +0.0000 -0.3933
+0.0000 -0.0000 -0.0000 -0.0000 -0.0000 +1.0000 +0.1883
T[R0][6]
+6.0000 +1.0000 +2.0000 +3.0000 +4.0000 +5.0000 +6.0000
+7.0000 +0.0000 +0.0000 +0.0000 +0.0000 +1.0000 -0.4564
+0.0000 +0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +0.0064
+0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +0.0000 -0.1281
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 +0.0000 -0.3933
+0.0000 +1.0000 -0.0000 -0.0000 -0.0000 -0.0000 +0.1883
Press return to continue.
Ab : The work of sort_c_mR();
+0.0000 +0.0000 +0.0000 +0.0000 +1.0000 -0.4564
+0.0000 +0.0000 +1.0000 +0.0000 +0.0000 +0.0064
+0.0000 +1.0000 +0.0000 +0.0000 +0.0000 -0.1281
+0.0000 +0.0000 +0.0000 +1.0000 +0.0000 -0.3933
+1.0000 -0.0000 -0.0000 -0.0000 -0.0000 +0.1883
Press return to continue.
gj_TP_mR(Ab) :
+1.00000 -0.00000 -0.00000 -0.00000 -0.00000 +0.18826
+0.00000 +1.00000 +0.00000 +0.00000 +0.00000 -0.12813
+0.00000 +0.00000 +1.00000 +0.00000 +0.00000 +0.00638
+0.00000 +0.00000 +0.00000 +1.00000 +0.00000 -0.39331
+0.00000 +0.00000 +0.00000 +0.00000 +1.00000 -0.45644
Press return to continue
Press X to stop