Calculs de normes matricielles

Calculs de normes matricielles - C - Programmation

Marsh Posté le 17-10-2007 à 15:47:38    

Salut tout le monde,
 
Je débute en programmation et j'aurais besoin d'un petit coup de pouce. Je dois calculer les normes matricielles (norme 1 et norme infinie). J'ai essayé de le sprogrammer, à la compilation ça passe, par contre je ne comprends pas pourquoi ça plante à l'exécution. Je vous mets ici mon code pour la norme 1 (pour la norme infinie il suffit juste de permuter les indices) :
 
main.c :

Code :
  1. #include <stdio.h>
  2. #include "functions.h"
  4. int main ()
  6. { int i,j,c,nrows,ncols;
  7.     double **mat;
  8. double res1,res3;
  10. printf("Choose the norm you want to calculate :\n 1. l1 norm\n 2. linf norm\n " );
  12. scanf("%d",&c);
  14. printf("Please input the number of rows and columns for the matrix. These numbers must be the same :\n" );
  16.     scanf("%d %d",&nrows,&ncols);
  18. mat = dmatrix(nrows,ncols); /* dynamic allocation for the matrix*/
  20. printf("Please input the %d*%d matrix elements :\n",nrows,ncols);
  22. /*nested loop*/
  24. for (i=0;i<nrows;i++)
  26.   {
  27.    for (j=0;j<ncols;j++)
  29.     {
  30.      scanf("%lf",&mat[i][j]);
  31.     }
  33.   }
  36. switch (c)
  38. {case 1: res1=norm1(mat,nrows,ncols);
  39.    printf("The result is %lf", res1);
  40.    break;
  43.      case 2: res3=norminf(mat,nrows,ncols);
  44.    printf("The result is %lf", res3);
  45.    break;
  47. }
  49. /* frees the memory for the matrix */
  51. for (i=0;i<nrows;i++)
  53.  {
  54.   free(mat[i]);
  55.  }
  57. free(mat);
  59. }


 
functions.c (contient les fonctions qui vont calculer les normes)
 

Code :
  1. #include <math.h>
  3. double norminf(double **mat ,int nrows, int ncols)
  4. {
  5. double sum, res;
  6. int i,j;
  7. res=0;
  8. sum=0;
  10. for(i=0;i<=nrows;i++)
  12.  {
  13.    for(j=0;j<=ncols;j++)
  15.   {
  16.     sum = sum + fabs(mat[i][j]);
  18.     if(sum>res)
  19.    {
  20.      res=sum;
  21.    }
  23.   }
  25.  }
  27. return(res);
  29. }
  32. double norm1(double **mat ,int nrows, int ncols)
  33. {
  34. double sum, res;
  35. int i,j;
  36. res=0;
  37. sum=0;
  39. for(j=0;j<=ncols;i++)
  41.  {   
  42.   for(i=0;i<=nrows;i++)
  44.   {
  45.    sum = sum + fabs(mat[i][j]);
  47.    if(sum>res)
  48.     {
  49.      res=sum;
  50.     }
  52.   }
  54.  }
  56. return(res);
  58. }


 
memory.c (gestion de la mémoire dynamique)

Code :
  1. #include <stdio.h>
  3. /* dynamic memory allocation function for a matrix */
  5. double ** dmatrix (int nrows, int ncols)
  7. {
  8. double **mat; /* matrix */
  9. int i;
  11. mat = (double **)malloc(nrows*sizeof(double*));
  13. for(i=0;i<nrows;i++)
  15.  {
  16.   mat[i] = (double *)malloc(ncols*sizeof(double));
  17.  }
  19. return mat;
  21. }


J'ai fait aussi un function.h pour les prototypes des fonctions. Est-ce que des programmeurs un peu plus expérimentés pourraient m'aider ? Merci bcp

Reply

Marsh Posté le 17-10-2007 à 15:47:38   

Reply

Marsh Posté le 17-10-2007 à 17:33:37    

Il faut avoir inférieur strictement au lieu d'inférieur ou égal dans les boucles de fonctions.c, car la matrice a été allouée avec les bornes excluses.

Reply

Sujets relatifs:

Leave a Replay

Make sure you enter the(*)required information where indicate.HTML code is not allowed