Instructions non effectué à partir d'un certain rang

Instructions non effectué à partir d'un certain rang - C - Programmation

Marsh Posté le 25-11-2011 à 22:14:51    

Bonjour,
j'ai un petit problème que j'ai repéré avec gdb mais je ne sais pas le résoudre. Pour vous expliquer, le programme marche correctement jusqu'à l'itération k=42 (incluse). Mais lorsque k passe à 43, les boucles sont effectuées mais aucunes instructions n'est réalisée.
Avez vous une raison ?
Merci !
Voila le code gdb :
 
260    x[k+1][i]=x[k][i]+a[k]*r[k][i];
(gdb)  
261    printf("x [%d] [%d] :%f\n",k,i,x[k][i]);
(gdb)  
x [43] [0] :0.989438
242   for (i=0;i<=n-1;i++)
(gdb)  
265    printf("\n" );
(gdb)  
 
266   k++;
(gdb)  
267   nbiter++;
(gdb)  
268  }while((k<=itermax-1)/*&&(norme2(r,k-1,n)/norme2(b,0,n)>e)*/);
(gdb)  
232   for (i=0;i<=n-1;i++)
(gdb)  
242   for (i=0;i<=n-1;i++)
(gdb)  
265    printf("\n" );
(gdb)  
 
266   k++;
(gdb)  
 
 
Voila le programme principale :
 

Code :
  1. void main (void)
  2. {
  3. double x[Max][Maxdim];
  4. double r[Max][Maxdim];
  5. double A[Maxdim][Maxdim];
  6. double b [1][Maxdim];
  7. double a [Maxdim];
  8. int i,n,w,z,k=0,nbiter=0,itermax,itermax2;
  9. double e;
  10. double Ar[Max][Maxdim];
  13. e=0.0000001;
  14. n=3;
  15. itermax=44;
  16. hilbert(A,b,n);             // hilbert dordre n avec sa fonction
  17.  for(w=0;w<n;w++)
  18.  {
  19.   for(z=0;z<n;z++)
  20.   {
  21.    x[w][z]=0;
  22.   }
  23.  }
  26.  do
  27. {
  28.  for (i=0;i<=n-1;i++)
  29.  {
  30.   b[0][0]=11.0/6;
  31.   b[0][1]=13.0/12;
  32.   b[0][2]=47.0/60;
  33.   r[k][i]=b[0][i]-indice_i_du_produit_matrice_vecteur(A,x,k,n,i);
  34.   printf("b[0][%d] : %f\n",i,b[0][i]);
  35.   printf("r[%d][%d] : %f\n",k,i,r[k][i]);
  36.  }
  37.  for (i=0;i<=n-1;i++)
  38.  {
  39.   produit_matrice_vecteur_2(A,r, Ar, k, n);
  41.   printf("Ar[%d][%d] : %f\n",k,i,Ar[k][i]);
  44.   a[k]=prodscal(r,r,k,n)/prodscal(r,Ar,k,n);
  45.   printf("prodscal(r,r,%d,n): %.10f\n",k,prodscal(r,r,k,n));
  46.   printf("prodscal(r,Ar,%d,n): %.10f\n",k,prodscal(r,Ar,k,n));
  47.   printf("a[%d]: %f\n",k,a[k]);
  48.   x[k+1][i]=x[k][i]+a[k]*r[k][i];
  49.   printf("x [%d] [%d] :%f\n",k,i,x[k][i]);
  52.  }
  53.   printf("\n" );
  54.  k++;
  55.  nbiter++;
  56. }while((k<=itermax-1)/*&&(norme2(r,k-1,n)/norme2(b,0,n)>e)*/);


Message édité par gilou le 26-11-2011 à 20:03:39
Reply

Marsh Posté le 25-11-2011 à 22:14:51   

Reply

Marsh Posté le 26-11-2011 à 16:02:48    

ton code ne compile pas. Poste un code qui compile.

Reply

Marsh Posté le 30-11-2011 à 19:02:13    

Je ne vois pas le probléme.
A tu essayé avec  while(k<itermax);

Reply

Sujets relatifs:

Leave a Replay

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