Discussion :

[salseropom]

Bonjour, après mon malloc j'ai une erreur

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
main.c 226 : malloc() failed

Voici mon code

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
 
double ** Particles=malloc(nbprocs*sizeof(*Particles));
if(Particles==NULL)
{
  MEMERROR;
  MPI_Finalize();
  exit(EXIT_FAILURE);
}
 
if(rank<restp1)
  {
    firstpart[rank] = (rank-1)*nbparp1;
    lastpart[rank]  = nbparp1*rank-1;  
  }
  else /* rank >= restp1 */
  {
    firstpart[rank] = nbpar*(rank-1)+rest;
    lastpart[rank]  = nbpar*rank+restm1; 
  }
  sizeParticles[rank]=neq*(lastpart[rank]-firstpart[rank]+1); /* ici  sizeParticles[rank]==780 ou 806 suivant la valeur de rank */
 
 
  for(i=1;i<nbprocs;++i)  /* nbprocs == 26 */
  {
    Particles[i]=malloc(sizeParticles[i]*sizeof(*Particles[i]));
    if(Particles[i]==NULL)
    {
      printf("i = %d\n",i);  /* echec quand i == 14 */
      MEMERROR;
      MPI_Finalize();
      exit(EXIT_FAILURE);
    }
  } /* end for(i=1;i<nbprocs;++i) */

Merci.

[Jean-Marc.Bourguet]

Il n'y a pas 226 lignes, ce n'est donc pas ton code complet.

Tu corromps vraisemblablement les structures de données de malloc.

[Fabllot]

Bonjour,
Es-tu certain que sizeParticles[i] est strictement positif ?

Sinon dans les malloc, j'ai plutôt l'habitude de me référer au type de l'objet et non à l'objet en lui-même, mais je ne sais pas si cela change quelque chose ...

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
double ** Particles =  (double*) malloc (nbprocs * sizeof (double*));
[...]
Particles[i] = malloc (sizeParticles[i] * sizeof (double));

[zooro]

Sinon dans les malloc, j'ai plutôt l'habitude de me référer au type de l'objet et non à l'objet en lui-même, mais je ne sais pas si cela change quelque chose ...

En utilisant la variable, et non son type, la maintenance est plus facile, puisque tu peux changer le type sans toucher au reste du code.

[Emmanuel Delahaye]

Bonjour,
Es-tu certain que sizeParticles[i] est strictement positif ?

Sinon dans les malloc, j'ai plutôt l'habitude de me référer au type de l'objet et non à l'objet en lui-même, mais je ne sais pas si cela change quelque chose ...

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
double ** Particles =  (double*) malloc (nbprocs * sizeof (double*));
[...]
Particles[i] = malloc (sizeParticles[i] * sizeof (double));

Mais que c'est laid...

http://emmanuel-delahaye.developpez....tes.htm#malloc

[salseropom]

Re,
j'ai identifié le problème :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
 
double ** Particles=malloc(nbprocs*sizeof(*Particles));
if(Particles==NULL)
{
  MEMERROR;
  MPI_Finalize();
  exit(EXIT_FAILURE);
}
 
if(rank<restp1)
  {
    firstpart[rank] = (rank-1)*nbparp1;
    lastpart[rank]  = nbparp1*rank-1;  
  }
  else /* rank >= restp1 */
  {
    firstpart[rank] = nbpar*(rank-1)+rest;
    lastpart[rank]  = nbpar*rank+restm1; 
  }
  sizeParticles[rank]=neq*(lastpart[rank]-firstpart[rank]+1); /* ici  sizeParticles[rank]==780 ou 806 suivant la valeur de rank */
 
 /* jusqu'ici tout va bien : sizeParticles[rank] me renvoie bien les bonnes valeurs. Aucun problème */
 
 
  for(i=1;i<nbprocs;++i)  /* nbprocs == 26 */
  {
    Particles[i]=malloc(sizeParticles[i]*sizeof(*Particles[i]));
    if(Particles[i]==NULL)
    {
      printf("i = %d\n",i);  /* echec quand i == 14 */
      MEMERROR;
      MPI_Finalize();
      exit(EXIT_FAILURE);
    }
  } /* end for(i=1;i<nbprocs;++i) */

ça bug lorsque i==14. Or 14 est exactement le nombre de noeuds de ma grappe de calcul. Je suis allé sur une autre grappe et ça bug à 5 où 5 est le nombre de noeuds de l'autre grappe. Donc je pense que le problème ne vient pas du malloc. Par contre, je ne sais pas d'où il pourrait venir... Avez-vous une idée ?

Merci encore.

[salseropom]

Ca y est, j'ai enfin trouvé mon bug !!

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
  for(i=1;i<nbprocs;++i)
  {  
    if(i<restp1)
    {
      firstpart[i] = (i-1)*nbparp1;
      lastpart[i]  = nbparp1*i-1;  
    }
    else /* i >= restp1 */
    {
      firstpart[i] = nbpar*(i-1)+rest;
      lastpart[i]  = nbpar*i+restm1; 
    }
 
    sizeParticles[i]=neq*(lastpart[i]-firstpart[i]+1);
  }

et là ça marche.