1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
import persalys
import openturns as ot
myStudy = persalys.Study("myStudy")
truncated_input = persalys.Input("X0", 0, ot.TruncatedDistribution(ot.Normal(), -1, 1))
user_defined_input = persalys.Input("X1", 2, ot.UserDefined([[1], [2], [3]], [1, 1, 2]))
hist_input = persalys.Input("X2", 1, ot.Histogram(0, [0.5, 1], [1, 0.5]))
normal_input = persalys.Input("X3", 1, ot.Normal(1, 2))
output = persalys.Output("Y0")
inputs = [truncated_input, user_defined_input, hist_input, normal_input]
outputs = [output]
formulas = ["X0+X1+X2+X3"]
model = persalys.SymbolicPhysicalModel("model", inputs, outputs, formulas)
myStudy.add(model)
script = myStudy.getPythonScript()
namespace = {}
exec(script, namespace)
for i in range(len(inputs)):
assert inputs[i].getDistribution() == namespace["inputs"][i].getDistribution()
|
