Discussion :

[elmoukrie]

Bonjour,

C'est pour réduire le temps des calculs d'une opération. On me conseille de paralléliser les opérations si cela est possible. Dans le code ci dessous, une matrice A (n,p) est générée. Chaque ligne (i) est un point (à six variables) pour lequel on calcule l'image Tmax(i). Les lignes sont indépendantes (heureusement) mais est ce que c'est possible de calculer en parallèle les images de toutes les lignes? On "soustraitra" le calcul chez un professionnel qui dispose de machines puissantes, mais il faut tout de même régler cette question de calcul en parallèle. Pour le code d'en bas, c'est sur la boucle for qu'il faut faire ce travail (je pense). On m'a parlé d'une commande forall au lieu de for, puis il faut que le code intègre d'autres commandes pour dire à la machine de gérer ce calcul en parallèle. Je n'ai que du brouillard de ces notions et j'espère que vous avez une idée. NB: je ne peux mettre la fonction hassane_25_essai

Merci d'avance.

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
 
 
n=5;  % n serait quelques centaines de milliers 8O
p=6; 
 
F_max= 0.002;
F_min= 0.001; 
ta_max=1200; 
ta_min= 600; 
U_max= 100;  
U_min= 1;   
Tj_max= 308.15;  
Tj_min=293.15; 
T0_max= 343.15;   
T0_min= 313.15;   
C_max=8;  
C_min=2.9;
 
An= lhsdesign(n,p);
 
 
 
Par_max=[F_max ta_max U_max Tj_max T0_max C_max];  
Par_min=[F_min ta_min U_min Tj_min T0_min C_min];  
 
A=An.*(Par_max-Par_min)+Par_min;    
 
TmaxA= []
 
for i=1:n
 
TmaxA(i)= hassane_25_essai(A(i,: ));
 
end
 
TmaxA

[Jerome Briot]

NB: je ne peux mettre la fonction hassane_25_essai

Dans ce cas, impossible de t'aider efficacement.

La parallélisation de calcul ne consiste pas seulement en l'ajout d'une simple commande (parfor dans ton cas).

Dans un premier temps, il faut s'assurer que le calcul est parallélisable… ce qui est impossible dans ton cas si tu ne nous montre pas la fonction concernée.

[elmoukrie]

Bonjour,

Je ne peux pas vous montrer la fonction et cela ne dépend pas de moi.

Mais pourquoi cela dépendrait de la fonction? La commande parfor est simple à utiliser d'après ce que j'ai vu sur le site de Matlab.

J'ai remplacé les boucles for que j'ai par parfor mais aucune amélioration n'est observée sur le temps de calcul.

Pourriez vous me dire plus?

[le fab]

parce que qu'il faut que chaque itération de ta boucle soit indépendante pour pouvoir envoyer des itérations séparées à différents cœurs de calcul
cela dit ils suffit d'essayer :
en remplacant for par parfor, et si tu constates une amélioration, c'est que ta boucle est parallelisable, sinon non
tu constates pas d'amélioration, ça veut dire qu'elle ne l'est pas

maintenant tu peux peut être la rendre parallélisable ...
mais sans le code on ne peut pas t'aider

[elmoukrie]

Chaque ligne de la matrice générée contient six valeurs de six variables qui seront intégrés dans la fonction appelé pour calculer l'image T_max (i). Les lignes sont indépendantes donc pour n lignes, je crois qu'il est possible de calculer ces n images séparément. Je ne vois pas en quoi montrer la fonction pourrait être utile! En tous cas, la fonction est en bas avec le petit code de calcul des images. ( Elle est en pièce jointe également)

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
 
 
function T_max = hassane_25_essai(Par)
 
clf
global k0_perh k_perh Ea_perh H_perh Keq_perh_30 Keq_perh Rperh k0_decomp_1 k_decomp_1 Ea_decomp_1 H_decomp_1 R_decomp_1 k0_decomp_2 k_decomp_2 Ea_decomp_2 H_decomp_2 R_decomp_2 k0_Ep_1 k_Ep_1 Ea_Ep_1 R_Ep_1 k0_Ep_2 Ea_Ep_2 k_Ep_2 R_Ep_2 k0_Ep_3 k_Ep_3 R_Ep_3 Ea_Ep_3 H_Ep H_RO k0_RO_W Ea_RO_W k_RO_W R_RO_W k0_RO_FA Ea_RO_FA k_RO_FA R_RO_FA k0_RO_PFA Ea_RO_PFA k_RO_PFA R_RO_PFA Temp_ref mj Cpj Qj UA R Cpr mr rho Tj0 Tj_degree Temp_degree qrx qexh Temp Vaq0 Vorg0 Vtot0 Taj FA_aj CpFA FA PFA W HP OA LA Int Ep RO Vol F Vaq tadd KFA K_FA W_initial CHP_initial 
 
 
% LES SIX VARIABLES D’ENTREE 
 
F= Par(1);  
tadd=Par(2); 
UA=Par(3);
Taj=Par(4);
Tj0=Par(5);
CHP_initial=Par(6);
 
k0_perh=0.15; 
Ea_perh=150000; 
H_perh=-5580; 
Keq_perh_30=0.96; 
k0_decomp_1=0.001; 
Ea_decomp_1=20000; 
H_decomp_1=-359000; 
k0_decomp_2=0.0009; 
Ea_decomp_2=20200; 
H_decomp_2=-163000; 
k0_Ep_1=0.00576; 
Ea_Ep_1=61600; 
k0_Ep_2=0.00437; 
Ea_Ep_2=15000; 
k0_Ep_3=0.004; 
Ea_Ep_3=15000; 
H_Ep=-230000; 
H_RO=-90000; 
k0_RO_W=0.000237; 
Ea_RO_W=150000; 
k0_RO_FA=0.00339; 
Ea_RO_FA=42100; 
k0_RO_PFA=0.0592; 
Ea_RO_PFA=70000; 
Temp_ref=66.85+273.15;
R=8.314; 
mj=0.3; 
Qj=15/60; 
Cpj=4000; 
Cpr= 2000; 
mr= 0.4; 
rho=1; 
Vaq0=0.14;   
Vorg0=0.12;              
Vtot0=Vaq0+Vorg0;  
FA_aj=24; 
CpFA=20; 
KFA=2.9E-6; 
K_FA=9;
W_initial=((1-(Vaq0*CHP_initial)*(34/(Vaq0*1000)))*1000)/18; 
 
tspan=[0:10:10000]; 
 
y0=[0 CHP_initial 0 W_initial 0.5 1.70 0.00 0.00 0.00 Tj0 Vtot0]; 
[t, y]=ode23s(@semibatch,tspan,y0); 
 
 
T_max = max(y(:,10));  %la variable sortie recherchée.
 
 
 
FA=y(:,1); 
HP=y(:,2);
PFA=y(:,3);
W=y(:,4);
OA=y(:,5);
LA=y(:,6);
Int=y(:,7);
Ep=y(:,8);
RO=y(:,9);
Temp=y(:,10);
Vol=y(:,11);
 
Vaq=Vol-Vorg0;
Tj_degree=Tj0-273.15;
Temp_degree=Temp-273.15;
 
k_perh=k0_perh*exp(-(Ea_perh/R)*((1./Temp)-(1/Temp_ref))); 
k_decomp_1=k0_decomp_1*exp(-(Ea_decomp_1/R)*((1./Temp)-(1/Temp_ref))); 
k_decomp_2=k0_decomp_2*exp(-(Ea_decomp_2/R)*((1./Temp)-(1/Temp_ref))); 
k_Ep_1=k0_Ep_1*exp(-(Ea_Ep_1/R)*((1./Temp)-(1/Temp_ref))); 
k_Ep_2=k0_Ep_2*exp(-(Ea_Ep_2/R)*((1./Temp)-(1/Temp_ref))); k_Ep_3=k0_Ep_3*exp(-(Ea_Ep_3/R)*((1./Temp)-(1/Temp_ref))); k_RO_W=k0_RO_W*exp(-(Ea_RO_W/R)*((1./Temp)-(1/Temp_ref))); k_RO_FA=k0_RO_FA*exp(-(Ea_RO_FA/R)*((1./Temp)-(1/Temp_ref))); k_RO_PFA=k0_RO_PFA*exp(-(Ea_RO_PFA/R)*((1./Temp)-(1/Temp_ref))); Keq_perh=Keq_perh_30*exp((-H_perh/R)*((1./Temp)-(1/303.15)));
 
Rperh=k_perh*(KFA^0.5).*((FA./W).^0.5).*(FA.*HP-(1./Keq_perh).*PFA.*W);
R_decomp_1=k_decomp_1.*PFA;
R_decomp_2=k_decomp_2.*PFA;
R_Ep_1=k_Ep_1.*PFA.*OA;
R_Ep_2=k_Ep_2.*PFA.*LA;
R_Ep_3=k_Ep_3.*PFA.*Int;
R_RO_W=k_RO_W*(KFA^0.5).*Ep.*((FA.*W).^0.5);
R_RO_FA=k_RO_FA*(KFA^0.5).*Ep.*FA.*((FA./W).^0.5);
R_RO_PFA=k_RO_PFA*(KFA^0.5).*Ep.*PFA.*((FA./W).^0.5);
 
qrx=-Vaq.*(Rperh*H_perh+R_decomp_1*H_decomp_1+R_decomp_2*H_decomp_2)-Vorg0*((R_Ep_1+R_Ep_2+R_Ep_3)*H_Ep+(R_RO_W+R_RO_FA+R_RO_PFA)*H_RO);
 
qexh=UA*(Tj0-Temp);
 
 
 
function dydt=semibatch(t,y)
 
global k0_perh k_perh Ea_perh H_perh Keq_perh_30 Keq_perh Rperh k0_decomp_1 k_decomp_1 Ea_decomp_1 H_decomp_1 R_decomp_1 k0_decomp_2 k_decomp_2 Ea_decomp_2 H_decomp_2 R_decomp_2 k0_Ep_1 k_Ep_1 Ea_Ep_1 R_Ep_1 k0_Ep_2 Ea_Ep_2 k_Ep_2 R_Ep_2 k0_Ep_3 k_Ep_3 R_Ep_3 Ea_Ep_3 H_Ep H_RO k0_RO_W Ea_RO_W k_RO_W R_RO_W k0_RO_FA Ea_RO_FA k_RO_FA R_RO_FA k0_RO_PFA Ea_RO_PFA k_RO_PFA R_RO_PFA Temp_ref mj Cpj Qj UA R Cpr mr rho Tj0 Tj_degree Temp_degree qrx qexh Temp Vaq0 Vorg0 Vtot0 Taj FA_aj CpFA FA PFA W HP OA LA Int Ep RO Vol F Vaq tadd KFA Vorg Vtot a K_FA mcp 
 
 
k_perh=k0_perh*exp(-(Ea_perh/R)*((1/y(10))-(1/Temp_ref))); 
k_decomp_1=k0_decomp_1*exp(-(Ea_decomp_1/R)*((1/y(10))-(1/Temp_ref))); k_decomp_2=k0_decomp_2*exp(-(Ea_decomp_2/R)*((1/y(10))-(1/Temp_ref))); 
 
k_Ep_1=k0_Ep_1*exp(-(Ea_Ep_1/R)*((1/y(10))-(1/Temp_ref))); 
k_Ep_2=k0_Ep_2*exp(-(Ea_Ep_2/R)*((1/y(10))-(1/Temp_ref))); k_Ep_3=k0_Ep_3*exp(-(Ea_Ep_3/R)*((1/y(10))-(1/Temp_ref))); k_RO_W=k0_RO_W*exp(-(Ea_RO_W/R)*((1/y(10))-(1/Temp_ref))); k_RO_FA=k0_RO_FA*exp(-(Ea_RO_FA/R)*((1/y(10))-(1/Temp_ref)));  k_RO_PFA=k0_RO_PFA*exp(-(Ea_RO_PFA/R)*((1/y(10))-(1/Temp_ref))); Keq_perh=Keq_perh_30*exp((-H_perh/R)*((1/y(10))-(1/303.15)));
 
 
Rperh=k_perh*(KFA^0.5)*((y(1)/y(4))^0.5)*(y(1)*y(2)-(1/Keq_perh)*y(3)*y(4));
R_decomp_1=k_decomp_1*y(3);
R_decomp_2=k_decomp_2*y(3);
R_Ep_1=k_Ep_1*y(3)*y(5);
R_Ep_2=k_Ep_2*y(3)*y(6);
R_Ep_3=k_Ep_3*y(3)*y(7);
R_RO_W=k_RO_W*(KFA^0.5)*y(8)*((y(1)*y(4))^0.5);
R_RO_FA=k_RO_FA*(KFA^0.5)*y(8)*y(1)*((y(1)/y(4))^0.5);
R_RO_PFA=k_RO_PFA*(KFA^0.5)*y(8)*y(3)*((y(1)/y(4))^0.5);
 
 
if t > tadd
    F=0;
end
 
 
Vaq=y(11)-Vorg0;
a=Vaq/y(11);
mcp=(Vaq*1.00+Vorg0*0.93)*Cpr;
 
 
dydt=[(1/(1+((1-a)/(a*K_FA))))*((F/Vaq)*(FA_aj-y(1))-Rperh+R_decomp_2+(1-a)*(R_Ep_1+R_Ep_2+R_Ep_3-R_RO_FA)/a);...
    -Rperh-(F*y(2)/Vaq); ((-F*y(3)/Vaq)+Rperh-R_decomp_2-R_decomp_1-(1-a)*(R_Ep_1+R_Ep_2+R_Ep_3+R_RO_PFA)/a);...
    Rperh+R_decomp_1-(F*y(4)/Vaq)-(1-a)*R_RO_W/a; -R_Ep_1; -R_Ep_2; R_Ep_2-R_Ep_3; (R_Ep_1+R_Ep_2+R_Ep_3)-R_RO_W-R_RO_FA-R_RO_PFA;...
    R_RO_W+R_RO_FA+R_RO_PFA; (1/mcp)*((-Vaq*(Rperh*H_perh+R_decomp_1*H_decomp_1+R_decomp_2*H_decomp_2)-Vorg0*((R_Ep_1+R_Ep_2+R_Ep_3)*H_Ep+(R_RO_W+R_RO_FA+R_RO_PFA)*H_RO))+UA*(Tj0-y(10))+FA_aj*F*CpFA*(Taj-y(10))); F];

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
 
 
n=5;  % n serait quelques centaines de milliers 8O
p=6; 
 
F_max= 0.002;
F_min= 0.001; 
ta_max=1200; 
ta_min= 600; 
U_max= 100;  
U_min= 1;   
Tj_max= 308.15;  
Tj_min=293.15; 
T0_max= 343.15;   
T0_min= 313.15;   
C_max=8;  
C_min=2.9;
 
 
 
An= lhsdesign(n,p);
 
 
 
Par_max=[F_max ta_max U_max Tj_max T0_max C_max];  
Par_min=[F_min ta_min U_min Tj_min T0_min C_min];  
 
A=An.*(Par_max-Par_min)+Par_min;    
 
T_max_A= []
 
for i=1:n
 
T_max_A(i)= hassane_25_essai(A(i,: ));
 
end
 
T_max_A

[le fab]

déjà il faut que tu supprime tes variables globales !!!
elles empêchent toute possibilité de paralleliser puisque théoriquement une itération peu modifier une variable globale pour l’itération suivante ...

[elmoukrie]

Bonjour,

Merci pour ton conseil! j'ai éliminé tous les variables globales et j'ai tout injecté dans la fonction dydt. je n'ai gardé évidement que les six variables qui changent.

A l'exécution, je vois effectivement un message qui dit que la tache est répartie sur un certain nombre de "cores" : deux dans mon cas. (alors que j'en ai 4).

LE TEMPS DE CALCUL EST NETTEMENT PETIT MAINTENANT.

Il est possible que je dise des bêtises mais comment savoir par exemple que cette parallélisation est optimisée, çàd qu' elle est faitE comme ça doit l'être et que tous les cores interviennent dans le calcul?

J'ai pas la toolbox statistique pour sobolset et lhsdesign sur le pc puissant qu'on m'a mis à disposition.Il parait que c'est une BONNE machine. Est ce qu'il y aurait une astuce?

[elmoukrie]

Bonjour,

Je n'ai toujours pas de réponse pour mon problème. Je vous prie de m'apporter des réponses si vous avez des idées. On m'a demandé d'enlever les paramètres globaux et c'est fait.

Auriez vous des propositions?

Merci