Discussion :

[[ZiP]]

Bonjour,

J'ai l'erreur suivante : Opération en virgule flottante incorrecte

Sur la ligne :

Code :

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

1
2
 
DutyCycle := (TH / (TH + TB)) * 100;

Voici mon programme :

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
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
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
 
unit UPrincipale;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, Math;
 
type
  TForm1 = class(TForm)
    Memo1: TMemo;
    Button1: TButton;
    procedure Button1Click(Sender: TObject);
  private
    { Déclarations privées }
  public
    { Déclarations publiques }
  end;
 
var
  Form1: TForm1;
 
implementation
 
{$R *.dfm}
 
procedure TForm1.Button1Click(Sender: TObject);
var
     R         : array of double;
     C         : array of double;
     i,j,k     : Integer;
     T,TH,TB   : double;
     DutyCycle : double;
begin
     SetLength(R,145);
     R[0] := 1;
     R[1] := 1.2;
     R[2] := 1.5;
     R[3] := 1.8;
     R[4] := 2.2;
     R[5] := 2.7;
     R[6] := 3.3;
     R[7] := 3.9;
     R[8] := 4.7;
     R[9] := 5.6;
     R[10] := 6.8;
     R[11] := 8.2;
     R[12] := 10;
     R[13] := 11;
     R[14] := 12;
     R[15] := 13;
     R[16] := 15;
     R[17] := 16;
     R[18] := 18;
     R[19] := 20;
     R[20] := 22;
     R[21] := 24;
     R[22] := 27;
     R[23] := 30;
     R[24] := 33;
     R[25] := 36;
     R[26] := 39;
     R[27] := 43;
     R[28] := 47;
     R[29] := 51;
     R[30] := 56;
     R[31] := 62;
     R[32] := 68;
     R[33] := 75;
     R[34] := 82;
     R[35] := 91;
     R[36] := 100;
     R[37] := 110;
     R[38] := 120;
     R[39] := 130;
     R[40] := 150;
     R[41] := 160;
     R[42] := 180;
     R[43] := 200;
     R[44] := 220;
     R[45] := 240;
     R[46] := 270;
     R[47] := 300;
     R[48] := 330;
     R[49] := 360;
     R[50] := 390;
     R[51] := 430;
     R[52] := 470;
     R[53] := 510;
     R[54] := 560;
     R[55] := 620;
     R[56] := 680;
     R[57] := 750;
     R[58] := 820;
     R[59] := 910;
     R[60] := 1000;
     R[61] := 1100;
     R[62] := 1200;
     R[63] := 1300;
     R[64] := 1500;
     R[65] := 1600;
     R[66] := 1800;
     R[67] := 2000;
     R[68] := 2200;
     R[69] := 2400;
     R[70] := 2700;
     R[71] := 3000;
     R[72] := 3300;
     R[73] := 3600;
     R[74] := 3900;
     R[75] := 4300;
     R[76] := 4700;
     R[77] := 5100;
     R[78] := 5600;
     R[79] := 6200;
     R[80] := 6800;
     R[81] := 7500;
     R[82] := 8200;
     R[83] := 9100;
     R[84] := 10000;
     R[85] := 11000;
     R[86] := 12000;
     R[87] := 13000;
     R[88] := 15000;
     R[89] := 16000;
     R[90] := 18000;
     R[91] := 20000;
     R[92] := 22000;
     R[93] := 24000;
     R[94] := 27000;
     R[95] := 30000;
     R[96] := 33000;
     R[97] := 36000;
     R[98] := 39000;
     R[99] := 43000;
     R[100] := 47000;
     R[101] := 51000;
     R[102] := 56000;
     R[103] := 62000;
     R[104] := 68000;
     R[105] := 75000;
     R[106] := 82000;
     R[107] := 91000;
     R[108] := 100000;
     R[109] := 110000;
     R[110] := 120000;
     R[111] := 130000;
     R[112] := 150000;
     R[113] := 160000;
     R[114] := 180000;
     R[115] := 200000;
     R[116] := 220000;
     R[117] := 240000;
     R[118] := 270000;
     R[119] := 300000;
     R[120] := 330000;
     R[121] := 360000;
     R[122] := 390000;
     R[123] := 430000;
     R[124] := 470000;
     R[125] := 510000;
     R[126] := 560000;
     R[127] := 620000;
     R[128] := 680000;
     R[129] := 750000;
     R[130] := 820000;
     R[131] := 910000;
     R[132] := 1000000;
     R[133] := 1200000;
     R[134] := 1500000;
     R[135] := 1800000;
     R[136] := 2200000;
     R[137] := 2700000;
     R[138] := 3300000;
     R[139] := 3900000;
     R[140] := 4700000;
     R[141] := 5600000;
     R[142] := 6800000;
     R[143] := 8200000;
     R[144] := 10000000;
 
     SetLength(C,23);
     C[0] := 1;
     C[1] := 1.5;
     C[2] := 2.2;
     C[3] := 2.7;
     C[4] := 3;
     C[5] := 3.3;
     C[6] := 3.9;
     C[7] := 4.7;
     C[8] := 6.8;
     C[9] := 8.2;
     C[10] := 10;
     C[11] := 12;
     C[12] := 15;
     C[13] := 22;
     C[14] := 33;
     C[15] := 47;
     C[16] := 68;
     C[17] := 100;
     C[18] := 150;
     C[19] := 220;
     C[20] := 330;
     C[21] := 470;
     C[22] := 680;
 
     for i:=0 to Length(C) do
     begin
          for j:=1 to Length(R) do
          begin
               for k:=1 to Length(R) do
               begin
                    TH := 0.693 * (R[k] + R[j]) * C[i] * 0.000001;
                    TB := 0.693 * R[k] * C[i] * 0.000001;
                    T := TH + TB;
                    DutyCycle := (TH / (TH + TB)) * 100;
                    if ((TH >= 1.05) and (TH <= 1.09)) and ((TB >= 1.05) and (TB <= 1.09)) and ((DutyCycle > 49) and (DutyCycle < 51)) then
                    begin
                         Memo1.Lines.Add('R1=' + FloatToStr(R[j]) + ', R2=' + FloatToStr(R[k]) + ', C1=' + FloatToStr(C[i]) + ', TH=' + FloatToStr(TH) + ', TB=' + FloatToStr(TB) + ', T=' + FloatToStr(T) + ', DutyCycle=' + FloatToStr(DutyCycle));
                         Memo1.Refresh;
                         Application.ProcessMessages;
                    end;
               end;
          end;
     end;
end;
 
end.

Savez-vous pourquoi et comment corriger ce problème ?

Merci,
ZiP

[[ZiP]]

J'avais un problème dans mes boucles...

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
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
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
 
unit UPrincipale;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, Math;
 
type
  TForm1 = class(TForm)
    Memo1: TMemo;
    Button1: TButton;
    procedure Button1Click(Sender: TObject);
  private
    { Déclarations privées }
  public
    { Déclarations publiques }
  end;
 
var
  Form1: TForm1;
 
implementation
 
{$R *.dfm}
 
procedure TForm1.Button1Click(Sender: TObject);
var
     R         : array of Extended;
     C         : array of Extended;
     i,j,k     : Integer;
     T,TH,TB   : Extended;
     DutyCycle : Extended;
begin
     SetLength(R,145);
     R[0] := 1;
     R[1] := 1.2;
     R[2] := 1.5;
     R[3] := 1.8;
     R[4] := 2.2;
     R[5] := 2.7;
     R[6] := 3.3;
     R[7] := 3.9;
     R[8] := 4.7;
     R[9] := 5.6;
     R[10] := 6.8;
     R[11] := 8.2;
     R[12] := 10;
     R[13] := 11;
     R[14] := 12;
     R[15] := 13;
     R[16] := 15;
     R[17] := 16;
     R[18] := 18;
     R[19] := 20;
     R[20] := 22;
     R[21] := 24;
     R[22] := 27;
     R[23] := 30;
     R[24] := 33;
     R[25] := 36;
     R[26] := 39;
     R[27] := 43;
     R[28] := 47;
     R[29] := 51;
     R[30] := 56;
     R[31] := 62;
     R[32] := 68;
     R[33] := 75;
     R[34] := 82;
     R[35] := 91;
     R[36] := 100;
     R[37] := 110;
     R[38] := 120;
     R[39] := 130;
     R[40] := 150;
     R[41] := 160;
     R[42] := 180;
     R[43] := 200;
     R[44] := 220;
     R[45] := 240;
     R[46] := 270;
     R[47] := 300;
     R[48] := 330;
     R[49] := 360;
     R[50] := 390;
     R[51] := 430;
     R[52] := 470;
     R[53] := 510;
     R[54] := 560;
     R[55] := 620;
     R[56] := 680;
     R[57] := 750;
     R[58] := 820;
     R[59] := 910;
     R[60] := 1000;
     R[61] := 1100;
     R[62] := 1200;
     R[63] := 1300;
     R[64] := 1500;
     R[65] := 1600;
     R[66] := 1800;
     R[67] := 2000;
     R[68] := 2200;
     R[69] := 2400;
     R[70] := 2700;
     R[71] := 3000;
     R[72] := 3300;
     R[73] := 3600;
     R[74] := 3900;
     R[75] := 4300;
     R[76] := 4700;
     R[77] := 5100;
     R[78] := 5600;
     R[79] := 6200;
     R[80] := 6800;
     R[81] := 7500;
     R[82] := 8200;
     R[83] := 9100;
     R[84] := 10000;
     R[85] := 11000;
     R[86] := 12000;
     R[87] := 13000;
     R[88] := 15000;
     R[89] := 16000;
     R[90] := 18000;
     R[91] := 20000;
     R[92] := 22000;
     R[93] := 24000;
     R[94] := 27000;
     R[95] := 30000;
     R[96] := 33000;
     R[97] := 36000;
     R[98] := 39000;
     R[99] := 43000;
     R[100] := 47000;
     R[101] := 51000;
     R[102] := 56000;
     R[103] := 62000;
     R[104] := 68000;
     R[105] := 75000;
     R[106] := 82000;
     R[107] := 91000;
     R[108] := 100000;
     R[109] := 110000;
     R[110] := 120000;
     R[111] := 130000;
     R[112] := 150000;
     R[113] := 160000;
     R[114] := 180000;
     R[115] := 200000;
     R[116] := 220000;
     R[117] := 240000;
     R[118] := 270000;
     R[119] := 300000;
     R[120] := 330000;
     R[121] := 360000;
     R[122] := 390000;
     R[123] := 430000;
     R[124] := 470000;
     R[125] := 510000;
     R[126] := 560000;
     R[127] := 620000;
     R[128] := 680000;
     R[129] := 750000;
     R[130] := 820000;
     R[131] := 910000;
     R[132] := 1000000;
     R[133] := 1200000;
     R[134] := 1500000;
     R[135] := 1800000;
     R[136] := 2200000;
     R[137] := 2700000;
     R[138] := 3300000;
     R[139] := 3900000;
     R[140] := 4700000;
     R[141] := 5600000;
     R[142] := 6800000;
     R[143] := 8200000;
     R[144] := 10000000;
 
     SetLength(C,23);
     C[0] := 1;
     C[1] := 1.5;
     C[2] := 2.2;
     C[3] := 2.7;
     C[4] := 3;
     C[5] := 3.3;
     C[6] := 3.9;
     C[7] := 4.7;
     C[8] := 6.8;
     C[9] := 8.2;
     C[10] := 10;
     C[11] := 12;
     C[12] := 15;
     C[13] := 22;
     C[14] := 33;
     C[15] := 47;
     C[16] := 68;
     C[17] := 100;
     C[18] := 150;
     C[19] := 220;
     C[20] := 330;
     C[21] := 470;
     C[22] := 680;
 
     for i:=0 to Length(C)-1 do
     begin
          for j:=0 to Length(R)-1 do
          begin
               for k:=0 to Length(R)-1 do
               begin
                    TH := 0.693 * (R[k] + R[j]) * C[i] * 0.000001;
                    TB := 0.693 * R[k] * C[i] * 0.000001;
                    T := TH + TB;
                    DutyCycle := (TH / (TH + TB)) * 100;
                    if ((TH >= 1.06) and (TH <= 1.08)) and ((TB >= 1.06) and (TB <= 1.08)) and ((DutyCycle >= 49) and (DutyCycle <= 51)) then
                    begin
                         Memo1.Lines.Add('R1=' + FloatToStr(R[j]) + 'Ω, R2=' + FloatToStr(R[k]) + 'Ω, C1=' + FloatToStr(C[i]) + 'µF, TH=' + FloatToStrF(TH,ffFixed,8,3) + 's, TB=' + FloatToStrF(TB,ffFixed,8,3) + 's, T=' + FloatToStrF(T,ffFixed,8,3) + ', Rapport Cyclique=' + FloatToStrF(DutyCycle,ffFixed,8,2) + '%');
                         Memo1.Refresh;
                         Application.ProcessMessages;
                    end;
               end;
          end;
     end;
end;
 
end.

C'est résolu.

ZiP

[Dr.Who]

je me permet une petite correction :

ta méthode à ces performances :
5 833 Millions de cycles, 76 705 boucles en 3.06 secondes
soit 25 066 lp/s

avec la correction on passe à :
147 Millions de cycles et torches le tout en 0.08 secondes
soit 958 812,5 lp/s (+933 746,5)
soit un gains de performances d'environ 3800%

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
 
var
  {
   on deporte stR et stC en dehors de la boucle, on les initialises
   directement à la declaration.
  }
  stR : array[0..144] of extended = (
    1, 1.2, 1.5, 1.8, 2.2, 2.7, 3.3, 3.9, 4.7, 5.6, 6.8, 8.2, 10, 11, 12,
    13, 15, 16, 18, 20, 22, 24, 27, 30, 33, 36, 39, 43, 47, 51, 56, 62, 68,
    75, 82, 91, 100, 110, 120, 130, 150, 160, 180, 200, 220, 240, 270, 300,
    330, 360, 390, 430, 470, 510, 560, 620, 680, 750, 820, 910, 1000, 1100,
    1200, 1300, 1500, 1600, 1800, 2000, 2200, 2400, 2700, 3000, 3300, 3600,
    3900, 4300, 4700, 5100, 5600, 6200, 6800, 7500, 8200, 9100, 10000,
    11000, 12000, 13000, 15000, 16000, 18000, 20000, 22000, 24000, 27000,
    30000, 33000, 36000, 39000, 43000, 47000, 51000, 56000, 62000, 68000,
    75000, 82000, 91000, 100000, 110000, 120000, 130000, 150000, 160000,
    180000, 200000, 220000, 240000, 270000, 300000, 330000, 360000, 390000,
    430000, 470000, 510000, 560000, 620000, 680000, 750000, 820000, 910000,
    1000000, 1200000, 1500000, 1800000, 2200000, 2700000, 3300000, 3900000,
    4700000, 5600000, 6800000, 8200000, 10000000
  );
 
  stC : array[0..22] of extended = (
    1, 1.5, 2.2, 2.7, 3, 3.3, 3.9, 4.7, 6.8, 8.2, 10, 12, 15, 22, 33, 47,
    68, 100, 150, 220, 330, 470, 680
  );
 
 
procedure TForm15.Button1Click(Sender: TObject);
var
  i,j,k     : Integer;
  hstC, hstR: Integer;
  T,TH,TB   : Extended;
  DutyCycle : Extended;
begin
  {
   on utilise high() et pas Length()-1
   on mets les valeurs en "cache" pour eviter les calls sur 
   high dans les for
  }
  hstC := high(stC);
  hstR := high(stR);
 
  // on bloque le refresh de memo1, rafraîchir c'est mourir!
  Memo1.Lines.BeginUpdate;
  try
    // on vide memo1
    Memo1.Clear;
 
    for i := 0 to hstC do // 145
      for j := 0 to hstR do // * 23
        for k := 0 to hstR do // * 23 = 76 705 itérations
        begin
          TH := 0.693 * (stR[k] + stR[j]) * stC[i] * 0.000001;
          TB := 0.693 * stR[k] * stC[i] * 0.000001;
          T := TH + TB;
          DutyCycle := (TH / (TH + TB)) * 100;
          if ((TH >= 1.06) and (TH <= 1.08)) and
             ((TB >= 1.06) and (TB <= 1.08)) and
             ((DutyCycle >= 49) and (DutyCycle <= 51)) then
          begin
            {
              format au lieu de nombreux FloatToStr
              evite 7 calls qui ralentissent la boucle
            }
            Memo1.Lines.Add(
              format('R1=%g &Omega;, R2=%g &Omega;, C1=%g µF, TH=%.3f s, TB=%.3f s, T=%.3f, Rapport Cyclique=%.2f%%',
                     [stR[j], stR[k], stC[i], TH, TB, T, DutyCycle])
            );
 
            {
            on appel pas processmessage ... trop lent, inutile. 
            on appel encore moins le refresh de memo1 ... inutile également
            }
 
          end;
        end;
  finally
    // Un appel à EndUpdate suffit à afficher la liste.
    Memo1.Lines.EndUpdate;
  end;
end;

voila et bon code.

[[ZiP]]

Bonjour,

Merci pour votre correction, j'ai appris plusieurs choses que je ne connaissais pas !

C'est même 23 * 145 * 145 soit 483575 itérations.

Cordialement,
ZiP

[Dr.Who]

ha oui en effet j'avais inversé ...


le gain reste le même quand même ...