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
|
class=SobolAnalysis name=aSobol physicalModel=aModel interestVariables=[y0] maximumCalls=1000 maximumCoefficientOfVariation=0.01 maximumElapsedTime=30 seed=2 maximum confidence interval length=-1 blockSize=4 replicationSize=250 confidenceLevel=0.95
result= class=SobolResult name=Unnamed outputNames=[y0] firstOrderIndices=class=Sample name=Unnamed implementation=class=SampleImplementation name=Unnamed size=1 dimension=2 description=[x0,x1] data=[[0.900967,0.0220335]] totalIndices=class=Sample name=Unnamed implementation=class=SampleImplementation name=Unnamed size=1 dimension=2 data=[[0.985481,0.0912614]] firstOrderIndicesInterval=[[0.877526, 0.924408]
[-0.108137, 0.152204]] totalIndicesInterval=[[0.844158, 1.1268]
[0.0665421, 0.115981]] callsNumber=1000 confidenceIntervalLength=0.207877
class=SobolAnalysis name=aSobol2 physicalModel=aModel interestVariables=[y0] maximumCalls=inf maximumCoefficientOfVariation=0.01 maximumElapsedTime=100000 seed=0 maximum confidence interval length=0.2 blockSize=1 replicationSize=1000 confidenceLevel=0.95
result= class=SobolResult name=Unnamed outputNames=[y0] firstOrderIndices=class=Sample name=Unnamed implementation=class=SampleImplementation name=Unnamed size=1 dimension=2 description=[x0,x1] data=[[0.892832,0.153005]] totalIndices=class=Sample name=Unnamed implementation=class=SampleImplementation name=Unnamed size=1 dimension=2 data=[[0.870144,0.107257]] firstOrderIndicesInterval=[[0.878084, 0.907581]
[0.0832421, 0.222768]] totalIndicesInterval=[[0.801019, 0.939269]
[0.0924361, 0.122078]] callsNumber=4000 confidenceIntervalLength=0.102617
class=SobolAnalysis name=aSobol3 physicalModel=aModel interestVariables=[y1] maximumCalls=1000 maximumCoefficientOfVariation=0.01 maximumElapsedTime=60 seed=2 maximum confidence interval length=0.01 blockSize=4 replicationSize=200 confidenceLevel=0.95
result= class=SobolResult name=Unnamed outputNames=[y1] firstOrderIndices=class=Sample name=Unnamed implementation=class=SampleImplementation name=Unnamed size=1 dimension=2 description=[x0,x1] data=[[0.90921,0.115756]] totalIndices=class=Sample name=Unnamed implementation=class=SampleImplementation name=Unnamed size=1 dimension=2 data=[[0.92307,0.0973193]] firstOrderIndicesInterval=[[0.886406, 0.932014]
[-0.00958092, 0.241092]] totalIndicesInterval=[[0.78735, 1.05879]
[0.0753105, 0.119328]] callsNumber=1600 confidenceIntervalLength=0.199635
#!/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)
aSobol = persalys.SobolAnalysis('aSobol', aModel)
interestVariables = ['y0']
aSobol.setInterestVariables(interestVariables)
aSobol.setMaximumCalls(1000)
aSobol.setMaximumConfidenceIntervalLength(-1)
aSobol.setMaximumElapsedTime(30)
aSobol.setBlockSize(4)
aSobol.setReplicationSize(250)
aSobol.setSeed(2)
aSobol.setConfidenceLevel(0.95)
myStudy.add(aSobol)
aSobol2 = persalys.SobolAnalysis('aSobol2', aModel)
interestVariables = ['y0']
aSobol2.setInterestVariables(interestVariables)
aSobol2.setMaximumConfidenceIntervalLength(0.2)
aSobol2.setMaximumElapsedTime(100000)
aSobol2.setBlockSize(1)
aSobol2.setReplicationSize(1000)
aSobol2.setSeed(0)
aSobol2.setConfidenceLevel(0.95)
myStudy.add(aSobol2)
aSobol3 = persalys.SobolAnalysis('aSobol3', aModel)
interestVariables = ['y1']
aSobol3.setInterestVariables(interestVariables)
aSobol3.setMaximumCalls(1000)
aSobol3.setMaximumConfidenceIntervalLength(0.01)
aSobol3.setMaximumElapsedTime(60)
aSobol3.setBlockSize(4)
aSobol3.setReplicationSize(200)
aSobol3.setSeed(2)
aSobol3.setConfidenceLevel(0.95)
myStudy.add(aSobol3)
|
