1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
-
-
+
+
-
+
-
-
+
+
-
-
+
+
|
%% CR for conical tank example - sensitivity
in "conical.cr";
OPERATOR sconical;
%% The actual port (1)
FOR ALL rho,g,V_0,V_0s,V,sV LET
sconical(rho,g,V_0,V_0s,effort,1,
FOR ALL COMP, rho,g,V_0,V_0s,V,sV LET
sconical(COMP, rho,g,V_0,V_0s,effort,1,
V,state,1,
sV,state,2)
= conical(rho,g,V_0,effort,1,V,state,1);
= conical(COMP, rho,g,V_0,effort,1,V,state,1);
%% The sensitivity port (2)
FOR ALL rho,g,V_0,V_0s,V,sV LET
sconical(rho,g,V_0,V_0s,effort,2,
FOR ALL COMP, rho,g,V_0,V_0s,V,sV LET
sconical(COMP, rho,g,V_0,V_0s,effort,2,
V,state,1,
sV,state,2)
= df(conical(rho,g,V_0,effort,1,V,state,1),V) * sV
+ df(conical(rho,g,V_0,effort,1,V,state,1),V_0)* V_0s;
= df(conical(COMP, rho,g,V_0,effort,1,V,state,1),V) * sV
+ df(conical(COMP, rho,g,V_0,effort,1,V,state,1),V_0)* V_0s;
END;;
|