Mathc matrices/Fichiers c : xi mr


Étude du code


Installer et compiler ce fichier dans votre répertoire de travail.

ai_mr.c
/* ------------------------------------ */
/*  Save as :   ai_mr.c                 */
/* ------------------------------------ */
#include "v_a.h"
/* ---------------------------------------------------------------------------
 Do    : Dynamically allocate a multidimensional array.

            (see : FAQ of the comp.lang.c group)

        You can keep the array's contents contiguous,

                 r0    r1    r2    ... rn
                 R_000|C_xxx|0_xxx|...|0_xxx


                R =  Number of rows.
                C =  Number of columns.

        The declaration of the sizes into the matrices, it is my work.

        So be careful.

        The first row and the first column are not used.

                      *********************

        The size of the row    of the matrix is into  A[R_SIZE][C0] = A[0][0]
        The size of the column of the matrix is into  A[C_SIZE][C0] = A[1][0]

                      *********************

        The first element of the matrix is z =(A[1][1],A[1][1])

        For a 10x10 matrix the last element is z = (A[10][10],A[10][10])

                      *********************

        If you want to duplicate a matrix A, you must use :

        double **T;

        T = i_RC_mR(  A[R_SIZE][C0],
                      A[C_SIZE][C0]);

        f_mR(T);

   -------------------------------------------------------------------------- */
   
/* ------------------------------------ */
/* Same as i_RC_mR() but work with the 
   functions.   
          
   ex :   i_mR(rsize_R(M),csize_R(M));  */
/* ------------------------------------ */ 
double **xi_mR(
int      r,
int      c
)
{
double **A;
int     ar;
int     ac;
int      i;

            if(r<R1||c<C1)
            {
            printf(" The size of the matrix must be positive integers.\n\n");
            printf(" double **i_mR(); \n\n");
            fflush(stdout);
            getchar();
            exit(EXIT_FAILURE);
            }

         ar = r + R1;
         ac = c + C1;

         A = malloc(ar * sizeof(*A));
             if(!A)
            {
            printf(" I was unable to allocate the memory you requested.\n\n");
            printf(" double **i_mR(); \n\n");
            printf(" **A = malloc(ar * sizeof(*A));\n\n");
            fflush(stdout);
            getchar();
            exit(EXIT_FAILURE);
            }

      A[0] = malloc(ar * ac * sizeof(**A) );
             if(!A[0])
            {
            printf(" I was unable to allocate the memory you requested.\n\n");
            printf(" double **i_mR();\n\n");
            printf(" A[0] = malloc(ar * ac * sizeof(**A) );\n\n");
            fflush(stdout);
            getchar();
            exit(EXIT_FAILURE);
            }

	for(i=R1; i<ar; i++) A[i] = A[0]+i*ac;

/* ----------- Give a value to the zero row and the zero column ------------- */

    A[R_SIZE][C0] = ar;
    A[C_SIZE][C0] = ac;

    for(r=R2; r<A[R_SIZE][C0]; r++)

        A[r][0] = 0.;

    for(c=C1; c<A[C_SIZE][C0]; c++)

        A[0][c] = c;

    m0_mR(A);

return(A);
}
/* ------------------------------------ */ 
/* ------------------------------------ */
void fun(int r,int c)
{
double **A     = r_mR(xi_mR(r,c),9);

  clrscrn();

  printf(" A[R%d,C%d]:          (p_mR();)\n",rsize_R(A),csize_R(A));
  p_mR(A,S4,P0,C8);

  printf("\n\n\n");
  printf(" A[R%d,C%d]:          (pall_mR();)\n",rsize_R(A),csize_R(A));
  pall_mR(A,S4,P0,C9);
 
  f_mR(A);
}
/* ------------------------------------ */
int main(void)
{
time_t t;

  srand(time(&t));
  
   do 
        fun(rp_I(R4),rp_I(C8));
        
    while(stop_w());
        
  return 0;
}
/* ------------------------------------ */
/* ------------------------------------ */


La fonction i_mR(); existant déjà dans la librairie je l'ai ici renommé xi_mR();

Nous connaissons la fonction pall_mR(); qui permet afficher la colonne zéro et la ligne zéro.


Étudions quelques sorties écran.

Exemple de sortie écran :

A[R4,C8]:          (p_mR();)

  -1   -3   +6   -3   -5   -3   -3   +4 
  -1   +8   -7   -9   +2   -1   +8   -5 
  -7   +6   -5   +6   +8   +2   -3   -1 
  -5   +8   -9   -5   -7   -1   -3   -5 




 A[R4,C8]:          (pall_mR();)

  +5   +1   +2   +3   +4   +5   +6   +7   +8 

  +9   -1   -3   +6   -3   -5   -3   -3   +4 
  +0   -1   +8   -7   -9   +2   -1   +8   -5 
  +0   -7   +6   -5   +6   +8   +2   -3   -1 
  +0   -5   +8   -9   -5   -7   -1   -3   -5 


 Press return to continue
 Press X      to stop


Nous allons étudier la fin de la fonction.

/* ----------- Give a value to the zero row and the zero column ------------- */

    A[R_SIZE][C0] = ar; /* On copie le nombre de  lignes  dans la case A[0;0]  */
    A[C_SIZE][C0] = ac; /* On copie le nombre de colonnes dans la case A[1;0]  */

    for(r=R2; r<A[R_SIZE][C0]; r++) /* On met des zéro dans la colonne zéro    */
                                    /* à partir de la ligne deux               */ 
        A[r][0] = 0.;

    for(c=1; c<A[C_SIZE][C0]; c++) /* On crée un index dans la ligne zéro      */

        A[0][c] = c;

    m0_mR(A);            /* On initialise la matrice comme une matrice zéro   */

return(A);               /* On retourne l'adresse de la matrice au programme   */
}


Passons maintenant au début de la fonction.

           if(r<R1||c<C1)            /* On vérifie que les tailles des lignes */ 
            {                        /* et des colonnes sont positives        */

            printf(" The size of the matrix must be positive integers.\n\n");
            printf(" double **i_mR(); \n\n");
            fflush(stdout);
            getchar();
            exit(EXIT_FAILURE);
            }

         ar = r + R1;                /* On rajoute  une  ligne, la  ligne  zéro */
         ac = c + C1;                /* On rajoute une colonne, la colonne zéro */


Maintenant vérifions si on a suffisamment d'espace pour le nombre de lignes

         A = malloc(ar * sizeof(**A));    /* Donne de l'espace pour ar lignes */
             if(!A)                       /* Si cela à échoué...             */ 
            {
            printf(" I was unable to allocate the memory you requested.\n\n");
            printf(" double **i_mR(); \n\n");
            printf(" **A = malloc(ar * sizeof(**A));\n\n");
            fflush(stdout);
            getchar();
            exit(EXIT_FAILURE);
            }

sizeof(**A) indique que l'on travaille avec des doubles. On crée donc un espace pour ar double.


Maintenant vérifions si on a suffisamment d'espace pour le nombre de colonnes

      A[0] = malloc(ar * ac * sizeof(**A) );  /* Donne de l'espace pour ar*ac colonnes */
             if(!A[0])                        /* Si cela à échoué...                   */ 
            {
            printf(" I was unable to allocate the memory you requested.\n\n");
            printf(" double **i_mR();\n\n");
            printf(" A[0] = malloc(ar * ac * sizeof(**A) );\n\n");
            fflush(stdout);
            getchar();
            exit(EXIT_FAILURE);
            }

	for(i=R1; i<ar; i++) A[i] = A[0]+i*ac; /* On associe à chaque ligne les ac colonnes */