je te rajouterai des que j'ai un peu de tmeps un petit source pour te montrer comment utiliser les matrice .. j'ia fias ç_a y a kk temps.. mais ça le merite de marcher et tres rapide pour des matrices 4x4
Le HEader
/*BONNIN Maleaume 2003,2004
mail: maleaume@msn.com
*/
#ifndef __MATRICE_H__
#define __MATRICE_H__
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define DIM 2
#ifndef __PI__
#define __PI__
#define PI 3.14159265358979
#endif
typedef struct Matrice{
int n;
int m;
double ** mat;
}matrice;
void mineur_matrice(matrice * min_A,matrice* A, int n,int m);
void inverse_matrice(matrice * A_1, matrice * A);
void trans_comatrice(matrice * CoA, matrice *A);
double determinant_matrice (matrice *m);
matrice * mineur_matrice(matrice* A, int n,int m);
void affiche_matrice(matrice* A);
void trans_matrice(matrice * C,matrice* A);
void mult_matrice(matrice * C,matrice* A, matrice *B);
void soust_matrice(matrice* C,matrice* A, matrice* B);
matrice* add_matrice(matrice* A, matrice* B);
matrice* alloc_matrice (int n,int m);
matrice* identite_matrice(int n);
matrice * mult_by_scal_matrice(matrice * A, double k);
void free_matrice(matrice *);
#endif
le source
/* BONNIN Maleaume 2003,2004
mail: maleaume@msn.com
*/
#include "matrice.h"
//--------------------------------------------------------------
// Fonction d'allocation d'une matrice (n,n)
// Remarque : on désalloue en cas d’échec en cours !
//--------------------------------------------------------------
matrice* alloc_matrice (int n,int m)
{
matrice * M = (matrice *)calloc(1,sizeof(matrice));
M->n = n;
M->m = m;
double **a;
a=(double **)calloc(n,sizeof(double *));
if (a!=NULL) {
for (int i=0; i<n; i++) {
a[i]=(double *)calloc(m,sizeof(double));
if (a[i]==NULL) {
for (int j=0; j<i; j++) free((void *)a[j]);
free((void *)a);
return NULL;
}
}
}
M->mat= a;
return M;
}
//--------------------------------------------------------------
// Fonction de désallocation d'une matrice (n,n)
//--------------------------------------------------------------
void free_matrice (matrice *M)
{
if (M->mat!=NULL) {
for (int i=0; i<M->n; i++) if (M->mat[i]!=NULL) free((void *)M->mat[i]);
free((void *)M->mat);
}
}
matrice* identite_matrice(int n){
matrice * C = alloc_matrice(n,n);
for(int i =0; i<n;i++)
C->mat[i][i] = 1;
return C;
}
matrice* add_matrice(matrice* A, matrice* B)
{
if(A->n != B->n || A->m!= B->m){
printf("erreur dans la dimension des matrices");
return NULL;
}
matrice* C= alloc_matrice(A->n,A->m);
for(int i = 0; i < C->n; i++)
for(int j = 0; j < C->m; j++)
C->mat[i][j] = A->mat[i][j] + B->mat[i][j];
return C;
}
void soust_matrice(matrice * C, matrice* A, matrice* B)
{
if(A->n != B->n || A->m!= B->m){
printf("erreur dans la dimension des matrices");
exit(0);
}
//matrice* C= alloc_matrice(A->n,A->m);
for(int i = 0; i < C->n; i++)
for(int j = 0; j < C->m; j++)
C->mat[i][j] = A->mat[i][j] - B->mat[i][j];
//return C;
}
matrice * mult_by_scal_matrice(matrice * A, double k){
if(k == 0.0){
matrice * C = alloc_matrice(A->n, A->m);
return C;
}
for(int i = 0; i<A->n;i++)
for(int j = 0; j < A->m;j++)
if(A->mat[i][j]!=0) A->mat[i][j]*=k;
return A;
}
void mult_matrice(matrice* C,matrice* A, matrice *B)
{
if(A->m != B->n)
{ printf("Erreur dans ladimensoin des matrices\n");
}
double somme;
// matrice * C = alloc_matrice(A->n,B->m);
for(int i = 0; i < A->n; i++)
for(int j = 0; j < B->m; j++){
somme = 0;
for(int k=0; k < A->m; k++)
somme += A->mat[i][k]*B->mat[k][j];
C->mat[i][j] = somme;
}
}
void trans_matrice(matrice * C,matrice* A)
{
//matrice* C= alloc_matrice(A->m,A->n);
for(int i = 0; i < A->n; i++)
for(int j = 0; j < A->m; j++)
C->mat[j][i] = A->mat[i][j];
//return C;
}
void affiche_matrice(matrice* A)
{
for(int i = 0; i < A->n; i++){
printf("\n");
for(int j = 0; j < A->m; j++)
printf("%f\t",A->mat[i][j]);
}
printf("\n\n");
}
double determinant_matrice (matrice *m)
{
if(m->n != m->n){
printf("La matrice n'est pas une matrice carrée\n");
return NULL;
}
matrice * sous_m = alloc_matrice(m->n-1, m->m-1);;
int signe=1;
double det=0;
if (m->n==1) return (m->mat[0][0]);
if(m->n==2)
return m->mat[0][0]*m->mat[1][1]-m->mat[1][0]*m->mat[0][1];
for (int j=0;j<m->n;j++){
// for(int j=0;j< m->m; j++){
mineur_matrice(sous_m,m,0,j);
signe=((j)%2 == 0)?1:-1;
det+=signe*m->mat[0][j]*determinant_matrice(sous_m);
}
free_matrice(sous_m);
free(sous_m);
return(det);
}
void mineur_matrice(matrice * min_A, matrice* A, int n,int m)
{
int ii = 0;
int jj = 0;
if(A->n == 1)
min_A->mat[0][0] = A->mat[0][0];
for(int i = 0; i < A->n; i++){
for(int j = 0; j < A->m; j++){
if(i != n && j!=m){
min_A->mat[ii][jj] = A->mat[i][j];
jj++;
}
}
if(i!=n) {ii++; jj=0;}
}
}
void trans_comatrice(matrice * CoA_t, matrice *A)
{
char signe = 1;
matrice* CoA = alloc_matrice(A->n,A->m);
matrice* min_A = alloc_matrice(A->n-1, A->m-1);
if(determinant_matrice(A)==0){
printf("La matrice est non inversible,sont déterminant est nul\n");
exit(0);
}
for(int i = 0; i < CoA->n; i++)
for(int j = 0; j < CoA->m; j++){
signe = ((i+j)%2 ==0)?1:-1;
mineur_matrice(min_A, A, i, j);
CoA->mat[i][j] = signe*determinant_matrice(min_A);
}
trans_matrice(CoA_t, CoA);
free_matrice(CoA);free(CoA);
free_matrice(min_A);free(min_A);
}
void inverse_matrice(matrice * A_1, matrice * A)
{
if(A->n != A->m){
printf("La matrice n'estpas une matrice carée, et ne peux donc etre inversée\n");
exit(0);
}
double det = determinant_matrice(A);
if(det == 0){
printf("Votre matrice est non inversible car de déterminant nul\n");
exit(0);
}
trans_comatrice(A_1,A);
for(int i = 0; i< A_1->n; i++)
for(int j = 0; j< A_1->m;j++)
A_1->mat[i][j]/=det;
}