Discussion :

[Emmanuel Delahaye]

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 :

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

void f (T const *p);

[crocodilex]

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.

[Franck.H]

Ok, c'est noté, c'est toujours bon a prendre comme conseil

[petdelascar]

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

[petdelascar]

J'ai tenté:

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
 
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 :

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
 
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);
        }
    }
...
}

[CGi]

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 :

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
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.

[petdelascar]

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 :

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
 
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;
       ...
    }
}
}

[Emmanuel Delahaye]

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)

[petdelascar]

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

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
 
/* 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 :

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
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.

[Emmanuel Delahaye]

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 :

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

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 :

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
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;
}

[petdelascar]

Voici le code compilable, merci encore, grandement:

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
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à.

[Emmanuel Delahaye]

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...

[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.

[Emmanuel Delahaye]

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.

[petdelascar]

'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.

[Emmanuel Delahaye]

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 !
<mode cynique = OFF>

[petdelascar]

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!