pointeur de pointeur

Version imprimable

Voir 40 message(s) de cette discussion en une page

05/05/2006, 19h55
Emmanuel Delahaye

Citation:

Envoyé par CSoldier

Je pense que tout dépend du role de ta fonction, si elle est succeptible de modifier des variables membres de la structure oui, l'argument de ta fonction devra être un pointeur sur une structure du meme type et pour l'appel tu passes donc en argument l'adresse de ta structure sinon si c'est juste pour une lecture seule, pas besoin de pointeur sur la structure.

Pour de très pretites structures (d'une taille de 1 ou 2 int), c'est acceptable (mais mon compilateur warn...).

Au-dessus, c'est non. (copie de données inutiles). Si c'est pour une lecture seule, on sait comment définir le paramètre :

Code:

void f (T const *p);
05/05/2006, 20h00
crocodilex

Citation:

Envoyé par Emmanuel Delahaye

Ca n'existe pas.

http://emmanuel-delahaye.developpez.....htm#param_tab

Merci, (j'ai encore appris quelque chose aujourd'hui) je dormirais mieux cette nuit.
05/05/2006, 20h00
Franck.H

Ok, c'est noté, c'est toujours bon a prendre comme conseil ;)
07/05/2006, 15h17
petdelascar

Merci monsieur super VIP, je pensais bien à un *tab[] mais je voyais pas trop comment faire.

J'ai tenté:

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
float a1=1, a2=2, a3=3, a4=4, a5=5, a6=6, a7=7, a8=8, a9=9;
 
float * test_a[9];
 
test_a [0] = &a1;
test_a [1] = &a2;
test_a [2] = &a3;
test_a [3] = &a4;
test_a [4] = &a5;
test_a [5] = &a6;
test_a [6] = &a7;
test_a [7] = &a8;
test_a [8] = &a9;
 
elmhes( test_a, sizeof  test_a / sizeof * test_a);

Mais le programme se bloque, je pense que c'est lié au fait que test_a est un tableau à 2 dimensions, si quelqu'un peut m'aider. Je redonne la fonction dont je cherche comment écrire le parramètre "a" <=> test_a:

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
 
void elmhes(float **a, int n)
{
int m,j,i;
float y,x;
 
for (m=2;m<n;m++) {
    x=0.0;
    i=m;
    for (j=m;j<=n;j++) {
        if (fabs(a[j][m-1]) > fabs(x)) {
        x=a[j][m-1];
        i=j;
        //printf("x=%f\n", x);
        }
    }
...
}

J' pas trop pigé les contrainte de ton problème
mais si tu veux passé un tableau à deux dimentions dont tu ne connais pas la taille à l'avance à une fonction tu peut utiliser la méthode suivante,
je te mes un exemple simple :
Code:

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 39 40 41 42 43 44 45 46 47 48 49 50 51 52 #include <stdlib.h> #include <stdio.h> void SetTab(int *tab, size_t x, size_t y) { size_t i,j; puts("- Remplissage du tableau."); for(j=0; j<y; j++) { for(i=0; i<x; i++) tab[i]=j*10+i; tab += x; } } void PrintTab(int tab[], size_t x, size_t y) { size_t i,j; puts("- Affichage du tableau :"); for(j=0; j<y; j++) { for(i=0; i<x; i++) printf("%d\n",tab[i]); puts("--------"); tab += x; } } void ClearTab(int *tab, size_t x, size_t y) { size_t i; x*=y; for(i=0; i<x; i++) tab[i]=0; puts("- Toutes les cellules ont ete mise a zero."); } int main() { int Tableau[3][3]; SetTab((int*)Tableau, 3, 3); PrintTab((int*)Tableau, 3, 3); ClearTab((int*)Tableau, 3, 3); PrintTab((int*)Tableau, 3, 3); #ifdef _WIN32 system("PAUSE"); /* Pour la console Windows */ #endif return 0; }
Le code à l'interieur des fonction le traite comme plusieurs tableau à une dimension.
et même dans la fonction ClearTab comme un unique tableau.

Citation:

J' pas trop pigé les contrainte de ton problème
mais si tu veux passé un tableau à deux dimentions dont tu ne connais pas la taille à l'avance à une fonction tu peut utiliser la méthode suivante,
je te mes un exemple simple :

Merci, mais mon problème est différent, j'ai une fonction d'une bibliothèque qui prend une matrice a en entrée (paramètré par **a), mais je ne comprend pas comment déffinir **a, il semble qu'il s'agisse d'un tableau 2 dimension, voici son code:
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 void elmhes(float **a, int n) { int m,j,i; float y,x; for (m=2;m<n;m++) { x=0.0; i=m; for (j=m;j<=n;j++) { if (fabs(a[j][m-1]) > fabs(x)) { x=a[j][m-1]; i=j; ... } } }

07/05/2006, 18h33
Emmanuel Delahaye

Citation:

Envoyé par petdelascar

Merci, mais mon problème est différent, j'ai une fonction d'une bibliothèque qui prend une matrice a en entrée (paramètré par **a), mais je ne comprend pas comment déffinir **a, il semble qu'il s'agisse d'un tableau 2 dimension, voici son code:

Un T** peut recevoir l'adresse d'un tableau de pointeurs sur T. Chaque pointeur du tableau peut lui même pointer sur un tableau d'élements de type T. Ca se comporte comme un tableau à deux dimensions, mais ce n'est pas linéaire. Peu importe.

Pour fabriquer un tel objet, voir dans la FAQ (tableau dynamique à 2 dimensions)

J'ai suivi vos directives, voici ce que j'ai fait:

Code:

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
 
/* METHODE 1 */
float **test_a;
int taille = 3, taille2 = 3;
 
/* Allocation de la 1er dimension */
test_a = malloc ( sizeof(*test_a)  *  taille);
/* Allocation des tableaux */
for (i=0; i<taille; i++)
{
   test_a[i] = malloc ( sizeof(**test_a) * taille2);
}
 
test_a [0][0] = 1;
test_a [0][1] = 2;
test_a [0][2] = 3;
test_a [1][0] = 4;
test_a [1][1] = 5;
test_a [1][2] = 6;
test_a [2][0] = 7;
test_a [2][1] = 8;
test_a [2][2] = 9;
 
 
elmhes(test_a, 3);

Mais ça bloque à l'execution dans
elmhes(test_a, 3);

voici pourtant elmhes qui fonctionne puisqu'extraite de Numerical recipes:

Code:

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
39
40
41
42
 
void elmhes(float **a, int n)
{
/*
Reduction to Hessenberg form by the elimination method. The real, nonsymmetric matrix
a[1..n][1..n] is replaced by an upper Hessenberg matrix with identical eigenvalues. Recommended,
but not required, is that this routine be preceded by balanc. On output, the
Hessenberg matrix is in elements a[i][j] with i = j+1. Elements with i > j+1 are to be
thought of as zero, but are returned with random values.
*/
 
int m,j,i;
float y,x;
 
for (m=2;m<n;m++) {
    x=0.0;
    i=m;
    for (j=m;j<=n;j++) {
        if (fabs(a[j][m-1]) > fabs(x)) {
        x=a[j][m-1];
        i=j;
        //printf("x=%f\n", x);
        }
    }
    if (i != m) {
        for (j=m-1;j<=n;j++) SWAP(a[i][j],a[m][j])
        for (j=1;j<=n;j++) SWAP(a[j][i],a[j][m])
    }
    if (x) {
        for (i=m+1;i<=n;i++) {
            if ((y=a[i][m-1]) != 0.0) {
                y /= x;
                a[i][m-1]=y;
                for (j=m;j<=n;j++)
                    a[i][j] -= y*a[m][j];
                for (j=1;j<=n;j++)
                    a[j][m] += y*a[j][i];
            }
        }
    }
}
}

merci, je pense être roche de la fin.

Citation:

Envoyé par petdelascar

J'ai suivi vos directives, voici ce que j'ai fait:

Merci de poster la version qui compile...

En attendant, j'ai fait des suppositions et j'ai instrumenté le code. Il y a au moins une tentative de débordement du tableau.

Code:

1
2
3
4
5
6
7
 
Assertion failed: j < n, file main.c, line 36
 
This application has requested the Runtime to terminate it in an unusual way.
Please contact the application's support team for more information.
 
Press ENTER to continue.

Ce code est incorrect.

Code:

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
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
 
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <assert.h>
 
#define SWAP(a, b)\
do\
{\
   float tmp = a;\
   a = b;\
   b = a;\
}\
while (0)
 
void elmhes(float **a, int n)
{
   /*
   Reduction to Hessenberg form by the elimination method. The real, nonsymmetric matrix
   a[1..n][1..n] is replaced by an upper Hessenberg matrix with identical eigenvalues. Recommended,
   but not required, is that this routine be preceded by balanc. On output, the
   Hessenberg matrix is in elements a[i][j] with i = j+1. Elements with i > j+1 are to be
   thought of as zero, but are returned with random values.
   */
 
   int m, j, i;
   float y, x;
 
   for (m = 2;m < n;m++)
   {
      x = 0.0;
      i = m;
      for (j = m;j <= n;j++)
      {
 
         assert(j >= 0);
         assert(j < n);
         assert(m - 1 >= 0);
         assert(m - 1 < n);
         if (fabs(a[j][m - 1]) > fabs(x))
         {
            x = a[j][m - 1];
            i = j;
            //printf("x=%f\n", x);
         }
      }
      if (i != m)
      {
         for (j = m - 1;j <= n;j++)
         {
            assert(i >= 0);
            assert(i < n);
            assert(j >= 0);
            assert(j < n);
            assert(m >= 0);
            assert(m < n);
            SWAP(a[i][j], a[m][j]);
         }
         for (j = 1;j <= n;j++)
         {
            assert(i >= 0);
            assert(i < n);
            assert(j >= 0);
            assert(j < n);
            assert(m >= 0);
            assert(m < n);
            SWAP(a[j][i], a[j][m]);
         }
      }
      if (x)
      {
         for (i = m + 1;i <= n;i++)
         {
            assert(i >= 0);
            assert(i < n);
            assert(m - 1 >= 0);
            assert(m - 1 < n);
            if ((y = a[i][m - 1]) != 0.0)
            {
               y /= x;
               assert(i >= 0);
               assert(i < n);
               assert(m - 1 >= 0);
               assert(m - 1 < n);
               a[i][m - 1] = y;
               for (j = m;j <= n;j++)
               {
                  assert(i >= 0);
                  assert(i < n);
                  assert(j >= 0);
                  assert(j < n);
                  assert(m >= 0);
                  assert(m < n);
                  a[i][j] -= y * a[m][j];
               }
               for (j = 1;j <= n;j++)
               {
                  assert(i >= 0);
                  assert(i < n);
                  assert(j >= 0);
                  assert(j < n);
                  assert(m >= 0);
                  assert(m < n);
                  a[j][m] += y * a[j][i];
               }
            }
         }
      }
   }
}
 
int main(void)
{
   float **test_a;
   int taille = 3, taille2 = 3;
   int i;
 
   /* Allocation de la 1er dimension */
   test_a = malloc ( sizeof(*test_a) * taille);
   /* Allocation des tableaux */
   for (i = 0; i < taille; i++)
   {
      test_a[i] = malloc ( sizeof(**test_a) * taille2);
   }
 
   test_a [0][0] = 1;
   test_a [0][1] = 2;
   test_a [0][2] = 3;
   test_a [1][0] = 4;
   test_a [1][1] = 5;
   test_a [1][2] = 6;
   test_a [2][0] = 7;
   test_a [2][1] = 8;
   test_a [2][2] = 9;
 
#if 1
 
   elmhes(test_a, 3);
#endif
 
   return 0;
}

Voici le code compilable, merci encore, grandement:

Code:

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
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
 
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
 
void elmhes(float **a, int n);
 
#define SWAP(g,h) {y=(g);(g)=(h);(h)=y;}
 
main ()
{
 
/* METHODE 1 */
float **test_a;
int taille = 3, taille2 = 3, i, j;
 
/* Allocation de la 1er dimension */
test_a = malloc ( sizeof(*test_a)  *  taille);
/* Allocation des tableaux */
for (i=0; i<taille; i++)
{
   test_a[i] = malloc ( sizeof(**test_a) * taille2);
}
 
test_a [0][0] = 1;
test_a [0][1] = 2;
test_a [0][2] = 3;
test_a [1][0] = 4;
test_a [1][1] = 5;
test_a [1][2] = 6;
test_a [2][0] = 7;
test_a [2][1] = 8;
test_a [2][2] = 9;
 
elmhes(test_a, 3);
 
/* Printing */
 
for (i = 0; i < 3; i++) {
    for (j = 0; j < 3; j++) {
        //printf ("xn(%d) = %f\n", i, xn[i]);
        //printf ("test_X(%d,%d) = %f\n", i, j, gsl_matrix_get (test_X, i, j));
        //sort = * test_a [i][i];
        printf ("test_a(%d,%d) = %f\n", i, j, test_a[i][j]);
    }
}
 
}
 
 
void elmhes(float **a, int n)
{
/*
Reduction to Hessenberg form by the elimination method. The real, nonsymmetric matrix
a[1..n][1..n] is replaced by an upper Hessenberg matrix with identical eigenvalues. Recommended,
but not required, is that this routine be preceded by balanc. On output, the
Hessenberg matrix is in elements a[i][j] with i = j+1. Elements with i > j+1 are to be
thought of as zero, but are returned with random values.
*/
 
int m,j,i;
float y,x;
 
for (m=2;m<n;m++) {
    x=0.0;
    i=m;
    for (j=m;j<=n;j++) {
        if (fabs(a[j][m-1]) > fabs(x)) {
        x=a[j][m-1];
        i=j;
        //printf("x=%f\n", x);
        }
    }
    if (i != m) {
        for (j=m-1;j<=n;j++) SWAP(a[i][j],a[m][j])
        for (j=1;j<=n;j++) SWAP(a[j][i],a[j][m])
    }
    if (x) {
        for (i=m+1;i<=n;i++) {
            if ((y=a[i][m-1]) != 0.0) {
                y /= x;
                a[i][m-1]=y;
                for (j=m;j<=n;j++)
                    a[i][j] -= y*a[m][j];
                for (j=1;j<=n;j++)
                    a[j][m] += y*a[j][i];
            }
        }
    }
}
}

En cachant l'appel de elmhes, on voit que l'erreur vient de là.

08/05/2006, 14h48
Emmanuel Delahaye

Citation:

Envoyé par petdelascar

En cachant l'appel de elmhes, on voit que l'erreur vient de là.

Ok, j'ai bien vu. Je t'ai démontré qu'il y avait des débordements de tableau. Le comportement est indéterminé.

Reprend mon code instrumenté (assert) et corrige jusqu'à ce qu'il aille jusqu'au bout. Ca devrait faire le ménage...
08/05/2006, 19h13
petdelascar

L'erreur venait bien d'un dépassement de tableau à cause d'une mauvaise implémentation de la fonction elmhes. n ne représente pas le nombre de lignes ou de colones de la matrice mais l'indice maximal du tableau correspondant, c'est à dire nombre de lignes - 1.

Merci beaucoup Emmanuel, et aux autres crocodiles.

;)
08/05/2006, 19h33
Emmanuel Delahaye

Citation:

Envoyé par petdelascar

L'erreur venait bien d'un dépassement de tableau à cause d'une mauvaise implémentation de la fonction elmhes. n ne représente pas le nombre de lignes ou de colones de la matrice mais l'indice maximal du tableau correspondant, c'est à dire nombre de lignes - 1.

J'espère que tu as compris le coté 'instrumentation' avec assert(). C'est un outil de mise au point très précieux.
08/05/2006, 22h10
petdelascar

Citation:

'espère que tu as compris le coté 'instrumentation' avec assert(). C'est un outil de mise au point très précieux.

Oui, je ne l'avais jamais vu celui là, c'est génial les forums on apprends pleins de trucs, c'est également un beau témoignage de l'altruisme.
08/05/2006, 22h19
Emmanuel Delahaye

Citation:

Envoyé par petdelascar

Oui, je ne l'avais jamais vu celui là, c'est génial les forums on apprends pleins de trucs, c'est également un beau témoignage de l'altruisme.

Arrête, je vais pleurer

"Tout le monde il est beau, tout le monde il est gentil..."

<mode cynique = ON>
Mon but est que tu écrives du code solide, pour m'éviter d'avoir un jour à le débugger, c'est tout ! :mouarf::mouarf::mouarf:
<mode cynique = OFF>
09/05/2006, 11h40
petdelascar

Citation:

Mon but est que tu écrives du code solide, pour m'éviter d'avoir un jour à le débugger, c'est tout !

C'est peine perdue, il y a trop de programmeurs sur terre, alors c'est de l'altruisme, aller, on y croit!

Voir 40 message(s) de cette discussion en une page