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
|
class=SRCAnalysis name=aSRC physicalModel=aModel interestVariables=[y0] simulationsNumber=1000 seed=2 blockSize=1
result= class=SRCResult name=Unnamed outputNames=[y0] indices=class=Sample name=Unnamed implementation=class=SampleImplementation name=Unnamed size=1 dimension=2 description=[x0,x1] data=[[0.90576,0.0939758]] indices intervals=[[0.881219, 0.9375]
[0.0849848, 0.105183]] signed indices=class=Sample name=Unnamed implementation=class=SampleImplementation name=Unnamed size=1 dimension=2 data=[[0.951714,0.306555]] signed indices intervals=[[0.938732, 0.968245]
[0.291522, 0.32432]] r2=class=Point name=Unnamed dimension=1 values=[0.999735] calls number=1000
class=SRCAnalysis name=aSRC2 physicalModel=aModel interestVariables=[y1] simulationsNumber=1000 seed=2 blockSize=1
result= class=SRCResult name=Unnamed outputNames=[y1] indices=class=Sample name=Unnamed implementation=class=SampleImplementation name=Unnamed size=1 dimension=2 description=[x0,x1] data=[[0.90576,0.0939758]] indices intervals=[[0.881219, 0.9375]
[0.0849848, 0.105183]] signed indices=class=Sample name=Unnamed implementation=class=SampleImplementation name=Unnamed size=1 dimension=2 data=[[0.951714,0.306555]] signed indices intervals=[[0.938732, 0.968245]
[0.291522, 0.32432]] r2=class=Point name=Unnamed dimension=1 values=[0.999735] calls number=1000
#!/usr/bin/env python
import openturns as ot
import persalys
myStudy = persalys.Study('myStudy')
persalys.Study.Add(myStudy)
dist_x0 = ot.Normal(0, 1)
x0 = persalys.Input('x0', 0, dist_x0, '')
dist_x1 = ot.Normal(0, 1)
x1 = persalys.Input('x1', 0, dist_x1, '')
x2 = persalys.Input('x2', 10, '')
fake_y0 = persalys.Output('fake_y0', '')
fake_y0.setIsSelected(False)
y0 = persalys.Output('y0', '')
y1 = persalys.Output('y1', '')
inputs = [x0, x1, x2]
outputs = [fake_y0, y0, y1]
formulas = ['x0', '10+3*x0+x1', '3*x0 + x1 + x2']
aModel = persalys.SymbolicPhysicalModel('aModel', inputs, outputs, formulas)
myStudy.add(aModel)
aSRC = persalys.SRCAnalysis('aSRC', aModel, 1000)
interestVariables = ['y0']
aSRC.setInterestVariables(interestVariables)
aSRC.setBlockSize(1)
aSRC.setSeed(2)
myStudy.add(aSRC)
aSRC2 = persalys.SRCAnalysis('aSRC2', aModel, 1000)
interestVariables = ['y1']
aSRC2.setInterestVariables(interestVariables)
aSRC2.setBlockSize(1)
aSRC2.setSeed(2)
myStudy.add(aSRC2)
|
