problem d'affectation illogique!

Version imprimable

bonjour!
s.v.p ,quelqu'un peu m’expliquer pourquoi le code suivant ne fonction pas !
(il s'agit d'un programme de simplexe )

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
for(int j=1;j<=varecart+1;j++)
for(int i=1;i<=(varbase+varecart+1);i++) 
      {
      if(j==ve){ mat[i][j]=mat[i][j]/mat[vs][ve];
 
        //  printf("%2.2f \t",mat[i][j]/mat[vs][ve]);
 
                  }
       else
                  {
                   mat[i][j]=mat[i][j]-mat[vs][j]*mat[i][ve]/mat[vs][ve];
 
       //printf("%2.2f \t",(mat[i][j]-mat[vs][j]*mat[i][ve]/mat[vs][ve])); 
                   }
       }

a vrai dire dans les printf les opérations fonctionnent normale mais l'affectation ne marche ps !!
voici le code de mon programme !!

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
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
 
 
 
float mat[15][15];
int varbase;
int varecart;
 
////////////////////////////////////////////////////////////////////////////////
void init()
{
  for(int i=0;i<=15;i++)
  for(int j=0;j<=15;j++)
  mat[i][j]=0;
}
////////////////////////////////////////////////////////////////////////////////
 
int object()
{
printf(" La Fonction Objectif : \n");
 
printf("saisie le nbr de varaibles de base:\n");
scanf("%d",&varbase);
 
for(int i=1;i<=varbase;i++)
{
printf("X%d=",i);
scanf("%f",&mat[i][1]);
}
return 0;   
}
 
////////////////////////////////////////////////////////////////////////////////
 
int contrainte()
{
  printf("saisie le nombnre de contraintes :\n");
  scanf("%d",&varecart);
  for(int j=2;j<=varecart+1;j++)
  { 
   printf("Saisie les variables de contrainte num %d :\n",j-1);
   for(int i=1;i<=varbase;i++)
    {
     printf("X%d=",i); 
     scanf("%f",&mat[i][j]);     
    }
  }
 
  for(int j=2;j<=varecart+1;j++) 
  {
  printf("solution de contrainte num %d:\n",j-1);
  scanf("%f",&mat[varbase+varecart+1][j]);
  }  
 
}
////////////////////////////////////////////////////////////////////////////////
int iden()
{
    //////////fct objectif
    for(int i=varbase+1;i<=varbase+varecart;i++)
    mat[i][1]=0;
    //mat[varbase+varecart+1][1]=1;
    ///////// matrice iden
    for(int j=2;j<=varecart+1;j++)
    for(int i=varbase+1;i<=varbase+varecart;i++)
    if(i-varbase==j-1)mat[i][j]=1;
    else mat[i][j]=0;
 return 0;   
}
////////////////////////////////////////////////////////////////////////////////
 
int affichsys()
{
 
    printf("le systeme :\n\n");
    printf("la fonction Objectif :\n MaxZ=\t");
    for(int i=1;i<=varbase;i++)
    printf("+%.fX%d",mat[i][1],i);
    printf("\n les contraintes :\n");
    for(int j=2;j<=varecart+1;j++)
    { 
      printf("\nla contrainte num %d\n",j-1);
      for(int i=1;i<=varbase;i++)
      printf("+%.fx%d ",mat[i][j],i);
 
    }       
    printf("\n");
    system("PAUSE");
 
    ///////////////////////////////////////  affichage en mode tableau
     printf("\n\n le tableau \n\n"); 
    for(int j=1;j<=varecart+1;j++)
    {
     printf("\n");
     for(int i=1;i<=varbase+varecart+1;i++)
     printf("%2.2f\t",mat[i][j]);
    } 
return 0;    
}
////////////////////////////////////////////////////////////////////////////////
 
int max_o()
{   
    int vs=0; 
    float max=0;
    for(int i=1;i<=varbase+varecart;i++)
    if (mat[i][1]>max)
        {
         vs=i;
         max=mat[i][1];
         }
 
return vs;
 
}
 
////////////////////////////////////////////////////////////////////////////////
 
int pivot()
{ 
 
  int vs=max_o(); 
 // printf("le var sortante x%d \n",vs); 
  float def_pivot=10000;
  int ve=0;
  for(int j=2;j<=varecart+1;j++)
  {
   if(mat[varbase+varecart+1][j]!=0&&mat[vs][j]!=0&&((mat[varbase+varecart+1][j]/mat[vs][j])<=def_pivot))
    {  
     def_pivot=mat[varbase+varecart+1][j]/mat[vs][j];
     ve=j; 
    }
  }     
    // printf("les pivot est %.f \n",mat[vs][ve]);  
 return ve;   
}
////////////////////////////////////////////////////////////////////////////////
 
int pivotage()
{
 
  int ve=pivot(); //la ligne pivot
  printf("la ligne pivot %d \n",ve);
  int vs=max_o();//la colomn pivot
  printf("la colone pivot%d\n",vs);
  printf("le pivot %2.2f \n",mat[vs][ve]);
  system("PAUSE");
 
  for(int j=1;j<=varecart+1;j++)
   { printf("\n");
    for(int i=1;i<=(varbase+varecart+1);i++) 
      {if(j==ve){ mat[i][j]=mat[i][j]/mat[vs][ve];
                    printf("%2.2f",mat[i][j]);        
        //  printf("%2.2f \t",mat[i][j]/mat[vs][ve]);
 
                 }
       else
                 {
       mat[i][j]=(mat[i][j]-mat[vs][j]*mat[i][ve]/mat[vs][ve]);
      printf("%2.2f",mat[i][j]);
       //printf("%2.2f \t",(mat[i][j]-mat[vs][j]*mat[i][ve]/mat[vs][ve]));
 
 
                 }
       }   
 }
     system("PAUSE");
 return 0;
}
 
////////////////////////////////////////////////////////////////////////////////
int main()
{
       object();
       contrainte();
       iden();
       affichsys();
       pivotage();
       //affichsys();
     // printf("le var sortante x%d \n",max_o()); 
     // printf("les pivot est %.f \n",mat[max_o()][pivot()]);       
 
 
 
system("PAUSE");
return 0;   
}

30/12/2011, 11h42
therwald
Pourrais-tu préciser en quoi cela ne fonctionne pas?
- segfault
- valeurs aberrantes
- ??
Les boucles de lecture/écriture de tableau avec des bornes calculées sont des sources fréquentes d'erreur par débordement...
A côté de ça les variables globales sont dangereuses car elles peuvent générer des interactions difficiles à suivre...mieux vaut utiliser des références transmises à l'appel des fonctions.
31/12/2011, 14h38
Neckara

Les tableaux commencent à 0 et non à 1

La première case d'un float test[10][10] est test[0][0].

La dernière colonne ou la dernière ligne est 9 : la ligne et la colonne 10 n'existe pas.

As-tu penser à vérifier que tu regardes bien des cases existantes?
Ex : test[14][2] n'existe pas et provoquera des segfault.

Sinon plutôt que d'utiliser un float test[10][10], déclare plutôt un (float *test) et fait un
test = (float *)malloc(sizeof(float)*ligne*colonne);

tu accèdera à chacun des élément par : *(test + colonne + ligne*nb_colonne)
Et n'oublie pas le free(test);
31/12/2011, 15h49
Joel F

Citation:

Envoyé par Neckara

Sinon plutôt que d'utiliser un float test[10][10], déclare plutôt un (float *test) et fait un
test = (float *)malloc(sizeof(float)*ligne*colonne);

tu accèdera à chacun des élément par : *(test + colonne + ligne*nb_colonne)
Et n'oublie pas le free(test);

Ca c'ets plutot du C, en C++, on preferera std::vector
31/12/2011, 18h32
Neckara

Les vecteurs sont des liste chaînées non?
Sinon, il y a les new/delete qui sont "équivalentes" à malloc/free
31/12/2011, 22h50
Joel F

Citation:

Envoyé par Neckara

Les vecteurs sont des liste chaînées non?

Non, c'est un bloc de memoire contigu; std::list est une liste doublement chainee. Ne jamais manipulé votre memoire manuellement, c'est le meilleur moyen de se vautrer.

Bonsoir,

Citation:

Envoyé par mdj-sofiane

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
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
#include<stdio.h>
#include<stdlib.h>
#include<string.h>



float mat[15][15];
int varbase;
int varecart;

////////////////////////////////////////////////////////////////////////////////
void init()
{
  for(int i=0;i<=15;i++)
  for(int j=0;j<=15;j++)
  mat[i][j]=0;
}
////////////////////////////////////////////////////////////////////////////////

int object()
{
printf(" La Fonction Objectif : \n");

printf("saisie le nbr de varaibles de base:\n");
scanf("%d",&varbase);

for(int i=1;i<=varbase;i++)
{
printf("X%d=",i);
scanf("%f",&mat[i][1]);
}
return 0;   
}

////////////////////////////////////////////////////////////////////////////////

int contrainte()
{
  printf("saisie le nombnre de contraintes :\n");
  scanf("%d",&varecart);
  for(int j=2;j<=varecart+1;j++)
  { 
   printf("Saisie les variables de contrainte num %d :\n",j-1);
   for(int i=1;i<=varbase;i++)
    {
     printf("X%d=",i); 
     scanf("%f",&mat[i][j]);     
    }
  }

  for(int j=2;j<=varecart+1;j++) 
  {
  printf("solution de contrainte num %d:\n",j-1);
  scanf("%f",&mat[varbase+varecart+1][j]);
  }  

}
////////////////////////////////////////////////////////////////////////////////
int iden()
{
    //////////fct objectif
    for(int i=varbase+1;i<=varbase+varecart;i++)
    mat[i][1]=0;
    //mat[varbase+varecart+1][1]=1;
    ///////// matrice iden
    for(int j=2;j<=varecart+1;j++)
    for(int i=varbase+1;i<=varbase+varecart;i++)
    if(i-varbase==j-1)mat[i][j]=1;
    else mat[i][j]=0;
 return 0;   
}
////////////////////////////////////////////////////////////////////////////////

int affichsys()
{
    
    printf("le systeme :\n\n");
    printf("la fonction Objectif :\n MaxZ=\t");
    for(int i=1;i<=varbase;i++)
    printf("+%.fX%d",mat[i][1],i);
    printf("\n les contraintes :\n");
    for(int j=2;j<=varecart+1;j++)
    { 
      printf("\nla contrainte num %d\n",j-1);
      for(int i=1;i<=varbase;i++)
      printf("+%.fx%d ",mat[i][j],i);
    
    }       
    printf("\n");
    system("PAUSE");

    ///////////////////////////////////////  affichage en mode tableau
     printf("\n\n le tableau \n\n"); 
    for(int j=1;j<=varecart+1;j++)
    {
     printf("\n");
     for(int i=1;i<=varbase+varecart+1;i++)
     printf("%2.2f\t",mat[i][j]);
    } 
return 0;    
}
////////////////////////////////////////////////////////////////////////////////

int max_o()
{   
    int vs=0; 
    float max=0;
    for(int i=1;i<=varbase+varecart;i++)
    if (mat[i][1]>max)
        {
         vs=i;
         max=mat[i][1];
         }
   
return vs;

}

////////////////////////////////////////////////////////////////////////////////

int pivot()
{ 
  
  int vs=max_o(); 
 // printf("le var sortante x%d \n",vs); 
  float def_pivot=10000;
  int ve=0;
  for(int j=2;j<=varecart+1;j++)
  {
   if(mat[varbase+varecart+1][j]!=0&&mat[vs][j]!=0&&((mat[varbase+varecart+1][j]/mat[vs][j])<=def_pivot))
    {  
     def_pivot=mat[varbase+varecart+1][j]/mat[vs][j];
     ve=j; 
    }
  }     
    // printf("les pivot est %.f \n",mat[vs][ve]);  
 return ve;   
}
////////////////////////////////////////////////////////////////////////////////
 
int pivotage()
{
 
  int ve=pivot(); //la ligne pivot
  printf("la ligne pivot %d \n",ve);
  int vs=max_o();//la colomn pivot
  printf("la colone pivot%d\n",vs);
  printf("le pivot %2.2f \n",mat[vs][ve]);
  system("PAUSE");
  
  for(int j=1;j<=varecart+1;j++)
   { printf("\n");
    for(int i=1;i<=(varbase+varecart+1);i++) 
      {if(j==ve){ mat[i][j]=mat[i][j]/mat[vs][ve];
                    printf("%2.2f",mat[i][j]);        
        //  printf("%2.2f \t",mat[i][j]/mat[vs][ve]);
                     
                 }
       else
                 {
       mat[i][j]=(mat[i][j]-mat[vs][j]*mat[i][ve]/mat[vs][ve]);
      printf("%2.2f",mat[i][j]);
       //printf("%2.2f \t",(mat[i][j]-mat[vs][j]*mat[i][ve]/mat[vs][ve]));
  
          
                 }
       }   
 }
     system("PAUSE");
 return 0;
}

////////////////////////////////////////////////////////////////////////////////
int main()
{
       object();
       contrainte();
       iden();
       affichsys();
       pivotage();
       //affichsys();
     // printf("le var sortante x%d \n",max_o()); 
     // printf("les pivot est %.f \n",mat[max_o()][pivot()]);       
      
      
      
system("PAUSE");
return 0;   
}

dans un tableau de 15 cases, les cases sont indexées de 0 à 14 inclus.