seg fault avec gmp

**kkouet'** · 26/05/2010, 15h17

alors voilà je suis nouvelle sur ce forum et j'ai une segfault que j'arrive pas à supprimer.
je m'explique, je fais un projet pour lequel je dois manipuler de très grands nombres (ex : 230979485048261039076034528529086839797) j'utilise donc la librairie gmp.
mon erreur est la suivante : j'ai un très grand entier (type mpz_t) que je dois convertir en binaire (dans mon cas, créer un tableau de 128 bits), ma fonction de conversion fonctionne parfaitement mais à la fin j'ai une segfault qui apparait alors que le résultat est correct (pour moi elle est là juste pour le fun

). j'ai donc utilisé gdb pour comprendre un peu mieux l'erreur et voilà le résultat :

Program received signal SIGSEGV, Segmentation fault.
0x001befa8 in ?? () from /lib/tls/i686/cmov/libc.so.6
(gdb) up
#1 0x001bf0cf in exit () from /lib/tls/i686/cmov/libc.so.6
(gdb) up
#2 0x001a6b5e in __libc_start_main () from /lib/tls/i686/cmov/libc.so.6
(gdb) up
#3 0x080489a1 in _start () at ../sysdeps/i386/elf/start.S:119
119 ../sysdeps/i386/elf/start.S: Aucun fichier ou dossier de ce type.
in ../sysdeps/i386/elf/start.S
Current language: auto
The current source language is "auto; currently asm".

j'avoue que j'ai jamais compris grand chose à gdb, je l'utilise pour situer la ligne de code fautive mais là il m'aide pas du tout !

voilà mon programme qui fait la conversion :

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
 
void gmp2bin(mpz_t n, int size, int* bin_tab){
        int i=size-1;
        mpz_t nb;
        mpz_init_set(nb,n);
 
        if(mpz_cmp_ui(nb,0)==0){
                while(i>=0){
                        bin_tab[i] = 0;
                        i--;
                }
        }
        else{
                while(mpz_cmp_ui(nb,0)>0){
                       if(mpz_even_p(nb) == 0){
                                bin_tab[i] = 1;
                                mpz_sub_ui(nb,nb,1);
                                mpz_divexact_ui(nb,nb,2);
                        }
                        else{
                                bin_tab[i] = 0;
                                mpz_divexact_ui(nb,nb,2);
                        }
                        i--;
                }
                while(i>=0){
                        bin_tab[i] = 0;
                        i--;
                }
        }
}

voilà j'espère que j'ai été assez claire et que quelqu'un pourra me venir en aide !

merci d'avance !

**Grulim** · 26/05/2010, 16h26

Peux-tu préciser l'OS (linux?), le compilo (gcc?), sa version, les options de compilation et aussi la version de gdb ?

**kkouet'** · 26/05/2010, 16h33

alors je suis sous ubuntu avec gcc version 4.4.1, avec comme options : -pg -g -Wall -W -pedantic -O0 et la version de gdb est la 7.0-ubuntu

(merci j'avais pas pensé à préciser tout ça)

**Grulim** · 26/05/2010, 16h42

un SEG FAULT peut résulter (entre autre) d'un dépassement de borne dans un tableau.
Comment est appelée la fonction gmp2bin et plus précisément, comment sont initialisés les paramètres size et bin_tab ?

**kkouet'** · 26/05/2010, 17h02

alors voilà c'était vraiment juste pour tester donc c'est pas bien compliqué.
comme ma fonction marchait malgré ma segfault j'ai continué en l'utilisant dans mon projet et là plus de segfault donc je pense que tu as raison, c'est au moment de la tester que je me plante.

même si ça marche je veux bien que tu me dise ou j'ai fait l'erreur parce que je pense que je pourrais la refaire ailleurs...

en tout cas merci de ton aide !

Code :

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

1
2
3
4
5
6
7
8
9
 
mpz_t res;
mpz_set_ui(res,17);
int* tab = malloc(20*sizeof(int));
 
gmp2bin(res,20,tab);
 
for(j=0;j<20;j++) printf("%i",tab[j]);
printf("\n");

**Grulim** · 26/05/2010, 17h19

Je ne connais pas la librairie gmp donc je ne me prononcerais pas sur son utilisation (correcte ou non)

.
Pour le reste, je vois une seule chose :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
                while(mpz_cmp_ui(nb,0)>0){
                       if(mpz_even_p(nb) == 0){
                                bin_tab[i] = 1;
                                mpz_sub_ui(nb,nb,1);
                                mpz_divexact_ui(nb,nb,2);
                        }
                        else{
                                bin_tab[i] = 0;
                                mpz_divexact_ui(nb,nb,2);
                        }
                        i--;
                }

Es-tu sure de sortir de cette boucle à temps ? cad avant que i < 0...

**Grulim** · 26/05/2010, 17h24

Si tu me permets de te donner quelques conseils amicaux :
1. Dans tes boucle while, je préfèrerais que tu utilises la pré-incrémentation/décrémentation plutôt que la post.
2. ton code peut être simplifié en utilisant un memset.
3. un tableau de int pour stocker des 0 et des 1, pourquoi pas un tableau de char ou d'unsigned char plutôt ?

**kkouet'** · 27/05/2010, 09h29

merci pour tes conseils, c'est toujours bon à prendre !
par contre je comprends pas pourquoi tu préfères que je mette le "i--", dans mon cas, au début, j'ai toujours eu l'habitude de le mettre à la fin et je trouve que c'est plus lisible. ensuite memset je connaissais pas (j'ai appris le C toute seule alors il y a plein de fonctions que je connais pas et j'ai souvent tendance à ré-inventer la roue

) et pour mon tableau je préfère les int, je trouve que c'est plus pratique à gérer. mais merci quand même pour tes suggestions !

**Grulim** · 27/05/2010, 14h46

Envoyé par kkouet'

merci pour tes conseils, c'est toujours bon à prendre !

A ton service

Envoyé par kkouet'

par contre je comprends pas pourquoi tu préfères que je mette le "i--", dans mon cas, au début, j'ai toujours eu l'habitude de le mettre à la fin et je trouve que c'est plus lisible

Dans ton cas, il n'y aura pas de différence, puisque tu utilises un type simple (int) mais pour des types plus compliqués, il se pourrait (note le conditionnel) que ton compilo ne sache pas optimiser ton code. C'est plus flagrant en C++.

Envoyé par kkouet'

ensuite memset je connaissais pas (j'ai appris le C toute seule alors il y a plein de fonctions que je connais pas et j'ai souvent tendance à ré-inventer la roue

)

Utilise les fonctions du Runtime C, tu verras ça améliorera la clarté de ton code (pour les autres) et tu arriveras au même résultat plus vite (si, si je t'assure

)

Envoyé par kkouet'

et pour mon tableau je préfère les int, je trouve que c'est plus pratique à gérer.

Ah bon, pourquoi ?

Envoyé par kkouet'

mais merci quand même pour tes suggestions !

Au plaisir.

Tu n'as pas répondu à ma question :

Es-tu sure de sortir de cette boucle à temps ? cad avant que i < 0...

**Grulim** · 27/05/2010, 14h48

voila le code un peu remanié:

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
36
37
38
void gmp2bin2(mpz_t n, int size, int* bin_tab){
    int i;
    mpz_t nb;
 
    memset(bin_tab, 0, size*sizeof(bin_tab[0]));
 
    if(mpz_cmp_ui(nb,0)!=0){
        mpz_init_set(nb,n);
        for(int i=size-1; i > 0 && mpz_cmp_ui(nb,0)>0; --i) {
           if(mpz_even_p(nb) == 0){
                bin_tab[i] = 1;
                mpz_sub_ui(nb,nb,1);
                mpz_divexact_ui(nb,nb,2);
            }
            else{
                mpz_divexact_ui(nb,nb,2);
            }
        }
    }
}
 
 
 
int main()
{
    mpz_t res;
    int* tab = malloc(20*sizeof(int));
    if(tab != NULL) {
        mpz_set_ui(res,17);
        gmp2bin(res,20,tab);
 
        for(j=0;j<20;j++) printf("%i",tab[j]);
        printf("\n");
 
        free(tab);
    }
    return 0;
}

**ironzorg** · 28/05/2010, 07h50

Code :

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

memset(bin_tab, 0, size*sizeof(bin_tab[0]));

Rien a voir avec le topic mais c'est un raccourci sympa a utiliser: pour recuperer la taille d'un tableau a une dimension, on peut connaitre sa taille grace a la macro si dessous (pour economiser un argument dans une fonction, en l'occurence `size').

Code :

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

1
2
 
#define T_SIZE(x) (sizeof(x) / sizeof(*x))

HTH.

gl · 28/05/2010, 10h19

Envoyé par Grulim

Dans ton cas, il n'y aura pas de différence, puisque tu utilises un type simple (int) mais pour des types plus compliqués, il se pourrait (note le conditionnel) que ton compilo ne sache pas optimiser ton code. C'est plus flagrant en C++.

C'est certes vrai en C++ (surcharge des opérateurs sur les classes utilisateur).

Mais en C, sur un compilateur un minimum efficace, cela fait longtemps que ce n'est plus vrai, puisque, pour faire simple, on ne peut appeler les opérateurs d'incrément et de décrément que sur des types simples.

Envoyé par ironzorg

Code :

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

memset(bin_tab, 0, size*sizeof(bin_tab[0]));

Rien a voir avec le topic mais c'est un raccourci sympa a utiliser: pour recuperer la taille d'un tableau a une dimension, on peut connaitre sa taille grace a la macro si dessous (pour economiser un argument dans une fonction, en l'occurence `size').

Code :

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

1
2
 
#define T_SIZE(x) (sizeof(x) / sizeof(*x))

C'est vrai pour un tableau statique. Dans le code de Grulim, bin_tab est un pointeur, ce n'est donc pas utilisable ici.

D'une manière générale, lorsque tu passes le tableau en paramètre d'une fonction, pour la fonction en question ce n'est qu'un pointeur et cette astuce n'est donc pas applicable, il est nécessaire de fournir la taille en paramètre de la fonction.

**Grulim** · 28/05/2010, 15h06

Envoyé par gl

Mais en C, sur un compilateur un minimum efficace, cela fait longtemps que ce n'est plus vrai, puisque, pour faire simple, on ne peut appeler les opérateurs d'incrément et de décrément que sur des types simples.

Pourrais-tu préciser ? Comme type non-simple, je pensais aux enums.
Est-ce les compilos efficaces qui signalent (interdisent?) cela, ou la norme (C99?) ?

Je travaille actuellement avec vieux compilo C (gcc 3.4.6 20060404), est-il concerné ?

NB: Ce sont de vrais questions, pas une vaine tentative de pinaillage

.