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
|
outs= [ Y0 ]
0 : [ 15.1833 ]
1 : [ 18.3833 ]
2 : [ 15.2912 ]
3 : [ 18.4912 ]
outs= [ Y0 ]
0 : [ -2.93473 ]
1 : [ 3.73435 ]
2 : [ 6.13775 ]
3 : [ 2.80071 ]
4 : [ -11.2465 ]
5 : [ -1.20628 ]
6 : [ 2.05204 ]
7 : [ -6.42253 ]
8 : [ 3.26745 ]
9 : [ -5.27323 ]
outs= [ Y0 ]
0 : [ 0 ]
1 : [ -4.77142 ]
2 : [ 4.77142 ]
3 : [ -2.86239 ]
4 : [ 10.1157 ]
5 : [ -8.88952 ]
6 : [ 1.63621 ]
7 : [ -4.68548 ]
8 : [ 7.56672 ]
9 : [ -11.2736 ]
outs= [ Y0 ]
0 : [ 3.72534 ]
1 : [ -1.06474 ]
2 : [ 7.66635 ]
3 : [ 2.30271 ]
4 : [ -5.13293 ]
5 : [ -14.7875 ]
6 : [ -2.70615 ]
7 : [ 5.56753 ]
8 : [ 14.1784 ]
9 : [ -7.54187 ]
outs= [ Y0 ]
0 : [ -1.46262 ]
1 : [ -4.13728 ]
2 : [ 0.470196 ]
3 : [ -7.44676 ]
4 : [ 9.94891 ]
5 : [ 0.487362 ]
6 : [ -0.81332 ]
7 : [ -1.40083 ]
8 : [ 7.6165 ]
9 : [ -4.50556 ]
outs= [ Y0 ]
0 : [ -0.173083 ]
1 : [ -3.33481 ]
2 : [ 1.16288 ]
3 : [ -2.67599 ]
4 : [ 7.9576 ]
5 : [ 8.19727 ]
6 : [ 4.53184 ]
7 : [ -1.96082 ]
8 : [ -11.6724 ]
9 : [ -0.4845 ]
outs= [ Y0 ]
0 : [ -1.35728 ]
1 : [ 1.89976 ]
2 : [ 0.555441 ]
3 : [ 13.583 ]
4 : [ -7.39928 ]
5 : [ -0.522775 ]
6 : [ 13.3877 ]
7 : [ -11.1484 ]
8 : [ 2.78968 ]
9 : [ -2.14925 ]
outs= [ Y0 ]
0 : [ -3.38162 ]
1 : [ 0.642078 ]
2 : [ -19.9116 ]
3 : [ 2.39837 ]
4 : [ -2.40688 ]
5 : [ -2.12083 ]
6 : [ -4.27945 ]
7 : [ -0.0333236 ]
8 : [ 2.95844 ]
9 : [ 14.3288 ]
outs= [ Y0 ]
0 : [ -3.38162 ]
1 : [ 0.642078 ]
2 : [ -19.9116 ]
3 : [ 2.39837 ]
4 : [ -2.40688 ]
5 : [ -2.12083 ]
6 : [ -4.27945 ]
7 : [ -0.0333236 ]
8 : [ 2.95844 ]
9 : [ 14.3288 ]
outs= [ Y0 ]
0 : [ -3.0804 ]
1 : [ -7.10749 ]
2 : [ 11.3192 ]
3 : [ 0.123019 ]
4 : [ -3.3712 ]
5 : [ -2.64171 ]
6 : [ 10.7842 ]
7 : [ 0.275507 ]
8 : [ 4.78088 ]
9 : [ 5.79249 ]
#!/usr/bin/env python
import openturns as ot
import persalys
anOTStudy = persalys.Study('anOTStudy')
persalys.Study.Add(anOTStudy)
dist_X0 = ot.Normal(0, 1)
X0 = persalys.Input('X0', 1, dist_X0, '')
dist_X1 = ot.Normal(0, 1)
X1 = persalys.Input('X1', 2, dist_X1, '')
Y0 = persalys.Output('Y0', '')
inputs = [X0, X1]
outputs = [Y0]
formulas = ['sin(X0)+8*X1']
aModelPhys = persalys.SymbolicPhysicalModel('aModelPhys', inputs, outputs, formulas)
anOTStudy.add(aModelPhys)
values = [[0.9, 1.1],
[1.8, 2.2]]
aDesign_1 = persalys.GridDesignOfExperiment('aDesign_1', aModelPhys, values)
aDesign_1.setBlockSize(1)
interestVariables = ['Y0']
aDesign_1.setInterestVariables(interestVariables)
anOTStudy.add(aDesign_1)
inputColumns = [0, 2]
aDesign_2 = persalys.ImportedDesignOfExperiment('aDesign_2', aModelPhys, 'normal.csv', inputColumns)
aDesign_2.setBlockSize(1)
interestVariables = ['Y0']
aDesign_2.setInterestVariables(interestVariables)
anOTStudy.add(aDesign_2)
aDesign_3 = persalys.ProbabilisticDesignOfExperiment('aDesign_3', aModelPhys, 10, 'QUASI_MONTE_CARLO')
aDesign_3.setSeed(0)
aDesign_3.setBlockSize(1)
interestVariables = ['Y0']
aDesign_3.setInterestVariables(interestVariables)
anOTStudy.add(aDesign_3)
aDesign_4 = persalys.FixedDesignOfExperiment('aDesign_4', aModelPhys)
inputSample = [
[-0.0791479, 0.475551],
[-1.79937, -0.0113437],
[0.409358, 0.908542],
[1.33205, 0.166384],
[0.585946, -0.710739],
[-1.26768, -1.72914],
[-0.52963, -0.275117],
[0.241872, 0.666002],
[-0.451024, 1.82679],
[1.03206, -1.05003]]
aDesign_4.setOriginalInputSample(inputSample)
aDesign_4.setBlockSize(1)
interestVariables = ['Y0']
aDesign_4.setInterestVariables(interestVariables)
anOTStudy.add(aDesign_4)
#!/usr/bin/env python
import openturns as ot
import persalys
anOTStudy2 = persalys.Study('anOTStudy2')
persalys.Study.Add(anOTStudy2)
dist_X0 = ot.Normal(0, 1)
X0 = persalys.Input('X0', 1, dist_X0, '')
dist_X1 = ot.Normal(0, 1)
X1 = persalys.Input('X1', 2, dist_X1, '')
dist_X2 = ot.Normal(0, 1)
X2 = persalys.Input('X2', 3, dist_X2, '')
Y0 = persalys.Output('Y0', '')
inputs = [X0, X1, X2]
outputs = [Y0]
formulas = ['sin(X0)+8*X1*X2']
aModelPhys = persalys.SymbolicPhysicalModel('aModelPhys', inputs, outputs, formulas)
anOTStudy2.add(aModelPhys)
aDesign_5 = persalys.ProbabilisticDesignOfExperiment('aDesign_5', aModelPhys, 10, 'SALHS')
aDesign_5.setSeed(0)
aDesign_5.setBlockSize(1)
interestVariables = ['Y0']
aDesign_5.setInterestVariables(interestVariables)
anOTStudy2.add(aDesign_5)
aDesign_7 = persalys.ProbabilisticDesignOfExperiment('aDesign_7', aModelPhys, 10, 'SALHS')
aDesign_7.setSeed(0)
aDesign_7.setBlockSize(1)
interestVariables = ['Y0']
aDesign_7.setInterestVariables(interestVariables)
anOTStudy2.add(aDesign_7)
aDesign_8 = persalys.ProbabilisticDesignOfExperiment('aDesign_8', aModelPhys, 10, 'MCLHS')
aDesign_8.setSeed(0)
aDesign_8.setBlockSize(1)
interestVariables = ['Y0']
aDesign_8.setInterestVariables(interestVariables)
anOTStudy2.add(aDesign_8)
aDesign_9 = persalys.ProbabilisticDesignOfExperiment('aDesign_9', aModelPhys, 10, 'MCLHS')
aDesign_9.setSeed(0)
aDesign_9.setBlockSize(1)
interestVariables = ['Y0']
aDesign_9.setInterestVariables(interestVariables)
anOTStudy2.add(aDesign_9)
aDesign_10 = persalys.ProbabilisticDesignOfExperiment('aDesign_10', aModelPhys, 10, 'MCLHS')
aDesign_10.setSeed(0)
aDesign_10.setBlockSize(1)
interestVariables = ['Y0']
aDesign_10.setInterestVariables(interestVariables)
anOTStudy2.add(aDesign_10)
|
