matrice.c

/* Putain de matrice des mes couilles */

#include <stdio.h>
#include <stdlib.h>
#include <string.h>


#define NMAX 100

typedef struct {
    int M[NMAX][NMAX];
    int N;
} TMatrice;


void afficherMatrice( TMatrice const m )
{
    int i, j;

    for( i = 0; i < m.N; ++i ) {
        for( j = 0; j < m.N; ++j )
            printf("%2d ", m.M[i][j]);

        printf("\n");
    }
}

void mulMat(
    TMatrice const m1 ,
    TMatrice const m2 ,
    TMatrice * m3 )
{
    int i, j, k, s;

    if ( m3 == NULL || m1.N != m2.N )
        return;

    for( i = 0; i < m1.N; ++i ) {
        for( j = 0; j < m1.N; ++j ) {
            for( s = 0, k = 0; k < m1.N; ++k )
                s += ( m1.M[i][k] * m2.M[k][j] );

            m3->M[i][j] = s;
        }
    }

    m3->N = m1.N;
}

int main()
{
    TMatrice m3;
    TMatrice m1 = {
        { { 1, 2, 1 },
          { 3, 2, 2 },
          { 0, 5, 1 } },
        3 // 3*3
    };
    TMatrice m2 = {
        { { 1, 3, 1 },
          { 2, 2, 0 },
          { 3, 1, 1 } },
        3 // 3*3
    };

    printf("M1:\n");
    afficherMatrice(m1);
    printf("M2:\n");
    afficherMatrice(m2);

    printf("\nM3 = M1 * M2:\n");
    mulMat(m1, m2, &m3);
    afficherMatrice(m3);

    return 0;
}