[C++] petit pb de seg. fault

petit pb de seg. fault [C++] - C++ - Programmation

Marsh Posté le 26-07-2005 à 17:26:38    

Bonjour,
 
je suis sur cette erreur depuis un certain temps et j'ai peur d'avoir perdu tout recul pour la résoudre.
 

Code :
  1. #include "fm_volume.h"
  3. using namespace std;
  5. int main ()
  6. {
  7.   FmVolume* firstVolume = new FmVolume("1exam10_06_02" );
  9.   cout<<"yop"<<endl;
  10.   return 1;
  11. }


 
la classe FmVolume étant définie dans fm_volume.h.
Mon probleme c'est à l'execution : superbe segmentation fault apres "yop" ... :pt1cable:  
Se peut-il qu'il s'agisse d'erreurs de la classe FmVolume qui seraient "descendues" ?
Merci pour vos suggestions...


Message édité par nifnef le 26-07-2005 à 18:17:00
Reply

Marsh Posté le 26-07-2005 à 17:26:38   

Reply

Marsh Posté le 26-07-2005 à 17:31:18    

# #include <cstdlib>
# #include <cstdio>
# #include <string>
# #include <fstream>
 
ça sert à quoi tout ça ?

Reply

Marsh Posté le 26-07-2005 à 17:38:28    

a pas grand chose effectivement...

Reply

Marsh Posté le 26-07-2005 à 18:03:29    

tu ne voudrais pas nous montrer ton code de fm_volume.h/cpp ?

Reply

Marsh Posté le 26-07-2005 à 18:07:08    

mouerf vous allez crier... mais si la seg fault est apres le yop, ca peut venir d'avant?
je vai le chercher...
le .h

Code :
  1. #ifndef __FM_VOLUME_H__
  2. #define __FM_VOLUME_H__
  4. #include <cstdlib>
  5. #include <cstdio>
  6. #include <string>
  7. #include <fstream>
  8. #include <iostream>
  10. #include <gsl/gsl_matrix.h>
  11. #include <gsl/gsl_vector.h>
  13. using namespace std;
  15. const float SIXTEEN_BITS_VALUE = 65535.0;
  16. const float TWELVE_BITS_VALUE_F = 32767.0;
  17. const int TWELVE_BITS_VALUE_I = 32767;
  18. const int HEIGHT_BITS_VALUE = 256;
  19. const int INIT_VALUE = -1;
  20. const string HDR_SUFF = ".hdr";
  21. const string IMG_SUFF = ".img";
  24. class FmVolume
  25. {
  26. private:
  27.   string fm_fname;
  28.   int fm_dims[4];
  29.   double fm_dim[3];
  30.   //gsl_vector* fm_pinfo;
  31.   //gsl_matrix* fm_mat;
  32.   string fm_descrip;
  33.   int fm_n;
  34.   gsl_matrix* fm_data[];
  35.   //fm_private;
  36.   int fm_error;
  38. public:
  39.   FmVolume(string filename);
  40.   ~FmVolume();
  42. private:
  43.   int loadIMG();
  44.   float read_header(FILE *fp_hdr, FILE *fp_img, int plane, int *flag);
  45.   int read_image(FILE *fp_img, unsigned short *data, int plane, int *flag, float scale);
  46. };
  48. #endif /* __FM_VOLUME_H__ */


le .cpp est plutot lourd...


Message édité par nifnef le 27-07-2005 à 15:19:11
Reply

Marsh Posté le 26-07-2005 à 18:10:51    

il y a une (au moins) seg fault dans le cpp mais je l'ai mise en commentaire et celle qui me bloque s'est alors dévoilée
si on peut voir ce qui se passe sans le cpp ca serait cool, il y a environ 400 lignes je me vois mal les poster.

Reply

Marsh Posté le 26-07-2005 à 18:32:18    

c'est quand tu veux pour donner du code ...
 
sinon RTFM, tu lances ton debugger et tu te débrouilles.

Reply

Marsh Posté le 26-07-2005 à 18:42:46    

Voici le cpp (arrangé et tronqué)

Code :
  1. #include "fm_volume.h"
  3. #include "HEADER.h"
  5. using namespace std;
  7. FmVolume::FmVolume(string filename)
  8. {
  9.   fm_fname = filename;
  10.   loadIMG();
  11. }
  14. FmVolume::~FmVolume()
  15. {
  16.   //gsl_vector_free(fm_pinfo); jamais alloué pour l'instant
  17.   //gsl_matrix_free(fm_mat); idem
  18.   for(int i=0; i<fm_dims[3]; i++)
  19.     gsl_matrix_free(fm_data[i]);
  20. }
  22. int
  23. FmVolume::loadIMG()
  24. {
  25.   int *flag = new int[2];
  26.   unsigned short* data;
  27.   int plane = INIT_VALUE;
  28.   float scale, thickness;
  29.   FILE* fp_hdr, *fp_img;
  30.  
  31.   if (fm_fname.size()== 0)
  32.     {
  33.       cerr << "No file selected."<<endl;
  34.       return EXIT_FAILURE;
  35.     }
  36.  
  37.   int hdrRead = (int)read_header(fp_hdr, fp_img, plane,flag);
  38.   if(hdrRead == 0.0)
  39.     return EXIT_FAILURE;
  41.   cout<<"Name : "<<fm_fname<<endl;
  42.   cout<<"Number of slices : "<<fm_dims[3]<<endl;
  43.   cout<<"Data type : "<<fm_dims[0]<<endl;
  45.  
  46. if(fm_error != 0)
  47.    {
  48.      cerr << "Error while loading " <<fm_fname<< endl;
  49.       return EXIT_FAILURE;
  50.    }
  53. for(int i=0; i<fm_dims[3]; i++)
  54.     fm_data[i] = gsl_matrix_alloc(fm_dims[1],fm_dims[2]);
  56.   data = new unsigned short[fm_dims[1]*fm_dims[2]];
  58. /*Recup des infos pour chaque plan*/
  60. gsl_matrix* temp = gsl_matrix_alloc(fm_dims[1],fm_dims[2]);
  62.   for (plane = 0; plane < fm_dims[3]; plane++)
  63.     {
  64.       scale = read_header(fp_hdr, fp_img, plane,flag);
  65.       if (scale > 0.0)
  66.         {
  67.   int imgRead = 0;//read_image(fp_hdr, fp_img, data,plane,flag,scale);
  68.    if (imgRead == EXIT_SUCCESS)
  69.     {
  70.       cout<<"plane "<<plane<<endl;
  71.       for (int i=0; i<fm_dims[1]; i++)
  72.  for (int j=0; j<fm_dims[2]; j++)
  73.       gsl_matrix_set(temp,i,j,data[j*fm_dims[1]+i]);
  74.      }
  75.     else
  76.      return EXIT_FAILURE;
  77.         }
  78.      fm_data[plane] = temp;
  79.     }
  80.   gsl_matrix_free(temp);
  82. cout<<"Slices width = "<<fm_dims[1]<<endl;
  83. cout<<"Slices height = "<<fm_dims[2]<<endl;
  85. thickness = fm_dim[2];
  86. cout<<"Slices thickness = "<<thickness<<endl;
  87. if (fp_hdr != NULL)
  88. fclose(fp_hdr);
  89. if (fp_img != NULL)
  90. fclose(fp_img);
  92. return EXIT_SUCCESS;
  93. }
  97. float
  98. FmVolume::read_header(FILE* fp_hdr, FILE* fp_img, int plane, int *flag)
  99. {
  100. struct hdr hdr;
  101.   string filename;
  102.   float scale;
  105.   if(plane<0)
  106.     {
  107.       fm_error=0;
  108.    
  109.        filename = fm_fname + HDR_SUFF;
  110.        fp_hdr = fopen(filename.c_str(),"rb" );
  112.        if (fp_hdr == NULL)
  113.         {
  114.           cerr<<"cannot open file " <<filename<<" for input."<<endl;
  115.   fm_error=1;
  116.           return 0.0;
  117.         }
  118.    
  119.        if(fread((char *)&hdr,1,sizeof(struct hdr),fp_hdr)!=sizeof(struct hdr))
  120.        {
  121.           cerr<<"file read error for file "<<filename<<endl;
  122.   fm_error=1;
  123.           return 0.0;
  124.   }
  126.       /*Copy required data from header struct*/
  127.       fm_dims[0] = hdr.bits;
  128.       fm_dims[1] = hdr.x_dim;
  129.       fm_dims[2] = hdr.y_dim;
  130.       fm_dims[3] = hdr.z_dim;
  131.      
  132.       fm_dim[0] = hdr.x_size;
  133.       fm_dim[1] = hdr.y_size;
  134.       fm_dim[2] = hdr.z_size;
  135.      
  136.       if(fm_dims[0]!=8 && fm_dims[0]!=16 && fm_dims[0]!=1)
  137.         {
  138.           cerr<<"Can't read data with "<<fm_dims[0]<<" bits/pixel ("<<filename<<" )"<<endl;
  139.   fm_error=1;
  140.           return 0.0;
  141.         }
  143.       flag[0]=hdr.glmin;
  144.       flag[1]=hdr.glmax;
  145.    
  146. /*Calculate global rescaling factor if image is 16 bit and compilation is 8 bit*/
  148. if(fm_dims[0]==16 && sizeof(unsigned short)==1)
  149.   {
  150.     if(hdr.glmin>=0)
  151.             {
  152.               if(hdr.glmax==0)
  153.  {
  154.    cerr<<"Header global maximum of zero is less than or equal to header global minimum ("<<filename<<" )\n";
  155.    fm_error=1;
  156.    return 0.0;
  157.  }
  158.               if(hdr.glmax>TWELVE_BITS_VALUE_I)
  159.  return SIXTEEN_BITS_VALUE/hdr.glmax;
  160.               else
  161.  return TWELVE_BITS_VALUE_F/hdr.glmax;
  162.             }
  163.     else
  164.     {
  165.       if(hdr.glmax<=-(TWELVE_BITS_VALUE_I + 1))
  166.     {
  167.       cerr<<"Header global maximum is out of possible range ("<<fm_fname<<" )"<<endl;
  168.       fm_error=1;
  169.       return 0.0;
  170.     }
  171.       return SIXTEEN_BITS_VALUE/(hdr.glmax+(TWELVE_BITS_VALUE_I + 1));
  172.             }
  173.   }
  174.       return 1.0;
  175.     }
  177.   else
  178.     {
  179.       fm_error=0;
  180.      
  181.       /* If initialization wasn't done--reject*/
  182.       if(fp_hdr==NULL)
  183. {
  184.   cerr<<"Coding error--failure to initialize read_header\n";
  185.   fp_img =NULL;
  186.   fm_error = 1;
  187.   return 0.0;
  188. }
  189.      
  190.       scale=1;
  191.      
  192.       /* If image file isn't already open, do so now and find plane*/
  193.       if(fp_img==NULL)
  194. {
  195.   filename = fm_fname+IMG_SUFF;
  196.   fp_img = fopen(filename.c_str(),"rb" );
  197.   if(fp_img == NULL)
  198.     {
  199.       cerr<<"cannot open file "<<filename<<" for input\n";
  200.       fm_error=1;
  201.             return 0.0;
  202.     }
  203.   if(plane>0)
  204.     {
  205.       if(fseek(fp_img,(long)plane*(fm_dims[0]/8)*fm_dims[1]*fm_dims[2],0)<0)
  206.  {
  207.    cerr<<"unable to seek to plane "<<plane<<" in file "<<filename<<"\n";
  208.    fm_error=1;
  209.    return 0.0;
  210.    }
  211.       return scale;
  212.     }
  213. }
  214.       /* Since image file is already open, we should be ready to read next plane*/
  216.       return scale;
  217.     }
  218. }


 
et le debugger me dit :
 
Program received signal SIGSEGV, Segmentation fault.
0x00000000 in ?? ()
 
... :pfff:


Message édité par nifnef le 28-07-2005 à 16:46:19
Reply

Marsh Posté le 26-07-2005 à 18:45:54    

en gros, je cherche à lire une image: dans le fichier .hdr on trouve les infos permettant de lire le .img qui contient les données brutes.

Reply

Marsh Posté le 27-07-2005 à 10:09:18    

nifnef a écrit :

Bonjour,
 
je suis sur cette erreur depuis un certain temps et j'ai peur d'avoir perdu tout recul pour la résoudre.
 

Code :
  1. #include "fm_volume.h"
  3. using namespace std;
  5. int main ()
  6. {
  7.   FmVolume* firstVolume = new FmVolume("1exam10_06_02" );
  9.   cout<<"yop"<<endl;
  10.   return 1;
  11. }


 
la classe FmVolume étant définie dans fm_volume.h.
Mon probleme c'est à l'execution : superbe segmentation fault apres "yop" ... :pt1cable:  
Se peut-il qu'il s'agisse d'erreurs de la classe FmVolume qui seraient "descendues" ?
Merci pour vos suggestions...


 
 
manquerait pas un

Code :
  1. delete firstvolume

des fois?

Reply

Marsh Posté le 27-07-2005 à 10:09:18   

Reply

Marsh Posté le 27-07-2005 à 12:54:40    

soit je met mon delete et ca plante dans le destructeur :  
 
Program received signal SIGSEGV, Segmentation fault.
gsl_matrix_free (m=0x0) at init_source.c:185
185       if (m->owner)
 
soit je ne le met pas :
 
Program received signal SIGSEGV, Segmentation fault.
0x40316f45 in _IO_cleanup () from /lib/libc.so.6
 
et je ne suis capable de résoudre aucun de ces pb

Reply

Marsh Posté le 27-07-2005 à 13:12:41    

ben t'as ta réponse. tu fais un gsm_matrix_free(0) ...

Reply

Marsh Posté le 27-07-2005 à 15:45:52    

Effectivement, j'essayais de detruire fm_mat qui n'etait pas alloué.
Mais ca a l'air de pas être le seul pb!
Bon, je crois que c'est pas clair pour moi les destructeurs ... j'ai lu a droite a gauche que je dois detruire les pointeurs?  
Donc dans mon cas, je détruis mon tableau de matrices (qui a été alloué et initialisé) et c'est tout?  
(ai-je raison de vouloir faire un free de ma matrice temporaire temp?)
 
si je fais delete[] fm_data apres, j'ai ca :
 
Program received signal SIGSEGV, Segmentation fault.
0x4031aaa4 in _int_free () from /lib/libc.so.6
 
ca m'enerve, j'ai l'impression que qq chose m'echappe

Reply

Marsh Posté le 27-07-2005 à 16:35:50    

je comprend pas trop les problèmes des gens avec l'allocation dynamique, y en a plein qui galère avec ca et pourtant c'est logique, tu alloue tu libère ca doit être automatique normalement (en tout cas pour moi) y a rien de plus à savoir.


---------------
http://www.blastmanu.info
Reply

Marsh Posté le 27-07-2005 à 16:58:22    

je comprendrai si ca marchait (non, je blague c'est pas garanti)... mais dans mon cas, tu peux me dire comment je dois libérer proprement?
j'ai essayé delete [] fm_data, j'ai aussi désalloué une part une les matrices du tableau (avec gsl_matrix_free) et ca marche pas.  
j'ai rajouté ca mais toujours pas d'amelioration ...

Code :
  1. delete fp_img;
  2.   delete fp_hdr;
  3.   delete [] flag;
  4.   delete [] data;


 
Program received signal SIGSEGV, Segmentation fault.
~basic_string (this=0xffffffff) at atomicity.h:40
40        __asm__ __volatile__ ("lock; xaddl %0,%2"
Current language:  auto; currently c++


Message édité par nifnef le 27-07-2005 à 17:20:32
Reply

Sujets relatifs:

Leave a Replay

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