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
|
obs= class=Observations name=obs physicalModel=aModelPhys observedInputs=[X3,X1] observedOutputs=[Y2,Y1] fileName=obs_normal2.csv inputColumns=[7,3] outputColumns=[6,2]
InvalidArgumentException occurred: True
InvalidArgumentException occurred: True
obs2= class=Observations name=obs2 physicalModel=aModelPhys observedInputs=[X3,X1] observedOutputs=[Y2,Y1]
#!/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, '')
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', '')
inputs = [X0, X1, X2, X3]
outputs = [Y0, Y1, Y2]
formulas = ['sin(X0)+8*X1', 'sin(X0)+8*X1', 'sin(X0)+8*X1']
aModelPhys = persalys.SymbolicPhysicalModel('aModelPhys', inputs, outputs, formulas)
anOTStudy.add(aModelPhys)
inputColumns = [4, 5, 1]
outputColumns = [3, 0]
inputNames = ['X3', 'X2', 'X1']
outputNames = ['Y1', 'Y0']
obs = persalys.Observations('obs', aModelPhys, 'obs_normal2.csv', inputColumns, outputColumns, inputNames, outputNames)
anOTStudy.add(obs)
inputSample = ot.Sample([
[1.20548, 1.43725],
[0.261018, -0.0907838],
[0.44549, -0.856712],
[1.78236, -0.241223],
[1.00434, 0.29995],
[-0.752817, 1.85579],
[-0.743622, 0.538071],
[-0.181004, 1.21381],
[-0.89182, 0.677958],
[-1.66086, -0.633066]]
)
inputSample.setDescription(['X3', 'X1'])
outputSample = ot.Sample([
[-0.438266, -0.355007],
[-0.470526, -1.31178],
[-0.560206, -1.03808],
[0.351418, -0.721533],
[1.36783, 0.539345],
[-0.382992, -0.671291],
[0.716353, 0.655027],
[0.377922, -0.353552],
[0.103474, -0.483642],
[-0.506925, -0.510764]]
)
outputSample.setDescription(['Y2', 'Y1'])
obs2 = persalys.Observations('obs2', aModelPhys, inputSample, outputSample)
anOTStudy.add(obs2)
|
