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
|
#! /usr/bin/env python
import openturns as ot
import openturns.testing
import persalys
import os
ot.RandomGenerator.SetSeed(0)
anOTStudy = persalys.Study("anOTStudy")
# Model
X0 = persalys.Input("X0", 1, ot.Normal())
X1 = persalys.Input("X1", 2)
X2 = persalys.Input("X2", 2)
X3 = persalys.Input("X3", 2)
Y0 = persalys.Output("Y0")
Y1 = persalys.Output("Y1")
Y2 = persalys.Output("Y2")
model = persalys.SymbolicPhysicalModel(
"aModelPhys", [X0, X1, X2, X3], [Y0, Y1, Y2], ["sin(X0)+8*X1"] * 3
)
anOTStudy.add(model)
# Observations ##
filename = "obs_normal2.csv"
sample = ot.Normal(8).getSample(10)
sample.exportToCSVFile(filename)
obs = persalys.Observations(
"obs", model, filename, [7, 3], [6, 2], ["X3", "X1"], ["Y2", "Y1"]
)
anOTStudy.add(obs)
print("obs=", obs)
openturns.testing.assert_almost_equal(
obs.getSample(), sample.getMarginal([3, 7, 2, 6]), 1e-16
)
# - reload file
sample2 = sample * 3
sample2.exportToCSVFile(filename)
obs.setFileName(filename)
openturns.testing.assert_almost_equal(
obs.getSample(), sample2.getMarginal([3, 7, 2, 6]), 1e-16
)
# - change columns
obs.setColumns([4, 5], ["X3", "X1"], [3, 0], ["Y2", "Y1"])
openturns.testing.assert_almost_equal(
obs.getSample(), sample2.getMarginal([5, 4, 0, 3]), 1e-16
)
try:
obs.setColumns([4, 5, 0], ["X3", "X1"], [3, 0], ["Y2", "Y1"])
except Exception as e:
print(
"InvalidArgumentException occurred: %s" % ("InvalidArgumentException" in str(e))
)
try:
obs.setColumns([4, 5, 0], ["X0", "X1", "X2"], [], [])
except Exception as e:
print(
"InvalidArgumentException occurred: %s"
% (
"InvalidArgumentException : Define observations for at least an output"
in str(e)
)
)
obs.setColumns([4, 5, 1], ["X0", "X1", "X2"], [3, 0], ["Y1", "Y2"])
assert obs.getSample().getDescription() == ["X0", "X1", "X2", "Y1", "Y2"]
# - change observed variables
obs.setColumns([4, 5, 1], ["X3", "X2", "X1"], [3, 0], ["Y1", "Y0"])
assert obs.getSample().getDescription() == ["X1", "X2", "X3", "Y0", "Y1"]
# Observations ##
inS = sample.getMarginal([7, 3])
inS.setDescription(["X3", "X1"])
outS = sample.getMarginal([6, 2])
outS.setDescription(["Y2", "Y1"])
obs2 = persalys.Observations("obs2", model, inS, outS)
anOTStudy.add(obs2)
print("obs2=", obs2)
openturns.testing.assert_almost_equal(
obs2.getSample(), sample.getMarginal([3, 7, 2, 6]), 1e-16
)
# script
script = anOTStudy.getPythonScript()
print(script)
exec(script)
os.remove(filename)
|
